Search

Finding a software engineering approach that allows for portability, rapid development, open collaboration, and performance for high performance computing on GPUs and CPUs is a challenge. We implement a portability scheme using the Numba compiler for Python in Monte Carlo / Dynamic Code (MC/DC), a new neutron transport application for rapid Monte Carlo methods development. Using this scheme, we have built MC/DC as a single source, single language, single compiler application that can run as a pure Python, compiled CPU, or compiled GPU solver. In GPU mode, we use Numba paired with an asynchronous GPU scheduler called Harmonize to increase GPU performance. We present performance results for a time-dependent problem on both the CPU and GPU and compare them to a production code., Comment: 11 pages, 5 figures

The techniques used to generate pseudo-random numbers for Monte Carlo (MC) applications bear many implications on the quality and speed of that programs work. As a random number generator (RNG) slows, the production of random numbers begins to dominate runtime. As RNG output grows in correlation, the final product becomes less reliable. These difficulties are further compounded by the need for reproducibility and parallelism. For reproducibility, the numbers generated to determine any outcome must be the same each time a simulation is run. However, the concurrency that comes with most parallelism introduces race conditions. To have both reproducibility and concurrency, separate RNG states must be tracked for each independently schedulable unit of simulation, forming independent random number streams. We propose an alternative to the stride-based parallel LCG seeding approach that scales more practically with increased concurrency and workload by generating seeds through hashing and allowing for repeated outputs. Data gathered from normality tests of tally results from simple MC transport benchmark calculations indicates that the proposed hash-based RNG does not significantly affect the tally result normality property as compared to the conventional stride-based RNG., Comment: 4 Pages, 3 Figures, Submitted to 2024 ANS Annual Conference

Airway-related quantitative imaging biomarkers are crucial for examination, diagnosis, and prognosis in pulmonary diseases. However, the manual delineation of airway trees remains prohibitively time-consuming. While significant efforts have been made towards enhancing airway modelling, current public-available datasets concentrate on lung diseases with moderate morphological variations. The intricate honeycombing patterns present in the lung tissues of fibrotic lung disease patients exacerbate the challenges, often leading to various prediction errors. To address this issue, the 'Airway-Informed Quantitative CT Imaging Biomarker for Fibrotic Lung Disease 2023' (AIIB23) competition was organized in conjunction with the official 2023 International Conference on Medical Image Computing and Computer Assisted Intervention (MICCAI). The airway structures were meticulously annotated by three experienced radiologists. Competitors were encouraged to develop automatic airway segmentation models with high robustness and generalization abilities, followed by exploring the most correlated QIB of mortality prediction. A training set of 120 high-resolution computerised tomography (HRCT) scans were publicly released with expert annotations and mortality status. The online validation set incorporated 52 HRCT scans from patients with fibrotic lung disease and the offline test set included 140 cases from fibrosis and COVID-19 patients. The results have shown that the capacity of extracting airway trees from patients with fibrotic lung disease could be enhanced by introducing voxel-wise weighted general union loss and continuity loss. In addition to the competitive image biomarkers for prognosis, a strong airway-derived biomarker (Hazard ratio>1.5, p<0.0001) was revealed for survival prognostication compared with existing clinical measurements, clinician assessment and AI-based biomarkers., Comment: 19 pages

Among the sample of Galactic transient X-ray binaries (SXTs) discovered to date, about 70 have been proposed as likely candidates to host a black hole. Yet, only 19 have been dynamically confirmed. Such a reliable confirmation requires phase-resolved spectroscopy of their companion stars, which is generally feasible when the system is in a quiescent state. However, since most of the SXT population lies in the galactic plane, which is strongly affected by interstellar extinction, their optical brightness during quiescence usually falls beyond the capabilities of the current instrumentation ($R\gtrsim22$). To overcome these limitations and thereby increase the number of confirmed Galactic black holes, a correlation between the full-width at half maximum (FWHM) of the H$\alpha$ line and the semi-amplitude of the donor's radial velocity curve ($K_2$) was presented in the past. Here, we extend the FWHM-$K_2$ correlation to the near-infrared (NIR), exploiting disc lines such as He I $\lambda$10830, Pa$\gamma$, and Br$\gamma$, in a sample of dynamically confirmed black-hole SXTs. We obtain $K_2 = 0.22(3) ~\textrm{FWHM}$, in good agreement with the optical correlation derived using H$\alpha$. The similarity of the two correlations seems to imply that the widths of H$\alpha$ and the NIR lines are consistent in quiescence. When combined with information on orbital periods, the NIR correlation allows us to constrain the mass of the compact object of systems in quiescence by using single-epoch spectroscopy. We anticipate that this new correlation will give access to highly reddened black-hole SXTs, which cannot be otherwise studied at optical wavelengths., Comment: Accepted for publication in A&A Letters

