We propose a novel approach for determining the orbital inclination of lowmass Xray binary systems by modelling the H and H line profiles emitted by the accretion disc with a Newtonian i.e. nonrelativistic version of DISKLINE. We applied the model to two sample sources Swift J1357.20933 and MAXI J1305704 which are both transient black hole systems and analyse two observations that were collected during a quiescent state and one observation of an outburst. The line profile is well described by the DISKLINE model although we had to add a Gaussian line to describe the deep inner core of the doublepeaked profile which the DISKLINE model was unable to reproduce. The H emission lines in the spectrum of Swift 1357.20933 and the H emission lines in that of MAXI J1305704 during the quiescent state are consistent with a scenario in which these lines originate from a disc ring between 9.6  5710SUP3SUP RSUBgSUB and 1.94  2010SUP4SUP RSUBgSUB respectively. We estimate an inclination angle of 81  5 degrees for Swift J1357.20933 and an angle of 73  4 degrees for MAXI J1305704. This is entirely consistent with the values reported in the literature. In agreement with the recent literature our analysis of the outburst spectrum of MAXI J1305704 revealed that the radius of the emission region deviates from expected values. It is larger than the orbital separation of the system. This outcome implies several potential scenarios including line profile contamination an alternative disc configuration that deviates from the Keplerian model or even the possibility of a circumbinary disc. We caution that these results were derived from a simplistic model that may not fully describe the complicated physics of accretion discs. Despite these limitations our results for the inclination angles are remarkably consistent with recent complementary studies and the proposed description of the emitting region remains entirely plausible.
['accretion', 'accretion disks', 'line: profiles', 'stars: black holes', 'X-rays: binaries', 'X-rays: individuals: Swift J1357.2–0933', 'X-rays: individuals: MAXI J1305–704', 'Astrophysics - High Energy Astrophysical Phenomena']
Xray view of emission lines in optical spectra Spectral analysis of the two lowmass Xray binary systems Swift J1357.20933 and MAXI J1305704
{'No Header': '8  \n9  \n10  \n11  \n12  \n13  \n14  \n15  \n16  \n17  \n18', 'X-ray view of emission lines in optical spectra: Spectral analysis of the two low-mass X-ray binary systems Swift J1357.2-0933 and MAXI J1305-704': 'A. Anitra 1 , C. Miceli 1 , 2 , 3 , T. Di Salvo 1 , R. Iaria 1 , N. Degenaar 6 , Jon M. Miller 7 , F. Barra 1 , W. Leone 1 , 4 , L. Burderi 5  \n- 1 Università degli Studi di Palermo, Dipartimento di Fisica e Chimica, via Archirafi 36, I-90123 Palermo, Italy\n- 2 INAF / IASF Palermo, via Ugo La Malfa 153, I-90146 Palermo, Italy\n- 3 IRAP, Universitè de Toulouse, CNRS, UPS, CNES, 9, avenue du Colonel Roche BP 44346 F-31028 Toulouse, Cedex 4,France\n- 4 Department of Physics, University of Trento, Via Sommarive 14, 38122 Povo (TN), Italy\n- 5 Dipartimento di Fisica, Università degli Studi di Cagliari, SP Monserrato-Sestu, KM 0.7, Monserrato, 09042 Italy\n- 6 Anton Pannekoek Institute for Astronomy, University of Amsterdam, Postbus 94249, 1090 GE Amsterdam, The Netherlands\n- 7 Department of Astronomy, The University of Michigan, 1085 South University Avenue, Ann Arbor, MI, 48109, USA  \nSeptember 19, 2024', 'ABSTRACT': 'We propose a novel approach for determining the orbital inclination of low-mass X-ray binary systems by modelling the H α and H β line profiles emitted by the accretion disc, with a Newtonian (i.e. non-relativistic) version of diskline . We applied the model to two sample sources, Swift J1357.2-0933 and MAXI J1305-704, which are both transient black hole systems, and analyse two observations that were collected during a quiescent state and one observation of an outburst. The line profile is well described by the diskline model, although we had to add a Gaussian line to describe the deep inner core of the double-peaked profile, which the diskline model was unable to reproduce. The H β emission lines in the spectrum of Swift J1357.2-0933 and the H α emission lines in that of MAXI J1305-704 during the quiescent state are consistent with a scenario in which these lines originate from a disc ring between (9 . 6 -57) × 10 3 , Rg and (1 . 94 -20) × 10 4 , Rg, respectively.  \nWe estimate an inclination angle of 81 ± 5 degrees for Swift J1357.2-0933 and an angle of 73 ± 4 degrees for MAXI J1305-704. This is entirely consistent with the values reported in the literature.  \nIn agreement with the recent literature, our analysis of the outburst spectrum of MAXI J1305-704 revealed that the radius of the emission region deviates from expected values. It is larger than the orbital separation of the system. This outcome implies several potential scenarios, including line profile contamination, an alternative disc configuration that deviates from the Keplerian model, or even the possibility of a circumbinary disc.  \nWe caution that these results were derived from a simplistic model that may not fully describe the complicated physics of accretion discs. Despite these limitations, our results for the inclination angles are remarkably consistent with recent complementary studies, and the proposed description of the emitting region remains entirely plausible.  \nKey words. accretion, accretion discs - stars: Black hole - stars: individual: Swift J1357.2- stars: individual: MAXI J1305-704 X-rays: binaries - X-rays: general - eclipses', '1. Introduction': 'Low-mass X-ray binaries (LMXBs) are systems that are composed of a compact object, such as a black hole (BH) or neutron star (NS), and a donor star with a mass lower than 1 M ⊙ . The compact object in these systems accretes matter from its companion via Roche-lobe overflow, which leads to the formation of an accretion disc around it.  \nThe study of the double-peaked emission line profiles that are often present in the optical spectra of these sources is a powerful tool that allows us to determine the orbital parameters of the system. By examining parameters such as the orbital period and radial velocity semi-amplitude we can dynamically confirm the nature of the compact objects in these systems.  \nIt has been acknowledged that the central core within the double-peaked lines that are emitted by an accretion disc become deeper with increasing inclination. This is particularly noticeable for inclinations greater than 67 degrees (Horne & Marsh 1986). Casares et al. (2022) discovered a direct correlation between the  \ndepth of the inner core in the double-peaked H α emission line and the inclination angle in a group of quiescent BH LMXBs. By applying this method, the authors were able to determine the inclination angle of di ff erent sources with a high confidence level.  \nWang et al. (2009) examined the optical spectrum of the millisecond radio pulsar binary SDSS J102347.6 + 003841. They focused on the analysis of double-peaked emission lines to determine the accretion disc geometry. By adopting two Lorentz functions to fit the lines and calculating the flux, they were able to estimate key properties such as the disc temperature range, the outer and inner emission radii, and the overall mass of the disc.  \nIn this work, we investigate whether an estimate of the system inclination angle can be derived by applying the diskline model for the first time to fit the double-peaked emission lines related to the H-series transition that appears in the optical spectra of two candidate transient BH systems. The diskline model can describe simple pure Keplerian fully symmetric disc profiles, where the emissivity is a pure power law. This may not be the case for optical lines that are emitted in the outer disc,  \n19  \n20  \n21  \n22  \n23  \n24  \n25  \n26  \n27  \n28  \n29  \n30  \n31  \n32  \n33  \n34  \n35  \n36  \n37  \n38  \n39  \n40  \n41  \n42  \n43  \n44  \n45  \n46  \n47  \n48  \n49  \n50  \n51  \n52  \n53  \n54  \n55  \n56  \n57  \n58  \n59  \n60  \n61  \n62  \n63  \n64  \n65  \n66  \n67  \n68  \n69  \n70  \n71  \n72  \n73  \n74  \n75  \n76  \n77  \n78  \n79  \n80  \n81  \n82  \n83  \n84  \n85  \n86  \n87  \n88  \n89  \n90  \n91  \n92  \n93  \nwhich is most a ff ected by asymmetries, non-uniform emissivity patterns, strong outer disc features such as hot spots and bulges, departures from Keplerian flow due to tidal interactions, and so on. However, we show that the line profiles observed during Xray quiescence of the BH candidates Swift J1357.2-0933 and Maxi J1305-704 appear to be simple enough to provide useful constraints on their inclination angle.', '2. X-ray transient system samples': 'To ensure the accuracy of the inclination measurements, we chose two transient BH X-ray binaries that have been intensely studied in the literature and for which several indications of a high-inclination angle have been reported.', '2.1. Swift J1357.2-0933': 'Swift J1357.2-0933 (hereafter J1357) is an LMXB X-ray transient that was originally detected during its outburst in 2011. The distance to J1357 is estimated to be greater than 2.29 kpc, which places the source within the thick Galactic disc (Mata Sánchez et al. 2015). Studies of the modulation of the H α double-peak emission line during 14 months of outbursts have led to an orbital period of 2 . 5673 ± 0 . 0006 h (Casares et al. 2022). An estimate of the compact object mass was provided by Mata Sánchez et al. (2015) based on the correlation between the full width at half maximum (FWHM) and the radial velocity of the donor star (K2) (Casares 2015). This revealed that this source is one of the most massive BHs in our galaxy (MBH > 9 . 3 M ⊙ ). The optical light curve of the source displays periodic dips with a consistent pattern during the outburst (Armas Padilla et al. 2014). These dips are thought to be linked to a toroidal structure within the inner accretion disc that gradually moves outward as the outburst proceeds (Corral-Santana et al. 2013).  \nThe frequent dipping episodes in the system suggest a highinclination angle i larger than 70 degrees (see Corral-Santana et al. 2013; Torres et al. 2015). Anitra et al. (2023) performed phase-resolved spectroscopy during a quiescent state. They focused on the analysis of the H β line profile. The authors estimated the systemic velocity and the radial velocity semiamplitude of the black hole and detected compelling evidence of a narrow and variable inner core within the double-peaked H β emission line profile. These features were observed in highinclination cataclysmic variables (Schoembs & Hartmann 1983) and are probably associated with an occultation of the inner-disc emission by the outer rim bulge.  \nHigh-inclination systems typically exhibit eclipses in their light curves (Anitra et al. 2021). The absence of eclipses in this system appears to contradict a high inclination (see Iaria et al. 2018; Armas Padilla et al. 2014). Corral-Santana et al. (2013) justified this absence by considering an inclination angle of 78 degrees and a low mass ratio, proposing that the radius of the donor star might be comparable to or even smaller than the outer disc rim. Based on these indicators of a high-inclination nature, J1357 is an ideal candidate for our analysis.', '2.2. MAXI J1305-704': 'MAXI J1305-704 (hereafter J1305) was proposed to be a highinclination (BH) X-ray binary because a dipping behaviour was detected during its outburst (Suwa et al. 2012, Shidatsu et al. 2013, Morihana et al. 2013, Kennea et al. 2012). The source has  \nan estimated distance of d = 7 . 5 l -1 . 4 + 1 . 8 kpc, which places it in the thick Galactic disc (Mata Sánchez et al. 2021).  \nJ1305 was studied in detail in the X-rays, but follow-up at other wavelengths has been scarce (see Shaw et al. 2017; Miller et al. 2014; Shidatsu et al. 2013). The first optical spectroscopic analysis of this source during its 2012 outburst has been performed by Miceli et al. (2024). It focused on the double-peaked H α emission line, which revealed no conclusive evidence of the outflow features within the system.  \nMata Sánchez et al. (2021) inferred an orbital period of 9 . 456 ± 0 . 096 h and a MBH = 8 . 9 + 1 . 6 -1 . 0 M ⊙ by analysing the orbital modulation in the optical light curve in quiescence. The same authors obtained a constraint on the orbital inclination of i = 72 + 5 -8 deg, which supports the high-inclination scenario. The number of observations during outburst and quiescence and the well-constrained values of the inclination angle and other orbital parameters (e.g. orbital period and mass of the two objects) make J1305 one of the best candidates for our analysis.', '3. Observations': 'We analysed di ff erent sets of observations for the two sample sources.  \nJ1357 was observed by the 10.4-meter Gran Telescopio Canarias (GTC) at the Observatorio del Roque de los Muchachos (ORM) on the island of La Palma (Spain), using the Optical System for Imaging and low-Intermediate-Resolution Integrated Spectroscopy (OSIRIS). Eight observations were collected during a quiescence state. They were focused on the wavelength range from 3950 to 5700 Å. The R2000B grism was set alongside a 0.8" slit during the campaign. This configuration allowed us to obtain a spectral resolution of R = 1903 and a dispersion of D = 0.86 Å / pix. The dispersion is evaluated at λ c = 4755 Å , and the resolution was obtained by the FWHM of the skylines in the background spectra. Data were collected on 5 March 2016 between 03:32:11.4 and 06:19:5. Each observation was 1235 s long for a total exposure of 2.75 h and covered a full orbit (Anitra et al. 2023).  \nJ1305 was observed both during an outburst and in a quiescent state. The initial observation took place during its discovery outburst in 2012 and was made with the 6.5-meter Magellan Clay telescope located at Las Campanas Observatory in Chile. This observation involved the use of the Low Dispersion Survey Spectrograph (LDSS-3) and the VPH-ALL grism setup. The observations covered two consecutive nights, 2 May (02:18:34 UTC to 02:56:41 UTC) and 3 May (00:29:36 UTC to 01:09:03 UTC) with a 0.75\' slit. The exposure time for each individual spectrum was set at 300 s for a total of six spectra per night (Miceli et al. 2024). The setup configuration allowed us to focus on a wavelength range between 4500 - 6950 Åwith a resolution power of R = 826.  \nThe source was observed during quiescence on 31 March 2016 with the Very Large Telescope Unit Telescope 1 (VLTUT1; Paranal Observatory, Chile) using the focal reducer / lowdispersion spectrograph 2 (FORS2, Appenzeller et al. 1998) in long-slit mode. The data sets were composed of 16 spectra, with an exposure time of 1800 s each. They were collected consecutively for a total exposure of ∼ 9 h, that is, they covered almost one orbital period (Mata Sánchez et al. 2021). The spectra cover a wavelength range between 5800 - 7300 Å with a spectral resolution of R ∼ 2140 and a dispersion of 0.76 Å / pix. This was obtained by measuring the FWHM of the skylines in the background spectra.  \n94  \n95  \n96  \n97  \n98  \n99  \n100  \n101  \n102  \n103  \n104  \n105  \n106  \n107  \n108  \n109  \n110  \n111  \n112  \n113  \n114  \n115  \n116  \n117  \n118  \n119  \n120  \n121  \n122  \n123  \n124  \n125  \n126  \n127  \n128  \n129  \n130  \n131  \n132  \n133  \n134  \n135  \n136  \n137  \n138  \n139  \n140  \n141  \n142  \n143  \n144  \n145  \n146  \n147  \n148  \n149  \n150  \n151  \n152  \n153  \n154  \n155  \n156  \n157  \n158  \n159  \n160  \n161  \n162  \n163  \n164  \n165  \n166  \n167  \n168  \n169  \n170  \n171  \n172  \n173  \n174  \n175  \n176  \n177  \n178  \n179  \n180  \n181  \n182  \n183  \n184  \n185  \n186  \n187  \n188  \n189  \n190  \n191  \n192  \n193  \n194  \n195  \n196  \n197  \n198  \n199  \n200  \n201  \n202  \n203  \n204  \n205  \n206  \n207  \n208  \n209  \n210  \n211  \nWe reduced the data using standard procedures based on the iraf 1 software, molly tasks, and python packages from astropy and pyastronomy (Astropy Collaboration et al. 2022). These tools allowed us to correct the observed spectra for bias and flats and to calibrate the data set. We note that cosmic rays contaminated the spectra of J1357, and we corrected for them using the L.A. Cosmic task (van Dokkum et al. 2012). We adopted the optimal extraction technique (Naylor 1998) to extract the twodimensional images. The quiescence spectra of J1305 might be contaminated by the companion emission. However, as claimed by Mata Sánchez et al. (2021), the contamination is probably lower than 10% of the total flux and is therefore negligible.', '4.1. Diskline model': "The emission lines of the Balmer series are associated with the atomic hydrogen transition. They appear as broad symmetrical double-horned lines when observed from accretion discs at highinclination angles. In these systems, the emission is mainly influenced by the Doppler shift, which is caused by the orbital motion of matter in the accretion disc around the compact object.  \nIn a disc ring at a certain radius RD, matter moves at the corresponding Keplerian velocity Vkep. For a distant observer, some of the matter within the disc moves toward the observer, while it simultaneously moves away from the observer on the far side of the disc. This dual motion causes the observer to perceive the emission line as simultaneously blue- and red-shifted, and it causes the distinctive symmetrical double-horn profile. Because the emission comes from di ff erent regions of the disc surface, the overall line profile is additionally broadened by the velocity distribution in the disc. As first discussed by Smak (1981) and later by Horne & Marsh (1986), the final line profile is significantly a ff ected by the inclination angle of the system relative to the line of sight because the component of the Keplerian velocity Vkep along the line of sight is V D = ± V Kep ( R D) sin i . Furthermore, the separation of the peaks is determined by the outer radius of the emitting region, but the wings are shaped by the emissivity profile.  \nThe X-ray spectrum-fitting software xspec (Arnaud 1996) provides several models for describing an emission line that is shaped by Doppler e ff ects from a Keplerian accretion disc. One of the most commonly adopted models for fitting double-horn emission lines is the diskline model (Fabian et al. 1989). It is usually adopted to describe the iron K α emission line that characterises the X-ray reflection features observed in LMXBs (see e.g. di Salvo et al. 2009).  \nThe shape of the broad iron fluorescence line observed at X-ray wavelengths, similar to the hydrogen Balmer series lines in optical spectra, is a ff ected by Newtonian Doppler shift. Its profile appears to be quite di ff erent, however. The Fe K α line originates in the innermost part of the disc, where the velocity of matter reaches relativistic values. Consequently, the relativistic beaming a ff ects the emission and intensifies the blue peak of the line while weakening the red peak. Furthermore, the line is gravitationally redshifted by the strong gravitational pull near the compact object, causing it to shift to lower energies in a way that depends on the distance from the compact object. As a result, the emission line appears to be broadened and asymmetric.  \nThe diskline model is composed of the following free parameters that shape the emission line profile: the inner (Rin) and outer radius (Rout) of the emission region in the disc (expressed in units of gravitational radii, Rg = GM / c 2 ), the orbital inclination of the system, and the index, β , which describes the dependence of the disc emissivity on the distance from the compact object (E ∝ R β ). This model was used several times to obtain constraints on the inclination angle of LMXBs (see e.g. Anitra et al. 2021; Iaria et al. 2009; Cackett et al. 2008) when relativistic e ff ects were dominant. Although this model is usually adopted to describe a (relativistic) emission line, it can be adapted to model the hydrogen line profiles as well. When the value of the inner and outer radius is high enough, meaning that the emission region of the line is far from the compact object, the relativistic beaming and the gravitational redshift indeed become negligible, and only the Newtonian Doppler shift dominates (see Fig. 1). Under these conditions, the model is equivalent to the optically thin solution presented in Horne & Marsh (1986). The outer radius of the disc is also most strongly a ff ected by asymmetries, non-uniform emission line emissivity patterns, and strong outer disc features. For instance, self-absorption by an optically thick disc atmosphere can contribute to deepening the inner core of the line. Simultaneously, disc precession may play a significant role in determining the orbit-averaged velocity shifts that are observed in the emission line centroid (Torres et al. 2002), while a hot spot can enhance one of the two peaks of the line (Marsh et al. 1994). For these reasons, it is important to exercise caution with respect to the constraints derived from the analysis.  \nIn order to isolate the emission line profile on which we focused our analysis, we normalised the spectra of J1357 and J1305 by dividing them by their polynomial best-fit functions. Specifically, we applied a third-order polynomial fit for the first source and a fifth-order polynomial fit for the second source. We thus adapted our optical data sets for an analysis in the xspec environment by creating a unit diagonal response matrix of the data size. We adhered to the following naming convention: spectra linked to J1357 are denoted src1-Q, and those associated with J1305 are labelled src2-O during outburst and src2-Q during quiescence.  \nWe find it convenient to express the constrained best-fit radii in units of the expected tidal radius of the system, that is, the outer edge of the disc truncated by the tidal torque of the companion star (Frank et al. 2002). The tidal radius depends on several binary parameters, but its value can be shown to be close to RT = 0 . 9 R1 (Frank et al. 2002), where R1 is the Roche-lobe radius of the compact object. This can be calculated using the Eggleton (1983) equation  \nR 1 a = 0 . 49 q -2 / 3 0 . 6 q -2 / 3 + ln (1 + q -1 / 3 ) , (1)  \nwhere q is the mass ratio of the two bodies M1 and M2, and a is the orbital separation of the system, which can be calculated using Kepler's third law. Moreover, adopting black hole masses of about 8.9 solar masses for J1305 (Mata Sánchez et al. 2021) and 9.3 solar masses for J1357 (Mata Sánchez et al. 2015), we can define the gravitational radius Rg for the two systems as 13 km and 14 km, respectively. Consequently, we obtain an orbital separation and tidal radius of a = (1 . 0259 ± 0 . 0002) × 10 6 Rg and RT = (5 . 9 ± 0 . 004) × 10 4 Rg for J1357, and a = (2 . 524 ± 0 . 002) × 10 5 Rg, and RT = (1 . 45 ± 0 . 06) × 10 5 Rg for J1305.  \n212  \n213  \n214  \n215  \n216  \n217  \n218  \n219  \n220  \n221  \n222  \n223  \n224  \n225  \n226  \n227  \n228  \n229  \n230  \n231  \n232  \n233  \n234  \n235  \n236  \n237  \n238  \n239  \n240  \n241  \n242  \n243  \n244  \n245  \n246  \n247  \n248  \n249  \n250  \n251  \n252  \n253  \n254  \n255  \n256  \n257  \n258  \n259  \n260  \n261  \n262  \n263  \n264  \n265  \n266  \n267  \n268  \n269  \n270  \n271  \n272  \n273  \n274  \n275  \n276  \n277  \n278  \n279  \n280  \n281  \n282  \n283  \n284  \n285  \n286  \n287  \n288  \n289  \n290  \n291  \n292  \n293  \n294  \n295  \n296  \nFig. 1. Series of examples showing the diskline model profiles with fixed inclinations ( i = 83 degrees), emissivity index ( β = 2.5), and energy centroid (E = 4857.522 Å) for all models. The inner and outer radius were systematically varied within ranges of (100 -10 4 ) Rg and (10 4 -6 × 10 4 ) Rg, respectively, to highlight the variations in the line profiles. The model normalisations are scaled for visual clarity.  \n<!-- image -->", '4.2. J1357 spectroscopy': 'The wavelength coverage of our spectra allowed us to focus on the H β line region with a su ffi ciently high signal-to-noise ratio (S / N) 2 .  \nWe fitted the eight spectra simultaneously using a power-law model to account for the normalised continuum (the associated parameters are not physically relevant), to which we added the diskline model to fit the H β line profile. The observations were obtained within a time interval of 2.75 h (one orbital period), and we therefore expect very little spectral variation in the di ff erent spectra. As a result, we constrained specific parameters, such as the inclination angle, emissivity index, and inner and outer radii, to be the same for all the spectra. The model, called model 1, achieved a poor fit to the data, with a χ 2 / d . o . f . = 4089.1 / 2811. Nonetheless, the values for the inner and outer radii of the line emission region within the disc, (9 . 6 ± 0 . 2) × 10 3 Rg ( ∼ 0 . 16 RT ) and 5 . 4 + 0 . 1 -0 . 2 × 10 4 Rg ( ∼ 0 . 92 RT ), respectively, agree with expectations. Additionally, the emissivity index characterises a region of the disc that is located farther from the compact object. However, as shown in Fig. 2 (second panel), residuals lie in the region between the two peaks.  \nTo improve the fit without modifying the diskline model, we attempted two di ff erent approaches. First, we masked out the core of the line and excluded these data from the fit. The resulting best-fit parameters are consistent with those presented previously, with a χ 2 / dof = 4071 . 72 / 2805. Second, we introduced a Gaussian absorption line for which the centroid and FWHM parameters were linked for all the spectra, but we al-  \nalized flux  \nNor  \nFig. 2. J1357 spectra collected during a quiescent phase and residuals in units of sigma with respect to the diskline model described in the test (top panel). The middle and bottom panel shows the residuals obtained by using only the diskline model ( model 1) and those obtained with the same model plus a Gaussian absorption line ( model 2), respectively.  \n<!-- image -->  \nlowed an independent variable depth to determine whether each spectrum significantly required this feature (see Sect. 5). This model, called model 2, achieved a better fit to the data with a χ 2 / dof = 3849 . 0 / 2794. We tested the improvement of the fit using the statistical test Ftest and obtained a probability of a chance improvement of ∼ 1 . 1 × 10 -28 . This means that including the Gaussian line in absorption improves the quality of the fit with a confidence level (c.l.) higher than 7 σ . The new model provides a more solid constraint on the inclination angle (i = 83 + 5 -3 degrees), and the H β emission line is consistent with being emitted by a ring in the disc with a radius between 9 . 6 + 0 . 2 -0 . 1 × 10 3 Rg ( ∼ 0 . 16 RT ) and 5 . 7 + 0 . 1 -0 . 2 × 10 4 Rg ( ∼ 0 . 91 RT ). However, for this high number of degrees of freedom, a reduced χ 2 of 1.38 is unacceptable because it implies a very low null-hypothesis probability. In order to obtain reliable error bars for the best-fit spectral parameters, we therefore scaled the errors by the q χ 2 red to achieve a reduced χ 2 of approximately 1 and repeated the fit. Using this approach, we obtained the same values for the inner and outer radii of the emitting region, while the fit was able to provide just a lower limit for the inclination angle of 80 degrees. Notably, this finding only holds true for errors calculated at the 90% c.l. Conversely, when errors are calculated at the 68% c.l, the inclination angle is constrained to be 83 . 8 + 5 . 4 -1 . 8 degrees (see Table A.1).', '4.3. J1305 spectroscopy': 'The wavelength coverage and spectral resolution provided by the telescope set-up allowed us to perform a high-resolution analysis focused on the H α emission line. In the spectra collected during the outburst, we identify H α and H β emission lines, corresponding to the Balmer series and He II 4686 Å as the Bowen blend (a mixture of N iii and C iii emission lines; see e.g. Steeghs  \n297  \n298  \n299  \n300  \n301  \n302  \n303  \n304  \n305  \n306  \n307  \n308  \n309  \n310  \n311  \n312  \n313  \n314  \n315  \n316  \n317  \n318  \n319  \n320  \n321  \n322  \n323  \n324  \n325  \n326  \n327  \n328  \n329  \n330  \n331  \n332  \n333  \n334  \n335  \n336  \n337  \n338  \n339  \n340  \n341  \n342  \n343  \n344  \n345  \n346  \n347  \n348  \n349  \n350  \n351  \n352  \n353  \n354  \n355  \n356  \n357  \n358  \n359  \n360  \n361  \n362  \n363  \n364  \n365  \n366  \n367  \n368  \n369  \n370  \n371  \nTable 1. Ratios of the blue and red peaks of the J1305 spectra collected during quiescence.  \n& Casares 2002). We decided to focus our analysis on the H α line, which is present in all the spectra in the two days and is the strongest isolated feature (Miceli et al. 2024). We separately analysed the two days of optical observations. For the first day, we used four of the six spectra excluding, in particular, src2O1 because of its strong absorption component and src2-O3 because it does not show a double-peak profile, but a flat top profile (which recalls outflow features; Miceli et al. (2024), see also Cúneo et al. 2020).  \nWe applied model 1 to the data and linked the parameters of each spectrum as described in the previous section (see Fig. 3). On the second day of observations, a distinct absorption component is evident in all spectra. It is located at wavelengths longer than those associated with the H α line (about 6563 Å) (see Fig. 4). To take this feature into account, we added to the previous model a Gaussian absorption line for which the centroid and FWHM parameters we linked for all the spectra, but we allowed an independently variable depth. For both days, the adopted models achieve a reasonably good fit to the data, with a χ 2 / dof = 139 . 4 / 153 (first day) and χ 2 / dof = 149 . 0 / 222 (second day). The best-fit parameters are reported in Table A.1. On both days, the inclination angle is 70 ± 4 and 71 ± 4 degrees, which is consistent with the angle reported in the literature ( i = 72 + 5 -8 degrees; Mata Sánchez et al. 2021). Furthermore, for the first day, we derived inner and outer radii values of 1 . 20 + 0 . 09 -0 . 07 × 10 5 Rg ( ∼ 0 . 83 RT ) and 1 . 96 + 0 . 32 -0 . 27 × 10 6 Rg ( ∼ 13 RT ), respectively. For the second day, these values were 1 . 66 + 0 . 10 -0 . 08 × 10 5 Rg ( ∼ 1 . 1 RT ) and 2 . 73 + 0 . 9 -0 . 6 × 10 6 Rg ( ∼ 19 RT ) . A more detailed discussion of these parameters is presented in the following section.', '4.3.1. Quiescence': 'We chose 10 of the 16 spectra for the quiescent phase that showed a more regular line profiles. Our preference lay with line profiles with a symmetrical double peak because the asymmetry in the intensity of the two peaks can be caused by various effects, for instance a hot spot (Shafter et al. 1986), and this might complicate the fitting. However, the exclusion of the 6 spectra does not compromise the credibility of the model. As a test, we individually analysed the excluded spectra and obtained parameter values that were entirely consistent with those reported in the main analysis (see Sect. 4.4). The symmetry of a profile can be estimated by examining the ratios of the two peaks, as reported in Table 1. Following the same procedure as applied to J1357, we first fitted the spectra applying model 1. However,  \nFig. 3. J1305 spectra collected during an outburst phase and residuals in units of sigma with respect to the diskline ( model 1), related to the first day of observation.  \n<!-- image -->  \nmalized flux  \nNor  \nFig. 4. J1305 spectra collected during an outburst phase and residuals in units of sigma with respect to the diskline model plus an absorption Gaussian line, related to the second day of observation.  \n<!-- image -->  \nthe latter did not reach a good fit to the data, with a χ 2 / dof of 6432.32 / 2650, because residuals were clearly present in the core region (see the first panel of Fig. 5). Therefore, we applied model 2 and linked the parameters of di ff erent spectra as described for the J1357 spectra. The model achieved a better fit to the data, ensuring a χ 2 / dof of 4161.1 / 2620. The χ 2 value had to be rejected here as well. Therefore, as discussed in the previous section, we multiplied the errors by the square root of the reduced χ 2 and repeated the fitting. We obtain a constraint on the inclination angle of 72 . 6 + 1 . 4 -1 . 3 degrees at 90% c.l., as well as best-fit values for the inner and outer radii of the emitting region associated with the H α emission line of 1 . 94 + 0 . 07 -0 . 09 × 10 4 Rg ( ∼ 0 . 1 RT ) and 2 . 01 ± 0 . 04 × 10 5 Rg ( ∼ 1 . 3 RT ), respectively.  \n372  \n373  \n374  \n375  \n376  \n377  \n378  \n379  \n380  \n381  \n382  \n383  \n384  \n385  \n386  \n387  \n388  \n389  \n390  \n391  \n392  \n393  \n394  \n395  \n396  \n397  \n398  \n399  \n400  \n401  \n402  \n403  \n404  \n405  \n406  \n407  \n408  \n409  \n410  \n411  \n412  \n413  \n414  \n415  \n416  \n417  \nFig. 5. J1305 spectra collected during a quiescent phase and residuals in units of sigma with respect to the diskline model described in the text (top panel). The middle and bottom panels show the residuals obtained with the diskline model alone ( model 1) and those obtained with the same model plus a Gaussian absorption line ( model 2).  \n<!-- image -->', '4.4. Analysis of the parameter degeneration': 'The line profiles can easily be a ff ected by external factors, such as hot spots or disc precession (e.g.Orosz et al. (1994); Mason et al. (2000); Marsh et al. (1994)). These potential pitfalls for a model might prompt concerns about the reliability of the approach. We therefore explore degeneracies and systematics quantitatively here to ensure that the previously derived constraints are statistically reliable.  \nWhen multiple parameter combinations yield similar or indistinguishable results during the model-fitting process, this indicates that the best-fit values of the individual parameters cannot be determined based on the data alone. To investigate this degeneracy, we employed the steppar command in xspec to perform a fit while stepping the value of the parameters through a given range. This approach allowed us to evaluate all possible pairings of the parameters of interest, measuring how the χ 2 value changes in relation to the best fit for each combination. This allowed us to identify any value combinations that might result in an equivalent spectral fit. We performed the steppar on data sets collected during quiescence from the two sources and focussed on parameters such as inclination, outer radius, and inner radius, which significantly a ff ect the line profile. In Fig. 6, we present contour plots obtained from the analysis.  \nThe contours of the regions are clearly defined and concentric around the best-fit value, indicating robustness in the parameter estimates. In addition, the contours lack extended overlapping or elongated areas that stretch along the diagonals, which would typically indicate degeneration. It is important to observe that even in regions with a diagonal trend (panels a and d), the contours delineate parameter values that lie within the error limits reported in the text. This shows that the same χ 2 values cannot be achieved with parameter pairs that deviate from those established by the best fits.  \nTable 2. Inclination values obtained by separately fitting the spectra related to the two sources in quiescence.  \nBecause the physics within an accretion disc is complex and the approach of the diskline model is simplistic, the high precision of the obtained inclination constraints (especially from the quiescence spectrum of J1305) may cast doubt on the reliability of the results. An insightful way to re-evaluating and discussing the error associated with the best-fit parameters, especially the inclination angle, involves fitting each spectrum individually and considering the dispersion of the resulting values.  \nWe therefore fitted the spectra acquired during quiescence for J1357 and J1305 separately, allowing all parameters to vary freely 3 . In Table 2 we list the inclination values for each fit. Moreover, we decided to also include the spectra that were rejected in the analysis presented in the text to avoid any bias in the estimation of the inclination angle error. We ensured through a steppar analysis throughout the entire range of possible inclination values that the fit did not yield a value corresponding to a local minimum of the χ 2 . Subsequently, we determined the best fit as the average of the obtained values, with the associated error calculated as the semi-dispersion because it represents the maximum dispersion around the average value. This yielded an inclination of 81 ± 5 degrees for J1357 and 73 ± 4 degrees for J1305, and these inclinations agree with the values derived from the analysis. Moreover, the error bars obtained with this method are consistently similar to those constrained by the diskline model.', '5. Discussion': 'We analysed GTC observations of J1357 collected during quiescence, along with VLT and Magellan observations of J1305, taken during quiescence and outburst phases, respectively. Our study focused on the double-peaked profiles of the H β and H α emission lines. With the aim of gaining insights into the system  \n418  \n419  \n420  \n421  \n422  \n423  \n424  \n425  \n426  \n427  \n428  \n429  \n430  \n431  \n432  \n433  \n434  \n435  \n436  \n437  \n438  \n439  \n440  \n441  \n442  \n443  \n444  \n445  \n446  \n447  \n448  \n449  \n450  \n451  \n452  \n453  \n454  \n455  \n456  \n457  \n458  \n459  \nFig. 6. Contour plots for inclination and the outer and inner radii of the spectra related to J1357 (first column) and J1305 (second column). The contours represent the 1 σ , 2 σ , and 3 σ confidence levels, and the cross marks the best-fit values obtained from the best fit.  \n<!-- image -->  \ngeometry, we employed the diskline model, which is complementary to traditional analysis approaches. As discussed in the previous section, the model including the diskline component along with a (Gaussian) absorption line provides a good fit to the data of both systems. This significantly improved the quality of the fit with respect to the model without the absorption line.  \nThe addition of this feature may be necessary due to the potential limitations in the assumptions made by the diskline model. It describes a Gaussian emission line profile, modified  \nby the Doppler e ff ect induced by a Keplerian velocity distribution, in an accretion disc. It is clear that the model does not take self-absorption or other features due to the optically thick nature of the disc into account, which was also noted by Orosz et al. (1994). This e ff ect is particularly strong for high-inclination angles, as appears to be the case of J1357 and J1305. This agrees with our results, for which an additional Gaussian absorption enabled a better description of the observed spectra.  \n460  \n461  \n462  \n463  \n464  \n465  \n466  \n467  \n468  \n469  \n470  \n471  \n472  \n473  \n474  \n475  \n476  \n477  \n478  \n479  \n480  \n481  \n482  \n483  \n484  \n485  \n486  \n487  \n488  \n489  \n490  \n491  \n492  \n493  \n494  \n495  \n496  \n497  \n498  \n499  \n500  \n501  \n502  \n503  \n504  \n505  \n506  \n507  \n508  \n509  \n510  \n511  \n512  \n513  \n514  \n515  \n516  \n517  \n518  \n519  \n520  \n521  \n522  \n523  \n524  \n525  \n526  \n527  \n528  \nUpon visual examination of the individual spectra, it becomes clear that the central core of the line displays narrow and variable absorption that changes in depth and even reaches (sometimes going below) the continuum average level. Features like this were mainly detected in cataclysmic variables that were observed at a high-inclination angle (with i > 75 degrees, Schoembs & Hartmann 1983), and they are thought to be caused by occultation of the inner regions of the accretion disc. Di ff erently from the e ff ects caused by self-absorption, narrow cores are variable with the orbital phase, and they can reach much deeper than the optically thick absorption e ff ects. They can reach even below the normalised continuum (Rayne & Whelan 1981). Narrow variable cores in the J1357 spectra were first claimed by Mata Sánchez et al. (2015), who noted a variation in the normalised flux of the H α core. This was later confirmed by Anitra et al. (2023), who analysed high-resolution data focused on the H β emission line. Mata Sánchez et al. (2021) performed dynamical studies of the same VLT data of J1305 as we presented here. They reported narrow cores in the spectra.  \nAs further proof of this, we calculated the significance of the absorption feature in each spectrum in units of sigma by comparing the normalisation of the line with its uncertainty at the 68% c.l. We found that the intensity of the line in J1357 src1Q2, src1-Q3, and src1-Q6 is ≤ 3 σ , meaning that during these orbital phases, the absorption core of the line is negligible, but it is stronger in src1-Q1, src1-Q7 and src1-Q8. In the quiescent spectra of J1305, all the absorption lines are statistically significant at more than 7 σ , except for the line associated with src2-Q15, where this feature is not necessary. This result is very similar to the result reported by Marsh et al. (1987) in the dwarf nova system Z Cha, which exhibited a core depth variability along the orbit in every line of the Balmer series in the spectrum.  \nAs reported in Sect. 4.3, spectra collected during the outburst do not show narrow variable cores. However, the spectra collected on the second day of observation display absorption features at higher wavelengths compared to the H β line, which may resemble an inverted P-Cygni profile, which is associated with inflows (Cúneo et al. 2020). Miceli et al. (2024) have extensively discussed the nature of this phenomenon and proposed that a broad and variable absorption component is observed in all spectra, which can a ff ect the shape of the line. Nevertheless, the origin of this feature remains the subject of ongoing debate.', '5.1.1. Inclination angle': 'The best-fit parameters of both sources provide a description of the geometry of the system in line with the hypothesis of a high inclination. We obtain an angle of 72 . 6 + 1 . 4 -1 . 3 degrees (73 ± 4 degrees considering the semi-dispersion around the average value) for J1305 during quiescence and angles of 70 ± 4 and 71 ± 4 degrees for the two days of outburst, respectively. These results are consistent with the previous constraint reported by Mata Sánchez et al. (2021), i = 72 + 5 -8 degrees.  \nRecently, Casares et al. (2022) studied the inner core depth in J1357 spectra and derived an inclination angle of 87 . 4 + 2 . 6 -5 . 6 degrees, which matches the estimate we find by applying the diskline model (i > 80 degrees at 90% c.l., i = 80 + 2 -6 degrees at 68 %c.l. and 81 ± 5 considering the average values). However, the lack of eclipses in both the X-ray and optical light curves of this source (Corral-Santana et al. 2013) contradicts the expected behaviour of an edge-on configuration. Corral-Santana et al. (2013) provided an explanation for this phenomenon, which may be due  \nFig. 7. Inclination angle obtained with Eq. 2 by varying the mass ratio q in a range of 0 to 0.5. The dashed black lines represent the lower and upper bounds of the inclination values determined for J1357 with model 2. Consequently, the data points highlighted with a red cross symbolise the lower and upper limits of the mass ratio q allowed for the system in order to preserve the lack of eclipses at this specific inclination.  \n<!-- image -->  \nto the low mass ratio of the system, implying that the radius of the donor star is either comparable to or smaller than that of the outer rim of the disc. Consequently, even in an edge-on configuration, the central disc region is not obstructed by eclipses of the companion star. It is possible to give quantitative credibility to this hypothesis by evaluating whether the mass ratio q given in the literature allows the occurrence or absence of eclipses at a particular inclination angle. As stated by Iaria et al. (2018), the angle θ between the line of sight and the equatorial plane can be described by the following equation:  \ntan θ =         R 2 2 -x 2 a 2 -GLYPH<16> R 2 2 -x 2 GLYPH<17>         1 / 2 , (2)  \nwhere R 2 is the Roche-lobe radius of the companion star, and x is the obscured region during the eclipse. Without an eclipse, x can be assumed to be zero, while the Roche-lobe radius can be calculated using Eq. 1 by substituting q -1 with q . By substituting R 2 into the previous equation, the angle θ is clearly only a function of the mass ratio q . In order to investigate the combinations of q and θ that preserve the condition without an eclipse, we varied the q parameter within the range of 0 to 0.5 and determined the corresponding values of θ . The results are shown in Fig. 7.  \nFor an inclination of 83 . 8 + 5 . 4 -1 . 8 degrees obtained at 68% c.l, we can establish a q ratio range between 0.0006 and 0.0275, which is consistent with the value reported by Mata Sánchez et al. (2015). Moreover, with the mass function for M 1, f(M1) = 11 . 0 ± 2 . 1 M ⊙ (Mata Sánchez et al. 2015), we can analytically estimate the mass of the companion star M 2. We establish a range that falls between 0 . 0006 ± 0 . 0001 M ⊙ and 0 . 32 ± 0 . 06 M ⊙ (note that Mata Sánchez et al. 2015 reported a lower limit of M 2 > 0 . 4 M ⊙ ). In other words, the inclination estimated for this system together with the limits on the mass ratio are in line  \n529  \n530  \n531  \n532  \n533  \n534  \n535  \n536  \n537  \n538  \n539  \n540  \n541  \n542  \n543  \n544  \n545  \n546  \n547  \n548  \n549  \n550  \n551  \n552  \n553  \n554  \n555  \n556  \n557  \n558  \n559  \n560  \n561  \n562  \n563  \n564  \n565  \n566  \n567  \n568  \n569  \n570  \n571  \n572  \n573  \n574  \n575  \n576  \n577  \n578  \n579  \n580  \n581  \n582  \n583  \n584  \n585  \n586  \n587  \n588  \n589  \n590  \n591  \n592  \n593  \n594  \n595  \n596  \n597  \n598  \n599  \n600  \n601  \n602  \n603  \n604  \n605  \n606  \n607  \n608  \n609  \n610  \n611  \n612  \n613  \n614  \n615  \nwith the literature, and the mass range for the companion star M 2 characterises it as a low-mass star that fills its Roche lobe.', '5.1.2. Inner and outer radius of the emitting region': 'We found that the H β line in J1357 is emitted from a region bounded between (0 . 16 -0 . 91) RT , while the H α emission line observed in the J1305 spectrum describes two di ff erent emitting regions for outburst and quiescence. During the quiescent phase, the emission line is located between (0 . 1 -1 . 3) RT , while during the outburst, the emitting region extends to radii between (0 . 83 -13) RT and (1 . 1 -19) RT .  \nTable A.1 shows that the outer radius of J1357 is smaller than the expected tidal radius, indicating that either the accretion disc does not extend to the tidal radius or that the H β emission line does not originate in the outermost region of the accretion disc. For J1305, the best-fit value for the outer radius obtained during quiescence is consistent with the expected geometry of the system (although the best-fit range includes values that exceed the tidal radius). On the other hand, the outer radius obtained for the outburst spectra is entirely out of scale. A value of about 1319 RT is higher by up to an order of magnitude than the orbital separation. This deviation is not due to a miscalculation or to a limitation of the application of diskline under certain circumstances. Miceli et al. (2024) reported a Doppler-shift velocity of the peaks during the outburst of approximately 400 km s -1 . Considering a Keplerian disc, we can apply the third Keplerian law to calculate the corresponding radius, and even in this case, the derived radius is four times the tidal radius. This result could be due to the e ff ects during the outburst on the line profile. It is possible that the velocity distribution of the emitting region is not Keplerian. In standard discs (Shakura & Sunyaev 1973), the in-flowing material reaches a stable equilibrium in which the radial component of the velocity is negligible compared to the azimuthal component. However, if there is a surge in viscosity within that region, the radial velocity may cease to be negligible, and the way in which matter accretes may deviate from the way in the Keplerian regime.  \nAnalternative explanation might be linked to the existence of a circumbinary disc around the binary system. Due to the highaccretion rate, outflows of matter can occur, and the ejected material may subsequently remain around the system, causing the formation of a viscous toroidal disc around the binary system (Chen & Podsiadlowski 2019). Although these discs are typically quite cold, the intense emission produced by the central source (Lx ∼ 10 37 erg s -1 for J1305 during outburst, see Miller et al. 2014) might irradiate these regions, thus providing the energy to trigger the atomic transitions that give rise to the H α emission line. This would explain an emitting region greater than the orbital separation itself. However, these are only hypothetical answers to a question that requires further analysis to be understood.', '5.2. Temperatures and hydrogen ionisation': 'Although our results seem to delineate an emitting region that is consistent with the geometry of the disc, at least during quiescence periods, it is crucial to consider the temperatures at these radii. It is indeed true that the emission lines H β and H α are associated with the energy transitions of neutral hydrogen. This means that when the temperatures are too high, the amount of neutral hydrogen will be negligible.  \nShakura & Sunyaev (1973) described the structure of the accretion disc around a BH under the assumption that the disc is optically thick and geometrically thin. The authors divided the disc into three ideal regions, depending on the predominant type of pressure and cross-section, and derived a system of equations for each of them through which the physics of the accretion disc can be completely described. We analysed the outer part of the disc, and the equation that describes the variation in the central temperature of the disc with the radius therefore is  \nT = 8 . 6 · 10 7 α -1 / 5 ˙ m 3 / 10 m -1 / 5 r -3 / 4 GLYPH<16> 1 -r -1 / 2 GLYPH<17> 3 / 10 , (3)  \nwhere m , ˙ m , andr are non-dimensional parameters that contain the dependence on the BH mass, the accretion rate, and the radius, respectively (see Shakura & Sunyaev (1973) for the explicit dependence), while α is the accretion disc viscosity parameter, which usually lies in the range 0.1- 0.4 (King et al. 2007). In this case, we assumed α = 0 . 1. As the density does not vary significantly across the vertical height ( z ) at these large radii, it can reasonably be assumed that the surface temperature of the disc is approximately equal to the central temperature (Shakura &Sunyaev 1973).  \nConsidering the values of Rin and Rout from our analysis, we estimated the temperature ranges for the H β and H α emitting regions during the quiescent phases of J1357 and J1305, as we have already discussed the reliability of the measurements during the outburst. We obtained temperature ranges of 2081 ± 255 K to 353 ± 40 K for J1305 and 2203 ± 317 K to 580 ± 81 K for J1357. This is consistent with a higher quiescent luminosity of LX ∼ 1 . 6 × 10 32 , erg , s -1 for J1305 (Mata Sánchez et al. 2021) compared to J1357, for which LX ∼ 1 . 3 × 10 31 , erg , s -1 (Armas Padilla et al. 2014).  \nThe degree of ionisation for the H-atom, that is, the fraction of ionised hydrogen with respect to neutral hydrogen at the estimated temperatures, can be computed using the Saha equation (Payne 1924),  \nlog H + H = log u + u + log 2 + 5 2 log T -5040 χ ion T -log Pe -0 . 48 . (4)  \nHere, Pe = nekT is the electron pressure, where ne is the electron density, which is typically assumed to be equal to 10 18 cm -3 (Sincell & Krolik 1998). χ ion is the ionization energy (13.6 eV for the H atom), and u is the partition function of the atom (u(H) = 2 and u(H + ) = 1 (see Böhm-Vitense 1992). Assuming the temperatures we inferred above, we obtain a range of H + / H between 10 -25 and 10 -114 for J1357, and between 10 -45 and 10 -186 for J1305. This implies that the amount of ionised hydrogen with respect to neutral hydrogen is negligible, and therefore, the radii estimated for the emitting regions of the H α and H β lines for the two systems appear to be physically plausible.', '6. Conclusion': 'We presented a new approach for analysing the emission lines found in the optical spectra of binary systems that involves the use of the diskline model in order to obtain constraints on the system geometry.  \nWe analysed two observations in quiescence of the X-ray binary BH candidates J1357 and J1305 and an observation collected during an outburst of the latter source. The best-fit parameters allowed us to provide a reasonable description of the geometry of these systems. The H β and H α emission lines in the quiescent spectra of J1357 and J1305 are emitted by a ring in the  \n616  \n617  \n618  \n619  \n620  \n621  \n622  \n623  \n624  \n625  \n626  \n627  \n628  \n629  \n630  \n631  \n632  \n633  \n634  \n635  \n636  \n637  \n638  \n639  \n640  \n641  \n642  \n643  \n644  \n645  \n646  \n647  \n648  \n649  \n650  \n651  \n652  \n653  \n654  \n655  \n656  \n657  \n658  \n659  \n660  \n661  \n662  \n663  \n664  \n665  \n666  \n667  \n668  \n669  \n670  \n671  \n672  \n673  \n674  \n675  \n676  \n677  \n678  \n679  \n680  \n681  \n682  \n683  \n684  \n685  \n686  \n687  \n688  \n689  \n690  \n691  \n692  \n693  \n694  \n695  \n696  \n697  \n698  \n699  \n700  \n701  \n702  \n703  \n704  \n705  \n706  \n707  \n708  \n709  \n710  \n711  \n712  \n713  \n714  \n715  \n716  \n717  \n718  \n719  \n720  \n721  \n722  \n723  \n724  \n725  \n726  \n727  \n728  \n729  \n730  \n731  \n732  \n733  \n734  \n735  \n736  \ndisc between (0 . 16 -0 . 91) RT and (0 . 1 -1 . 3) RT , respectively. The analysis of the outburst spectra yielded an emission region that is not aligned with the expected system geometry, suggesting a range for the outer radius from the tidal radius to values exceeding the orbital separation of the binary system. We can put forward some hypotheses to account for this behaviour, including a non-Keplerian flow in the outer disc or the existence of a circumbinary disc.  \nOur analysis reveals that the inclination angles of the two systems closely match the expected values. This confirms their high-inclination nature. The proposed method requires further investigation because the physics of emission lines from accretion discs is far more complex, and the diskline model provides a simplified description of the emission profile modified by Doppler e ff ects in a Keplerian flow, without any consideration of self-absorption e ff ects in the optically thick atmosphere of accretion discs. However, the application of the diskline model to the disc emission lines in the optical band can be a powerful tool for providing a reliable geometrical description of these sources, allowing us to give precise estimates of the inclination angle, and also to gather information on the expected temperature and ionisation level of the emitting region. Further analyses similar to the analysis we presented here will be useful to provide additional evidence regarding the reliability of this method and to fully explore its capabilities. For example, applying the method to an eclipsing source would represent an optimal approach to test its validity. It would also be interesting to ascertain whether di ff erent emission lines within the same spectrum yield di ff erent insights. In theory, access to a diverse array of emission lines would allow us to delineate the emission from distinct rings of the disc.  \nIt is clear, therefore, that obtaining credible constraints on the parameters of the accretion disc is a challenging task. The application of this method might provide a compatibility check which together with the results that have been obtained with other methods might shed some light on the disc geometry.  \nAcknowledgements. The authors acknowledge the financial contribution PRININAF 2019 with the project "Probing the geometry of accretion: from theory to observations" (PI: Belloni).', 'References': "Anitra, A., Di Salvo, T., Iaria, R., et al. 2021, A&A, 654, A160 Anitra, A., Mata Sánchez, D., Muñoz-Darias, T., et al. 2023, A&A, 679, A145 Appenzeller, I., Fricke, K., Fürtig, W., et al. 1998, The Messenger, 94, 1  \nArmas Padilla, M., Wijnands, R., Altamirano, D., et al. 2014, MNRAS, 439, 3908 Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference Series, Vol. 101, Astronomical Data Analysis Software and Systems V, ed. G. H. Jacoby & J. Barnes, 17 Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167 Böhm-Vitense, E. 1992, Introduction to stellar astrophysics. Volume 3. Stellar structure and evolution., Vol. 3 Cackett, E. M., Miller, J. M., Bhattacharyya, S., et al. 2008, ApJ, 674, 415 Casares, J. 2015, ApJ, 808, 80 Casares, J., Muñoz-Darias, T., Torres, M. A. P., et al. 2022, MNRAS, 516, 2023 Chen, W.-C. & Podsiadlowski, P. 2019, ApJ, 876, L11 Corral-Santana, J. M., Casares, J., Shahbaz, T., et al. 2013, in Revista Mexicana de Astronomia y Astrofisica Conference Series, Vol. 42, Revista Mexicana de Astronomia y Astrofisica Conference Series, 3-4 Cúneo, V. A., Muñoz-Darias, T., Sánchez-Sierras, J., et al. 2020, MNRAS, 498, 25 di Salvo, T., D'Aí, A., Iaria, R., et al. 2009, MNRAS, 398, 2022 Eggleton, P. P. 1983, ApJ, 268, 368 Fabian, A. C., Rees, M. J., Stella, L., & White, N. E. 1989, MNRAS, 238, 729 Frank, J., King, A., & Raine, D. J. 2002, Accretion Power in Astrophysics: Third Edition  \nArticle number, page 10 of 11  \nHorne, K. & Marsh, T. R. 1986, MNRAS, 218, 761 Iaria, R., D'Aí, A., di Salvo, T., et al. 2009, A&A, 505, 1143 Iaria, R., Gambino, A. F., Di Salvo, T., et al. 2018, MNRAS, 473, 3490 Kennea, J., Miller, J. M., Beardmore, A., Degenaar, N., & Reynolds, M. T. 2012,  \nThe Astronomer's Telegram, 4071, 1 King, A. R., Pringle, J. E., & Livio, M. 2007, MNRAS, 376, 1740 Marsh, T. R., Horne, K., & Shipman, H. L. 1987, MNRAS, 225, 551 Marsh, T. R., Robinson, E. L., & Wood, J. H. 1994, MNRAS, 266, 137 Mason, E., Skidmore, W., Howell, S. B., et al. 2000, MNRAS, 318, 440 Mata Sánchez, D., Muñoz-Darias, T., Casares, J., Corral-Santana, J. M., & Shahbaz, T. 2015, MNRAS, 454, 2199  \nMata Sánchez, D., Rau, A., Álvarez Hernández, A., et al. 2021, MNRAS, 506, 581  \nMiceli, C., Mata Sánchez, D., Anitra, A., et al. 2024, A&A, 684, A67  \n750  \nMiller, J. M., Raymond, J., Kallman, T. R., et al. 2014, ApJ, 788, 53  \n751  \nMorihana, K., Sugizaki, M., Nakahira, S., et al. 2013, PASJ, 65, L10  \n752  \nNaylor, T. 1998, MNRAS, 296, 339  \n753  \nOrosz, J. A., Bailyn, C. D., Remillard, R. A., McClintock, J. E., & Foltz, C. B.  \n754  \n1994, ApJ, 436, 848  \n755  \nPayne, C. H. 1924, Nature, 113, 783  \n756  \n- Rayne, M. W. & Whelan, J. A. J. 1981, MNRAS, 196, 73  \n757  \nSchoembs, R. & Hartmann, K. 1983, A&A, 128, 37  \n758  \nShafter, A. W., Szkody, P., & Thorstensen, J. R. 1986, ApJ, 308, 765  \n759  \nShakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337  \n760  \n- Shaw, A. W., Charles, P. A., Casares, J., & Steeghs, D. 2017, in 7 years of MAXI: monitoring X-ray Transients, ed. M. Serino, M. Shidatsu, W. Iwakiri, &T. Mihara, 45\n- Shidatsu, M., Ueda, Y., Nakahira, S., et al. 2013, The Astrophysical Journal, 779, 26  \nSincell, M. W. & Krolik, J. H. 1998, Advances in Space Research, 21, 113 Smak, J. 1981, Acta Astron., 31, 395  \nSteeghs, D. & Casares, J. 2002, ApJ, 568, 273  \n768  \n- Suwa, F., Negoro, H., Nakahira, S., et al. 2012, The Astronomer's Telegram, 4035, 1\n- Torres, M. A. P., Callanan, P. J., Garcia, M. R., et al. 2002, ApJ, 569, 423  \n771  \nTorres, M. A. P., Jonker, P. G., Miller-Jones, J. C. A., et al. 2015, MNRAS, 450,  \n772  \n4292  \n773  \nvan Dokkum, P. G., Bloom, J., & Tewes, M. 2012, L.A.Cosmic: Laplacian Cosmic Ray Identification, Astrophysics Source Code Library, record ascl:1207.005  \nWang, Z., Archibald, A. M., Thorstensen, J. R., et al. 2009, ApJ, 703, 2017  \n777  \n774  \n775  \n776  \n769  \n770  \n761  \n762  \n763  \n764  \n765  \n766  \n767  \n737  \n738  \n739  \n740  \n741  \n742  \n743  \n744  \n745  \n746  \n747  \n748  \n749  \n778", 'T able A.1. Best-fit v alues for the parameters of the tw o models described in the te xt that include the diskline component. Uncertainties are at the': '90%  \nc.l.'}