@article{ToalaBowmanVanReethetal.2022, author = {Toal{\´a}, Jes{\´u}s Alberto and Bowman, Dominic and Van Reeth, Timothy and Todt, Helge Tobias and Dsilva, Karan and Shenar, Tomer and Koenigsberger, Gloria Suzanne and Estrada-Dorado, Sandino and Oskinova, Lida and Hamann, Wolf-Rainer}, title = {Multiple variability time-scales of the early nitrogen-rich Wolf-Rayet star WR 7}, series = {Monthly notices of the Royal Astronomical Society}, volume = {514}, journal = {Monthly notices of the Royal Astronomical Society}, number = {2}, publisher = {Oxford University Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stac1455}, pages = {2269 -- 2277}, year = {2022}, abstract = {We present the analysis of the optical variability of the early, nitrogen-rich Wolf-Rayet (WR) star WR 7. The analysis of multisector Transiting Exoplanet Survey Satellite (TESS) light curves and high-resolution spectroscopic observations confirm multiperiodic variability that is modulated on time-scales of years. We detect a dominant period of 2.6433 +/- 0.0005 d in the TESS sectors 33 and 34 light curves in addition to the previously reported high-frequency features from sector 7. We discuss the plausible mechanisms that may be responsible for such variability in WR 7, including pulsations, binarity, co-rotating interaction regions (CIRs), and clumpy winds. Given the lack of strong evidence for the presence of a stellar or compact companion, we suggest that WR 7 may pulsate in quasi-coherent modes in addition to wind variability likely caused by CIRs on top of stochastic low-frequency variability. WR 7 is certainly a worthy target for future monitoring in both spectroscopy and photometry to sample both the short (less than or similar to 1 d) and long (greater than or similar to 1000 d) variability time-scales.}, language = {en} } @article{AlmeidaSanaTayloretal.2017, author = {Almeida, Leonardo A. and Sana, H. and Taylor, W. and Barb{\´a}, Rodolfo and Bonanos, Alceste Z. and Crowther, Paul and Damineli, Augusto and de Koter, A. and de Mink, Selma E. and Evans, C. J. and Gieles, Mark and Grin, Nathan J. and H{\´e}nault-Brunet, V. and Langer, Norbert and Lennon, D. and Lockwood, Sean and Ma{\´i}z Apell{\´a}niz, Jes{\´u}s and Moffat, A. F. J. and Neijssel, C. and Norman, C. and Ram{\´i}rez-Agudelo, O. H. and Richardson, N. D. and Schootemeijer, Abel and Shenar, Tomer and Soszyński, Igor and Tramper, Frank and Vink, J. S.}, title = {The tarantula massive binary monitoring}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {598}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201629844}, pages = {36}, year = {2017}, abstract = {Context: Massive binaries play a crucial role in the Universe. Knowing the distributions of their orbital parameters is important for a wide range of topics from stellar feedback to binary evolution channels and from the distribution of supernova types to gravitational wave progenitors, yet no direct measurements exist outside the Milky Way. Aims: The Tarantula Massive Binary Monitoring project was designed to help fill this gap by obtaining multi-epoch radial velocity (RV) monitoring of 102 massive binaries in the 30 Doradus region. Methods: In this paper we analyze 32 FLAMES/GIRAFFE observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single-lined (SB1) and 31 double-lined (SB2) spectroscopic binaries. Results: Overall, the binary fraction and orbital properties across the 30 Doradus region are found to be similar to existing Galactic samples. This indicates that within these domains environmental effects are of second order in shaping the properties of massive binary systems. A small difference is found in the distribution of orbital periods, which is slightly flatter (in log space) in 30 Doradus than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, orbital periods in 30 Doradus can be as short as 1.1 d, somewhat shorter than seen in Galactic samples. Equal mass binaries (q> 0.95) in 30 Doradus are all found outside NGC 2070, the central association that surrounds R136a, the very young and massive cluster at 30 Doradus's core. Most of the differences, albeit small, are compatible with expectations from binary evolution. One outstanding exception, however, is the fact that earlier spectral types (O2-O7) tend to have shorter orbital periods than later spectral types (O9.2-O9.7). Conclusions: Our results point to a relative universality of the incidence rate of massive binaries and their orbital properties in the metallicity range from solar (Z⊙) to about half solar. This provides the first direct constraints on massive binary properties in massive star-forming galaxies at the Universe's peak of star formation at redshifts z ~ 1 to 2 which are estimated to have Z ~ 0.5 Z⊙.}, language = {en} } @article{ShenarRichardsonSablowskietal.2017, author = {Shenar, Tomer and Richardson, N. D. and Sablowski, Daniel P. and Hainich, Rainer and Sana, H. and Moffat, A. F. J. and Todt, Helge Tobias and Hamann, Wolf-Rainer and Oskinova, Lida and Sander, Andreas Alexander Christoph and Tramper, Frank and Langer, Norbert and Bonanos, Alceste Z. and de Mink, Selma E. and Gr{\"a}fener, G. and Crowther, Paul and Vink, J. S. and Almeida, Leonardo A. and de Koter, A. and Barb{\´a}, Rodolfo and Herrero, A. and Ulaczyk, Krzysztof}, title = {The tarantula massive binary monitoring}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {598}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201629621}, pages = {16}, year = {2017}, abstract = {We present the first SB2 orbital solution and disentanglement of the massive Wolf-Rayet binary R145 (P = 159 d) located in the Large Magellanic Cloud. The primary was claimed to have a stellar mass greater than 300 M-circle dot, making it a candidate for being the most massive star known to date. While the primary is a known late-type, H-rich Wolf-Rayet star (WN6h), the secondary has so far not been unambiguously detected. Using moderate-resolution spectra, we are able to derive accurate radial velocities for both components. By performing simultaneous orbital and polarimetric analyses, we derive the complete set of orbital parameters, including the inclination. The spectra are disentangled and spectroscopically analyzed, and an analysis of the wind-wind collision zone is conducted. The disentangled spectra and our models are consistent with a WN6h type for the primary and suggest that the secondary is an O3.5 If*/WN7 type star. We derive a high eccentricity of e = 0 : 78 and minimum masses of M-1 sin(3) i approximate to M-2 sin(3) i = 13 +/- 2 M-circle dot, with q = M-2/M-1 = 1.01 +/- 0.07. An analysis of emission excess stemming from a wind-wind collision yields an inclination similar to that obtained from polarimetry (i = 39 +/- 6 degrees). Our analysis thus implies M-1 = 53(-20)(+40) and M2 = 54(-20)(+40) M-circle dot, excluding M-1 > 300 M-circle dot. A detailed comparison with evolution tracks calculated for single and binary stars together with the high eccentricity suggests that the components of the system underwent quasi-homogeneous evolution and avoided mass-transfer. This scenario would suggest current masses of approximate to 80 M-circle dot and initial masses of M-i,M-1 approximate to 10(5) and M-i,M-2 approximate to 90 M-circle dot, consistent with the upper limits of our derived orbital masses, and would imply an age of approximate to 2.2 Myr.}, language = {en} } @article{ShenarOskinovaJaervinenetal.2017, author = {Shenar, Tomer and Oskinova, Lida and Jaervinen, S. P. and Luckas, P. and Hainich, Rainer and Todt, Helge Tobias and Hubrig, Swetlana and Sander, Andreas Alexander Christoph and Ilyin, Ilya and Hamann, Wolf-Rainer}, title = {Constraining the weak-wind problem}, series = {Contributions Of The Astronomical Observatory Skalnate Pleso}, volume = {48}, journal = {Contributions Of The Astronomical Observatory Skalnate Pleso}, number = {1}, publisher = {Astronomick{\´y} {\´U}stav SAV}, address = {Tatransk{\´a} Lomnica}, issn = {1335-1842}, doi = {10.1051/0004-6361/201731291}, pages = {139 -- 143}, year = {2017}, abstract = {Mass-loss rates of massive, late type main sequence stars are much weaker than currently predicted, but their true values are very difficult to measure. We suggest that confined stellar winds of magnetic stars can be exploited to constrain the true mass-loss rates M of massive main sequence stars. We acquired UV, X-ray, and optical amateur data of HD 54879 (09.7 V), one of a few O-type stars with a detected atmospheric magnetic field (B-d greater than or similar to 2 kG). We analyze these data with the Potsdam Wolf-Rayet (PoWR) and XSPEC codes. We can roughly estimate the mass-loss rate the star would have in the absence of a magnetic field as log M-B=0 approximate to -9.0 M-circle dot yr(-1). Since the wind is partially trapped within the Alfven radius rA greater than or similar to 12 R-*,, the true mass-loss rate of HD 54879 is log M less than or similar to -10.2 M-circle dot yr(-1). Moreover, we find that the microturbulent, macroturbulent, and projected rotational velocities are lower than previously suggested (< 4 km s(-1)). An initial mass of 16 M-circle dot and an age of 5 Myr are inferred. We derive a mean X-ray emitting temperature of log T-x = 6.7 K and an X-ray luminosity of log L-x = 32 erg s(-1). The latter implies a significant X-ray excess (log L-x/L-Bol approximate to - 6.0), most likely stemming from collisions at the magnetic equator. A tentative period of P approximate to 5 yr is derived from variability of the Ha line. Our study confirms that strongly magnetized stars lose little or no mass, and supplies important constraints on the weak-wind problem of massive main sequence stars.}, language = {en} } @article{HainichOskinovaShenaretal.2018, author = {Hainich, Rainer and Oskinova, Lida and Shenar, Tomer and Marchant Campos, Pablo and Eldridge, J. J. and Sander, Andreas Alexander Christoph and Hamann, Wolf-Rainer and Langer, Norbert and Todt, Helge Tobias}, title = {Observational properties of massive black hole binary progenitors}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {609}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201731449}, pages = {62}, year = {2018}, abstract = {Context: The first directly detected gravitational waves (GW 150914) were emitted by two coalescing black holes (BHs) with masses of ≈ 36 M⊙ and ≈ 29 M⊙. Several scenarios have been proposed to put this detection into an astrophysical context. The evolution of an isolated massive binary system is among commonly considered models. Aims: Various groups have performed detailed binary-evolution calculations that lead to BH merger events. However, the question remains open as to whether binary systems with the predicted properties really exist. The aim of this paper is to help observers to close this gap by providing spectral characteristics of massive binary BH progenitors during a phase where at least one of the companions is still non-degenerate. Methods: Stellar evolution models predict fundamental stellar parameters. Using these as input for our stellar atmosphere code (Potsdam Wolf-Rayet), we compute a set of models for selected evolutionary stages of massive merging BH progenitors at different metallicities. Results: The synthetic spectra obtained from our atmosphere calculations reveal that progenitors of massive BH merger events start their lives as O2-3V stars that evolve to early-type blue supergiants before they undergo core-collapse during the Wolf-Rayet phase. When the primary has collapsed, the remaining system will appear as a wind-fed high-mass X-ray binary. Based on our atmosphere models, we provide feedback parameters, broad band magnitudes, and spectral templates that should help to identify such binaries in the future. Conclusions: While the predicted parameter space for massive BH binary progenitors is partly realized in nature, none of the known massive binaries match our synthetic spectra of massive BH binary progenitors exactly. Comparisons of empirically determined mass-loss rates with those assumed by evolution calculations reveal significant differences. The consideration of the empirical mass-loss rates in evolution calculations will possibly entail a shift of the maximum in the predicted binary-BH merger rate to higher metallicities, that is, more candidates should be expected in our cosmic neighborhood than previously assumed.}, language = {en} } @article{RamiaramanantsoaMoffatHarmonetal.2018, author = {Ramiaramanantsoa, Tahina and Moffat, Anthony F. J. and Harmon, Robert and Ignace, R. and St-Louis, Nicole and Vanbeveren, Dany and Shenar, Tomer and Pablo, Herbert and Richardson, Noel D. and Howarth, Ian D. and Stevens, Ian R. and Piaulet, Caroline and St-Jean, Lucas and Eversberg, Thomas and Pigulski, Andrzej and Popowicz, Adam and Kuschnig, Rainer and Zoclonska, Elzbieta and Buysschaert, Bram and Handler, Gerald and Weiss, Werner W. and Wade, Gregg A. and Rucinski, Slavek M. and Zwintz, Konstanze and Luckas, Paul and Heathcote, Bernard and Cacella, Paulo and Powles, Jonathan and Locke, Malcolm and Bohlsen, Terry and Chen{\´e}, Andr{\´e}-Nicolas and Miszalski, Brent and Waldron, Wayne L. and Kotze, Marissa M. and Kotze, Enrico J. and B{\"o}hm, Torsten}, title = {BRITE-Constellation high-precision time-dependent photometry of the early O-type supergiant zeta Puppis unveils the photospheric drivers of its small- and large-scale wind structures}, series = {Monthly notices of the Royal Astronomical Society}, volume = {473}, journal = {Monthly notices of the Royal Astronomical Society}, number = {4}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stx2671}, pages = {5532 -- 5569}, year = {2018}, abstract = {From 5.5 months of dual-band optical photometric monitoring at the 1 mmag level, BRITE-Constellation has revealed two simultaneous types of variability in the O4I(n)fp star ζ Puppis: one single periodic non-sinusoidal component superimposed on a stochastic component. The monoperiodic component is the 1.78-d signal previously detected by Coriolis/Solar Mass Ejection Imager, but this time along with a prominent first harmonic. The shape of this signal changes over time, a behaviour that is incompatible with stellar oscillations but consistent with rotational modulation arising from evolving bright surface inhomogeneities. By means of a constrained non-linear light-curve inversion algorithm, we mapped the locations of the bright surface spots and traced their evolution. Our simultaneous ground-based multisite spectroscopic monitoring of the star unveiled cyclical modulation of its He ii λ4686 wind emission line with the 1.78-d rotation period, showing signatures of corotating interaction regions that turn out to be driven by the bright photospheric spots observed by BRITE. Traces of wind clumps are also observed in the He ii λ4686 line and are correlated with the amplitudes of the stochastic component of the light variations probed by BRITE at the photosphere, suggesting that the BRITE observations additionally unveiled the photospheric drivers of wind clumps in ζ Pup and that the clumping phenomenon starts at the very base of the wind. The origins of both the bright surface inhomogeneities and the stochastic light variations remain unknown, but a subsurface convective zone might play an important role in the generation of these two types of photospheric variability.}, language = {en} } @article{SanderFuerstKretschmaretal.2018, author = {Sander, Andreas Alexander Christoph and F{\"u}rst, F. and Kretschmar, P. and Oskinova, Lida and Todt, Helge Tobias and Hainich, Rainer and Shenar, Tomer and Hamann, Wolf-Rainer}, title = {Coupling hydrodynamics with comoving frame radiative transfer}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {610}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201731575}, pages = {19}, year = {2018}, abstract = {Aims. To gain a realistic picture of the donor star in Vela X-1, we constructed a hydrodynamically consistent atmosphere model describing the wind stratification while properly reproducing the observed donor spectrum. To investigate how X-ray illumination affects the stellar wind, we calculated additional models for different X-ray luminosity regimes. Methods. We used the recently updated version of the Potsdam Wolf-Rayet code to consistently solve the hydrodynamic equation together with the statistical equations and the radiative transfer. Results. The wind flow in Vela X-1 is driven by ions from various elements, with Fe III and S III leading in the outer wind. The model-predicted mass-loss rate is in line with earlier empirical studies. The mass-loss rate is almost unaffected by the presence of the accreting NS in the wind. The terminal wind velocity is confirmed at u(infinity) approximate to 600 km s(-1). On the other hand, the wind velocity in the inner region where the NS is located is only approximate to 100 km s(-1), which is not expected on the basis of a standard beta-velocity law. In models with an enhanced level of X-rays, the velocity field in the outer wind can be altered. If the X-ray flux is too high, the acceleration breaks down because the ionization increases. Conclusions. Accounting for radiation hydrodynamics, our Vela X-1 donor atmosphere model reveals a low wind speed at the NS location, and it provides quantitative information on wind driving in this important HMXB.}, language = {en} } @article{RamachandranHamannHainichetal.2018, author = {Ramachandran, Varsha and Hamann, Wolf-Rainer and Hainich, Rainer and Oskinova, Lida and Shenar, Tomer and Sander, Andreas Alexander Christoph and Todt, Helge Tobias and Gallagher, John S.}, title = {Stellar population of the superbubble N206 in the LMC II. Parameters of the OB and WR stars, and the total massive star feedback}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {615}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201832816}, pages = {72}, year = {2018}, abstract = {Context. Clusters or associations of early-type stars are often associated with a "superbubble" of hot gas. The formation of such superbubbles is caused by the feedback from massive stars. The complex N206 in the Large Magellanic Cloud (LMC) exhibits a superbubble and a rich massive star population. Aims. Our goal is to perform quantitative spectral analyses of all massive stars associated with the N206 superbubble in order to determine their stellar and wind parameters. We compare the superbubble energy budget to the stellar energy input and discuss the star formation history of the region. Results. We present the stellar and wind parameters of the OB stars and the two Wolf-Rayet (WR) binaries in the N206 complex. Twelve percent of the sample show Oe/Be type emission lines, although most of them appear to rotate far below critical. We found eight runaway stars based on their radial velocity. The wind-momentum luminosity relation of our OB sample is consistent with the expectations. The Hertzsprung-Russell diagram (HRD) of the OB stars reveals a large age spread (1-30 Myr), suggesting different episodes of star formation in the complex. The youngest stars are concentrated in the inner part of the complex, while the older OB stars are scattered over outer regions. We derived the present day mass function for the entire N206 complex as well as for the cluster NGC2018. The total ionizing photon flux produced by all massive stars in the N206 complex is Q(0) approximate to 5 x 10(50) s(-1), and the mechanical luminosity of their stellar winds amounts to L-mec = 1.7 x 10(38) erg s(-1). Three very massive Of stars are found to dominate the feedback among 164 OB stars in the sample. The two WR winds alone release about as much mechanical luminosity as the whole OB star sample. The cumulative mechanical feedback from all massive stellar winds is comparable to the combined mechanical energy of the supernova explosions that likely occurred in the complex. Accounting also for the WR wind and supernovae, the mechanical input over the last five Myr is approximate to 2.3 x 10(52) erg. Conclusions. The N206 complex in the LMC has undergone star formation episodes since more than 30 Myr ago. From the spectral analyses of its massive star population, we derive a current star formation rate of 2.2 x 10(-3) M-circle dot yr(-1). From the combined input of mechanical energy from all stellar winds, only a minor fraction is emitted in the form of X-rays. The corresponding input accumulated over a long time also exceeds the current energy content of the complex by more than a factor of five. The morphology of the complex suggests a leakage of hot gas from the superbubble.}, language = {en} } @article{ShenarHainichTodtetal.2018, author = {Shenar, Tomer and Hainich, Rainer and Todt, Helge Tobias and Moffat, Anthony F. J. and Sander, Andreas Alexander Christoph and Oskinova, Lida and Ramachandran, Varsha and Munoz, M. and Pablo, H. and Sana, Hugues and Hamann, Wolf-Rainer}, title = {The shortest-period Wolf-Rayet binary in the small magellanic cloud}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {616}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201833006}, pages = {15}, year = {2018}, abstract = {Context. SMC AB6 is the shortest-period (P = 6.5 d) Wolf-Rayet (WR) binary in the Small Magellanic Cloud. This binary is therefore a key system in the study of binary interaction and formation of WR stars at low metallicity. The WR component in AB6 was previously found to be very luminous (log L = 6.3 [L-circle dot]) compared to its reported orbital mass (approximate to 8 M-circle dot), placing it significantly above the Eddington limit. Aims. Through spectroscopy and orbital analysis of newly acquired optical data taken with the Ultraviolet and Visual Echelle Spectrograph (UVES), we aim to understand the peculiar results reported for this system and explore its evolutionary history. Methods. We measured radial velocities via cross-correlation and performed a spectral analysis using the Potsdam Wolf-Rayet model atmosphere code. The evolution of the system was analyzed using the Binary Population and Spectral Synthesis evolution code. Results. AB6 contains at least four stars. The 6.5 d period WR binary comprises the WR primary (WN3:h, star A) and a rather rapidly rotating (v(eq) = 265 km s(-1)) early O-type companion (O5.5 V, star B). Static N III and N IV emission lines and absorption signatures in He lines suggest the presence of an early-type emission line star (O5.5 I(f), star C). Finally, narrow absorption lines portraying a long-term radial velocity variation show the existence of a fourth star (O7.5 V, star D). Star D appears to form a second 140 d period binary together with a fifth stellar member, which is a B-type dwarf or a black hole. It is not clear that these additional components are bound to the WR binary. We derive a mass ratio of M-O/M-WR = 2.2 +/- 0.1. The WR star is found to be less luminous than previously thought (log L = 5.9 [L-circle dot]) and, adopting M-O = 41 M-circle dot for star B, more massive (M-WR = 18 M-circle dot). Correspondingly, the WR star does not exceed the Eddington limit. We derive the initial masses of M-i,M-WR = 60 M-circle dot and M-i,M-O = 40 M-circle dot and an age of 3.9 Myr for the system. The WR binary likely experienced nonconservative mass transfer in the past supported by the relatively rapid rotation of star B. Conclusions. Our study shows that AB6 is a multiple - probably quintuple - system. This finding resolves the previously reported puzzle of the WR primary exceeding the Eddington limit and suggests that the WR star exchanged mass with its companion in the past.}, language = {en} } @article{RamiaramanantsoaRatnasingamShenaretal.2018, author = {Ramiaramanantsoa, Tahina and Ratnasingam, Rathish and Shenar, Tomer and Moffat, Anthony F. J. and Rogers, Tamara M. and Popowicz, Adam and Kuschnig, Rainer and Pigulski, Andrzej and Handler, Gerald and Wade, Gregg A. and Zwintz, Konstanze and Weiss, Werner W.}, title = {A BRITE view on the massive O-type supergiant V973 Scorpii}, series = {Monthly notices of the Royal Astronomical Society}, volume = {480}, journal = {Monthly notices of the Royal Astronomical Society}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/sty1897}, pages = {972 -- 986}, year = {2018}, abstract = {Stochastically triggered photospheric light variations reaching similar to 40 mmag peak-to-valley amplitudes have been detected in the O8 Iaf supergiant V973 Scorpii as the outcome of 2 months of high-precision time-resolved photometric observations with the BRIght Target Explorer (BRITE) nanosatellites. The amplitude spectrum of the time series photometry exhibits a pronounced broad bump in the low-frequency regime (less than or similar to 0.9 d(-1)) where several prominent frequencies are detected. A time-frequency analysis of the observations reveals typical mode lifetimes of the order of 5-10 d. The overall features of the observed brightness amplitude spectrum of V973 Sco match well with those extrapolated from two-dimensional hydrodynamical simulations of convectively driven internal gravity waves randomly excited from deep in the convective cores of massive stars. An alternative or additional possible source of excitation from a sub-surface convection zone needs to be explored in future theoretical investigations.}, language = {en} }