We present evidence for a 0.278(8) d (=6.7 h) orbital period in the X-ray transient GRS 1716-249 (=N Oph 93), based on a superhump modulation detected during the 1995 mini-outburst plus ellipsoidal variability in quiescence. With a quiescent magnitude of r=23.19+-0.15 N Oph 93 is too faint to warrant a full dynamical study through dedicated time-resolved spectroscopy. Instead, we apply the FWHM-K2 correlation to the disc Halpha emission line detected in Gran Telescopio Canarias spectra and obtain K2=521+-52 km/s. This leads to a mass function f(M)=4.1+-1.2 Msun, thus indicating the presence of a black hole in this historic X-ray transient. Furthermore, from the depth of the Halpha trough and the quiescent light curve we constrain the binary inclination to i=61+-15 deg, while the detection of superhumps sets an upper limit to the donor to compact star mass ratio q=M2/M1<=0.25. Our de-reddened (r-i) colour is consistent with a ~K6 main sequence star that fills its Roche lobe in a 0.278 d orbit. Using all this information we derive a compact object mass M1=6.4+3.2-2.0 Msun at 68 per cent confidence. We also constrain the distance to GRS 1716-249 to 6.9+-1.1 kpc, placing the binary ~0.8 kpc above the Galactic Plane, in support of a large natal kick., Comment: Accepted for publication in MNRAS, 12 pages, 9 figures, 2 Tables

Accreting white dwarfs are known to show signatures of wind-type outflows in the ultraviolet. At optical wavelengths, however, wind detections have only been reported for a few sources. We present GTC-10.4m optical spectroscopy of four accreting white dwarfs (BZ Cam, V751 Cyg, MV Lyr, and V425 Cas) observed during luminous epochs, when their optical emission is expected to be dominated by the accretion disc. We focused the analysis on four emission lines: H$\alpha$ and He I $\lambda$5876, $\lambda$6678, $\lambda$7065. Line profiles are complex and variable on short (minutes) and long (days to weeks) time scales, with transient absorption and emission components. Among them, we detect strong blue-shifted absorptions at $\gtrsim 1000$ km s$^{-1}$. These high-velocity components, present only in the blue wing of the emission lines, are observed in all four sources and could be associated with accretion disc winds. For MV Lyr and V425 Cas, these would represent the first detection of optical outflows in these objects, while in the case of BZ Cam and V751 Cyg, the presence of outflows has been previously reported. This study suggests that, in addition to ultraviolet winds, optical outflows might be also common in accreting white dwarfs. We discuss the observational properties of these winds and their possible similarity to those detected in accreting black holes and neutrons stars., Comment: Accepted for publication in A&A

We report on the first simultaneous high-time resolution X-ray and infrared (IR) observations of a neutron star low mass X-ray binary in its hard state. We performed $\approx 2\,$h of simultaneous observations of 4U 1728-34 using HAWK-I@VLT, XMM-Newton and NuSTAR. The source displayed significant X-ray and IR variability down to sub-second timescales. By measuring the cross-correlation function between the infrared and X-ray lightcurves, we discovered a significant correlation with an infrared lead of $\approx 30-40\,$ms with respect to the X-rays. We analysed the X-ray energy dependence of the lag, finding a marginal increase towards higher energies. Given the sign of the lag, we interpret this as possible evidence of Comptonization from external seed photons. We discuss the origin of the IR seed photons in terms of cyclo-synchrotron radiation from an extended hot flow. Finally, we also observed the IR counterpart of a type-I X-ray burst, with a delay of $\approx7.2\,$s. Although some additional effects may be at play, by assuming that this lag is due to light travel time between the central object and the companion star, we find that 4U 1728-34 must have an orbital period longer than $3\,$h and an inclination higher than 8$^\circ$., Comment: Accepted for publication in MNRAS, 11 pages 7 figures

