Showing 1–50 of 98 results for author: Eriksson, K

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

The ACROBAT 2022 Challenge: Automatic Registration Of Breast Cancer Tissue

Authors: Philippe Weitz, Masi Valkonen, Leslie Solorzano, Circe Carr, Kimmo Kartasalo, Constance Boissin, Sonja Koivukoski, Aino Kuusela, Dusan Rasic, Yanbo Feng, Sandra Sinius Pouplier, Abhinav Sharma, Kajsa Ledesma Eriksson, Stephanie Robertson, Christian Marzahl, Chandler D. Gatenbee, Alexander R. A. Anderson, Marek Wodzinski, Artur Jurgas, Niccolò Marini, Manfredo Atzori, Henning Müller, Daniel Budelmann, Nick Weiss, Stefan Heldmann , et al. (16 additional authors not shown)

Abstract: The alignment of tissue between histopathological whole-slide-images (WSI) is crucial for research and clinical applications. Advances in computing, deep learning, and availability of large WSI datasets have revolutionised WSI analysis. Therefore, the current state-of-the-art in WSI registration is unclear. To address this, we conducted the ACROBAT challenge, based on the largest WSI registration… ▽ More The alignment of tissue between histopathological whole-slide-images (WSI) is crucial for research and clinical applications. Advances in computing, deep learning, and availability of large WSI datasets have revolutionised WSI analysis. Therefore, the current state-of-the-art in WSI registration is unclear. To address this, we conducted the ACROBAT challenge, based on the largest WSI registration dataset to date, including 4,212 WSIs from 1,152 breast cancer patients. The challenge objective was to align WSIs of tissue that was stained with routine diagnostic immunohistochemistry to its H&E-stained counterpart. We compare the performance of eight WSI registration algorithms, including an investigation of the impact of different WSI properties and clinical covariates. We find that conceptually distinct WSI registration methods can lead to highly accurate registration performances and identify covariates that impact performances across methods. These results establish the current state-of-the-art in WSI registration and guide researchers in selecting and developing methods. △ Less

Submitted 29 May, 2023; originally announced May 2023.

The carbon star DY Persei may be a cool R Coronae Borealis variable

Authors: D. A. Garcia-Hernandez, N. Kameswara Rao, D. L. Lambert, K. Eriksson, A. B. S. Reddy, T. Masseron