We prove the large deviation principle for several entropy and cross entropy estimators based on return times and waiting times on shift spaces over finite alphabets. We consider shift-invariant probability measures satisfying some decoupling conditions which imply no form of mixing nor ergodicity. We establish precise relations between the rate functions of the different estimators, and between these rate functions and the corresponding pressures, one of which is the R\'enyi entropy function. For the most commonly used definition of return times, the large-deviation rate function is proved to be nonconvex, except in marginal cases. The results apply in particular to irreducible Markov chains, equilibrium measures for Bowen-regular potentials, g-measures, invariant Gibbs states for absolutely summable interactions in statistical mechanics, and also to probability measures which may be far from Gibbsian, including some hidden Markov models and repeated quantum measurement processes., Comment: Minor corrections, additions and clarifications. Updated references

Ultracompact X-ray binaries (UCXBs) are a distinctive but elusive family of low-mass X-ray binaries (LMXBs) characterised by their tight orbits and degenerate donor stars. Here we present UltraCompCAT, the first online and comprehensive catalogue of UCXBs. The initial version of UltraCompCAT comprises 49 sources, including 20 'confirmed' UCXBs (those with a measured orbital period shorter than 80 min) and 25 systems that we label as 'candidate' based on their multi-wavelength phenomenology. For completeness, we also include four LMXBs with orbital periods in the range of 80 to 120 min, since they might be related (e.g. close progenitors) or even part of the UCXB population that evolved towards longer periods. We discuss the orbital period and Galactic distribution of the catalogue's sample. We provide evidence for the presence of at least two separate groups of UCXBs. One formed by persistent systems with orbital periods shorter than 30 min and a second group of transient objects (70 per cent) with periods in the range of 40 to 60 min. We show that the former group is dominated by sources formed in globular clusters, while the latter accounts for the (known) UCXB population in the Galactic field. We discuss the possible evolutionary channels for both groups., Comment: Accepted in Astronomy & Astrophysics (see online catalogue here https://research.iac.es/proyecto/compactos/UltraCompCAT/ )

We present simultaneous optical and near-infrared spectra obtained during the 2021 outburst of the black hole transient 4U 1543-47. The X-ray hardness-intensity diagram and the comparison with similar systems reveal a luminous outburst, probably reaching the Eddington luminosity, as well as a long-lasting excursion to the so-called ultra-luminous state. VLT/X-shooter spectra were taken in two epochs 14 days apart during the early and brightest part of the outburst, while the source was in this ultra-luminous accretion state. The data show strong H and HeI emission lines, as well as high-excitation HeII and OIII transitions. Most lines are single-peaked in both spectra, except for the OIII lines that exhibit evident double-peaked profiles during the second epoch. The Balmer lines are embedded in broad absorption wings that we believe are mainly produced by the contribution of the A2V donor to the optical flux, which we estimate to be in the range of 11 to 14 per cent in the $r$ band during our observations. Although no conspicuous outflow features are found, we observe some wind-related line profiles, particularly in the near-infrared. Such lines include broad emission line wings and skewed red profiles, suggesting the presence of a cold (i.e. low ionisation) outflow with similar observational properties to those found in other low-inclination black hole transients., Comment: Accepted for publication in A&A

Accretion disks around compact objects are expected to enter an unstable phase at high luminosity. One instability may occur when the radiation pressure generated by accretion modifies the disk viscosity, resulting in the cyclic depletion and refilling of the inner disk on short timescales. Such a scenario, however, has only been quantitatively verified for a single stellar-mass black hole. Although there are hints of these cycles in a few isolated cases, their apparent absence in the variable emission of most bright accreting neutron stars and black holes has been a lingering puzzle. Here we report the presence of the same multiwavelength instability around an accreting neutron star. Moreover, we show that the variability across the electromagnetic spectrum-from radio to X-ray-of both black holes and neutron stars at high accretion rates can be explained consistently if the accretion disks are unstable, producing relativistic ejections during transitions that deplete or refill the inner disk. Such new association allows us to identify the main physical components responsible for the fast multiwavelength variability of highly accreting compact objects., Comment: Published in Nature. 26 pages, 10 figures. DOI: 10.1038/s41586-022-05648-3

We present a new method to derive binary inclinations in quiescent black hole (BH) X-ray transients (XRTs), based on the depth of the trough (T) from double-peaked Ha emission profiles arising in accretion discs. We find that the inclination angle (i) is linearly correlated with T in phase-averaged spectra with sufficient orbital coverage (>~50 per cent) and spectral resolution, following i (deg)=93.5 x T +23.7. The correlation is caused by a combination of line opacity and local broadening, where a leading (excess broadening) component scales with the de-projected velocity of the outer disc. Interestingly, such scaling allows to estimate the fundamental ratio M1/Porb by simply resolving the intrinsic width of the double-peak profile. We apply the T-i correlation to derive binary inclinations for GRO J0422+32 and Swift J1357-0933, two BH XRTs where strong flickering activity has hindered determining their values through ellipsoidal fits to photometric light curves. Remarkably, the inclination derived for GRO J0422+32 (i=55.6+-4.1) implies a BH mass of 2.7+0.7-0.5 Msun thus placing it within the gap that separates BHs from neutron stars. This result proves that low-mass BHs exist in nature and strongly suggests that the so-called "mass gap" is mainly produced by low number statistics and possibly observational biases. On the other hand, we find that Swift J1357-0933 contains a 10.9+1.7-1.6 Msun BH seen nearly edge on (i=87.4+2.6-5.6 deg). Such extreme inclination, however, should be treated with caution since it relies on extrapolating the T-i correlation beyond i>~75 deg, where it has not yet been tested., Comment: Accepted for publication in MNRAS, 21 pages, 15 figures, 6 Tables; references and typos corrected

MAXI J1348-630 is a low mass X-ray binary discovered in 2019 during a bright outburst. During this event, the system sampled both hard and soft states following the standard evolution. We present multi-epoch optical and near-infrared spectroscopy obtained with X-shooter at the Very Large Telescope. Our dataset includes spectra taken during the brightest phases of the outburst as well as the decay towards quiescence. We study the evolution of the main emission lines, paying special attention to the presence of features commonly associated with accretion disc winds, such as blue-shifted absorptions, broad emission line wings and flat-top profiles. We find broad emission line wings in H-alpha during the hard-to-soft transition and blue-shifted absorption troughs at ~-500 km/s in H-beta, HeI-5876, H-alpha and Pa-beta during the bright soft-intermediate state. In addition, flat-top profiles are seen throughout the outburst. We interpret these observables as signatures of a cold (i.e. optical to infrared) accretion disc wind present in the system. We discuss the properties of the wind and compare them with those seen in other X-ray transients. In particular, the wind velocity that we observe is low when compared to those of other systems, which might be a direct consequence of the relatively low binary inclination, as suggested by several observables. This study strengthen the hypothesis that cold winds are a common feature in low mass X-ray binaries and that they can also be detected in low inclination objects via high-quality optical and infrared spectroscopy., Comment: Accepted for publication in A&A

All disc-accreting astrophysical objects produce powerful outflows. In binaries containing neutron stars (NS) or black holes, accretion often takes place during violent outbursts. The main disc wind signatures during these eruptions are blue-shifted X-ray absorption lines, which are preferentially seen in disc-dominated "soft states". By contrast,optical wind-formed lines have recently been detected in "hard states", when a hot corona dominates the luminosity. The relationship between these signatures is unknown, and no erupting system has revealed wind-formed lines between the X-ray and optical bands yet, despite the many strong resonance transitions in this ultraviolet (UV) region. Here, we show that the transient NS binary Swift J1858.6-0814 exhibits wind-formed, blue-shifted absorption associated with C IV, N V and He II in time-resolved UV spectroscopy during a luminous hard state. This represents the first evidence for a warm, moderately ionized outflow component in this state. Simultaneously observed optical lines also display transient blue-shifted absorption. Decomposing the UV data into constant and variable components, the blue-shifted absorption is associated with the former. This implies that the outflow is not connect to the luminous flares in the data. The joint presence of UV and optical wind features reveals a multi-phase and/or stratified outflow from the outer disc. This type of persistent mass loss across all accretion states has been predicted by radiation-hydrodynamic simulations and helps to explain the shorter-than-expected outbursts duration., Comment: Published in Nature. Submitted: 9 July 2021

In the precision medicine era, there is a growing need for precision radiotherapy where the planned radiation dose needs to be optimally determined by considering a myriad of patient-specific information in order to ensure treatment efficacy. Existing artificial-intelligence (AI) methods can recommend radiation dose prescriptions within the scope of this available information. However, treating physicians may not fully entrust the AI's recommended prescriptions due to known limitations or when the AI recommendation may go beyond physicians' current knowledge. This paper lays out a systematic method to integrate expert human knowledge with AI recommendations for optimizing clinical decision making. Towards this goal, Gaussian process (GP) models are integrated with deep neural networks (DNNs) to quantify the uncertainty of the treatment outcomes given by physicians and AI recommendations, respectively, which are further used as a guideline to educate clinical physicians and improve AI models performance. The proposed method is demonstrated in a comprehensive dataset where patient-specific information and treatment outcomes are prospectively collected during radiotherapy of $67$ non-small cell lung cancer patients and retrospectively analyzed.