Abstract: Optical and near-IR photometry suggests that the carbon star DY Persei exhibits fadings similar to those of R Coronae Borealis (RCB) variables. Photometric surveys of the Galaxy and Magellanic Clouds uncovered new DY Per variables with infrared photometry identifying them with cool carbon stars, perhaps, with an unusual tendency to shed mass. In an attempt to resolve DY Per's identity crisis -- a… ▽ More Optical and near-IR photometry suggests that the carbon star DY Persei exhibits fadings similar to those of R Coronae Borealis (RCB) variables. Photometric surveys of the Galaxy and Magellanic Clouds uncovered new DY Per variables with infrared photometry identifying them with cool carbon stars, perhaps, with an unusual tendency to shed mass. In an attempt to resolve DY Per's identity crisis -- a cool carbon giant or a cool RCB variable? -- we analyze a high-resolution H&K band spectrum of DY Per. The CO first-overtone bands in the K-band of DY Per show a high abundance of 18O such that 16O/18O = 4 +- 1, a ratio sharply at odds with published results for `regular' cool carbon giants with 16O/18O ~ 1000 but this exceptionally low ratio is characteristic of RCB-variables and HdC stars. This similarity suggests that DY Per indeed may be a cool RCB variable. Current opinion considers RCB-variables to result from merger of a He onto a CO white dwarf; observed abundances of these H-deficient stars including the exceptionally low 16O/18O ratios are in fair accord with predicted compositions for white dwarf merger products. A H-deficiency for DY Per is not directly observable but is suggested from the strength of a HF line and an assumption that F may be overabundant, as observed and predicted for RCB stars. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ (16 pages and 4 figures)

Synthetic photometry for carbon-rich giants. V. Effects of grain-size-dependent dust opacities

Authors: Kjell Eriksson, Susanne Höfner, Bernhard Aringer

Abstract: The properties and the evolution of asymptotic giant branch (AGB) stars are strongly influenced by their mass loss through a stellar wind. This is believed to be caused by radiation pressure due to the absorption and scattering of the stellar radiation by the dust grains formed in the atmosphere. The optical properties of dust are often estimated using the small particle limit (SPL) approximation,… ▽ More The properties and the evolution of asymptotic giant branch (AGB) stars are strongly influenced by their mass loss through a stellar wind. This is believed to be caused by radiation pressure due to the absorption and scattering of the stellar radiation by the dust grains formed in the atmosphere. The optical properties of dust are often estimated using the small particle limit (SPL) approximation, and it has been used frequently in modelling AGB stellar winds when performing radiation-hydrodynamics (RHD) simulations. We aim to investigate the effects of replacing the SPL approximation by detailed Mie calculations of the size-dependent opacities for grains of amorphous carbon forming in C-rich AGB star atmospheres. We performed RHD simulations for a large grid of carbon star atmosphere+wind models with different effective temperatures, luminosities, stellar masses, carbon excesses, and pulsation properties. Also, a posteriori radiative transfer calculations for many radial structures (snapshots) of these models were done, resulting in spectra and filter magnitudes. We find, when giving up the SPL approximation, the wind models become more strongly variable and more dominated by gusts, although the average mass-loss rates and outflow speeds do not change significantly; the increased radiative pressure on the dust throughout its formation zone does, however, result in smaller grains and lower condensation fractions (and thus higher gas-to-dust ratios). The photometric K magnitudes are generally brighter, but at V the effects of using size-dependent dust opacities are more complex: brighter for low mass-loss rates and dimmer for massive stellar winds. Given the large effects on spectra and photometric properties, it is necessary to use the detailed dust optical data instead of the simple SPL approximation in stellar atmosphere+wind modelling where dust is formed. △ Less

Submitted 17 March, 2023; originally announced March 2023.

Comments: 14 pages, 24 figures. Accepted for publication in Astronomy \& Astrophysics

Journal ref: A&A 673, A21 (2023)

ACROBAT -- a multi-stain breast cancer histological whole-slide-image data set from routine diagnostics for computational pathology

Authors: Philippe Weitz, Masi Valkonen, Leslie Solorzano, Circe Carr, Kimmo Kartasalo, Constance Boissin, Sonja Koivukoski, Aino Kuusela, Dusan Rasic, Yanbo Feng, Sandra Kristiane Sinius Pouplier, Abhinav Sharma, Kajsa Ledesma Eriksson, Leena Latonen, Anne-Vibeke Laenkholm, Johan Hartman, Pekka Ruusuvuori, Mattias Rantalainen

Abstract: The analysis of FFPE tissue sections stained with haematoxylin and eosin (H&E) or immunohistochemistry (IHC) is an essential part of the pathologic assessment of surgically resected breast cancer specimens. IHC staining has been broadly adopted into diagnostic guidelines and routine workflows to manually assess status and scoring of several established biomarkers, including ER, PGR, HER2 and KI67.… ▽ More The analysis of FFPE tissue sections stained with haematoxylin and eosin (H&E) or immunohistochemistry (IHC) is an essential part of the pathologic assessment of surgically resected breast cancer specimens. IHC staining has been broadly adopted into diagnostic guidelines and routine workflows to manually assess status and scoring of several established biomarkers, including ER, PGR, HER2 and KI67. However, this is a task that can also be facilitated by computational pathology image analysis methods. The research in computational pathology has recently made numerous substantial advances, often based on publicly available whole slide image (WSI) data sets. However, the field is still considerably limited by the sparsity of public data sets. In particular, there are no large, high quality publicly available data sets with WSIs of matching IHC and H&E-stained tissue sections. Here, we publish the currently largest publicly available data set of WSIs of tissue sections from surgical resection specimens from female primary breast cancer patients with matched WSIs of corresponding H&E and IHC-stained tissue, consisting of 4,212 WSIs from 1,153 patients. The primary purpose of the data set was to facilitate the ACROBAT WSI registration challenge, aiming at accurately aligning H&E and IHC images. For research in the area of image registration, automatic quantitative feedback on registration algorithm performance remains available through the ACROBAT challenge website, based on more than 37,000 manually annotated landmark pairs from 13 annotators. Beyond registration, this data set has the potential to enable many different avenues of computational pathology research, including stain-guided learning, virtual staining, unsupervised pre-training, artefact detection and stain-independent models. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Satellites in the Ti~1$s$ core level spectra of SrTiO$_3$ and TiO$_2$

Authors: Atsushi Hariki, Keisuke Higashi, Tatsuya Yamaguchi, Jiebin Li, Curran Kalha, Manfred Mascheck, Susanna K. Eriksson, Tomas Wiell, Frank M. F. de Groot, Anna Regoutz

Abstract: Satellites in core level spectra of photoelectron spectroscopy (PES) can provide crucial information on the electronic structure and chemical bonding in materials, particular in transition metal oxides. This paper explores satellites of the Ti 1$s$ and 2$p$ core level spectra of SrTiO$_3$ and TiO$_2$. Conventionally, soft x-ray PES (SXPS) probes the Ti 2$p$ core level; however, it is not ideal to… ▽ More Satellites in core level spectra of photoelectron spectroscopy (PES) can provide crucial information on the electronic structure and chemical bonding in materials, particular in transition metal oxides. This paper explores satellites of the Ti 1$s$ and 2$p$ core level spectra of SrTiO$_3$ and TiO$_2$. Conventionally, soft x-ray PES (SXPS) probes the Ti 2$p$ core level; however, it is not ideal to fully capture satellite features due to its inherent spin-orbit-splitting (SOS). Here, hard x-ray PES(HAXPES) provides access to the Ti 1$s$ spectrum instead, which allows us to study intrinsic charge responses upon core-hole creation without the complication from SOS and with favorable intrinsic linewidths. The experimental spectra are theoretically analyzed by two impurity models, including an Anderson impurity model (AIM) built on local density approximation (LDA) and dynamical mean-field theory (DMFT), and a conventional TiO$_6$ cluster model. The theoretical results emphasize the importance of explicit inclusion of higher-order Ti-O charge-transfer processes beyond the nearest-neighboring Ti-O bond to simulate the core level spectra of SrTiO$_3$ and TiO$_2$. The AIM approach with continuous bath orbitals provided by LDA+DMFT represents the experimental spectra well. Crucially, with the aid of the LDA+DMFT method, this paper provides a robust prescription of how to use the computationally cheap cluster model in fitting analyses of core level spectra. △ Less

Submitted 6 January, 2023; v1 submitted 17 June, 2022; originally announced June 2022.

Journal ref: Phys. Rev. B 106, 205138, 2022

Quantum diffusion as a process of linear quantum dynamics

Authors: Karl-Erik Eriksson

Abstract: Quantum diffusion, as developed in the 1990s, could explain how a system, subject to measurement, goes into an eigenstate of the measured observable. Here it is shown that quantum diffusion theory can be interpreted as a result within linear relativistic quantum mechanics. Thus, in contrast to what is widely believed, quantum measurement can be analyzed within the theory of quantum mechanics itsel… ▽ More Quantum diffusion, as developed in the 1990s, could explain how a system, subject to measurement, goes into an eigenstate of the measured observable. Here it is shown that quantum diffusion theory can be interpreted as a result within linear relativistic quantum mechanics. Thus, in contrast to what is widely believed, quantum measurement can be analyzed within the theory of quantum mechanics itself. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: 11 pages

Dynamic atmospheres and winds of cool luminous giants, II. Gradual Fe enrichment of wind-driving silicate grains

Authors: Susanne Höfner, Sara Bladh, Bernhard Aringer, Kjell Eriksson

Abstract: The winds observed around AGB stars are generally attributed to radiation pressure on dust formed in the dynamical atmospheres of these long-period variables. The composition of wind-driving grains is affected by a feedback between their optical properties and the resulting heating due to stellar radiation. We explore the gradual Fe enrichment of wind-driving silicate grains in M-type AGB stars to… ▽ More The winds observed around AGB stars are generally attributed to radiation pressure on dust formed in the dynamical atmospheres of these long-period variables. The composition of wind-driving grains is affected by a feedback between their optical properties and the resulting heating due to stellar radiation. We explore the gradual Fe enrichment of wind-driving silicate grains in M-type AGB stars to derive typical values for Fe/Mg and test the effects on wind properties and synthetic spectra. We present new DARWIN models that allow for the growth of silicate grains with a variable Fe/Mg ratio and predict mass-loss rates, wind velocities, and grain properties. Synthetic spectra and other observables are computed with the COMA code. The self-regulating feedback between grain composition and radiative heating, in combination with quickly falling densities in the stellar wind, leads to low values of Fe/Mg, typically a few percent. Nevertheless, the models show distinct silicate features around 10 and 18 microns. Fe enrichment affects visual and near-IR photometry moderately, and the new DARWIN models agree well with observations in (J-K) vs. (V-K) and Spitzer color-color diagrams. The enrichment of the silicate dust with Fe is a secondary process, taking place in the stellar wind on the surface of large Fe-free grains that have initiated the outflow. Therefore, the mass-loss rates are basically unaffected, while the wind velocities tend to be slightly higher than in corresponding models with Fe-free silicate dust. The gradual Fe enrichment of silicate grains in the inner wind region should produce signatures observable in mid-IR spectro-interferometrical measurements. Mass-loss rates derived from existing DARWIN models, based on Fe-free silicates, can be applied to stellar evolution models since the mass-loss rates are not significantly affected by the inclusion of Fe in the silicate grains. △ Less

Submitted 29 October, 2021; originally announced October 2021.

Comments: Accepted for publication in Astronomy & Astrophysics (15 pages, 14 figures)

Journal ref: A&A 657, A109 (2022)

Give quantum mechanics a chance: use relativistic quantum mechanics to analyze measurement!

Authors: Karl-Erik Eriksson

Abstract: At the time of publication of H. Everett's Relative-State Formulation (1957) and DeWitt's Many-Worlds Interpretation (1970), quantum mechanics was available in a more modern and adequate version than the one used by these authors. We show that with the more modern quantum theory, quantum measurement could have been analyzed along more conventional lines in a one-world cosmology. Bell criticized th… ▽ More At the time of publication of H. Everett's Relative-State Formulation (1957) and DeWitt's Many-Worlds Interpretation (1970), quantum mechanics was available in a more modern and adequate version than the one used by these authors. We show that with the more modern quantum theory, quantum measurement could have been analyzed along more conventional lines in a one-world cosmology. Bell criticized the Everett-DeWitt theory quite sharply in 1987 but this seems not to have affected the acceptance of the old quantum mechanics as the framework for analysis of measurement. △ Less

Submitted 14 July, 2021; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: 28 pages, 2 figures

Effects of nitridation on SiC/SiO$_2$ structures studied by hard X-ray photoelectron spectroscopy

Authors: Judith Berens, Sebastian Bichelmaier, Nathalie K. Fernando, Pardeep K. Thakur, Tien-Lin Lee, Manfred Mascheck, Tomas Wiell, Susanna K. Eriksson, J. Matthias Kahk, Johannes Lischner, Manesh V. Mistry, Thomas Aichinger, Gregor Pobegen, Anna Regoutz

Abstract: SiC is set to enable a new era in power electronics impacting a wide range of energy technologies, from electric vehicles to renewable energy. Its physical characteristics outperform silicon in many aspects, including band gap, breakdown field, and thermal conductivity. The main challenge for further development of SiC-based power semiconductor devices is the quality of the interface between SiC a… ▽ More SiC is set to enable a new era in power electronics impacting a wide range of energy technologies, from electric vehicles to renewable energy. Its physical characteristics outperform silicon in many aspects, including band gap, breakdown field, and thermal conductivity. The main challenge for further development of SiC-based power semiconductor devices is the quality of the interface between SiC and its native dielectric SiO$_2$. High temperature nitridation processes can improve the interface quality and ultimately the device performance immensely, but the underlying chemical processes are still poorly understood. Here, we present an energy-dependent hard X-ray photoelectron spectroscopy (HAXPES) study probing non-destructively SiC and SiO$_2$ and their interface in device stacks treated in varying atmospheres. We successfully combine laboratory- and synchrotron-based HAXPES to provide unique insights into the chemistry of interface defects and their passivation through nitridation processes. △ Less

Submitted 21 April, 2020; v1 submitted 13 December, 2019; originally announced December 2019.

Constraining the thermally-pulsing asymptotic giant branch phase with resolved stellar populations in the Small Magellanic Cloud

Authors: Giada Pastorelli, Paola Marigo, Léo Girardi, Yang Chen, Stefano Rubele, Michele Trabucchi, Bernhard Aringer, Sara Bladh, Alessandro Bressan, Josefina Montalbán, Martha L. Boyer, Julianne J. Dalcanton, Kjell Eriksson, Martin A. T. Groenewegen, Susanne Höfner, Thomas Lebzelter, Ambra Nanni, Philip Rosenfield, Peter R. Wood, Maria-Rosa L. Cioni

Abstract: The thermally-pulsing asymptotic giant branch (TP-AGB) experienced by low- and intermediate-mass stars is one of the most uncertain phases of stellar evolution and the models need to be calibrated with the aid of observations. To this purpose, we couple high-quality observations of resolved stars in the Small Magellanic Cloud (SMC) with detailed stellar population synthesis simulations computed wi… ▽ More The thermally-pulsing asymptotic giant branch (TP-AGB) experienced by low- and intermediate-mass stars is one of the most uncertain phases of stellar evolution and the models need to be calibrated with the aid of observations. To this purpose, we couple high-quality observations of resolved stars in the Small Magellanic Cloud (SMC) with detailed stellar population synthesis simulations computed with the TRILEGAL code. The strength of our approach relies on the detailed spatially-resolved star formation history of the SMC, derived from the deep near-infrared photometry of the VISTA survey of the Magellanic Clouds, as well as on the capability to quickly and accurately explore a wide variety of parameters and effects with the COLIBRI code for the TP-AGB evolution. Adopting a well-characterized set of observations -- star counts and luminosity functions -- we set up a calibration cycle along which we iteratively change a few key parameters of the TP-AGB models until we eventually reach a good fit to the observations. Our work leads to identify two best-fitting models that mainly differ in the efficiencies of the third dredge-up and mass loss in TP-AGB stars with initial masses larger than about 3 M$_{\odot}$. On the basis of these calibrated models we provide a full characterization of the TP-AGB stellar population in the SMC in terms of stellar parameters (initial masses, C/O ratios, carbon excess, mass-loss rates). Extensive tables of isochrones including these improved models are publicly available. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: accepted for publication in MNRAS. 29 pages, 24 figures, 2 appendices. One online appendix available as ancillary file on this page

Carbon star wind models at solar and sub-solar metallicities: a comparative study. I. Mass loss and the properties of dust-driven winds

Authors: S. Bladh, K. Eriksson, P. Marigo, S. Liljegren, B. Aringer

Abstract: The atmospheres and winds of C-type AGB stars are modeled with the 1D spherically symmetric radiation-hydrodynamical code DARWIN. To explore the metallicity-dependence of mass loss we calculate model grids at three different chemical abundances. Since carbon may be dredged up during the thermal pulses as AGB stars evolve, we keep the carbon abundance as a free parameter. The models in these three… ▽ More The atmospheres and winds of C-type AGB stars are modeled with the 1D spherically symmetric radiation-hydrodynamical code DARWIN. To explore the metallicity-dependence of mass loss we calculate model grids at three different chemical abundances. Since carbon may be dredged up during the thermal pulses as AGB stars evolve, we keep the carbon abundance as a free parameter. The models in these three different grids all have a current mass of one solar mass; effective temperatures of 2600K, 2800K, 3000K, or 3200K; and stellar luminosities equal to log(L/Lsun)=3.70, 3.85, or 4.00. The models show that mass loss in carbon stars is facilitated by high luminosities, low effective temperatures, and a high carbon excess (C-O) at both solar and subsolar metallicities. Similar combinations of effective temperature, luminosity, and carbon excess produce outflows at both solar and subsolar metallicities. There are no large systematic differences in the mass-loss rates and wind velocities produced by these wind models with respect to metallicity, nor any systematic difference concerning the distribution of grain sizes or how much carbon is condensed into dust. DARWIN models at subsolar metallicity have approximately 15% lower mass-loss rates compared to DARWIN models at solar metallicity with the same stellar parameters and carbon excess. For both solar and subsolar environments typical grain sizes range between 0.1 and 0.5 micron, the degree of condensed carbon varies between 5% and 40%, and the gas-to-dust ratios between 500 and 10000. Furthermore, stellar evolution models can use the mass-loss rates calculated from DARWIN models at solar metallicity when modeling the AGB phase at subsolar metallicities if carbon excess is used as the critical abundance parameter instead of the C/O ratio. △ Less

Submitted 14 February, 2019; originally announced February 2019.

Comments: 14 pages, 13 figures, accepted to A&A 28/01/2019

Journal ref: A&A 623, A119 (2019)

Scattering theory of the bifurcation in quantum measurement

Authors: Karl-Erik Eriksson, Kristian Lindgren

Abstract: We model quantum measurement of a two-level system $μ$. Previous obstacles for understanding the measurement process are removed by basing the analysis of the interaction between $μ$ and the measurement device on quantum field theory. We show how microscopic details of the measurement device can influence the transition to a final state. A statistical analysis of the ensemble of initial states rev… ▽ More We model quantum measurement of a two-level system $μ$. Previous obstacles for understanding the measurement process are removed by basing the analysis of the interaction between $μ$ and the measurement device on quantum field theory. We show how microscopic details of the measurement device can influence the transition to a final state. A statistical analysis of the ensemble of initial states reveals that those initial states that are efficient in leading to a transition to a final state, result in either of the expected eigenstates for $μ$, with probabilities that agree with the Born rule. △ Less

Submitted 14 May, 2019; v1 submitted 4 January, 2019; originally announced January 2019.

Comments: Change of title and minor revisions of main text and supplemental material; main text 8 pages, 2 figures; supplemental material 9 pages, 2 figures

On the Interplay between Robustness and Dynamic Planning for Adaptive Radiation Therapy

Authors: Michelle Böck, Kjell Eriksson, Anders Forsgren

Abstract: Interfractional geometric uncertainties can lead to deviations of the actual delivered dose from the prescribed dose distribution. To better handle these uncertainties during treatment, the authors propose a dynamic framework for robust adaptive radiation therapy in which a variety of robust adaptive treatment strategies are introduced and evaluated. This variety is a result of optimization variab… ▽ More Interfractional geometric uncertainties can lead to deviations of the actual delivered dose from the prescribed dose distribution. To better handle these uncertainties during treatment, the authors propose a dynamic framework for robust adaptive radiation therapy in which a variety of robust adaptive treatment strategies are introduced and evaluated. This variety is a result of optimization variables with various degrees of freedom within robust optimization models that vary in their grade of conservativeness. We formulate an adaptive framework with dependence on time and/or uncertainty or neither. The framework is evaluated with either expected-value-, worst-case-, or conditional-value-at-risk-(CVaR)-optimization. The goal of this study is to identify the mathematical properties of the proposed robust adaptive strategies in order to restore the planned dose distribution. We apply a result from convex analysis to show that in our framework, the time-dependent problem has a time-independent optimal solution. This implies that the time-dependent problem can be solved by dynamically solving the corresponding time-independent problem. This is especially relevant in a clinical context, since it implies that adaptive radiation therapy can be implemented in clinical practice using commercially available treatment planning systems. We demonstrate that time-and-scenario-dependent variables are most compatible with worst-case-optimization and that including dose prediction potentially improves treatment compared to non-adaptive planning. Contrarily, time-and-scenario-independent variables appear most compatible with expected-value-optimization. Non-adaptive robust strategies may provide sufficient target coverage, comparable to robust adaptive strategies, if the occurring uncertainties follow the same distribution as those included in the robust model. △ Less

Submitted 1 November, 2018; originally announced November 2018.

Journal ref: Biomedical Physics & Engineering Express 2019

Quantum Field Theory with Classical Sources - Linearized Quantum Gravity

Authors: Bo-Sture K. Skagerstam, Karl-Erik Eriksson, Per K. Rekdal

Abstract: In a previous work and in terms of an exact quantum-mechanical framework, $\hbar$-independent causal and retarded expectation values of the second-quantized electro-magnetic fields in the Coulomb gauge were derived in the presence of a conserved classical electric current. The classical $\hbar$-independent Maxwell's equations then naturally emerged. In the present work, we extend these considerati… ▽ More In a previous work and in terms of an exact quantum-mechanical framework, $\hbar$-independent causal and retarded expectation values of the second-quantized electro-magnetic fields in the Coulomb gauge were derived in the presence of a conserved classical electric current. The classical $\hbar$-independent Maxwell's equations then naturally emerged. In the present work, we extend these considerations to linear gravitational quantum deviations around a flat Minkowski space-time in a Coulomb-like gauge. The emergence of the classical causal and properly retarded linearized classical theory of general relativity with a conserved classical energy-momentum tensor is then outlined. The quantum-mechanical framework also provides for a simple approach to classical quadrupole gravitational radiation of Einstein and microscopic spontaneous graviton emission and/or absorption processes. △ Less

Submitted 19 June, 2018; originally announced June 2018.

Gaia Data Release 2: Observational Hertzsprung-Russell diagrams

Authors: Gaia Collaboration, C. Babusiaux, F. van Leeuwen, M. A. Barstow, C. Jordi, A. Vallenari, D. Bossini, A. Bressan, T. Cantat-Gaudin, M. van Leeuwen, A. G. A. Brown, T. Prusti, J. H. J. de Bruijne, C. A. L. Bailer-Jones, M. Biermann, D. W. Evans, L. Eyer, F. Jansen, S. A. Klioner, U. Lammers, L. Lindegren, X. Luri, F. Mignard, C. Panem, D. Pourbaix , et al. (428 additional authors not shown)

Abstract: We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the select… ▽ More We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies. △ Less

Submitted 13 August, 2018; v1 submitted 25 April, 2018; originally announced April 2018.

Comments: Published in the A&A Gaia Data Release 2 special issue. Tables 2 and A.4 corrected. Tables available at http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/616/A10

Journal ref: A&A 616, A10 (2018)

Increased accuracy of planning tools for optimization of dynamic multileaf collimator delivery of radiotherapy through reformulated objective functions

Authors: Lovisa Engberg, Kjell Eriksson, Anders Forsgren

Abstract: The purpose of this study is to examine in a clinical setting a novel formulation of objective functions for intensity-modulated radiotherapy treatment plan multicriteria optimization (MCO) that we suggested in a recent study. The proposed objective functions are extended with dynamic multileaf collimator (DMLC) delivery constraints from the literature, and a tailored interior point method is desc… ▽ More The purpose of this study is to examine in a clinical setting a novel formulation of objective functions for intensity-modulated radiotherapy treatment plan multicriteria optimization (MCO) that we suggested in a recent study. The proposed objective functions are extended with dynamic multileaf collimator (DMLC) delivery constraints from the literature, and a tailored interior point method is described to efficiently solve the resulting optimization formulation. In a numerical planning study involving three patient cases, DMLC plans Pareto optimal to the MCO formulation with the proposed objective functions are generated. Evaluated based on pre-defined plan quality indices, these DMLC plans are compared to conventionally generated DMLC plans. Comparable or superior plan quality is observed. Supported by these results, the proposed objective functions are argued to have a potential to streamline the planning process, since they are designed to overcome the methodological shortcomings associated with the conventional penalty-based objective functions assumed to cause the current need for time-consuming trial-and-error parameter tuning. In particular, the increased accuracy of the planning tools imposed by the proposed objective functions has the potential to make the planning process less complicated. These conclusions position the proposed formulation as an alternative to existing methods for automated planning. △ Less

Submitted 23 October, 2018; v1 submitted 2 February, 2018; originally announced February 2018.

Journal ref: Physics in Medicine & Biology 2018, volume 63, number 12, page 125012

Causality in Quantum Field Theory with Classical Sources - Quantum Electrodynamics

Authors: Bo-Sture K. Skagerstam, Karl-Erik Eriksson, Per K. Rekdal

Abstract: In an exact quantum-mechanical framework, we show that expectation values of the second-quantized electro-magnetic fields in the Coulomb gauge, and in the presence of classical sources, automatically lead to causal and retarded electro-magnetic field strengths. The classical $\hbar$-independent Maxwell's equations naturally emerge from this fundamental quantum-mechanical approach in terms of expec… ▽ More In an exact quantum-mechanical framework, we show that expectation values of the second-quantized electro-magnetic fields in the Coulomb gauge, and in the presence of classical sources, automatically lead to causal and retarded electro-magnetic field strengths. The classical $\hbar$-independent Maxwell's equations naturally emerge from this fundamental quantum-mechanical approach in terms of expectation values of quantum fields, and are therefore also consistent with the special theory of relativity. The fundamental difference between interference phenomena due to the linear nature of the classical Maxwell theory as, e.g., in classical optics, and interference effects of quantum states is clarified. The framework outlined also provides for a simple approach to, e.g., spontaneous photon emission and/or absorption processes as well as to the classical Vavilov-Cherenkov radiation. The inherent and necessary quantum fluctuations, limiting a precise space-time knowledge of expectation values of the quantum fields considered, are, finally, recalled. △ Less

Submitted 30 March, 2019; v1 submitted 30 January, 2018; originally announced January 2018.

Comments: 27 pages, one figure, 59 references. Updated references and correction of typograhical errors

Bifurcation in Quantum Measurement

Authors: Karl-Erik Eriksson, Martin Cederwall, Kristian Lindgren, Erik Sjöqvist

Abstract: We present a generic model of (non-destructive) quantum measurement. Being formulated within reversible quantum mechanics, the model illustrates a mechanism of a measurement process --- a transition of the measured system to an eigenstate of the measured observable. The model consists of a two-level system $μ$ interacting with a larger system $A$, consisting of smaller subsystems. The interaction… ▽ More We present a generic model of (non-destructive) quantum measurement. Being formulated within reversible quantum mechanics, the model illustrates a mechanism of a measurement process --- a transition of the measured system to an eigenstate of the measured observable. The model consists of a two-level system $μ$ interacting with a larger system $A$, consisting of smaller subsystems. The interaction is modelled as a scattering process. Restricting the states of $A$ to product states leads to a bifurcation process: In the limit of a large system $A$, the initial states of $A$ that are efficient in leading to a final state are divided into two separated subsets. For each of these subsets, $μ$ ends up in one of the eigenstates of the measured observable. The probabilities obtained in this branching confirm the Born rule. △ Less

Submitted 3 November, 2017; v1 submitted 4 August, 2017; originally announced August 2017.

Comments: A revised version that includes a more general presentation of the model (in Sect. 4) and a larger revision of the Introduction

Pulsation-induced atmospheric dynamics in M-type AGB stars. Effects on wind properties, photometric variations and near-IR CO line profiles

Authors: S. Liljegren, S. Höfner, K. Eriksson, W. Nowotny

Abstract: Wind-driving in asymptotic giant branch (AGB) stars is commonly attributed to a two-step process. First, matter in the stellar atmosphere is levitated by shock waves, induced by stellar pulsation, and second, this matter is accelerated by radiation pressure on dust, resulting in a wind. In dynamical atmosphere and wind models the effects of the stellar pulsation are often simulated by a simplistic… ▽ More Wind-driving in asymptotic giant branch (AGB) stars is commonly attributed to a two-step process. First, matter in the stellar atmosphere is levitated by shock waves, induced by stellar pulsation, and second, this matter is accelerated by radiation pressure on dust, resulting in a wind. In dynamical atmosphere and wind models the effects of the stellar pulsation are often simulated by a simplistic prescription at the inner boundary. We test a sample of dynamical models for M-type AGB stars, for which we kept the stellar parameters fixed to values characteristic of a typical Mira variable but varied the inner boundary condition. The aim was to evaluate the effect on the resulting atmosphere structure and wind properties. The results of the models are compared to observed mass-loss rates and wind velocities, photometry, and radial velocity curves, and to results from 1D radial pulsation models. Dynamical atmosphere models are calculated, using the DARWIN code for different combinations of photospheric velocities and luminosity variations. The inner boundary is changed by introducing an offset between maximum expansion of the stellar surface and the luminosity and/or by using an asymmetric shape for the luminosity variation. Models that resulted in realistic wind velocities and mass-loss rates, when compared to observations, also produced realistic photometric variations. For the models to also reproduce the characteristic radial velocity curve present in Mira stars (derived from CO $Δv = 3$ lines), an overall phase shift of 0.2 between the maxima of the luminosity and radial variation had to be introduced. We find that a group of models with different boundary conditions (29 models, including the model with standard boundary conditions) results in realistic velocities and mass-loss rates, and in photometric variations. △ Less

Submitted 26 June, 2017; originally announced June 2017.

Journal ref: A&A 606, A6 (2017)

An exponential limit shape of random $q$-proportion Bulgarian solitaire

Authors: Kimmo Eriksson, Markus Jonsson abd Jonas Sjöstrand

Abstract: We introduce \emph{$p_n$-random $q_n$-proportion Bulgarian solitaire} ($0<p_n,q_n\le 1$), played on $n$ cards distributed in piles. In each pile, a number of cards equal to the proportion $q_n$ of the pile size rounded upward to the nearest integer are candidates to be picked. Each candidate card is picked with probability $p_n$, independently of other candidate cards. This generalizes Popov's ran… ▽ More We introduce \emph{$p_n$-random $q_n$-proportion Bulgarian solitaire} ($0<p_n,q_n\le 1$), played on $n$ cards distributed in piles. In each pile, a number of cards equal to the proportion $q_n$ of the pile size rounded upward to the nearest integer are candidates to be picked. Each candidate card is picked with probability $p_n$, independently of other candidate cards. This generalizes Popov's random Bulgarian solitaire, in which there is a single candidate card in each pile. Popov showed that a triangular limit shape is obtained for a fixed $p$ as $n$ tends to infinity. Here we let both $p_n$ and $q_n$ vary with $n$. We show that under the conditions $q_n^2 p_n n/{\log n}\rightarrow \infty$ and $p_n q_n \rightarrow 0$ as $n\to\infty$, the $p_n$-random $q_n$-proportion Bulgarian solitaire has an exponential limit shape. △ Less

Submitted 21 March, 2017; originally announced March 2017.

Limit shapes of stable configurations of a generalized Bulgarian solitaire

Authors: Kimmo Eriksson, Markus Jonsson, Jonas Sjöstrand

Abstract: Bulgarian solitaire is played on $n$ cards divided into several piles; a move consists of picking one card from each pile to form a new pile. In a recent generalization, $σ$-Bulgarian solitaire, the number of cards you pick from a pile is some function $σ$ of the pile size, such that you pick $σ(h)\le h$ cards from a pile of size $h$. Here we consider a special class of such functions. Let us call… ▽ More Bulgarian solitaire is played on $n$ cards divided into several piles; a move consists of picking one card from each pile to form a new pile. In a recent generalization, $σ$-Bulgarian solitaire, the number of cards you pick from a pile is some function $σ$ of the pile size, such that you pick $σ(h)\le h$ cards from a pile of size $h$. Here we consider a special class of such functions. Let us call $σ$ well-behaved if $σ(1)=1$ and if both $σ(h)$ and $h-σ(h)$ are non-decreasing functions of $h$. Well-behaved $σ$-Bulgarian solitaire has a geometric interpretation in terms of layers at certain levels being picked in each move. It also satisfies that if a stable configuration of $n$ cards exists it is unique. Moreover, if piles are sorted in order of decreasing size ($λ_1 \ge λ_2\ge \dots$) then a configuration is convex if and only if it is a stable configuration of some well-behaved $σ$-Bulgarian solitaire. If sorted configurations are represented by Young diagrams and scaled down to have unit height and unit area, the stable configurations corresponding to an infinite sequence of well-behaved functions ($σ_1, σ_2, \dots$) may tend to a limit shape $φ$. We show that every convex $φ$ with certain properties can arise as the limit shape of some sequence of well-behaved $σ_n$. For the special case when $σ_n(h)=\lceil q_n h \rceil$ for $0 < q_n \le 1$, these limit shapes are triangular (in case $q_n^2 n\rightarrow 0$), or exponential (in case $q_n^2 n\rightarrow \infty$), or interpolating between these shapes (in case $q_n^2 n\rightarrow C>0$). △ Less

Submitted 21 March, 2017; originally announced March 2017.

Aperture synthesis imaging of the carbon AGB star R Sculptoris: Detection of a complex structure and a dominating spot on the stellar disk

Authors: M. Wittkowski, K. -H. Hofmann, S. Höfner, J. B. Le Bouquin, W. Nowotny, C. Paladini, J. Young, J. -P. Berger, M. Brunner, I. de Gregorio-Monsalvo, K. Eriksson, J. Hron, E. M. L. Humphreys, M. Lindqvist, M. Maercker, S. Mohamed, H. Olofsson, S. Ramstedt, G. Weigelt

Abstract: We present near-infrared interferometry of the carbon-rich asymptotic giant branch (AGB) star R Sculptoris. The visibility data indicate a broadly circular resolved stellar disk with a complex substructure. The observed AMBER squared visibility values show drops at the positions of CO and CN bands, indicating that these lines form in extended layers above the photosphere. The AMBER visibility va… ▽ More We present near-infrared interferometry of the carbon-rich asymptotic giant branch (AGB) star R Sculptoris. The visibility data indicate a broadly circular resolved stellar disk with a complex substructure. The observed AMBER squared visibility values show drops at the positions of CO and CN bands, indicating that these lines form in extended layers above the photosphere. The AMBER visibility values are best fit by a model without a wind. The PIONIER data are consistent with the same model. We obtain a Rosseland angular diameter of 8.9+-0.3 mas, corresponding to a Rosseland radius of 355+-55 Rsun, an effective temperature of 2640+-80 K, and a luminosity of log L/Lsun=3.74+-0.18. These parameters match evolutionary tracks of initial mass 1.5+-0.5 Msun and current mass 1.3+-0.7 Msun. The reconstructed PIONIER images exhibit a complex structure within the stellar disk including a dominant bright spot located at the western part of the stellar disk. The spot has an H-band peak intensity of 40% to 60% above the average intensity of the limb-darkening-corrected stellar disk. The contrast between the minimum and maximum intensity on the stellar disk is about 1:2.5. Our observations are broadly consistent with predictions by dynamic atmosphere and wind models, although models with wind appear to have a circumstellar envelope that is too extended compared to our observations. The detected complex structure within the stellar disk is most likely caused by giant convection cells, resulting in large-scale shock fronts, and their effects on clumpy molecule and dust formation seen against the photosphere at distances of 2-3 stellar radii. △ Less

Submitted 8 February, 2017; originally announced February 2017.

Comments: 13 pages, 11 figures, 4 pages appendix, accepted for publication in Astronomy and Astrophysics (A&A)

Journal ref: A&A 601, A3 (2017)

The adventure of carbon stars - Observations and modelling of a set of C-rich AGB stars

Authors: Gioia Rau, Josef Hron, Claudia Paladini, Benard Aringer, Kjell Eriksson, Paola Marigo, Walter Nowotny, Rebekka Grellmann

Abstract: Modeling stellar atmospheres is a complex, intriguing task in modern astronomy. A systematic comparison of models with multi-technique observations is the only efficient way to constrain them. Aims: We performed a self-consistent modeling of the atmospheres of six C-rich AGB stars: R Lep, R Vol, Y Pav, AQ Sgr, U Hya and X TrA, with the aim of enlarging the knowledge of the dynamic processes occurr… ▽ More Modeling stellar atmospheres is a complex, intriguing task in modern astronomy. A systematic comparison of models with multi-technique observations is the only efficient way to constrain them. Aims: We performed a self-consistent modeling of the atmospheres of six C-rich AGB stars: R Lep, R Vol, Y Pav, AQ Sgr, U Hya and X TrA, with the aim of enlarging the knowledge of the dynamic processes occurring in their atmospheres. Methods. We used VLTI/MIDI interferometric observations, in combination with spectro-photometric data, and compared them with self-consistent dynamic models atmospheres. Results: The models can reproduce SED data well at wavelengths longwards of 1μm, and the interferometric observations between 8μm and 10μm. Discrepancies could be due to a combination of data- and model-related effects. The models best fitting the Miras are significantly extended, with a prominent shell-like structure, while the models best fitting the non-Miras are more compact, showing lower average mass-loss. The mass loss is of episodic or multi-periodic nature, but causes the visual amplitudes to be notably larger than the observed ones. Stellar parameters were derived from the model fitting: T_eff, L_bol, M, C/O, $\dot{\text{M}}$. Our findings agree well with literature values within the uncertainties. T_eff and L_bol are also in good agreement with the T derived from the angular diameter θ_{(V-K)} and L_bol from the SED fitting, except for AQ Sgr. Finally, θ_Ross and θ_(V-K) agree with each other better for the Miras than for the non-Miras,which is probably connected to the episodic nature of the latter models. We also located the stars in the H-R diagram, comparing them with evolutionary tracks. We found that the main derived properties from the model fitting are in good agreement with TP-AGB evolutionary calculations for carbon stars (COLIBRI code). △ Less

Submitted 16 January, 2017; originally announced January 2017.

Comments: 22 pages. Accepted for publication in A&A

Journal ref: A&A 600, A92 (2017)

Abundances of disk and bulge giants from hi-res optical spectra: I. O, Mg, Ca, and Ti in the Solar neighborhood and Kepler field samples

Authors: H. Jönsson, N. Ryde, T. Nordlander, A. Pehlivan, H. Hartman, P. Jönsson, K. Eriksson

Abstract: The galactic bulge is a significant part of our galaxy, but it is hard to observe, being both distant and covered by dust in the disk. Therefore there do not exist many hi-res optical spectra of bulge stars with large wavelength coverage, whose determined abundances can be compared with nearby, similarly analyzed stellar samples. We aim to determine the, for chemical evolution models, so important… ▽ More The galactic bulge is a significant part of our galaxy, but it is hard to observe, being both distant and covered by dust in the disk. Therefore there do not exist many hi-res optical spectra of bulge stars with large wavelength coverage, whose determined abundances can be compared with nearby, similarly analyzed stellar samples. We aim to determine the, for chemical evolution models, so important alpha elements of a sample of bulge giants using hi-res optical spectra with large wavelength coverage. The abundances found will be compared to similarly derived abundances from similar spectra of similar stars in the local thin and thick disks. In this first paper we focus on the Solar neighborhood reference sample. We use spectral synthesis to derive the stellar parameters as well as the elemental abundances of both the local as well as the bulge samples of giants. Special care is taken to benchmark our method of determining stellar parameters against independent measurements of effective temperatures from angular diameter measurements and surface gravities from asteroseismology. In this first paper we present the method used to determine the stellar parameters as well as the elemental abundances, evaluate them, and present the results for our local disk sample of 291 giants. When comparing our determined spectroscopic temperatures to those derived from angular diameter measurements, we reproduce these with a systematic difference of +10 K and a standard deviation of 53 K. The spectroscopic gravities are reproducing the ones determined from asteroseismology with a systematic offset of +0.10 dex and a standard deviation of 0.12 dex. When it comes to the abundance trends, our sample of local disk giants is closely following that of other works analyzing solar neighborhood dwarfs, showing that the much brighter giant stars are as good abundance probes as the often used dwarfs. △ Less

Submitted 16 November, 2016; originally announced November 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 598, A100 (2017)

Towards Robust Adaptive Radiation Therapy Strategies

Authors: Michelle Böck, Anders Forsgren, Kjell Eriksson, Björn Hårdemark

Abstract: The authors propose robust adaptive strategies based on stochastic minimax optimization for a series of simulated treatments on a one-dimensional patient phantom. The plan applied during the first fractions should be able to handle anticipated systematic and random errors. At scheduled fractions, the impact of the measured errors on the delivered dose distribution is evaluated. For a patient recei… ▽ More The authors propose robust adaptive strategies based on stochastic minimax optimization for a series of simulated treatments on a one-dimensional patient phantom. The plan applied during the first fractions should be able to handle anticipated systematic and random errors. At scheduled fractions, the impact of the measured errors on the delivered dose distribution is evaluated. For a patient receiving a dose that does not satisfy clinical goals, the plan is reoptimized based on these individually measured errors. The adapted plan is then applied during subsequent fractions until a new scheduled adaptation becomes necessary. In the first adaptive strategy, the measured systematic and random error scenarios and their assigned probabilities are updated to guide the robust reoptimization. In the second strategy, the grade of conservativeness is adapted in response to the measured dose delivery errors. In the third strategy, the uncertainty margins around the target are recalculated based on the measured errors. The simulated systematic and random errors are either similar to the anticipated errors or unpredictably larger, in order to critically evaluate the performance of the adaptive strategies. The simulations show that robustly optimized treatment plans provide sufficient treatment quality for those treatment error scenarios similar to the anticipated error scenarios. Combining robust planning with adaptation leads to improved OAR protection. In case of unpredictably larger treatment errors, the first strategy performs best at improving treatment quality in terms of target coverage and OAR protection. According to the simulations, these robust adaptive treatment strategies are able to identify candidates for an adaptive treatment, thus giving the opportunity to provide individualized plans, and improve their treatment quality through adaptation. △ Less

Submitted 4 November, 2016; originally announced November 2016.

Journal ref: Medical Physics 2017

Explicit optimization of plan quality measures in intensity-modulated radiation therapy treatment planning

Authors: Lovisa Engberg, Kjell Eriksson, Anders Forsgren, Björn Hårdemark

Abstract: Conventional planning objectives in optimization of intensity-modulated radiotherapy treatment (IMRT) plans are designed to minimize the violation of dose-volume histogram (DVH) thresholds using penalty functions. Although successful in guiding the DVH curve towards these thresholds, conventional planning objectives offer limited control of the individual points on the DVH curve (doses-at-volume)… ▽ More Conventional planning objectives in optimization of intensity-modulated radiotherapy treatment (IMRT) plans are designed to minimize the violation of dose-volume histogram (DVH) thresholds using penalty functions. Although successful in guiding the DVH curve towards these thresholds, conventional planning objectives offer limited control of the individual points on the DVH curve (doses-at-volume) used to evaluate plan quality. In this study, we abandon the usual penalty-function framework and propose planning objectives that more explicitly relate to DVH statistics. The proposed planning objectives are based on mean-tail-dose, resulting in convex optimization. We also demonstrate how to adapt a standard optimization method to the proposed formulation in order to obtain a substantial reduction in computational cost. We investigate the potential of the proposed planning objectives as tools for optimizing DVH statistics through juxtaposition with the conventional planning objectives on two patient cases. Sets of treatment plans with differently balanced planning objectives are generated using either the proposed or the conventional approach. Dominance in the sense of better distributed doses-at-volume is observed in plans optimized within the proposed framework, indicating that the DVH statistics are better optimized and more efficiently balanced using the proposed planning objectives. △ Less

Submitted 25 January, 2017; v1 submitted 5 July, 2016; originally announced July 2016.

Journal ref: Medical Physics 2017, volume 44, pages 2045-2053

Dust-driven winds of AGB stars: The critical interplay of atmospheric shocks and luminosity variations

Authors: S. Liljegren, S. Höfner, W. Nowotny, K. Eriksson

Abstract: Winds of AGB stars are thought to be driven by a combination of pulsation-induced shock waves and radiation pressure on dust. In dynamic atmosphere and wind models, the stellar pulsation is often simulated by prescribing a simple sinusoidal variation in velocity and luminosity at the inner boundary of the model atmosphere. We experiment with different forms of the luminosity variation in order to… ▽ More Winds of AGB stars are thought to be driven by a combination of pulsation-induced shock waves and radiation pressure on dust. In dynamic atmosphere and wind models, the stellar pulsation is often simulated by prescribing a simple sinusoidal variation in velocity and luminosity at the inner boundary of the model atmosphere. We experiment with different forms of the luminosity variation in order to assess the effects on the wind velocity and mass-loss rate, when progressing from the simple sinusoidal recipe towards more realistic descriptions. Using state-of-the-art dynamical models of C-rich AGB stars, a range of different asymmetric shapes of the luminosity variation and a range of phase shifts of the luminosity variation relative to the radial variation are tested. These tests are performed on two stellar atmosphere models. The first model has dust condensation and, as a consequence, a stellar wind is triggered, while the second model lacks both dust and wind. The first model with dust and stellar wind is very sensitive to moderate changes in the luminosity variation. There is a complex relationship between the luminosity minimum, and dust condensation: changing the phase corresponding to minimum luminosity can either increase or decrease mass-loss rate and wind velocity. The luminosity maximum dominates the radiative pressure on the dust, which in turn, is important for driving the wind. These effects of changed luminosity variation are coupled with the dust formation. In contrast there is very little change to the structure of the model without dust. Changing the luminosity variation, both by introducing a phase shift and by modifying the shape, influences wind velocity and the mass-loss rate. To improve wind models it would probably be desirable to extract boundary conditions from 3D dynamical interior models or stellar pulsation models. △ Less

Submitted 22 March, 2016; originally announced March 2016.

Comments: 11 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 589, A130 (2016)

Markov chains on graded posets: Compatibility of up-directed and down-directed transition probabilities

Authors: Kimmo Eriksson, Markus Jonsson, Jonas Sjöstrand

Abstract: We consider two types of discrete-time Markov chains where the state space is a graded poset and the transitions are taken along the covering relations in the poset. The first type of Markov chain goes only in one direction, either up or down in the poset (an \emph{up chain} or \emph{down chain}). The second type toggles between two adjacent rank levels (an \emph{up-and-down chain}). We introduc… ▽ More We consider two types of discrete-time Markov chains where the state space is a graded poset and the transitions are taken along the covering relations in the poset. The first type of Markov chain goes only in one direction, either up or down in the poset (an \emph{up chain} or \emph{down chain}). The second type toggles between two adjacent rank levels (an \emph{up-and-down chain}). We introduce two compatibility concepts between the up-directed transition probabilities (an \emph{up rule}) and the down-directed (a \emph{down rule}), and we relate these to compatibility between up-and-down chains. This framework is used to prove a conjecture about a limit shape for a process on Young's lattice. Finally, we settle the questions whether the reverse of an up chain is a down chain for some down rule and whether there exists an up or down chain at all if the rank function is not bounded. △ Less

Submitted 11 December, 2015; originally announced December 2015.

Exploring wind-driving dust species in cool luminous giants III. Wind models for M-type AGB stars: dynamic and photometric properties

Authors: S. Bladh, S. Höfner, B. Aringer, K. Eriksson

Abstract: Stellar winds observed in asymptotic giant branch (AGB) stars are usually attributed to a combination of stellar pulsations and radiation pressure on dust. Shock waves triggered by pulsations propagate through the atmosphere, compressing the gas and lifting it to cooler regions, which create favourable conditions for grain growth. If sufficient radiative acceleration is exerted on the newly formed… ▽ More Stellar winds observed in asymptotic giant branch (AGB) stars are usually attributed to a combination of stellar pulsations and radiation pressure on dust. Shock waves triggered by pulsations propagate through the atmosphere, compressing the gas and lifting it to cooler regions, which create favourable conditions for grain growth. If sufficient radiative acceleration is exerted on the newly formed grains through absorption or scattering of stellar photons, an outflow can be triggered. Strong candidates for wind-driving dust species in M-type AGB stars are magnesium silicates (Mg$_2$SiO$_4$ and MgSiO$_3$). Such grains can form close to the stellar surface, they consist of abundant materials and, if they grow to sizes comparable to the wavelength of the stellar flux maximum, they experience strong acceleration by photon scattering. We use a frequency-dependent radiation-hydrodynamics code with a detailed description for the growth of Mg$_2$SiO$_4$ grains to calculate the first extensive set of time-dependent wind models for M-type AGB stars. The resulting wind properties, visual and near-IR photometry and mid-IR spectra are compared with observations.We show that the models can produce outflows for a wide range of stellar parameters. We also demonstrate that they reproduce observed mass-loss rates and wind velocities, as well as visual and near-IR photometry. However, the current models do not show the characteristic silicate features at 10 and 18 $μ$m as a result of the cool temperature of Mg$_2$SiO$_4$ grains in the wind. Including a small amount of Fe in the grains further out in the circumstellar envelope will increase the grain temperature and result in pronounced silicate features, without significantly affecting the photometry in the visual and near-IR wavelength regions. △ Less

Submitted 30 January, 2015; originally announced February 2015.

Comments: 11 pages, 14 figures

Systematic trend of water vapour absorption in red giant atmospheres revealed by high resolution TEXES 12 micron spectra

Authors: N. Ryde, J. Lambert, M. Farzone, M. J. Richter, E. Josselin, G. M. Harper, K. Eriksson, T. K. Greathouse

Abstract: The structures of the outer atmospheres of red giants are very complex. The notion of large optically thick molecular spheres around the stars (MOLspheres) has been invoked in order to explain e.g. spectro-interferometric observations. However, high-resolution spectra in the mid-IR do not easily fit into this picture. They rule out any large sphere of water vapour in LTE surrounding red giants. Ou… ▽ More The structures of the outer atmospheres of red giants are very complex. The notion of large optically thick molecular spheres around the stars (MOLspheres) has been invoked in order to explain e.g. spectro-interferometric observations. However, high-resolution spectra in the mid-IR do not easily fit into this picture. They rule out any large sphere of water vapour in LTE surrounding red giants. Our aim here is to investigate high-resolution, mid-infrared spectra for a range of red giants, from early-K to mid M. We have recorded 12 microns spectra of 10 well-studied bright red giants, with TEXES on the IRTF. We find that all giants in our study cooler than 4300 K, spanning a range of effective temperatures, show water absorption lines stronger than expected. The strengths of the lines vary smoothly with spectral type. We identify several spectral features in the wavelength region that undoubtedly are formed in the photosphere. From a study of water-line ratios of the stars, we find that the excitation temperatures, in the line-forming regions, are several hundred Kelvin lower than expected from a classical photospheric model. This could either be due to an actually lower temperature structure in the outer regions of the photospheres caused by, for example, extra cooling, or due to non-LTE level populations, affecting the source function and line opacities. We have demonstrated that these diagnostically interesting water lines are a general feature of red giants across spectral types, and we argue for a general explanation of their formation rather than explanations requiring specific properties. Since the water lines are neither weak (filled in by emission) nor appear in emission, as predicted by LTE MOLsphere models in their simplest forms, the evidence for the existence of such large optically-thick, molecular spheres enshrouding the stars is weakened. (abbreviated) △ Less

Submitted 15 October, 2014; originally announced October 2014.

Comments: Accepted for publication in A&A

Journal ref: A&A 573, A28 (2015)

Evolved stars with complex atmospheres - the high spectral resolution, mid-IR view

Authors: N. Ryde, J. Lambert, M. J. Richter, E. Josselin, G. M. Harper, K. Eriksson, A. Boogert, C. DeWitt, T. Encrenaz, T. Greathouse, D. Jaffe, K. Kulas, M. McKelvey, J. Najita, W. Vacca

Abstract: The physical structures of the outer atmospheres of red giants are not known. They are certainly complex and a range of recent observations are showing that we need to embrace to non-classical atmosphere models to interpret these regions. This region's properties is of importance, not the least, for the understanding of the mass-loss mechanism for these stars, which is not still understood. Here,… ▽ More The physical structures of the outer atmospheres of red giants are not known. They are certainly complex and a range of recent observations are showing that we need to embrace to non-classical atmosphere models to interpret these regions. This region's properties is of importance, not the least, for the understanding of the mass-loss mechanism for these stars, which is not still understood. Here, we present observational constraints of the outer regions of red giants, based on mid-IR, high spectral resolution spectra. We also discuss possible non-LTE effects and highlight a new non-LTE code that will be used to analyse the spectra of these atmospheric layers. We conclude by mentioning our new SOFIA/EXES observations of red giants at 6 microns, where the vibration-rotation lines of water vapour can be detected and spectrally resolved for the first time. △ Less

Submitted 26 August, 2014; originally announced August 2014.

Comments: Talk presented at the "Why Galaxies Care About AGB Stars III" conference in Vienna, July 2014. It will be published in the Astronomical Society of the Pacific Conference Series

The Hamburg/ESO R-process Enhanced Star survey (HERES) X. HE 2252-4225, one more r-process enhanced and actinide-boost halo star

Authors: L. Mashonkina, N. Christlieb, K. Eriksson

Abstract: We report on a detailed abundance analysis of the r-process enhanced giant star, HE 2252-4225 ([Fe/H] = -2.63, [r/Fe] = 0.80). Determination of stellar parameters and element abundances was based on analysis of high-quality VLT/UVES spectra. The surface gravity was calculated from the NLTE ionisation balance between Fe I and Fe II. Accurate abundances were determined for a total of 38 elements, in… ▽ More We report on a detailed abundance analysis of the r-process enhanced giant star, HE 2252-4225 ([Fe/H] = -2.63, [r/Fe] = 0.80). Determination of stellar parameters and element abundances was based on analysis of high-quality VLT/UVES spectra. The surface gravity was calculated from the NLTE ionisation balance between Fe I and Fe II. Accurate abundances were determined for a total of 38 elements, including 22 neutron-capture elements beyond Sr and up to Th. This object is deficient in carbon, as expected for a giant star with Teff < 4800 K. The stellar Na-Zn abundances are well fitted by the yields of a single supernova of 14.4 Msun. For the neutron-capture elements in the Sr-Ru, Ba-Yb, and Os-Ir regions, the abundance pattern of HE 2252-4225 is in excellent agreement with the average abundance pattern of the strongly r-process enhanced stars CS 22892-052, CS 31082-001, HE 1219-0312, and HE 1523-091. This suggests a common origin of the first, second, and third r-process peak elements in HE 2252-4225 in the classical r-process. We tested the solar r-process pattern based on the most recent s-process calculations of Bisterzo et al. (2014) and found that elements in the range from Ba to Ir match it very well. No firm conclusion can be drawn about the relationship between the fisrt neutron-capture peak elements, Sr to Ru, in HE 2252-4225 and the solar r-process, due to the uncertainty in the solar r-process. The investigated star has an anomalously high Th/Eu abundance ratio, so that radioactive age dating results in a stellar age of tau = 1.5+-1.5 Gyr that is not expected for a very metal-poor halo star. △ Less

Submitted 21 July, 2014; originally announced July 2014.

Comments: 20 pages, 6 tables, 9 figures, accepted for publication in A&A

Journal ref: A&A 569, A43 (2014)

Synthetic photometry for carbon-rich giants. IV. An extensive grid of dynamic atmosphere and wind models

Authors: Kjell Eriksson, Walter Nowotny, Susanne Höfner, Bernhard Aringer, Astrid Wachter

Abstract: The evolution and spectral properties of stars on the AGB are significantly affected by mass loss through dusty stellar winds. Dynamic atmosphere and wind models are an essential tool for studying these evolved stars, both individually and as members of stellar populations, to understand their contribution to the integrated light and chemical evolution of galaxies. This paper is part of a series… ▽ More The evolution and spectral properties of stars on the AGB are significantly affected by mass loss through dusty stellar winds. Dynamic atmosphere and wind models are an essential tool for studying these evolved stars, both individually and as members of stellar populations, to understand their contribution to the integrated light and chemical evolution of galaxies. This paper is part of a series testing state-of-the-art atmosphere and wind models of carbon stars against observations, and making them available for use in various theoretical and observational studies. We have computed low-resolution spectra and photometry (in the wavelength range 0.35-25 mu) for a grid of 540 dynamic models with stellar parameters typical of solar-metallicity C-rich AGB stars and with a range of pulsation amplitudes. The models cover the dynamic atmosphere and dusty outflow (if present), assuming spherical symmetry, and taking opacities of gas-phase species and dust grains consistently into account. To characterize the time-dependent dynamic and photometric behaviour of the models in a concise way we defined a number of classes for models with and without winds. Comparisons with observed data in general show a quite good agreement for example regarding mass-loss rates vs. (J-K) colours or K magnitudes vs. (J-K) colours. Some exceptions from the good overall agreement, however, are found and attributed to the range of input parameters (e.g. relatively high carbon excesses) or intrinsic model assumptions (e.g. small particle limit for grain opacities). While current results indicate that some changes in model assumptions and parameter ranges should be made in the future to bring certain synthetic observables into better agreement with observations, it seems unlikely that these pending improvements will significantly affect the mass-loss rates of the models. △ Less

Submitted 29 April, 2014; originally announced April 2014.

Comments: 28 pages, 15 figures. Table B.1, an 11-page table, is only available at CDS

The Chemical Evolution of Fluorine in the Bulge - High-resolution K-band spectra of giants in three fields

Authors: H. Jönsson, N. Ryde, G. M. Harper, K. Cunha, M. Schultheis, K. Eriksson, C. Kobayashi, V. V. Smith, M. Zoccali

Abstract: Possible main formation sites of F in the Universe include AGB stars, the ν-process in Type II SNe, and/or W-R stars. The importance of the W-R stars has theoretically been questioned and they are probably not needed in the modelling of the chemical evolution of F in the solar neighborhood. It has, however, been suggested that W-R stars are indeed needed to explain the chemical evolution of F in t… ▽ More Possible main formation sites of F in the Universe include AGB stars, the ν-process in Type II SNe, and/or W-R stars. The importance of the W-R stars has theoretically been questioned and they are probably not needed in the modelling of the chemical evolution of F in the solar neighborhood. It has, however, been suggested that W-R stars are indeed needed to explain the chemical evolution of F in the Bulge. The molecular spectral data of the often used HF-molecule has not been presented in a complete and consistent way and has recently been debated in the literature. In this article we determine the [F/O] vs. [O/H] trend in the Bulge to investigate the possible contribution from W-R stars. Additionally, we present here a HF line list for the K- and L-bands (including the often used 23358.33 Å line) and an accompanying partition function. The F abundances were determined using spectral fitting from hi-res NIR spectra of eight K giants recorded by the spectrograph CRIRES. We have also re-analyzed five previously published Bulge giants using our new HF molecular data. We find that the F-O abundance in the Bulge probably cannot be explained with chemical evolution models including only AGB-stars and the ν-process in SNe Type II, i.e. a significant amount of F production in W-R stars is likely needed to explain the F abundance in the Bulge. Concerning the HF line list, we find that a possible reason for the inconsistencies in the literature, with two different excitation energies being used, is two different definitions of the zero-point energy for the HF molecule and therefore also two accompanying different dissociation energies. Both line lists are correct, as long as the corresponding consistent partition function is used in the spectral synthesis. However, we suspect this has not been the case in several earlier works leading to F abundances 0.3 dex too high. △ Less

Submitted 11 March, 2014; originally announced March 2014.

The Limit Shape of a Stochastic Bulgarian Solitaire

Authors: Kimmo Eriksson, Markus Jonsson, Jonas Sjöstrand

Abstract: We consider a stochastic version of Bulgarian solitaire: A number of cards are distributed in piles; in every round a new pile is formed by cards from the old piles, and each card is picked independently with a fixed probability. This game corresponds to a multi-square birth-and-death process on Young diagrams of integer partitions. We prove that this process converges in a strong sense to an expo… ▽ More We consider a stochastic version of Bulgarian solitaire: A number of cards are distributed in piles; in every round a new pile is formed by cards from the old piles, and each card is picked independently with a fixed probability. This game corresponds to a multi-square birth-and-death process on Young diagrams of integer partitions. We prove that this process converges in a strong sense to an exponential limit shape as the number of cards tends to infinity. Furthermore, we bound the probability of deviation from the limit shape and relate this to the number of rounds played in the solitaire. △ Less

Submitted 11 December, 2015; v1 submitted 10 September, 2013; originally announced September 2013.

Comments: Splitting up Theorem 1 into Theorem 1 and Theorem 2 to clarify the average shape result. Results unchanged

MSC Class: 60J10; 05A17

Fundamental parameters of 16 late-type stars derived from their angular diameter measured with VLTI/AMBER

Authors: P. Cruzalèbes, A. Jorissen, Y. Rabbia, S. Sacuto, A. Chiavassa, E. Pasquato, B. Plez, K. Eriksson, A. Spang, O. Chesneau

Abstract: Thanks to their large angular dimension and brightness, red giants and supergiants are privileged targets for optical long-baseline interferometers. Sixteen red giants and supergiants have been observed with the VLTI/AMBER facility over a two-years period, at medium spectral resolution (R=1500) in the K band. The limb-darkened angular diameters are derived from fits of stellar atmospheric models o… ▽ More Thanks to their large angular dimension and brightness, red giants and supergiants are privileged targets for optical long-baseline interferometers. Sixteen red giants and supergiants have been observed with the VLTI/AMBER facility over a two-years period, at medium spectral resolution (R=1500) in the K band. The limb-darkened angular diameters are derived from fits of stellar atmospheric models on the visibility and the triple product data. The angular diameters do not show any significant temporal variation, except for one target: TX Psc, which shows a variation of 4% using visibility data. For the eight targets previously measured by Long-Baseline Interferometry (LBI) in the same spectral range, the difference between our diameters and the literature values is less than 5%, except for TX Psc, which shows a difference of 11%. For the 8 other targets, the present angular diameters are the first measured from LBI. Angular diameters are then used to determine several fundamental stellar parameters, and to locate these targets in the Hertzsprung-Russell Diagram (HRD). Except for the enigmatic Tc-poor low-mass carbon star W Ori, the location of Tc-rich stars in the HRD matches remarkably well the thermally-pulsating AGB, as it is predicted by the stellar-evolution models. For pulsating stars with periods available, we compute the pulsation constant and locate the stars along the various sequences in the Period -- Luminosity diagram. We confirm the increase in mass along the pulsation sequences, as predicted by the theory, except for W Ori which, despite being less massive, appears to have a longer period than T Cet along the first-overtone sequence. △ Less

Submitted 13 June, 2013; originally announced June 2013.

Comments: 15 pages, 9 figures, 6 tables

Synthetic photometry for carbon-rich giants. III. Tracing the sequence of mass-losing galactic C-type Miras

Authors: W. Nowotny, B. Aringer, S. Hoefner, K. Eriksson

Abstract: With the help of model calculations we aim at reproducing the observational photometric findings for a large sample of well-characterised galactic C-type Mira variables losing mass at different rates. We used dynamic model atmospheres, describing the outer layers of C-rich Miras, which are severly affected by dynamic effects. Based on the resulting structures and under the assumptions of chemical… ▽ More With the help of model calculations we aim at reproducing the observational photometric findings for a large sample of well-characterised galactic C-type Mira variables losing mass at different rates. We used dynamic model atmospheres, describing the outer layers of C-rich Miras, which are severly affected by dynamic effects. Based on the resulting structures and under the assumptions of chemical equilibrium as well as LTE, we computed synthetic spectra and synthetic broad-band photometry (Johnson-Cousins-Glass). A set of five representative models with different stellar parameters describes a sequence from less to more evolved objects with steadily increasing mass-loss rates. This allowed us to study the significant influence of circumstellar dust on the spectral energy distributions and the (amplitudes of) lightcurves in different filters. We tested the photometric properties (mean NIR magnitudes, colours, and amplitudes) and other characteristics of the models (mass-loss rates, periods, and bolometric corrections) by comparing these with the corresponding observational data adopted from the literature. Using different kinds of diagrams we illustrate where the models are located in a supposed evolutionary sequence defined by observed C-type Mira samples. Based on comparisons of galactic targets with empirical relations derived for C stars in the Large Magellanic Cloud we discuss the relevance of metallicity and excess carbon (C--O) for the development of dust-driven winds. Having investigated the dynamic model atmospheres from different (mainly photometric) perspectives, we conclude that our modelling approach (meaning the combination of numerical method and a suitable choice of model parameters) is able to describe C-rich long-period variables over a wide range of mass-loss rates, i.e., from moderately pulsating objects without any dusty wind to highly dust-enshrouded Carbon Miras. △ Less

Submitted 15 February, 2013; originally announced February 2013.

Comments: Accepted for publication in A&A

Conjugacy of Coxeter elements

Authors: Henrik Eriksson, Kimmo Eriksson

Abstract: For a Coxeter group (W,S), a permutation of the set S is called a Coxeter word and the group element represented by the product is called a Coxeter element. Moving the first letter to the end of the word is called a rotation and two Coxeter elements are rotation equivalent if their words can be transformed into each other through a sequence of rotations and legal commutations. We prove that Coxe… ▽ More For a Coxeter group (W,S), a permutation of the set S is called a Coxeter word and the group element represented by the product is called a Coxeter element. Moving the first letter to the end of the word is called a rotation and two Coxeter elements are rotation equivalent if their words can be transformed into each other through a sequence of rotations and legal commutations. We prove that Coxeter elements are conjugate if and only if they are rotation equivalent. This was known for some special cases but not for Coxeter groups in general. △ Less

Submitted 12 February, 2013; originally announced February 2013.

MSC Class: 20F55; 68R15

Journal ref: Electr. J. Combinatorics 16(2), 9p, 2009

The wind of the M-type AGB star RT Virginis probed by VLTI/MIDI

Authors: Stéphane Sacuto, Sofia Ramstedt, Susanne Höfner, Hans Olofsson, Sara Bladh, Kjell Eriksson, Bernhard Aringer, Daniela Klotz, Matthias Maercker

Abstract: We study the circumstellar environment of the M-type AGB star RT Vir using mid-infrared high spatial resolution observations from the ESO-VLTI focal instrument MIDI. The aim of this study is to provide observational constraints on theoretical prediction that the winds of M-type AGB objects can be driven by photon scattering on iron-free silicate grains located in the close environment (about 2 to… ▽ More We study the circumstellar environment of the M-type AGB star RT Vir using mid-infrared high spatial resolution observations from the ESO-VLTI focal instrument MIDI. The aim of this study is to provide observational constraints on theoretical prediction that the winds of M-type AGB objects can be driven by photon scattering on iron-free silicate grains located in the close environment (about 2 to 3 stellar radii) of the star. We interpreted spectro-interferometric data, first using wavelength-dependent geometric models. We then used a self-consistent dynamic model atmosphere containing a time-dependent description of grain growth for pure forsterite dust particles to reproduce the photometric, spectrometric, and interferometric measurements of RT Vir. Since the hydrodynamic computation needs stellar parameters as input, a considerable effort was first made to determine these parameters. MIDI differential phases reveal the presence of an asymmetry in the stellar vicinity. Results from the geometrical modeling give us clues to the presence of aluminum and silicate dust in the close circumstellar environment (< ~5 stellar radii). Comparison between spectro-interferometric data and a self-consistent dust-driven wind model reveals that silicate dust has to be present in the region between 2 to 3 stellar radii to reproduce the 59 and 63 m baseline visibility measurements around 9.8 micron. This gives additional observational evidence in favor of winds driven by photon scattering on iron-free silicate grains located in the close vicinity of an M-type star. However, other sources of opacity are clearly missing to reproduce the 10-13 micron visibility measurements for all baselines. This study is a first attempt to understand the wind mechanism of M-type AGB stars by comparing photometric, spectrometric, and interferometric measurements with state-of-the-art, self-consistent dust-driven wind models. The agreement of the dynamic model atmosphere with interferometric measurements in the 8-10 micron spectral region gives additional observational evidence that the winds of M-type stars can be driven by photon scattering on iron-free silicate grains. Finally, a larger statistical study and progress in advanced self-consistent 3D modeling are still required to solve the remaining problems. △ Less

Submitted 24 January, 2013; originally announced January 2013.

Comments: Astronomy and Astrophysics (2013)

Comparative Modelling of the Spectra of Cool Giants

Authors: T. Lebzelter, U. Heiter, C. Abia, K. Eriksson, M. Ireland, H. Neilson, W. Nowotny, J. Maldonado, T. Merle, R. Peterson, B. Plez, C. I. Short, G. M. Wahlgren, C. Worley, B. Aringer, S. Bladh, P. de Laverny, A. Goswami, A. Mora, R. P. Norris, A. Recio-Blanco, M. Scholz, F. Thevenin, T. Tsuji, G. Kordopatis , et al. (2 additional authors not shown)

Abstract: Our ability to extract information from the spectra of stars depends on reliable models of stellar atmospheres and appropriate techniques for spectral synthesis. Various model codes and strategies for the analysis of stellar spectra are available today. We aim to compare the results of deriving stellar parameters using different atmosphere models and different analysis strategies. The focus is set… ▽ More Our ability to extract information from the spectra of stars depends on reliable models of stellar atmospheres and appropriate techniques for spectral synthesis. Various model codes and strategies for the analysis of stellar spectra are available today. We aim to compare the results of deriving stellar parameters using different atmosphere models and different analysis strategies. The focus is set on high-resolution spectroscopy of cool giant stars. Spectra representing four cool giant stars were made available to various groups and individuals working in the area of spectral synthesis, asking them to derive stellar parameters from the data provided. The results were discussed at a workshop in Vienna in 2010. Most of the major codes currently used in the astronomical community for analyses of stellar spectra were included in this experiment. We present the results from the different groups, as well as an additional experiment comparing the synthetic spectra produced by various codes for a given set of stellar parameters. Similarities and differences of the results are discussed. Several valid approaches to analyze a given spectrum of a star result in quite a wide range of solutions. The main causes for the differences in parameters derived by different groups seem to lie in the physical input data and in the details of the analysis method. This clearly shows how far from a definitive abundance analysis we still are. △ Less

Submitted 12 September, 2012; originally announced September 2012.

Comments: accepted for publication in A&A. This version includes also the online tables. Reference spectra will later be available via the CDS

Explicit Time-Stepping for Stiff ODEs

Authors: Kenneth Eriksson, Claes Johnson, Anders Logg

Abstract: We present a new strategy for solving stiff ODEs with explicit methods. By adaptively taking a small number of stabilizing small explicit time steps when necessary, a stiff ODE system can be stabilized enough to allow for time steps much larger than what is indicated by classical stability analysis. For many stiff problems the cost of the stabilizing small time steps is small, so the improvement i… ▽ More We present a new strategy for solving stiff ODEs with explicit methods. By adaptively taking a small number of stabilizing small explicit time steps when necessary, a stiff ODE system can be stabilized enough to allow for time steps much larger than what is indicated by classical stability analysis. For many stiff problems the cost of the stabilizing small time steps is small, so the improvement is large. We illustrate the technique on a number of well-known stiff test problems. △ Less

Submitted 12 May, 2012; originally announced May 2012.

MSC Class: 65L05; 65L07; 65L20; 65L50; 65L60; 65L70; 65M12

Journal ref: SIAM Journal on Scientific Computing 25(4), pp. 1142--1157 (2003)

Determination of rotation periods in solar-like stars with irregular sampling: the Gaia case

Authors: E. Distefano, A. C. Lanzafame, A. F. Lanza, S. Messina, A. J. Korn, K. Eriksson, J. Cuypers

Abstract: We present a study on the determination of rotation periods (P) of solar-like stars from the photometric irregular time-sampling of the ESA Gaia mission, currently scheduled for launch in 2013, taking into account its dependence on ecliptic coordinates. We examine the case of solar-twins as well as thousands of synthetic time-series of solar-like stars rotating faster than the Sun. In the case of… ▽ More We present a study on the determination of rotation periods (P) of solar-like stars from the photometric irregular time-sampling of the ESA Gaia mission, currently scheduled for launch in 2013, taking into account its dependence on ecliptic coordinates. We examine the case of solar-twins as well as thousands of synthetic time-series of solar-like stars rotating faster than the Sun. In the case of solar twins we assume that the Gaia unfiltered photometric passband G will mimic the variability of the total solar irradiance (TSI) as measured by the VIRGO experiment. For stars rotating faster than the Sun, light-curves are simulated using synthetic spectra for the quiet atmosphere, the spots, and the faculae combined by applying semi-empirical relationships relating the level of photospheric magnetic activity to the stellar rotation and the Gaia instrumental response. The capabilities of the Deeming, Lomb-Scargle, and Phase Dispersion Minimisation methods in recovering the correct rotation periods are tested and compared. The false alarm probability (FAP) is computed using Monte Carlo simulations and compared with analytical formulae. The Gaia scanning law makes the rate of correct detection of rotation periods strongly dependent on the ecliptic latitude (beta). We find that for P ~ 1 d, the rate of correct detection increases with ecliptic latitude from 20-30 per cent at beta ~ 0° to a peak of 70 per cent at beta=45°, then it abruptly falls below 10 per cent at beta > 45°. For P > 5 d, the rate of correct detection is quite low and for solar twins is only 5 per cent on average. △ Less

Submitted 3 February, 2012; originally announced February 2012.

Comments: 12 pages, 18 figures, accepted by MNRAS

Analysis of old very metal rich stars in the solar neighbourhood

Authors: M. Trevisan, B. Barbuy, K. Eriksson, B. Gustafsson, M. Grenon, L. Pompéia

Abstract: A sample of mostly old metal-rich dwarf and turn-off stars with high eccentricity and low maximum height above the Galactic plane has been identified. From their kinematics, it was suggested that the inner disk is their most probable birthplace. Their chemical imprints may therefore reveal important information about the formation and evolution of the still poorly understood inner disk. To probe t… ▽ More A sample of mostly old metal-rich dwarf and turn-off stars with high eccentricity and low maximum height above the Galactic plane has been identified. From their kinematics, it was suggested that the inner disk is their most probable birthplace. Their chemical imprints may therefore reveal important information about the formation and evolution of the still poorly understood inner disk. To probe the formation history of these stellar populations, a detailed analysis of a sample of very metal-rich stars is carried out. We derive the metallicities, abundances of α elements, ages, and Galactic orbits. The analysis of 71 metal-rich stars is based on optical high-resolution échelle spectra obtained with the FEROS spectrograph at the ESO 1.52-m Telescope at La Silla, Chile. The metallicities and abundances of C, O, Mg, Si, Ca, and Ti were derived based on LTE detailed analysis, employing the MARCS model atmospheres. We confirm the high metallicity of these stars reaching up to [Fe I/H]~0.58, and the sample of metal-rich dwarfs can be kinematically subclassified in samples of thick disk, thin disk, and intermediate stellar populations. All sample stars show solar α-Fe ratios, and most of them are old and still quite metal rich. The orbits suggest that the thin disk, thick disk and intermediate populations were formed at Galactocentric distances of ~8 kpc, ~6 kpc, and ~7 kpc, respectively. The mean maximum height of the thick disk subsample of Z_max~380 pc, is lower than for typical thick disk stars. A comparison of α-element abundances of the sample stars with bulge stars shows that the oxygen is compatible with a bulge or inner thick disk origin. Our results suggest that models of radial mixing and dynamical effects of the bar and bar/spiral arms might explain the presence of these old metal-rich dwarf stars in the solar neighbourhood. △ Less

Submitted 28 September, 2011; originally announced September 2011.

Comments: 28 pages, 27 figures, Accepted for publication in Astronomy & Astrophysics

Sulphur abundances in halo giants from the [S I] line at 1082 nm and the S I triplet around 1045 nm

Authors: H. Jönsson, N. Ryde, P. E. Nissen, R. Collet, K. Eriksson, M. Asplund, B. Gustafsson

Abstract: It is still debated whether or not the Galactic chemical evolution of sulphur in the halo followed the constant or flat trend with [Fe/H], ascribed to the result of explosive nucleosynthesis in type II SNe. The aim of this study is to try to clarify this situation by measuring the sulphur abundance in a sample of halo giants using two diagnostics; the S I triplet around 1045 nm and the [S I] line… ▽ More It is still debated whether or not the Galactic chemical evolution of sulphur in the halo followed the constant or flat trend with [Fe/H], ascribed to the result of explosive nucleosynthesis in type II SNe. The aim of this study is to try to clarify this situation by measuring the sulphur abundance in a sample of halo giants using two diagnostics; the S I triplet around 1045 nm and the [S I] line at 1082 nm. The latter of the two is not believed to be sensitive to non-LTE effects. We can thereby minimize the uncertainties in the diagnostic used and estimate the usefulness of the triplet in sulphur determination in halo K giants. We will also be able to compare our sulphur abundance differences from the two diagnostics with the expected non-LTE effects in the 1045 nm triplet previously calculated by others. High-resolution near-infrared spectra of ten K giants were recorded using the spectrometer CRIRES mounted on VLT. Two standard settings were used; one covering the S I triplet and one covering the [S I] line. The sulphur abundances were determined individually with equivalent widths and synthetic spectra for the two diagnostics using tailored 1D model atmospheres and relying on non-LTE corrections from the litterature. Effects of convective inhomogeneities in the stellar atmospheres are investigated. We corroborate the flat trend in the [S/Fe] vs. [Fe/H] plot for halo stars found in other works and cannot find a scatter nor a rise in [S/Fe] obtained in some other previous studies. We find the sulphur abundances deduced from the non-LTE corrected triplet somewhat lower than the abundances from the [S I] line, possibly indicating too large non-LTE corrections. Considering 3D modeling, however, they might instead be too small. Further we show that the [S I] line is possible to use as a sulphur diagnostic down to [Fe/H] = -2.3 in giants. △ Less

Submitted 12 April, 2011; originally announced April 2011.

Comments: Accepted for publication in Astronomy & Astrophysics

Ices in the Quiescent IC 5146 Dense Cloud

Authors: J. E. Chiar, Y. J. Pendleton, L. J. Allamandola, A. C. A. Boogert, K. Ennico, T. P. Greene, T. R. Geballe, J. V. Keane, C. J. Lada, R. E. Mason, T. L. Roellig, S. A. Sandford, A. G. G. M. Tielens, M. W. Werner, D. C. B. Whittet, L. Decin, K. Eriksson

Abstract: This paper presents spectra in the 2 to 20 micron range of quiescent cloud material located in the IC 5146 cloud complex. The spectra were obtained with NASA's Infrared Telescope Facility (IRTF) SpeX instrument and the Spitzer Space Telescope's Infrared Spectrometer. We use these spectra to investigate dust and ice absorption features in pristine regions of the cloud that are unaltered by embedded… ▽ More This paper presents spectra in the 2 to 20 micron range of quiescent cloud material located in the IC 5146 cloud complex. The spectra were obtained with NASA's Infrared Telescope Facility (IRTF) SpeX instrument and the Spitzer Space Telescope's Infrared Spectrometer. We use these spectra to investigate dust and ice absorption features in pristine regions of the cloud that are unaltered by embedded stars. We find that the H2O-ice threshold extinction is 4.03+/-0.05 mag. Once foreground extinction is taken into account, however, the threshold drops to 3.2 mag, equivalent to that found for the Taurus dark cloud, generally assumed to be the touchstone quiescent cloud against which all other dense cloud and embedded young stellar object observations are compared. Substructure in the trough of the silicate band for two sources is attributed to CH3OH and NH3 in the ices, present at the ~2% and ~5% levels, respectively, relative to H2O-ice. The correlation of the silicate feature with the E(J-K) color excess is found to follow a much shallower slope relative to lines of sight that probe diffuse clouds, supporting the previous results by Chiar et al. (2007). △ Less

Submitted 12 February, 2011; originally announced February 2011.

Comments: 13 pages, 13 figures with multiple parts, accepted for publication in Astrophysical Journal, Feb. 2011

A grid of MARCS model atmospheres for S stars

Authors: Sophie Van Eck, Pieter Neyskens, Bertrand Plez, Alain Jorissen, Bengt Edvardsson, Kjell Eriksson, Bengt Gustafsson, Uffe-Grae Jorgensen, Ake Nordlund

Abstract: S-type stars are late-type giants whose atmosphere is enriched in carbon and s-process elements because of either extrinsic pollution by a binary companion or intrinsic nucleosynthesis and dredge-up on the thermally-pulsing AGB. A large grid of S-star model atmospheres has been computed covering the range 2700 < Teff < 4000 K with 0.5 < C/O < 0.99. ZrO and TiO band strength indices as well as VJHK… ▽ More S-type stars are late-type giants whose atmosphere is enriched in carbon and s-process elements because of either extrinsic pollution by a binary companion or intrinsic nucleosynthesis and dredge-up on the thermally-pulsing AGB. A large grid of S-star model atmospheres has been computed covering the range 2700 < Teff < 4000 K with 0.5 < C/O < 0.99. ZrO and TiO band strength indices as well as VJHKL photometry are needed to disentangle Teff, C/O and [s/Fe]. A "best-model finding tool" was developed using a set of well-chosen indices and checked against photometry as well as low- and high-resolution spectroscopy. It is found that applying M-star model atmospheres (i.e., with a solar C/O ratio) to S stars can lead to errors on Teff up to 400K. We constrain the parameter space occupied by S stars of the vast sample of Henize stars in terms of Teff, [C/O] and [s/Fe]. △ Less

Submitted 9 November, 2010; originally announced November 2010.

Comments: 5 pages, 4 figures, "Why Galaxies Care About AGB Stars II" Conference, Vienna Aug. 2010, to appear in ASP Conference Series

An abundance study of red-giant-branch stars in the Hercules dwarf spheroidal galaxy

Authors: D. Adén, K. Eriksson, S. Feltzing, E. K. Grebel, A. Koch, M. I. Wilkinson

Abstract: Using high-resolution spectroscopy, we provide a determination of [Fe/H] and [Ca/H] for confirmed red-giant branch member stars of the Hercules dwarf spheroidal galaxy. Based on this we explore the ages of the prevailing stellar populations in Hercules, and the enrichment history from supernovae. Additionally, we provide a new simple metallicity calibration for Stromgren photometry for metal-poor,… ▽ More Using high-resolution spectroscopy, we provide a determination of [Fe/H] and [Ca/H] for confirmed red-giant branch member stars of the Hercules dwarf spheroidal galaxy. Based on this we explore the ages of the prevailing stellar populations in Hercules, and the enrichment history from supernovae. Additionally, we provide a new simple metallicity calibration for Stromgren photometry for metal-poor, red giant branch stars. We find that the red-giant branch stars of the Hercules dSph galaxy are more metal-poor than estimated in our previous study that was based on photometry alone. Additionally, we find an abundance trend such that [Ca/Fe] is higher for more metal-poor stars, and lower for more metal-rich stars, with a spread of about 0.8 dex. The [Ca/Fe] trend suggests an early rapid chemical enrichment through supernovae of type II, followed by a phase of slow star formation dominated by enrichment through supernovae of type Ia. A comparison with isochrones indicates that the red giants in Hercules are older than 10 Gyr. △ Less

Submitted 27 October, 2010; originally announced October 2010.

Comments: 12 pages, 11 figures. Accepted for publication in A&A

Fluorine Abundances in Galactic AGB Stars

Authors: C. Abia, K. Cunha, S. Cristallo., P. de Laverny, I. Dominguez, K. Eriksson, L. Gialanella, K. Hinkle, G. Imbriani, A Recio-Blanco, V. V. Smith, O. Straniero, R. Wahlin

Abstract: An analysis of the fluorine abundance in Galactic AGB carbon stars (24 N-type, 5 SC-type and 5 J-type) is presented. This study uses the state- of-the-art carbon rich atmosphere models and improved atomic and molecular line lists in the 2.3 μm region. F abundances significantly lower are obtained in comparison to previous study in the literature. The main reason of this difference is due to mo… ▽ More An analysis of the fluorine abundance in Galactic AGB carbon stars (24 N-type, 5 SC-type and 5 J-type) is presented. This study uses the state- of-the-art carbon rich atmosphere models and improved atomic and molecular line lists in the 2.3 μm region. F abundances significantly lower are obtained in comparison to previous study in the literature. The main reason of this difference is due to molecular blends. In the case of carbon stars of SC-type, differences in the model atmospheres are also relevant. The new F enhancements are now in agreement with the most recent theoretical nucleosynthesis models in low- mass AGB stars, solving the long standing problem of F in Galactic AGB stars. Nevertheless, some SC-type carbon stars still show larger F abundances than predicted by stellar models. The possibility that these stars are of larger mass is briefly discussed. △ Less

Submitted 26 April, 2010; originally announced April 2010.

Comments: 16 pages, 3 figures, Accepted in ApJ Letters

Oxygen isotopic ratios in cool R Coronae Borealis stars

Authors: D. Anibal Garcia-Hernandez, David L. Lambert, N. Kameswara Rao, Ken H. Hinkle, Kjell Eriksson

Abstract: We investigate the relationship between R Coronae Borealis (RCB) stars and hydrogen-deficient carbon (HdC) stars by measuring precise 16O/18O ratios for five cool RCB stars. The 16O/18O ratios are derived by spectrum synthesis from high-resolution (R=50,000) K-band spectra. Lower limits to the 16O/17O and 14N/15}N ratios as well as Na and S abundances (when possible) are also given. RCB stars in o… ▽ More We investigate the relationship between R Coronae Borealis (RCB) stars and hydrogen-deficient carbon (HdC) stars by measuring precise 16O/18O ratios for five cool RCB stars. The 16O/18O ratios are derived by spectrum synthesis from high-resolution (R=50,000) K-band spectra. Lower limits to the 16O/17O and 14N/15}N ratios as well as Na and S abundances (when possible) are also given. RCB stars in our sample generally display less 18O than HdC stars - the derived 16O/18O ratios range from 3 to 20. The only exception is the RCB star WX CrA, which seems to be a HdC-like star with 16O/18O=0.3. Our result of a higher 16O/18O ratio for the RCB stars must be accounted for by a theory of the formation and evolution of HdC and RCB stars. We speculate that a late dredge-up of products of He-burning, principally 12C and 16O, may convert a 18O-rich HdC star into a 18O-poor RCB star as the H-deficient star begins its final evolution from a cool supergiant to the top of the white dwarf cooling track. △ Less

Submitted 15 March, 2010; originally announced March 2010.

Comments: Accepted for publication in the Astrophysical Journal; 21 pages, 3 Tables, 11 Figures