We present twelve epochs of optical spectroscopy taken across the discovery outburst of the black hole candidate MAXI J1803-298 with the GTC and VLT telescopes. The source followed a standard outburst evolution with hard and soft states. The system displays a triangular shape in the hardness intensity diagram, consistent with that seen in high inclination black hole transients and the previously reported detection of X-ray dips. The two epochs observed during the initial hard state exhibited asymmetric emission line profiles, including a P-Cygni profile simultaneously detected in H-alpha and He I 6678, which indicates the presence of an optical wind in the system. The remaining spectra, obtained during the transition to the soft state and the subsequent decay, are instead characterized by narrower, double peaked emission lines embedded into broad absorption components. One epoch (intermediate state) also includes near-infrared coverage, revealing complex line profiles in the Paschen and Bracket series, which suggests that the outflow is still present during the outburst decay through the soft state. The growing list of low-mass X-ray binaries with optical and near-infrared outflow signatures indicates that these are common features. Furthermore, the lowest luminosity spectrum exhibits an H-alpha full-width-at-half-maximum of 1570 +- 100 km/s. This, together with previous constraints on the binary parameters, allows us to favor a compact object mass of ~ 3-10 Msun, further supporting its black hole nature., Comment: 10 pages, 5 figures. Accepted for publication in ApJ Letters on the 23rd of January 2022

Black hole low-mass X-ray binaries (BH LMXBs) evolve in a similar way during outburst. Based on the X-ray spectrum and variability, this evolution can be divided into three canonical states: low/hard, intermediate and high/soft state. BH LMXBs evolve from the low/hard to the high/soft state through the intermediate state in some outbursts (here called "full outbursts"). However, in other cases, BH LMXBs undergo outbursts in which the source never reaches the high/soft state, here called "Failed-Transition outburst" (FT outbursts). From a sample of 56 BH LMXBs undergoing 128 outbursts, we find that $\sim$36% of these BH LMXBs experienced at least one FT outburst, and that FT outbursts represent $\sim$33% of the outbursts of the sample, showing that these are common events. We compare all the available X-ray data of full and FT outbursts of BH LMXBs from RXTE/PCA, Swift/BAT and MAXI and find that FT and full outbursts cannot be distinguished from their X-ray light curves, HIDs or X-ray variability during the initial 10-60 days after the outburst onset. This suggests that both types of outbursts are driven by the same physical process. We also compare the optical and infrared (O/IR) data of FT and full outbursts of GX 339-4. We found that this system is generally brighter in O/IR bands before an FT outburst, suggesting that the O/IR flux points to the physical process that later leads to a full or an FT outburst. We discuss our results in the context of models that describe the onset and evolution of outbursts in accreting X-ray binaries., Comment: 30 pages, 20 figures, 2 tables. Accepted for publication in MNRAS main journal

We illustrate the mathematical theory of entropy production in repeated quantum measurement processes developed in a previous work by studying examples of quantum instruments displaying various interesting phenomena and singularities. We emphasize the role of the thermodynamic formalism, and give many examples of quantum instruments whose resulting probability measures on the space of infinite sequences of outcomes (shift space) do not have the (weak) Gibbs property. We also discuss physically relevant examples where the entropy production rate satisfies a large deviation principle but fails to obey the central limit theorem and the fluctuation-dissipation theorem. Throughout the analysis, we explore the connections with other, a priori unrelated topics like functions of Markov chains, hidden Markov models, matrix products and number theory.

We studied the outburst evolution and timing properties of the recently discovered X-ray transient MAXI J1348-630 as observed with NICER. We produced the fundamental diagrams commonly used to trace the spectral evolution, and power density spectra to study the fast X-ray variability. The main outburst evolution of MAXI J1348-630 is similar to that commonly observed in black hole transients. The source evolved from the hard state, through hard- and soft-intermediate states, into the soft state in the outburst rise, and back to the hard state in reverse during the outburst decay. At the end of the outburst, MAXI J1348-630 underwent two reflares with peak fluxes ~1 and ~2 orders of magnitude fainter than the main outburst, respectively. During the reflares, the source remained in the hard state only, without undergoing any state transitions, which is similar to the so-called "failed outbursts". Different types of quasi-periodic oscillations (QPOs) are observed at different phases of the outburst. Based on our spectral-timing results, we conclude that MAXI J1348-630 is a black hole candidate., Comment: 12 pages, 8 figures, 1 table, accepted for publication in MNRAS