Showing 1–50 of 52 results for author: Sturm, A

JWST ice band profiles reveal mixed ice compositions in the HH 48 NE disk

Authors: Jennifer B. Bergner, J. A. Sturm, Elettra L. Piacentino, M. K. McClure, Karin I. Oberg, A. C. A. Boogert, E. Dartois, M. N. Drozdovskaya, H. J. Fraser, Daniel Harsono, Sergio Ioppolo, Charles J. Law, Dariusz C. Lis, Brett A. McGuire, Gary J. Melnick, Jennifer A. Noble, M. E. Palumbo, Yvonne J. Pendleton, Giulia Perotti, Danna Qasim, W. R. M. Rocha, E. F. van Dishoeck

Abstract: Planet formation is strongly influenced by the composition and distribution of volatiles within protoplanetary disks. With JWST, it is now possible to obtain direct observational constraints on disk ices, as recently demonstrated by the detection of ice absorption features towards the edge-on HH 48 NE disk as part of the Ice Age Early Release Science program. Here, we introduce a new radiative tra… ▽ More Planet formation is strongly influenced by the composition and distribution of volatiles within protoplanetary disks. With JWST, it is now possible to obtain direct observational constraints on disk ices, as recently demonstrated by the detection of ice absorption features towards the edge-on HH 48 NE disk as part of the Ice Age Early Release Science program. Here, we introduce a new radiative transfer modeling framework designed to retrieve the composition and mixing status of disk ices using their band profiles, and apply it to interpret the H2O, CO2, and CO ice bands observed towards the HH 48 NE disk. We show that the ices are largely present as mixtures, with strong evidence for CO trapping in both H2O and CO2 ice. The HH 48 NE disk ice composition (pure vs. polar vs. apolar fractions) is markedly different from earlier protostellar stages, implying thermal and/or chemical reprocessing during the formation or evolution of the disk. We infer low ice-phase C/O ratios around 0.1 throughout the disk, and also demonstrate that the mixing and entrapment of disk ices can dramatically affect the radial dependence of the C/O ratio. It is therefore imperative that realistic disk ice compositions are considered when comparing planetary compositions with potential formation scenarios, which will fortunately be possible for an increasing number of disks with JWST. △ Less

Submitted 12 September, 2024; originally announced September 2024.

Comments: Accepted to ApJ. 24 pages, 15 figures

A min-max random game on a graph that is not a tree

Authors: Natalia Cardona-Tobón, Anja Sturm, Jan M. Swart

Abstract: We study a random game in which two players in turn play a fixed number of moves. For each move, there are two possible choices. To each possible outcome of the game we assign a winner in an i.i.d. fashion with a fixed parameter p. In the case where all different game histories lead to different outcomes, a classical result due to Pearl (1980) says that in the limit when the number of moves is lar… ▽ More We study a random game in which two players in turn play a fixed number of moves. For each move, there are two possible choices. To each possible outcome of the game we assign a winner in an i.i.d. fashion with a fixed parameter p. In the case where all different game histories lead to different outcomes, a classical result due to Pearl (1980) says that in the limit when the number of moves is large, there is a sharp threshold in the parameter p that separates the regimes in which either player has with high probability a winning strategy. We are interested in a modification of this game where the outcome is determined by the exact sequence of moves played by the first player and by the number of times the second player has played each of the two possible moves. We show that also in this case, there is a sharp threshold in the parameter p that separates the regimes in which either player has with high probability a winning strategy. Since in the modified game, different game histories can lead to the same outcome, the graph associated with the game is no longer a tree which means independence is lost. As a result, the analysis becomes more complicated and open problems remain. △ Less

Submitted 4 September, 2024; originally announced September 2024.

Comments: 31 pages, 7 figures

MSC Class: Primary: 82C26; Secondary: 60K35; 91A15; 91A50

A JWST/MIRI analysis of the ice distribution and PAH emission in the protoplanetary disk HH 48 NE

Authors: J. A. Sturm, M. K. McClure, D. Harsono, J. B. Bergner, E. Dartois, A. C. A. Boogert, M. A. Cordiner, M. N. Drozdovskaya, S. Ioppolo, C. J. Law, D. C. Lis, B. A. McGuire, G. J. Melnick, J. A. Noble, K. I. Öberg, M. E. Palumbo, Y. J. Pendleton, G. Perotti, W. R. M. Rocha, R. G. Urso, E. F. van Dishoeck

Abstract: Ice-coated dust grains provide the main reservoir of volatiles that play an important role in planet formation processes and may become incorporated into planetary atmospheres. However, due to observational challenges, the ice abundance distribution in protoplanetary disks is not well constrained. We present JWST/MIRI observations of the edge-on disk HH 48 NE carried out as part of the IRS program… ▽ More Ice-coated dust grains provide the main reservoir of volatiles that play an important role in planet formation processes and may become incorporated into planetary atmospheres. However, due to observational challenges, the ice abundance distribution in protoplanetary disks is not well constrained. We present JWST/MIRI observations of the edge-on disk HH 48 NE carried out as part of the IRS program Ice Age. We detect CO$_2$, NH$_3$, H$_2$O and tentatively CH$_4$ and NH$_4^+$. Radiative transfer models suggest that ice absorption features are produced predominantly in the 50-100 au region of the disk. The CO$_2$ feature at 15 micron probes a region closer to the midplane (z/r = 0.1-0.15) than the corresponding feature at 4.3 micron (z/r = 0.2-0.6), but all observations trace regions significantly above the midplane reservoirs where we expect the bulk of the ice mass to be located. Ices must reach a high scale height (z/r ~ 0.6; corresponding to modeled dust extinction Av ~ 0.1), in order to be consistent with the observed vertical distribution of the peak ice optical depths. The weakness of the CO$_2$ feature at 15 micron relative to the 4.3 micron feature and the red emission wing of the 4.3 micron CO$_2$ feature are both consistent with ices being located at high elevation in the disk. The retrieved NH$_3$ abundance and the upper limit on the CH$_3$OH abundance relative to H$_2$O are significantly lower than those in the interstellar medium (ISM), but consistent with cometary observations. Full wavelength coverage is required to properly study the abundance distribution of ices in disks. To explain the presence of ices at high disk altitudes, we propose two possible scenarios: a disk wind that entrains sufficient amounts of dust, thus blocking part of the stellar UV radiation, or vertical mixing that cycles enough ices into the upper disk layers to balance ice photodesorption. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 16 pages, 11 figures, accepted in A&A

Journal ref: A&A 689, A92 (2024)

Constraints on the gas-phase C/O ratio of DR Tau's outer disk from CS, SO, and C$_2$H observations

Authors: Jane Huang, Edwin A. Bergin, Romane Le Gal, Sean M. Andrews, Jaehan Bae, Luke Keyte, J. A. Sturm

Abstract: Millimeter wavelength observations of Class II protoplanetary disks often display strong emission from hydrocarbons and high CS/SO values, providing evidence that the gas-phase C/O ratio commonly exceeds 1 in their outer regions. We present new NOEMA observations of CS $5-4$, SO $7_6-6_5$ and $5_6-4_5$, C$_2$H $N=3-2$, HCN $3-2$, HCO$^+$ $3-2$, and H$^{13}$CO$^+$ $3-2$ in the DR Tau protoplanetary… ▽ More Millimeter wavelength observations of Class II protoplanetary disks often display strong emission from hydrocarbons and high CS/SO values, providing evidence that the gas-phase C/O ratio commonly exceeds 1 in their outer regions. We present new NOEMA observations of CS $5-4$, SO $7_6-6_5$ and $5_6-4_5$, C$_2$H $N=3-2$, HCN $3-2$, HCO$^+$ $3-2$, and H$^{13}$CO$^+$ $3-2$ in the DR Tau protoplanetary disk at a resolution of $\sim0.4''$ (80 au). Estimates for the disk-averaged CS/SO ratio range from $\sim0.4-0.5$, the lowest value reported thus far for a T Tauri disk. At a projected separation of $\sim180$ au northeast of the star, the SO moment maps exhibit a clump that has no counterpart in the other lines, and the CS/SO value decreases to $<0.2$ at its location. Thermochemical models calculated with DALI indicate that DR Tau's low CS/SO ratio and faint C$_2$H emission can be explained by a gas-phase C/O ratio that is $<1$ at the disk radii traced by NOEMA. Comparisons of DR Tau's SO emission to maps of extended structures traced by $^{13}$CO suggest that late infall may contribute to driving down the gas-phase C/O ratio of its disk. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: Accepted by ApJ

The contact process on dynamical random trees with degree dependence

Authors: Natalia Cardona-Tobón, Marcel Ortgiese, Marco Seiler, Anja Sturm

Abstract: The contact process is a simple model for the spread of an infection in a structured population. We investigate the case when the underlying structure evolves dynamically as a degree-dependent dynamical percolation model. Starting with a connected locally finite base graph we initially declare edges independently open with a probability that is allowed to depend on the degree of the adjacent verti… ▽ More The contact process is a simple model for the spread of an infection in a structured population. We investigate the case when the underlying structure evolves dynamically as a degree-dependent dynamical percolation model. Starting with a connected locally finite base graph we initially declare edges independently open with a probability that is allowed to depend on the degree of the adjacent vertices and closed otherwise. Edges are independently updated with a rate depending on the degrees and then are again declared open and closed with the same probabilities. We are interested in the contact process, where infections are only allowed to spread via open edges. Our aim is to analyze the impact of the update speed and the probability for edges to be open on the existence of a phase transition. For a general connected locally finite graph, our first result gives sufficient conditions for the critical value for survival to be strictly positive. Furthermore, in the setting of Bienaymé-Galton-Watson trees, we show that the process survives strongly with positive probability for any infection rate if the offspring distribution has a stretched exponential tail with an exponent depending on the percolation probability and the update speed. In particular, if the offspring distribution follows a power law and the connection probability is given by a product kernel and the update speed exhibits polynomial behaviour, we provide a complete characterisation of the phase transition. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 56 pages, 4 figures

MSC Class: Primary 60K35; Secondary 05C80; 82C22

Quantum Circuit Learning on NISQ Hardware

Authors: Niclas Schillo, Andreas Sturm

Abstract: Current quantum computers are small and error-prone systems for which the term noisy intermediate-scale quantum (NISQ) has become established. Since large scale, fault-tolerant quantum computers are not expected to be available in the near future, the task of finding NISQ suitable algorithms has received a lot of attention in recent years. The most prominent candidates in this context are variatio… ▽ More Current quantum computers are small and error-prone systems for which the term noisy intermediate-scale quantum (NISQ) has become established. Since large scale, fault-tolerant quantum computers are not expected to be available in the near future, the task of finding NISQ suitable algorithms has received a lot of attention in recent years. The most prominent candidates in this context are variational quantum algorithms. Due to their hybrid quantum-classical architecture they require fewer qubits and quantum gates so that they can cope with the limitations of NISQ computers. An important class of variational quantum algorithms is the quantum circuit learning (QCL) framework. Consisting of a data encoding and a trainable, parametrized layer, these schemes implement a quantum model function that can be fitted to the problem at hand. For instance, in combination with the parameter shift rule to compute derivatives, they can be used to solve differential equations. QCL and related algorithms have been widely studied in the literature. However, numerical experiments are usually limited to simulators and results from real quantum computers are scarce. In this paper we close this gap by executing QCL circuits on a superconducting IBM quantum processor in conjunction with an analysis of the hardware errors. We show that exemplary QCL circuits with up to three qubits are executable on the IBM quantum computer. For this purpose, multiple functions are learned and an exemplary differential equation is solved on the quantum computer. Moreover, we present how the QCL framework can be used to learn different quantum model functions in parallel, which can be applied to solve coupled differential equations in an efficient way. △ Less

Submitted 3 May, 2024; originally announced May 2024.

A Central Limit Theorem for Functions on Weighted Sparse Inhomogeneous Random Graphs

Authors: Anja Sturm, Moritz Wemheuer

Abstract: We prove a central limit theorem for a certain class of functions on sparse rank-one inhomogeneous random graphs endowed with additional i.i.d. edge and vertex weights. Our proof of the central limit theorem uses a perturbative form of Stein's method and relies on a careful analysis of the local structure of the underlying sparse inhomogeneous random graphs (as the number of vertices in the graph… ▽ More We prove a central limit theorem for a certain class of functions on sparse rank-one inhomogeneous random graphs endowed with additional i.i.d. edge and vertex weights. Our proof of the central limit theorem uses a perturbative form of Stein's method and relies on a careful analysis of the local structure of the underlying sparse inhomogeneous random graphs (as the number of vertices in the graph tends to infinity), which may be of independent interest, as well as a local approximation property of the function, which is satisfied for a number of combinatorial optimisation problems. These results extend recent work by Cao (2021) for Erdős-Rényi random graphs and additional i.i.d. weights only on the edges. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: 64 pages, 4 figures

Unlocking Quantum Optimization: A Use Case Study on NISQ Systems

Authors: Andreas Sturm, Bharadwaj Mummaneni, Leon Rullkötter

Abstract: The major advances in quantum computing over the last few decades have sparked great interest in applying it to solve the most challenging computational problems in a wide variety of areas. One of the most pronounced domains here are optimization problems and a number of algorithmic approaches have been proposed for their solution. For the current noisy intermediate-scale quantum (NISQ) computers… ▽ More The major advances in quantum computing over the last few decades have sparked great interest in applying it to solve the most challenging computational problems in a wide variety of areas. One of the most pronounced domains here are optimization problems and a number of algorithmic approaches have been proposed for their solution. For the current noisy intermediate-scale quantum (NISQ) computers the quantum approximate optimization algorithm (QAOA), the variational quantum eigensolver (VQE), and quantum annealing (QA) are the central algorithms for this problem class. The two former can be executed on digital gate-model quantum computers, whereas the latter requires a quantum annealer. Across all hardware architectures and manufactures, the quantum computers available today share the property of being too error-prone to reliably execute involved quantum circuits as they typically arise from quantum optimization algorithms. In order to characterize the limits of existing quantum computers, many component and system level benchmarks have been proposed. However, owing to the complex nature of the errors in quantum systems these benchmark fail to provide predictive power beyond simple quantum circuits and small examples. Application oriented benchmarks have been proposed to remedy this problem, but both, results from real quantum systems as well as use cases beyond constructed academic examples, remain very rare. This paper addresses precisely this gap by considering two industrial relevant use cases: one in the realm of optimizing charging schedules for electric vehicles, the other concerned with the optimization of truck routes. Our central contribution are systematic series of examples derived from these uses cases that we execute on different processors of the gate-based quantum computers of IBM as well as on the quantum annealer of D-Wave. △ Less

Submitted 10 April, 2024; originally announced April 2024.

JWST MIRI MRS Images Disk Winds, Water, and CO in an Edge-On Protoplanetary Disk

Authors: Nicole Arulanantham, M. K. McClure, Klaus Pontoppidan, Tracy L. Beck, J. A. Sturm, D. Harsono, A. C. A. Boogert, M. Cordiner, E. Dartois, M. N. Drozdovskaya, C. Espaillat, G. J. Melnick, J. A. Noble, M. E. Palumbo, Y. J. Pendleton, H. Terada, E. F. van Dishoeck

Abstract: We present JWST MIRI MRS observations of the edge-on protoplanetary disk around the young sub-solar mass star Tau 042021, acquired as part of the Cycle 1 GO program "Mapping Inclined Disk Astrochemical Signatures (MIDAS)." These data resolve the mid-IR spatial distributions of H$_2$, revealing X-shaped emission extending to ~200 au above the disk midplane with a semi-opening angle of $35 \pm 5$ de… ▽ More We present JWST MIRI MRS observations of the edge-on protoplanetary disk around the young sub-solar mass star Tau 042021, acquired as part of the Cycle 1 GO program "Mapping Inclined Disk Astrochemical Signatures (MIDAS)." These data resolve the mid-IR spatial distributions of H$_2$, revealing X-shaped emission extending to ~200 au above the disk midplane with a semi-opening angle of $35 \pm 5$ degrees. We do not velocity-resolve the gas in the spectral images, but the measured semi-opening angle of the H$_2$ is consistent with an MHD wind origin. A collimated, bipolar jet is seen in forbidden emission lines from [Ne II], [Ne III], [Ni II], [Fe II], [Ar II], and [S III]. Extended H$_2$O and CO emission lines are also detected, reaching diameters between ~90 and 190 au, respectively. Hot molecular emission is not expected at such radii, and we interpret its extended spatial distribution as scattering of inner disk molecular emission by dust grains in the outer disk surface. H I recombination lines, characteristic of inner disk accretion shocks, are similarly extended, and are likely also scattered light from the innermost star-disk interface. Finally, we detect extended PAH emission at 11.3 microns co-spatial with the scattered light continuum, making this the first low-mass T Tauri star around which extended PAHs have been confirmed, to our knowledge. MIRI MRS line images of edge-on disks provide an unprecedented window into the outflow, accretion, and scattering processes within protoplanetary disks, allowing us to constrain the disk lifetimes and accretion and mass loss mechanisms. △ Less

Submitted 20 March, 2024; v1 submitted 19 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJL on March 13th, 2024

Sub-Riemannian Random Walks: From Connections to Retractions

Authors: Michael Herrmann, Pit Neumann, Simon Schwarz, Anja Sturm, Max Wardetzky

Abstract: We study random walks on sub-Riemannian manifolds using the framework of retractions, i.e., approximations of normal geodesics. We show that such walks converge to the correct horizontal Brownian motion if normal geodesics are approximated to at least second order. In particular, we (i) provide conditions for convergence of geodesic random walks defined with respect to normal, compatible, and part… ▽ More We study random walks on sub-Riemannian manifolds using the framework of retractions, i.e., approximations of normal geodesics. We show that such walks converge to the correct horizontal Brownian motion if normal geodesics are approximated to at least second order. In particular, we (i) provide conditions for convergence of geodesic random walks defined with respect to normal, compatible, and partial connections and (ii) provide examples of computationally efficient retractions, e.g., for simulating anisotropic Brownian motion on Riemannian manifolds. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 21 pages

JWST Detects Neon Line Variability in a Protoplanetary Disk

Authors: C. C. Espaillat, T. Thanathibodee, C. V. Pittman, J. A. Sturm, M. K. McClure, N. Calvet, F. M. Walter, R. Franco-Hernandez, J. Muzerolle Page

Abstract: We report the first detection of variability in the mid-infrared neon line emission of a protoplanetary disk by comparing a JWST MIRI MRS spectrum of SZ Cha taken in 2023 with a Spitzer IRS SH spectrum of this object from 2008. We measure the [Ne III]-to-[Ne II] line flux ratio, which is a diagnostic of the high-energy radiation field, to distinguish between the dominance of EUV- or X-ray-driven d… ▽ More We report the first detection of variability in the mid-infrared neon line emission of a protoplanetary disk by comparing a JWST MIRI MRS spectrum of SZ Cha taken in 2023 with a Spitzer IRS SH spectrum of this object from 2008. We measure the [Ne III]-to-[Ne II] line flux ratio, which is a diagnostic of the high-energy radiation field, to distinguish between the dominance of EUV- or X-ray-driven disk photoevaporation. We find that the [Ne III]-to-[Ne II] line flux ratio changes significantly from $\sim1.4$ in 2008 to $\sim0.2$ in 2023. This points to a switch from EUV-dominated to X-ray-dominated photoevaporation of the disk. We present contemporaneous ground-based optical spectra of the Halpha emission line that show the presence of a strong wind in 2023. We propose that this strong wind prevents EUV radiation from reaching the disk surface while the X-rays permeate the wind and irradiate the disk. We speculate that at the time of the Spitzer observations, the wind was suppressed and EUV radiation reached the disk. These observations confirm that the MIR neon emission lines are sensitive to changes in high-energy radiation reaching the disk surface. This highlights the [Ne III]-to-[Ne II] line flux ratio as a tool to gauge the efficiency of disk photoevaporation in order to provide constraints on the planet-formation timescale. However, multiwavelength observations are crucial to interpret the observations and properly consider the star-disk connection. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: Accepted to ApJL

A JWST inventory of protoplanetary disk ices: The edge-on protoplanetary disk HH 48 NE, seen with the Ice Age ERS program

Authors: J. A. Sturm, M. K. McClure, T. L. Beck, D. Harsono, J. B. Bergner, E. Dartois, A. C. A. Boogert, J. E. Chiar, M. A. Cordiner, M. N. Drozdovskaya, S. Ioppolo, C. J. Law, H. Linnartz, D. C. Lis, G. J. Melnick, B. A. McGuire, J. A. Noble, K. I. Öberg, M. E. Palumbo, Y. J. Pendleton, G. Perotti, K. M. Pontoppidan, D. Qasim, W. R. M. Rocha, H. Terada , et al. (2 additional authors not shown)

Abstract: Ices are the main carriers of volatiles in protoplanetary disks and are crucial to our understanding of the chemistry that ultimately sets the organic composition of planets. The ERS program Ice Age on the JWST follows the ice evolution through all stages of star and planet formation. JWST/NIRSpec observations of the edge-on Class II protoplanetary disk HH~48~NE reveal spatially resolved absorptio… ▽ More Ices are the main carriers of volatiles in protoplanetary disks and are crucial to our understanding of the chemistry that ultimately sets the organic composition of planets. The ERS program Ice Age on the JWST follows the ice evolution through all stages of star and planet formation. JWST/NIRSpec observations of the edge-on Class II protoplanetary disk HH~48~NE reveal spatially resolved absorption features of the major ice components H$_2$O, CO$_2$, CO, and multiple weaker signatures from less abundant ices NH$_3$, OCN$^-$, and OCS. Isotopologue $^{13}$CO$_2$ ice has been detected for the first time in a protoplanetary disk. Since multiple complex light paths contribute to the observed flux, the ice absorption features are filled in by ice-free scattered light. The $^{12}$CO$_2$/$^{13}$CO$_2$ ratio of 14 implies that the $^{12}$CO$_2$ feature is saturated, without the flux approaching 0, indicative of a very high CO$_2$ column density on the line of sight, and a corresponding abundance with respect to hydrogen that is higher than ISM values by a factor of at least a few. Observations of rare isotopologues are crucial, as we show that the $^{13}$CO$_2$ observation allows us to determine the column density of CO$_2$ to be at an order of magnitude higher than the lower limit directly inferred from the observed optical depth. Radial variations in ice abundance, e.g., snowlines, are significantly modified since all observed photons have passed through the full radial extent of the disk. CO ice is observed at perplexing heights in the disk, extending to the top of the CO-emitting gas layer. We argue that the most likely interpretation is that we observe some CO ice at high temperatures, trapped in less volatile ices like H$_2$O and CO$_2$. Future radiative transfer models will be required to constrain the implications on our current understanding of disk physics and chemistry. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 16 pages, 8 figures, accepted for publication in A&A

The edge-on protoplanetary disk HH 48 NE II. Modeling ices and silicates

Authors: J. A. Sturm, M. K. McClure, J. B. Bergner, D. Harsono, E. Dartois, M. N. Drozdovskaya, S. Ioppolo, K. I. Öberg, C. J. Law, M. E. Palumbo, Y. J. Pendleton, W. R. M. Rocha, H. Terada, R. G. Urso

Abstract: The abundance and distribution of ice in protoplanetary disks (PPD) is critical to understand the linkage between the composition of circumstellar matter and the composition of exoplanets. Edge-on PPDs are a useful tool to constrain such ice composition and its location in the disk, as ice spectral signatures can be observed in absorption against the continuum emission arising from the warmer cent… ▽ More The abundance and distribution of ice in protoplanetary disks (PPD) is critical to understand the linkage between the composition of circumstellar matter and the composition of exoplanets. Edge-on PPDs are a useful tool to constrain such ice composition and its location in the disk, as ice spectral signatures can be observed in absorption against the continuum emission arising from the warmer central disk regions. The aim of this work is to model ice absorption features in PPDs and determine how well the abundance of the main ice species across the disk can be determined within the uncertainty of the physical parameter space. The edge-on PPD around HH 48 NE, a target of the JWST ERS program IceAge, is used as a reference system. We use RADMC-3D to raytrace the mid-infrared continuum. Using a constant parameterized ice abundance, ice opacities are added to the dust opacity in regions wherever the disk is cold enough for the main carbon, oxygen and nitrogen carriers to freeze out. The global abundance of the main ice carriers in HH 48 NE can be determined within a factor of 3, when taking the uncertainty of the physical parameters into account. Ice features in PPDs can be saturated at an optical depth <1, due to local saturation. Spatially observed ice optical depths cannot be directly related to column densities due to radiative transfer effects. Vertical snowlines will not be a clear transition due to the radially increasing height of the snowsurface, but their location may be constrained from observations using radiative transfer modeling. Radial snowlines are not really accesible. Not only the ice abundance, but also inclination, settling, grain size distribution and disk mass have strong impact on the observed ice absorption features in disks. Relative changes in ice abundance can be inferred from observations only if the source structure is well constrained △ Less

Submitted 8 May, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 18 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

The edge-on protoplanetary disk HH 48 NE I. Modeling the geometry and stellar parameters

Authors: J. A. Sturm, M. K. McClure, C. J. Law, D. Harsono, J. B. Bergner, E. Dartois, M. N. Drozdovskaya, S. Ioppolo, K. I. Öberg, M. E. Palumbo, Y. J. Pendleton, W. R. M. Rocha, H. Terada, R. G. Urso

Abstract: Context. Observations of edge-on disks are an important tool for constraining general protoplanetary disk properties that cannot be determined in any other way. However, most radiative transfer models cannot simultaneously reproduce the spectral energy distributions (SEDs) and resolved scattered light and submillimeter observations of these systems, due to the differences in geometry and dust prop… ▽ More Context. Observations of edge-on disks are an important tool for constraining general protoplanetary disk properties that cannot be determined in any other way. However, most radiative transfer models cannot simultaneously reproduce the spectral energy distributions (SEDs) and resolved scattered light and submillimeter observations of these systems, due to the differences in geometry and dust properties at different wavelengths. Aims. We simultaneously constrain the geometry of the edge-on protoplanetary disk HH 48 NE and the characteristics of the host star. HH 48 NE is part of the JWST early release science program Ice Age. This work serves as a stepping stone towards a better understanding of the disk physical structure and icy chemistry in this particular source. This kind of modeling lays the groundwork for studying other edge-on sources to be observed with the JWST. Methods. We fit a parameterized dust model to HH 48 NE by coupling the radiative transfer code RADMC-3D and an MCMC framework. The dust structure was fitted independently to a compiled SED, a scattered light image at 0.8 $μ$m and an ALMA dust continuum observation at 890 $μ$m. Results. We find that 90% of the dust mass in HH 48 NE is settled to the disk midplane, less than in average disks, and that the atmospheric layers of the disk contain exclusively large grains (0.3-10 $μ$m). The exclusion of small grains in the upper atmosphere likely has important consequences for the chemistry due to the deep penetration of high-energy photons. The addition of a relatively large cavity (ca. 50 au in radius) is necessary to explain the strong mid-infrared emission, and to fit the scattered light and continuum observations simultaneously. △ Less

Submitted 8 May, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 16 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 677, A17 (2023)

Theory and Implementation of the Quantum Approximate Optimization Algorithm: A Comprehensible Introduction and Case Study Using Qiskit and IBM Quantum Computers

Authors: Andreas Sturm

Abstract: The present tutorial aims to provide a comprehensible and easily accessible introduction into the theory and implementation of the famous Quantum Approximate Optimization Algorithm (QAOA). We lay our focus on practical aspects and step-by-step guide through the realization of a proof of concept quantum application based on a real-world use case. In every step we first explain the underlying theory… ▽ More The present tutorial aims to provide a comprehensible and easily accessible introduction into the theory and implementation of the famous Quantum Approximate Optimization Algorithm (QAOA). We lay our focus on practical aspects and step-by-step guide through the realization of a proof of concept quantum application based on a real-world use case. In every step we first explain the underlying theory and subsequently provide the implementation using IBM's Qiskit. In this way we provide a thorough understanding of the mathematical modelling and the (quantum) algorithms as well as the equally important knowledge how to properly write the code implementing those theoretical concepts. As another central aspect of this tutorial we provide extensive experiments on the 27 qubits state-of-the-art quantum computer ibmq_ehningen. From the discussion of these experiments we gain an overview on the current status of quantum computers and deduce which problem sizes can meaningfully be executed on today's hardware. △ Less

Submitted 23 January, 2023; originally announced January 2023.

An Ice Age JWST inventory of dense molecular cloud ices

Authors: M. K. McClure, W. R. M. Rocha, K. M. Pontoppidan, N. Crouzet, L. E. U. Chu, E. Dartois, T. Lamberts, J. A. Noble, Y. J. Pendleton, G. Perotti, D. Qasim, M. G. Rachid, Z. L. Smith, Fengwu Sun, Tracy L Beck, A. C. A. Boogert, W. A. Brown, P. Caselli, S. B. Charnley, Herma M. Cuppen, H. Dickinson, M. N. Drozdovskaya, E. Egami, J. Erkal, H. Fraser , et al. (17 additional authors not shown)

Abstract: Icy grain mantles are the main reservoir of the volatile elements that link chemical processes in dark, interstellar clouds with the formation of planets and composition of their atmospheres. The initial ice composition is set in the cold, dense parts of molecular clouds, prior to the onset of star formation. With the exquisite sensitivity of JWST, this critical stage of ice evolution is now acces… ▽ More Icy grain mantles are the main reservoir of the volatile elements that link chemical processes in dark, interstellar clouds with the formation of planets and composition of their atmospheres. The initial ice composition is set in the cold, dense parts of molecular clouds, prior to the onset of star formation. With the exquisite sensitivity of JWST, this critical stage of ice evolution is now accessible for detailed study. Here we show the first results of the Early Release Science program "Ice Age" that reveal the rich composition of these dense cloud ices. Weak ices, including, $^{13}$CO$_2$, OCN$^-$, $^{13}$CO, OCS, and COMs functional groups are now detected along two pre-stellar lines of sight. The $^{12}$CO$_2$ ice profile indicates modest growth of the icy grains. Column densities of the major and minor ice species indicate that ices contribute between 2 and 19% of the bulk budgets of the key C, O, N, and S elements. Our results suggest that the formation of simple and complex molecules could begin early in a water-ice rich environment. △ Less

Submitted 22 January, 2023; originally announced January 2023.

Comments: To appear in Nature Astronomy on January 23rd, 2023. 33 pages, 16 figures, 3 tables; includes extended and supplemental data sections. Part of the JWST Ice Age Early Release Science program's science enabling products. Enhanced spectra downloadable on Zenodo at the following DOI: 10.5281/zenodo.7501239

Heat kernel asymptotics for scaling limits of isoradial graphs

Authors: Simon Schwarz, Anja Sturm, Max Wardetzky

Abstract: We consider the asymptotics of the discrete heat kernel on isoradial graphs for the case where the time and the edge lengths tend to zero simultaneously. Depending on the asymptotic ratio between time and edge lengths, we show that two different regimes arise: (i) a Gaussian regime and (ii) a Poissonian regime, which resemble the short-time asymptotics of the heat kernel on (i) Euclidean spaces an… ▽ More We consider the asymptotics of the discrete heat kernel on isoradial graphs for the case where the time and the edge lengths tend to zero simultaneously. Depending on the asymptotic ratio between time and edge lengths, we show that two different regimes arise: (i) a Gaussian regime and (ii) a Poissonian regime, which resemble the short-time asymptotics of the heat kernel on (i) Euclidean spaces and (ii) graphs, respectively. △ Less

Submitted 23 January, 2024; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: v2: revised version, 19 pages

Contact process on a dynamical long range percolation

Authors: Marco Seiler, Anja Sturm

Abstract: In this paper we introduce a contact process on a dynamical long range percolation (CPDLP) defined on a complete graph $(V,\mathcal{E})$. A dynamical long range percolation is a Feller process defined on the edge set $\mathcal{E}$, which assigns to each edge the state of being open or closed independently. The state of an edge $e$ is updated at rate $v_e$ and is open after the update with probabil… ▽ More In this paper we introduce a contact process on a dynamical long range percolation (CPDLP) defined on a complete graph $(V,\mathcal{E})$. A dynamical long range percolation is a Feller process defined on the edge set $\mathcal{E}$, which assigns to each edge the state of being open or closed independently. The state of an edge $e$ is updated at rate $v_e$ and is open after the update with probability $p_e$ and closed otherwise. The contact process is then defined on top of this evolving random environment using only open edges for infection while recovery is independent of the background. First, we conclude that an upper invariant law exists and that the phase transitions of survival and non-triviality of the upper invariant coincide. We then formulate a comparison with a contact process with a specific infection kernel which acts as a lower bound. Thus, we obtain an upper bound for the critical infection rate. We also show that if the probability that an edge is open is low for all edges then the CPDLP enters an immunization phase, i.e. it will not survive regardless of the value of the infection rate. Furthermore, we show that on $V=\mathbb{Z}$ and under suitable conditions on the rates of the dynamical long range percolation the CPDLP will almost surely die out if the update speed converges to zero for any given infection rate $λ$. △ Less

Submitted 23 November, 2023; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: 31 pages, 11 figures, (revised version) Some minor change for clarification and corrections of typos. Article will soon appear in EJP

MSC Class: 60K35; 05C80; 82C22

Disentangling protoplanetary disk gas mass and carbon depletion through combined atomic and molecular tracers

Authors: J. A. Sturm, A. S. Booth, M. K. McClure, M. Leemker, E. F. van Dishoeck

Abstract: The total disk gas mass and elemental C, N, O composition of protoplanetary disks are crucial ingredients for our understanding of planet formation. Measuring the gas mass is complicated, since H$_2$ cannot be detected in the cold bulk of the disk and the elemental abundances with respect to hydrogen are degenerate with gas mass in all disk models. We present new NOEMA observations of CO, $^{13}$C… ▽ More The total disk gas mass and elemental C, N, O composition of protoplanetary disks are crucial ingredients for our understanding of planet formation. Measuring the gas mass is complicated, since H$_2$ cannot be detected in the cold bulk of the disk and the elemental abundances with respect to hydrogen are degenerate with gas mass in all disk models. We present new NOEMA observations of CO, $^{13}$CO, C$^{18}$O and optically thin C$^{17}$O $J$=2-1 lines, and use additional high angular resolution Atacama Large Millimeter Array millimeter continuum and CO data to construct a representative model of LkCa 15. The transitions that constrain the gas mass and carbon abundance most are C$^{17}$O 2-1, N${_2}$H$^+$ 3-2 and HD 1-0. Using these three molecules we find that the gas mass in the LkCa 15 disk is $M_\mathrm{g}=0.01 ^{+0.01}_{-0.004} M_{\odot}$, a factor of six lower than estimated before. The carbon abundance is C/H = ($3 \pm 1.5) \times10^{-5}$, implying a moderate depletion of elemental carbon by a factor of 3-9. All other analyzed transitions also agree with these numbers, within a modeling uncertainty of a factor of two. Using the resolved \ce{C2H} image we find a C/O ratio of $\sim$1, which is consistent with literature values of H$_2$O depletion in this disk. The lack of severe carbon depletion in the LkCa 15 disk is consistent with the young age of the disk, but contrasts with the higher depletions seen in older cold transition disks. Combining optically thin CO isotopologue lines with N$_2$H$^+$ is promising to break the degeneracy between gas mass and CO abundance. The moderate level of depletion for this source with a cold, but young disk, suggests that long carbon transformation timescales contribute to the evolutionary trend seen in the level of carbon depletion among disk populations, rather than evolving temperature effects and presence of dust traps alone. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 16 pages, 8 figures. Accepted for publication in A&A

Possible discovery of Calvera's supernova remnant

Authors: M. Arias, A. Botteon, C. G. Bassa, S. van der Jagt, R. J. van Weeren, S. P. O'Sullivan, Q. Bosschaart, R. S. Dullaart, M. J. Hardcastle, J. W. T. Hessels, T. Shimwell, M. M. Slob, J. A. Sturm, C. Tasse, N. C. M. A. Theijssen, J. Vink

Abstract: We report the discovery of a ring of low surface brightness radio emission around the Calvera pulsar, a high Galactic latitude, isolated neutron star, in the LOFAR Two-metre Sky Survey (LoTSS). It is centered at $α=14\mathrm{h}11\mathrm{m}12.6\mathrm{s}$, $δ=+79^\mathrm{o}23'15"$, has inner and outer radii of $14.2'$ and $28.4'$, and an integrated flux density at 144 MHz of $1.08\pm0.15$ Jy. The r… ▽ More We report the discovery of a ring of low surface brightness radio emission around the Calvera pulsar, a high Galactic latitude, isolated neutron star, in the LOFAR Two-metre Sky Survey (LoTSS). It is centered at $α=14\mathrm{h}11\mathrm{m}12.6\mathrm{s}$, $δ=+79^\mathrm{o}23'15"$, has inner and outer radii of $14.2'$ and $28.4'$, and an integrated flux density at 144 MHz of $1.08\pm0.15$ Jy. The ring center is offset by $4.9'$ from the location of the Calvera pulsar. H$α$ observations with the Isaac Newton Telescope show no coincident optical emission, but do show a small ($\sim20"$) optical structure internal to the ring. We consider three possible interpretations for the ring: that it is an H~II region, a supernova remnant (SNR), or an Odd Radio Circle (ORC). The positional coincidence of the ring, the pulsar, and an X-ray-emitting non-equilibrium ionisation plasma previously detected, lead us to prefer the SNR interpretation. If the source is indeed a SNR and its association with the Calvera pulsar is confirmed, then Calvera's SNR, or G118.4+37.0, will be one of few SNRs in the Galactic halo. △ Less

Submitted 28 July, 2022; originally announced July 2022.

Comments: Accepted in A&A

Rates of Convergence for Chains of Expansive Markov Operators

Authors: Neal Hermer, D. Russell Luke, Anja Sturm

Abstract: We provide conditions that guarantee local rates of convergence in distribution of iterated random functions that are not nonexpansive mappings in locally compact Hadamard spaces. Our results are applied to stochastic instances of common algorithms in optimization, stochastic tomography for X-FEL imaging, and a stochastic algorithm for the computation of Fréchet means in model spaces for phylogene… ▽ More We provide conditions that guarantee local rates of convergence in distribution of iterated random functions that are not nonexpansive mappings in locally compact Hadamard spaces. Our results are applied to stochastic instances of common algorithms in optimization, stochastic tomography for X-FEL imaging, and a stochastic algorithm for the computation of Fréchet means in model spaces for phylogenetic trees. △ Less

Submitted 31 October, 2023; v1 submitted 10 June, 2022; originally announced June 2022.

Comments: 29 pages, 75 references. Accepted to Transactions of Mathematics and its Applications. This is a more narrowly focused version of arXiv:2007.06479. arXiv admin note: text overlap with arXiv:2205.15897

MSC Class: 60J05; 46N10; 46N30; 65C40; 49J55

Nonexpansive Markov Operators and Random Function Iterations for Stochastic Fixed Point Problems

Authors: Neal Hermer, D. Russell Luke, Anja Sturm

Abstract: We study the convergence of random function iterations for finding an invariant measure of the corresponding Markov operator. We call the problem of finding such an invariant measure the stochastic fixed point problem. This generalizes earlier work studying the stochastic feasibility problem, namely, to find points that are, with probability 1, fixed points of the random functions. When no such po… ▽ More We study the convergence of random function iterations for finding an invariant measure of the corresponding Markov operator. We call the problem of finding such an invariant measure the stochastic fixed point problem. This generalizes earlier work studying the stochastic feasibility problem, namely, to find points that are, with probability 1, fixed points of the random functions. When no such points exist, the stochastic feasibility problem is called inconsistent, but still under certain assumptions, the more general stochastic fixed point problem has a solution and the random function iterations converge to an invariant measure for the corresponding Markov operator. We show how common structures in deterministic fixed point theory can be exploited to establish existence of invariant measures and convergence in distribution of the Markov chain. This framework specializes to many applications of current interest including, for instance, stochastic algorithms for large-scale distributed computation, and deterministic iterative procedures with computational error. The theory developed in this study provides a solid basis for describing the convergence of simple computational methods without the assumption of infinite precision arithmetic or vanishing computational errors. △ Less

Submitted 9 March, 2023; v1 submitted 31 May, 2022; originally announced May 2022.

Comments: 40 pages, 61 references, 2 figures. This is the revision of a narrowed version of arXiv:2007.06479 which could not be published in its entirety

MSC Class: 60J05; 46N10; 46N30; 65C40; 49J55

Journal ref: Journal of Convex Analysis 30 (2023), No. 4, 1073--1114

Gas temperature structure across transition disk cavities

Authors: M. Leemker, A. S. Booth, E. F. van Dishoeck, A. F. Pérez-Sánchez, J. Szulágyi, A. D. Bosman, S. Bruderer, S. Facchini, M. R. Hogerheijde, T. Paneque-Carreño, J. A. Sturm

Abstract: [Abridged] Most disks observed at high angular resolution show substructures. Knowledge about the gas surface density and temperature is essential to understand these. The aim of this work is to constrain the gas temperature and surface density in two transition disks: LkCa15 and HD 169142. We use new ALMA observations of the $^{13}$CO $J=6-5$ transition together with archival $J=2-1$ data of… ▽ More [Abridged] Most disks observed at high angular resolution show substructures. Knowledge about the gas surface density and temperature is essential to understand these. The aim of this work is to constrain the gas temperature and surface density in two transition disks: LkCa15 and HD 169142. We use new ALMA observations of the $^{13}$CO $J=6-5$ transition together with archival $J=2-1$ data of $^{12}$CO, $^{13}$CO and C$^{18}$O to observationally constrain the gas temperature and surface density. Furthermore, we use the thermochemical code DALI to model the temperature and density structure of a typical transition disk. The $6-5/2-1$ line ratio in LkCa15 constrains the gas temperature in the emitting layers inside the dust cavity to be up to 65 K, warmer than in the outer disk at 20-30 K. For the HD 169142, the peak brightness temperature constrains the gas in the dust cavity of HD 169142 to be 170 K, whereas that in the outer disk is only 100 K. Models also show that a more luminous central star, a lower abundance of PAHs and the absence of a dusty inner disk increase the temperature of the emitting layers and hence the line ratio in the gas cavity. The gas column density in the LkCa15 dust cavity drops by a factor >2 compared to the outer disk, with an additional drop of an order of magnitude inside the gas cavity at 10 AU. In the case of HD 169142, the gas column density drops by a factor of 200$-$500 inside the gas cavity, which could be due to a massive companion of several M$_{\mathrm{J}}$. The broad dust-depleted gas region from 10-68 AU for LkCa15 may imply several lower mass planets. This work demonstrates that knowledge of the gas temperature is important to determine the gas surface density and thus whether planets, and if so what kind of planets, are the most likely carving the dust cavities. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: Accepted for publication in Astronomy and astrophysics

Journal ref: A&A 663, A23 (2022)

Contact process in an evolving random environment

Authors: Marco Seiler, Anja Sturm

Abstract: In this paper we introduce a contact process in an evolving random environment (CPERE) on a connected and transitive graph with bounded degree, where we assume that this environment is described through an ergodic spin systems with finite range. We show that under a certain growth condition the phase transition of survival is independent of the initial configuration of the process. We study the in… ▽ More In this paper we introduce a contact process in an evolving random environment (CPERE) on a connected and transitive graph with bounded degree, where we assume that this environment is described through an ergodic spin systems with finite range. We show that under a certain growth condition the phase transition of survival is independent of the initial configuration of the process. We study the invariant laws of the CPERE and show that under aforementioned growth condition the phase transition for survival coincides with the phase transition of non-triviality of the upper invariant law. Furthermore, we prove continuity properties for the survival probability and derive equivalent conditions for complete convergence, in an analogous way as for the classical contact process. We then focus on the special case, where the evolving random environment is described through a dynamical percolation. We show that the contact process on a dynamical percolation on the $d$-dimensional integer dies out a.s. at criticality and complete convergence holds for all parameter choices. In the end we derive some comparison results between a dynamical percolation and ergodic spin systems with finite range such that we get bounds on the survival probability of a contact process in an evolving random environment and we determine in this case that complete convergence holds in a certain parameter regime. △ Less

Submitted 15 September, 2023; v1 submitted 30 March, 2022; originally announced March 2022.

Comments: 53 pages, 19 figures, (revised version) Some gaps in the proof of Proposition 2.6 were filled, which resulted in a stronger statement. Manuskript was accepted in EJP

MSC Class: 60K35; 05C80; 82C22

Efficient Random Walks on Riemannian Manifolds

Authors: Simon Schwarz, Michael Herrmann, Anja Sturm, Max Wardetzky

Abstract: According to a version of Donsker's theorem, geodesic random walks on Riemannian manifolds converge to the respective Brownian motion. From a computational perspective, however, evaluating geodesics can be quite costly. We therefore introduce approximate geodesic random walks based on the concept of retractions. We show that these approximate walks converge in distribution to the correct Brownian… ▽ More According to a version of Donsker's theorem, geodesic random walks on Riemannian manifolds converge to the respective Brownian motion. From a computational perspective, however, evaluating geodesics can be quite costly. We therefore introduce approximate geodesic random walks based on the concept of retractions. We show that these approximate walks converge in distribution to the correct Brownian motion as long as the geodesic equation is approximated up to second order. As a result we obtain an efficient algorithm for sampling Brownian motion on compact Riemannian manifolds. △ Less

Submitted 3 December, 2023; v1 submitted 2 February, 2022; originally announced February 2022.

Comments: 14 pages; v3: published version

MSC Class: 65C30; 60H35; 58J65

Journal ref: Foundations of Computational Mathematics (2023)

Tracing pebble drift and trapping using radial carbon depletion profiles in protoplanetary disks

Authors: J. A. Sturm, M. K. McClure, D. Harsono, S. Facchini, F. Long, M. Kama, E. A. Bergin, E. F. van Dishoeck

Abstract: The composition of planets may be largely determined by the chemical processing and accretion of icy pebbles in protoplanetary disks. Recent observations of protoplanetary disks hint at wide-spread depletion of gaseous carbon. The missing volatile carbon is likely frozen in CO and/or CO$_2$ ice on grains and locked into the disk through pebble trapping in pressure bumps or planetesimals. We presen… ▽ More The composition of planets may be largely determined by the chemical processing and accretion of icy pebbles in protoplanetary disks. Recent observations of protoplanetary disks hint at wide-spread depletion of gaseous carbon. The missing volatile carbon is likely frozen in CO and/or CO$_2$ ice on grains and locked into the disk through pebble trapping in pressure bumps or planetesimals. We present the results of the first successful ACA (Atacama Compact Array) [C I] $J$ = 1-0 mini-survey of seven protoplanetary disks. Using tailored azimuthally symmetric DALI (Dust And LInes) thermo-chemical disk models, supported by the [C I] $J$ = 1-0 and resolved CO isotopologue data, we determine the system-averaged elemental volatile carbon abundance in the outer disk of three sources. Six out of seven sources are detected in [C I] $J$ = 1-0 with ACA, four of which show a distinct disk component. Based on the modeling we find severe cold gaseous carbon depletion in the outer disk of DL Tau and moderate depletion in the outer disks of DR Tau and DO Tau. Combining the outer and inner disk carbon abundances, we demonstrate definitive evidence for radial drift in the disk of DL Tau, where the existence of multiple dust rings points to either short lived or leaky dust traps. We find dust locking in the compact and smooth disks of DO Tau and DR Tau, hinting at unresolved dust substructure. Comparing our results with stars of different ages and luminosities, we identify an observational evolutionary trend in gaseous carbon depletion that is consistent with dynamical models of CO depletion processes. Transport efficiency of solids in protoplanetary disks can significantly differ from what we expect based on the current resolved substructure in the continuum observations. This has important implications for our understanding of the impact of radial drift and pebble accretion on planetary compositions. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 23 pages, 19 figures, accepted to A&A

Journal ref: A&A 660, A126 (2022)

Laser-equipped gas reaction chamber for probing environmentally sensitive materials at near atomic scale

Authors: Heena Khanchandani, Ayman A. El-Zoka, Se-Ho Kim, Uwe Tezins, Dirk Vogel, Andreas Sturm, Dierk Raabe, Baptiste Gault, Leigh Stephenson

Abstract: Numerous metallurgical and materials science applications depend on quantitative atomic-scale characterizations of environmentally-sensitive materials and their transient states. Studying the effect upon materials subjected to thermochemical treatments in specific gaseous atmospheres is of central importance for specifically studying a material's resistance to certain oxidative or hydrogen environ… ▽ More Numerous metallurgical and materials science applications depend on quantitative atomic-scale characterizations of environmentally-sensitive materials and their transient states. Studying the effect upon materials subjected to thermochemical treatments in specific gaseous atmospheres is of central importance for specifically studying a material's resistance to certain oxidative or hydrogen environments. It is also important for investigating catalytic materials, direct reduction of an oxide, particular surface science reactions or nanoparticle fabrication routes. This manuscript realizes such experiments upon a thermochemical reaction chamber called the Reacthub and allows for transferring treated materials under cryogenic & ultrahigh vacuum workflow conditions for characterisation by either atom probe or scanning Xe+/electron microscopies. Two examples are discussed in the present study. One protocol was in the deuterium gas charging (25 kPa D2 at 200 °C) of a high-manganese twinning-induced-plasticity steel and characterization of the ingress and trapping of hydrogen at various features (grain boundaries in particular) in efforts to relate this to the steel's hydrogen embrittlement susceptibility. Deuterium was successfully detected after gas charging but most contrast originated from the complex ion FeOD+ signal and the feature may be an artefact. The second example considered the direct deuterium reduction (5 kPa D2 at 700 °C) of a single crystal wuestite sample, demonstrating that under a standard thermochemical treatment causes rapid reduction upon the nanoscale. Further studies are required for complete confidence about these phenomena, but these experiments successfully demonstrate that how an ex-situ thermochemical treatment can be realised that captures environmentally-sensitive transient states that can be analysed by atomic-scale by atom probe microscope. △ Less

Submitted 21 October, 2021; v1 submitted 26 July, 2021; originally announced July 2021.

A partial information decomposition for discrete and continuous variables

Authors: Kyle Schick-Poland, Abdullah Makkeh, Aaron J. Gutknecht, Patricia Wollstadt, Anja Sturm, Michael Wibral

Abstract: Conceptually, partial information decomposition (PID) is concerned with separating the information contributions several sources hold about a certain target by decomposing the corresponding joint mutual information into contributions such as synergistic, redundant, or unique information. Despite PID conceptually being defined for any type of random variables, so far, PID could only be quantified f… ▽ More Conceptually, partial information decomposition (PID) is concerned with separating the information contributions several sources hold about a certain target by decomposing the corresponding joint mutual information into contributions such as synergistic, redundant, or unique information. Despite PID conceptually being defined for any type of random variables, so far, PID could only be quantified for the joint mutual information of discrete systems. Recently, a quantification for PID in continuous settings for two or three source variables was introduced. Nonetheless, no ansatz has managed to both quantify PID for more than three variables and cover general measure-theoretic random variables, such as mixed discrete-continuous, or continuous random variables yet. In this work we will propose an information quantity, defining the terms of a PID, which is well-defined for any number or type of source or target random variable. This proposed quantity is tightly related to a recently developed local shared information quantity for discrete random variables based on the idea of shared exclusions. Further, we prove that this newly proposed information-measure fulfills various desirable properties, such as satisfying a set of local PID axioms, invariance under invertible transformations, differentiability with respect to the underlying probability density, and admitting a target chain rule. △ Less

Submitted 24 June, 2021; v1 submitted 23 June, 2021; originally announced June 2021.

Comments: 24 pages, 1 figure, Reason for replacement: Updated funding information

MSC Class: 94A15 (Primary) 28C15; 60B05; 28A50 (Secondary)

Dust dynamics in planet-driven spirals

Authors: Ardjan Sturm, Giovanni Rosotti, Carsten Dominik

Abstract: Context. Protoplanetary disks are known to host spiral features that are observed in scattered light, ALMA continuum and more recently in CO gas emission and gas dynamics. It is however unknown if spirals in gas and dust trace the same morphology. Aims. We study the morphology and amplitude of dusty spirals as function of Stokes number and the underlying mechanisms causing a difference from gas… ▽ More Context. Protoplanetary disks are known to host spiral features that are observed in scattered light, ALMA continuum and more recently in CO gas emission and gas dynamics. It is however unknown if spirals in gas and dust trace the same morphology. Aims. We study the morphology and amplitude of dusty spirals as function of Stokes number and the underlying mechanisms causing a difference from gas spirals. We then construct a model to relate the deviation from Keplerian rotation in the gas to a perturbation in surface density of gas and dust. Methods. We use FARGO-3D with dust implementation to numerically study the spirals, after which the results are interpreted using a semi-analytical model. This model is tested on observational data to predict the perturbation of the spiral in gas dynamics based on the continuum data. Results. We find that the pitch angle of a spiral does not differ significantly between gas and dust. The amplitude of the dust spiral decreases with Stokes number (St) and starts to fade out at a typical St > 0.1 as the dust becomes decoupled from the gas. The semianalytical model provides an accurate and fast representation of the spiral in surface density of the dust from the gas. We find a spiral in the TW Hya velocity residual map, never seen before, which is a feature in the vertical velocity and has a kink at the continuum gap, yielding strong evidence for a planet at 99 au. Conclusions. We built a model that gives an estimate of the underlying dynamics of dust in a spiral, which can serve as a proof of planetary origin of spirals and can be a probe for the Stokes number in the disk. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 11 pages, 8 figures, appendices, accepted to A&A 14 September 2020

Journal ref: A&A 643, A92 (2020)

Random Function Iterations for Stochastic Fixed Point Problems

Authors: Neal Hermer, D. Russell Luke, Anja Sturm

Abstract: We study the convergence of random function iterations for finding an invariant measure of the corresponding Markov operator. We call the problem of finding such an invariant measure the stochastic fixed point problem. This generalizes earlier work studying the stochastic feasibility problem}, namely, to find points that are, with probability 1, fixed points of the random functions [Hermer, Luke,… ▽ More We study the convergence of random function iterations for finding an invariant measure of the corresponding Markov operator. We call the problem of finding such an invariant measure the stochastic fixed point problem. This generalizes earlier work studying the stochastic feasibility problem}, namely, to find points that are, with probability 1, fixed points of the random functions [Hermer, Luke, Sturm, 2019]. When no such points exist, the stochastic feasibility problem is called inconsistent, but still under certain assumptions, the more general stochastic fixed point problem has a solution and the random function iterations converge to an invariant measure for the corresponding Markov operator. There are two major types of convergence: almost sure convergence of the iterates to a fixed point in the case of stochastic feasibility, and convergence in distribution more generally. We show how common structures in deterministic fixed point theory can be exploited to establish existence of invariant measures and convergence of the Markov chain. We show that weaker assumptions than are usually encountered in the analysis of Markov chains guarantee linear/geometric convergence. This framework specializes to many applications of current interest including, for instance, stochastic algorithms for large-scale distributed computation, and deterministic iterative procedures with computational error. The theory developed in this study provides a solid basis for describing the convergence of simple computational methods without the assumption of infinite precision arithmetic or vanishing computational errors. △ Less

Submitted 23 March, 2022; v1 submitted 13 July, 2020; originally announced July 2020.

Comments: 54 pages, 96 references

MSC Class: 60J05; 46N10; 46N30; 65C40; 49J55

A low-mass stellar companion to the young variable star RZ Psc

Authors: Grant M. Kennedy, Christian Ginski, Matthew A. Kenworthy, Myriam Benisty, Thomas Henning, Rob G. van Holstein, Quentin Kral, François Ménard, Julien Milli, Luis Henry Quiroga-Nuñez, Christian Rab, Tomas Stolker, Ardjan Sturm

Abstract: RZ Psc is a young Sun-like star with a bright and warm infrared excess that is occasionally dimmed significantly by circumstellar dust structures. Optical depth arguments suggest that the dimming events do not probe a typical sight line through the circumstellar dust, and are instead caused by structures that appear above an optically thick mid-plane. This system may therefore be similar to system… ▽ More RZ Psc is a young Sun-like star with a bright and warm infrared excess that is occasionally dimmed significantly by circumstellar dust structures. Optical depth arguments suggest that the dimming events do not probe a typical sight line through the circumstellar dust, and are instead caused by structures that appear above an optically thick mid-plane. This system may therefore be similar to systems where an outer disk is shadowed by material closer to the star. Here we report the discovery that RZ Psc hosts a 0.12$M_\odot$ companion at a projected separation of 23au. We conclude that the disk must orbit the primary star. While we do not detect orbital motion, comparison of the angle of linear polarization of the primary with the companion's on-sky position angle provides circumstantial evidence that the companion and disc may not share the same orbital plane. Whether the companion severely disrupts the disc, truncates it, or has little effect at all, will require further observations of both the companion and disc. △ Less

Submitted 28 May, 2020; originally announced May 2020.

Comments: MNRASL in press

Recursive tree processes and the mean-field limit of stochastic flows

Authors: Tibor Mach, Anja Sturm, Jan M. Swart

Abstract: Interacting particle systems can often be constructed from a graphical representation, by applying local maps at the times of associated Poisson processes. This leads to a natural coupling of systems started in different initial states. We consider interacting particle systems on the complete graph in the mean-field limit, i.e., as the number of vertices tends to infinity. We are not only interest… ▽ More Interacting particle systems can often be constructed from a graphical representation, by applying local maps at the times of associated Poisson processes. This leads to a natural coupling of systems started in different initial states. We consider interacting particle systems on the complete graph in the mean-field limit, i.e., as the number of vertices tends to infinity. We are not only interested in the mean-field limit of a single process, but mainly in how several coupled processes behave in the limit. This turns out to be closely related to recursive tree processes as studied by Aldous and Bandyopadyay in discrete time. We here develop an analogue theory for recursive tree processes in continuous time. We illustrate the abstract theory on an example of a particle system with cooperative branching. This yields an interesting new example of a recursive tree process that is not endogenous. △ Less

Submitted 18 March, 2020; v1 submitted 27 December, 2018; originally announced December 2018.

Comments: 68 pages, 4 figures, revised version with changed numbering of theorems

MSC Class: 82C22; 60J25; 60J80; 60K35

Trait-dependent branching particle systems with competition and multiple offspring

Authors: Gabriel Berzunza, Anja Sturm, Anita Winter

Abstract: In this work we model the dynamics of a population that evolves as a continuous time branching process with a trait structure and ecological interactions in form of mutations and competition between individuals. We generalize existing microscopic models by allowing individuals to have multiple offspring at a reproduction event. Furthermore, we allow the reproduction law to be influenced both by th… ▽ More In this work we model the dynamics of a population that evolves as a continuous time branching process with a trait structure and ecological interactions in form of mutations and competition between individuals. We generalize existing microscopic models by allowing individuals to have multiple offspring at a reproduction event. Furthermore, we allow the reproduction law to be influenced both by the trait type of the parent as well as by the mutant trait type. We look for tractable large population approximations. More precisely, under some natural assumption on the branching and mutation mechanisms, we establish a superprocess limit as solution of a well-posed martingale problem. Standard approaches do not apply in our case due to the lack of the branching property, which is a consequence of the dependency created by the competition between individuals. For showing uniqueness we therefore had to develop a generalization of Dawson's Girsanov Theorem that may be of independent interest. △ Less

Submitted 15 October, 2020; v1 submitted 28 August, 2018; originally announced August 2018.

Comments: 52

MSC Class: 60J80; 60J68; 60K35

Random Function Iterations for Consistent Stochastic Feasibility

Authors: Neal Hermer, D. Russell Luke, Anja Sturm

Abstract: We study the convergence of stochastic fixed point iterations in the consistent case (in the sense of Butnariu and Flåm (1995)) in several different settings, under decreasingly restrictive regularity assumptions of the fixed point mappings. The iterations are Markov chains and, for the purposes of this study, convergence is understood in very restrictive terms. We show that sufficient conditions… ▽ More We study the convergence of stochastic fixed point iterations in the consistent case (in the sense of Butnariu and Flåm (1995)) in several different settings, under decreasingly restrictive regularity assumptions of the fixed point mappings. The iterations are Markov chains and, for the purposes of this study, convergence is understood in very restrictive terms. We show that sufficient conditions for geometric (linear) convergence in expectation of stochastic projection algorithms presented in Nedić (2011), are in fact necessary for geometric (linear) convergence in expectation more generally of iterated random functions. △ Less

Submitted 21 September, 2018; v1 submitted 16 August, 2018; originally announced August 2018.

Comments: 29 pages, 4 figures

MSC Class: 60J05; 52A22; 49J55 (Primary) 49J53; 65K05 (Secondary)

Journal ref: Numerical Functional Analysis and Optimization 40(4):386--420 (2019)

Frontiers in Pigment Cell and Melanoma Research

Authors: Fabian V. Filipp, Stanca Birlea, Marcus W. Bosenberg, Douglas Brash, Pamela B. Cassidy, Suzie Chen, John August D'Orazio, Mayumi Fujita, Boon-Kee Goh, Meenhard Herlyn, Arup K. Indra, Lionel Larue, Sancy A. Leachman, Caroline Le Poole, Feng Liu-Smith, Prashiela Manga, Lluis Montoliu, David A. Norris, Yiqun Shellman, Keiran S. M. Smalley, Richard A. Spritz, Richard A. Sturm, Susan M. Swetter, Tamara Terzian, Kazumasa Wakamatsu , et al. (2 additional authors not shown)

Abstract: We identify emerging frontiers in clinical and basic research of melanocyte biology and its associated biomedical disciplines. We describe challenges and opportunities in clinical and basic research of normal and diseased melanocytes that impact current approaches to research in melanoma and the dermatological sciences. We focus on four themes: (1) clinical melanoma research, (2) basic melanoma re… ▽ More We identify emerging frontiers in clinical and basic research of melanocyte biology and its associated biomedical disciplines. We describe challenges and opportunities in clinical and basic research of normal and diseased melanocytes that impact current approaches to research in melanoma and the dermatological sciences. We focus on four themes: (1) clinical melanoma research, (2) basic melanoma research, (3) clinical dermatology, and (4) basic pigment cell research, with the goal of outlining current highlights, challenges, and frontiers associated with pigmentation and melanocyte biology. Significantly, this document encapsulates important advances in melanocyte and melanoma research including emerging frontiers in melanoma immunotherapy, medical and surgical oncology, dermatology, vitiligo, albinism, genomics and systems biology, epidemiology, pigment biophysics and chemistry, and evolution. △ Less

Submitted 24 August, 2018; v1 submitted 22 July, 2018; originally announced August 2018.

The Laplace project: an integrated suite for correlative atom probe tomography and electron microscopy under cryogenic and UHV conditions

Authors: Leigh T. Stephenson, Agnieszka Szczepaniak, Isabelle Mouton, Kristiane Ann Kathrin Rusitzka, Andrew J. Breen, Uwe Tezins, Andreas Sturm, Dirk Vogel, Yanhong Chang, Paraskevas Kontis, Alexander Rosenthal, Jeffrey D. Shepard, Urs Maier, Thomas F. Kelly, Dierk Raabe, Baptiste Gault

Abstract: We present sample transfer instrumentation and integrated protocols for the preparation and correlative characterization of environmentally-sensitive materials by both atom probe tomography and electron microscopy. Ultra-high vacuum cryogenic suitcases allow specimen transfer between preparation, processing and several imaging platforms without exposure to atmospheric contamination. For expedient… ▽ More We present sample transfer instrumentation and integrated protocols for the preparation and correlative characterization of environmentally-sensitive materials by both atom probe tomography and electron microscopy. Ultra-high vacuum cryogenic suitcases allow specimen transfer between preparation, processing and several imaging platforms without exposure to atmospheric contamination. For expedient transfers, we installed a fast-docking station equipped with a cryogenic pump upon three systems; two atom probes, a scanning electron microscope / Xe-plasma focused ion beam and a N$_2$-atmosphere glovebox. We also installed a plasma FIB with a solid-state cooling stage to reduce beam damage and contamination, through reducing chemical activity and with the cryogenic components as passive cryogenic traps. We demonstrate the efficacy of the new laboratory protocols by the successful preparation and transfer of two highly contamination- and temperature-sensitive samples - water and ice. Analysing pure magnesium atom probe data, we show that surface oxidation can be effectively suppressed using an entirely cryogenic protocol (during specimen preparation and during transfer). Starting with the cryogenically-cooled plasma FIB, we also prepared and transferred frozen ice samples while avoiding significant melting or sublimation, suggesting that we may be able to measure the nanostructure of other normally-liquid or soft materials. Isolated cryogenic protocols within the N$_2$ glove box demonstrate the absence of ice condensation suggesting that environmental control can commence from fabrication until atom probe analysis. △ Less

Submitted 28 May, 2018; originally announced May 2018.

Comments: 13 pages, 5 figures, to be submitted to PLOS ONE

The Algebraic Approach to Duality: An Introduction

Authors: Anja Sturm, Jan M. Swart, Florian Völlering

Abstract: This survey article gives an elementary introduction to the algebraic approach to Markov process duality, as opposed to the pathwise approach. In the algebraic approach, a Markov generator is written as the sum of products of simpler operators, which each have a dual with respect to some duality function. We discuss at length the recent suggestion by Giardinà, Redig, and others, that it may be a g… ▽ More This survey article gives an elementary introduction to the algebraic approach to Markov process duality, as opposed to the pathwise approach. In the algebraic approach, a Markov generator is written as the sum of products of simpler operators, which each have a dual with respect to some duality function. We discuss at length the recent suggestion by Giardinà, Redig, and others, that it may be a good idea to choose these simpler operators in such a way that they form an irreducible representation of some known Lie algebra. In particular, we collect the necessary background on representations of Lie algebras that is crucial for this approach. We also discuss older work by Lloyd and Sudbury on duality functions of product form and the relation between intertwining and duality. △ Less

Submitted 20 February, 2018; originally announced February 2018.

Comments: 69 pages. Lecture notes for a learning session at the program Genealogies of Interacting Particle Systems, July-August 2017, Institute for Mathematical Sciences, Singapore

MSC Class: Primary: 82C22; Secondary: 60K35; 17B10; 22E46

A new characterization of endogeny

Authors: Tibor Mach, Anja Sturm, Jan M. Swart

Abstract: Aldous and Bandyopadhyay have shown that each solution to a recursive distributional equation (RDE) gives rise to recursive tree process (RTP), which is a sort of Markov chain in which time has a tree-like structure and in which the state of each vertex is a random function of its descendants. If the state at the root is measurable with respect to the sigma field generated by the random functions… ▽ More Aldous and Bandyopadhyay have shown that each solution to a recursive distributional equation (RDE) gives rise to recursive tree process (RTP), which is a sort of Markov chain in which time has a tree-like structure and in which the state of each vertex is a random function of its descendants. If the state at the root is measurable with respect to the sigma field generated by the random functions attached to all vertices, then the RTP is said to be endogenous. For RTPs defined by continuous maps, Aldous and Bandyopadhyay showed that endogeny is equivalent to bivariate uniqueness, and they asked if the continuity hypothesis can be removed. We introduce a higher-level RDE that through its $n$-th moment measures contains all $n$-variate RDEs. We show that this higher-level RDE has minimal and maximal fixed points with respect to the convex order, and that these coincide if and only if the corresponding RTP is endogenous. As a side result, this allows us to answer the question of Aldous and Bandyopadhyay positively. △ Less

Submitted 1 October, 2018; v1 submitted 16 January, 2018; originally announced January 2018.

Comments: 17 pages. Corrected formula (6.12), part of which was accidentally deleted in the previous version

MSC Class: 60K35 (Primary); 60J05; 82C22; 60J80 (Secondary)

Coalescent results for diploid exchangeable population models

Authors: Matthias Birkner, Huili Liu, Anja Sturm

Abstract: We consider diploid bi-parental analogues of Cannings models: in a population of fixed size $N$ the next generation is composed of $V_{i,j}$ offspring from parents $i$ and $j$, where $V=(V_{i,j})_{1\le i\neq j \le N}$ is a (jointly) exchangeable (symmetric) array. Every individual carries two chromosome copies, each of which is inherited from one of its parents. We obtain general conditions, formu… ▽ More We consider diploid bi-parental analogues of Cannings models: in a population of fixed size $N$ the next generation is composed of $V_{i,j}$ offspring from parents $i$ and $j$, where $V=(V_{i,j})_{1\le i\neq j \le N}$ is a (jointly) exchangeable (symmetric) array. Every individual carries two chromosome copies, each of which is inherited from one of its parents. We obtain general conditions, formulated in terms of the vector of the total number of offspring to each individual, for the convergence of the properly scaled ancestral process for an $n$-sample of genes towards a ($Ξ$-)coalescent. This complements Möhle and Sagitov's (2001) result for the haploid case and sharpens the profile of Möhle and Sagitov's (2003) study of the diploid case, which focused on fixed couples, where each row of $V$ has at most one non-zero entry. We apply the convergence result to several examples, in particular to two diploid variations of Schweinsberg's (2003) model, leading to Beta-coalescents with two-fold and with four-fold mergers, respectively. △ Less

Submitted 27 March, 2018; v1 submitted 8 September, 2017; originally announced September 2017.

Comments: 41 pages, 1 figure

Pathwise duals of monotone and additive Markov processes

Authors: Anja Sturm, Jan M. Swart

Abstract: This paper develops a systematic treatment of monotonicity-based pathwise dualities for Markov processes taking values in partially ordered sets. We show that every Markov process that takes values in a finite partially ordered set and whose generator can be represented in monotone maps has a pathwise dual process. In the special setting of attractive spin systems this has been discovered earlier… ▽ More This paper develops a systematic treatment of monotonicity-based pathwise dualities for Markov processes taking values in partially ordered sets. We show that every Markov process that takes values in a finite partially ordered set and whose generator can be represented in monotone maps has a pathwise dual process. In the special setting of attractive spin systems this has been discovered earlier by Gray. We show that the dual simplifies a lot when the state space is a lattice (in the order-theoretic meaning of the word) and all monotone maps satisfy an additivity condition. This leads to a unified treatment of several well-known dualities, including Siegmund's dual for processes with a totally ordered state space, duality of additive spin systems, and a duality due to Krone for the two-stage contact process, and allows for the construction of new dualities as well. We show that the well-known representation of additive spin systems in terms of open paths in a graphical representation can be generalized to additive Markov processes taking values in general lattices, but for the process and its dual to be representable on the same underlying space, we need to assume that the lattice is distributive. In the final section, we show how our results can be generalized from finite state spaces to interacting particle systems with finite local state spaces. △ Less

Submitted 25 October, 2016; v1 submitted 21 October, 2015; originally announced October 2015.

Comments: 43 pages, 2 figures, final version, more extensive literature review, additional index and list of notation

MSC Class: 82C22 (Primary); 60J27; 06A06; 06B10 (Secondary)

Multi-type spatial branching models for local self-regulation I: Construction and an exponential duality

Authors: Andreas Greven, Anja Sturm, Anita Winter, Iljana Zähle

Abstract: We consider a spatial multi-type branching model in which individuals migrate in geographic space according to random walks and reproduce according to a state-dependent branching mechanism which can be sub-, super- or critical depending on the local intensity of individuals of the different types. The model is a Lotka-Volterra type model with a spatial component and is related to two models studie… ▽ More We consider a spatial multi-type branching model in which individuals migrate in geographic space according to random walks and reproduce according to a state-dependent branching mechanism which can be sub-, super- or critical depending on the local intensity of individuals of the different types. The model is a Lotka-Volterra type model with a spatial component and is related to two models studied in \cite{BlathEtheridgeMeredith2007} as well as to earlier work in \cite{Etheridge2004} and in \cite{NeuhauserPacala1999}. Our main focus is on the diffusion limit of small mass, locally many individuals and rapid reproduction. This system differs from spatial critical branching systems since it is not density preserving and the densities for large times do not depend on the initial distribution but mainly on the carrying capacities. We prove existence of the infinite particle model and the system of interacting diffusions as solutions of martingale problems or systems of stochastic equations. In the exchangeable case in which the parameters are not type dependent we show uniqueness of the solutions. For that purpose we establish a new exponential duality. △ Less

Submitted 14 September, 2015; originally announced September 2015.

MSC Class: 60K35

Limit laws of the empirical Wasserstein distance: Gaussian distributions

Authors: Thomas Rippl, Axel Munk, Anja Sturm

Abstract: We derive central limit theorems for the Wasserstein distance between the empirical distributions of Gaussian samples. The cases are distinguished whether the underlying laws are the same or different. Results are based on the (quadratic) Frechet differentiability of the Wasserstein distance in the Gaussian case. Extensions to elliptically symmetric distributions are discussed as well as several a… ▽ More We derive central limit theorems for the Wasserstein distance between the empirical distributions of Gaussian samples. The cases are distinguished whether the underlying laws are the same or different. Results are based on the (quadratic) Frechet differentiability of the Wasserstein distance in the Gaussian case. Extensions to elliptically symmetric distributions are discussed as well as several applications such as bootstrap and statistical testing. △ Less

Submitted 19 February, 2016; v1 submitted 15 July, 2015; originally announced July 2015.

Comments: 35 pages, 1 figure

MSC Class: 62E20; 60F05 (Primary); 58C20; 90B06 (Secondary)

On qualitative robustness of the Lotka--Nagaev estimator for the offspring mean of a supercritical Galton--Watson process

Authors: Dominic Schuhmacher, Anja Sturm, Henryk Zähle

Abstract: We characterize the sets of offspring laws on which the Lotka--Nagaev estimator for the mean of a supercritical Galton--Watson process is qualitatively robust. These are exactly the locally uniformly integrating sets of offspring laws, which may be quite large. If the corresponding global property is assumed instead, we obtain uniform robustness as well. We illustrate both results with a number… ▽ More We characterize the sets of offspring laws on which the Lotka--Nagaev estimator for the mean of a supercritical Galton--Watson process is qualitatively robust. These are exactly the locally uniformly integrating sets of offspring laws, which may be quite large. If the corresponding global property is assumed instead, we obtain uniform robustness as well. We illustrate both results with a number of concrete examples. As a by-product of the proof we obtain that the Lotka--Nagaev estimator is [locally] uniformly weakly consistent on the respective sets of offspring laws, conditionally on non-extinction. △ Less

Submitted 13 February, 2015; v1 submitted 15 September, 2014; originally announced September 2014.

MSC Class: 60J80; 62G05; 62G35

A particle system with cooperative branching and coalescence

Authors: Anja Sturm, Jan M. Swart

Abstract: In this paper, we introduce a one-dimensional model of particles performing independent random walks, where only pairs of particles can produce offspring ("cooperative branching"), and particles that land on an occupied site merge with the particle present on that site ("coalescence"). We show that the system undergoes a phase transition as the branching rate is increased. For small branching rate… ▽ More In this paper, we introduce a one-dimensional model of particles performing independent random walks, where only pairs of particles can produce offspring ("cooperative branching"), and particles that land on an occupied site merge with the particle present on that site ("coalescence"). We show that the system undergoes a phase transition as the branching rate is increased. For small branching rates, the upper invariant law is trivial, and the process started with finitely many particles a.s. ends up with a single particle. Both statements are not true for high branching rates. An interesting feature of the process is that the spectral gap is zero even for low branching rates. Indeed, if the branching rate is small enough, then we show that for the process started in the fully occupied state, the particle density decays as one over the square root of time, and the same is true for the decay of the probability that the process still has more than one particle at a later time if it started with two particles. △ Less

Submitted 28 May, 2015; v1 submitted 2 November, 2013; originally announced November 2013.

Comments: Published at http://dx.doi.org/10.1214/14-AAP1032 in the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AAP-AAP1032

Journal ref: Annals of Applied Probability 2015, Vol. 25, No. 3, 1616-1649

New results on pathwise uniqueness for the heat equation with colored noise

Authors: Thomas Rippl, Anja Sturm

Abstract: We consider strong uniqueness and thus also existence of strong solutions for the stochastic heat equation with a multiplicative colored noise term. Here, the noise is white in time and colored in q dimensional space ($q \geq 1$) with a singular correlation kernel. The noise coefficient is Hölder continuous in the solution. We discuss improvements of the sufficient conditions obtained in Mytnik, P… ▽ More We consider strong uniqueness and thus also existence of strong solutions for the stochastic heat equation with a multiplicative colored noise term. Here, the noise is white in time and colored in q dimensional space ($q \geq 1$) with a singular correlation kernel. The noise coefficient is Hölder continuous in the solution. We discuss improvements of the sufficient conditions obtained in Mytnik, Perkins and Sturm (2006) that relate the Hölder coefficient with the singularity of the correlation kernel of the noise. For this we use new ideas of Mytnik and Perkins (2011) who treat the case of strong uniqueness for the stochastic heat equation with multiplicative white noise in one dimension. Our main result on pathwise uniqueness confirms a conjecture that was put forward in their paper. △ Less

Submitted 18 December, 2012; originally announced December 2012.

Comments: arXiv admin note: text overlap with arXiv:0809.0248 by other authors

MSC Class: 60H15; 60K35; 60K37; 60J80; 60F05

On spatial coalescents with multiple mergers in two dimensions

Authors: Benjamin Heuer, Anja Sturm

Abstract: We consider the genealogy of a sample of individuals taken from a spatially structured population when the variance of the offspring distribution is relatively large. The space is structured into discrete sites of a graph G. If the population size at each site is large, spatial coalescents with multiple mergers, so called spatial Lambda-coalescents, for which ancestral lines migrate in space and c… ▽ More We consider the genealogy of a sample of individuals taken from a spatially structured population when the variance of the offspring distribution is relatively large. The space is structured into discrete sites of a graph G. If the population size at each site is large, spatial coalescents with multiple mergers, so called spatial Lambda-coalescents, for which ancestral lines migrate in space and coalesce according to some Lambda-coalescent mechanism, are shown to be appropriate approximations to the genealogy of a sample of individuals. We then consider as the graph G the two dimensional torus with side length 2L+1 and show that as L tends to infinity, and time is rescaled appropriately, the partition structure of spatial Lambda-coalescents of individuals sampled far enough apart converges to the partition structure of a non-spatial Kingman coalescent. From a biological point of view this means that in certain circumstances both the spatial structure as well as larger variances of the underlying offspring distribution are harder to detect from the sample. However, supplemental simulations show that for moderately large L the different structure is still evident. △ Less

Submitted 24 September, 2012; originally announced September 2012.

Comments: 29 pages, 6 figures

MSC Class: 60J25; 60K35; 92D25; 92D10

Subcritical contact processes seen from a typical infected site

Authors: Anja Sturm, Jan M. Swart

Abstract: What is the long-time behavior of the law of a contact process started with a single infected site, distributed according to counting measure on the lattice? This question is related to the configuration as seen from a typical infected site and gives rise to the definition of so-called eigenmeasures, which are possibly infinite measures on the set of nonempty configurations that are preserved unde… ▽ More What is the long-time behavior of the law of a contact process started with a single infected site, distributed according to counting measure on the lattice? This question is related to the configuration as seen from a typical infected site and gives rise to the definition of so-called eigenmeasures, which are possibly infinite measures on the set of nonempty configurations that are preserved under the dynamics up to a multiplicative constant. In this paper, we study eigenmeasures of contact processes on general countable groups in the subcritical regime. We prove that in this regime, the process has a unique spatially homogeneous eigenmeasure. As an application, we show that the exponential growth rate is continuously differentiable and strictly decreasing as a function of the recovery rate, and we give a formula for the derivative in terms of the eigenmeasures of the contact process and its dual. △ Less

Submitted 29 May, 2013; v1 submitted 21 October, 2011; originally announced October 2011.

Comments: Changed the organization of the proofs somewhat to more clearly make a link to classical results about quasi-invariant laws. 44 pages

MSC Class: 82C22 (Primary) 60K35; 82B43 (Secondary)

Tightness of voter model interfaces

Authors: Anja Sturm, Jan M. Swart

Abstract: Consider a long-range, one-dimensional voter model started with all zeroes on the negative integers and all ones on the positive integers. If the process obtained by identifying states that are translations of each other is positively recurrent, then it is said that the voter model exhibits interface tightness. In 1995, Cox and Durrett proved that one-dimensional voter models exhibit interface t… ▽ More Consider a long-range, one-dimensional voter model started with all zeroes on the negative integers and all ones on the positive integers. If the process obtained by identifying states that are translations of each other is positively recurrent, then it is said that the voter model exhibits interface tightness. In 1995, Cox and Durrett proved that one-dimensional voter models exhibit interface tightness if their infection rates have a finite third moment. Recently, Belhaouari, Mountford, and Valle have improved this by showing that a finite second moment suffices. The present paper gives a new short proof of this fact. We also prove interface tightness for a long range swapping voter model, which has a mixture of long range voter model and exclusion process dynamics. △ Less

Submitted 11 July, 2007; v1 submitted 29 June, 2007; originally announced June 2007.

Comments: Corrected a serious error in Lemma 3. Some additional insignificant corrections. 9 pages

MSC Class: 82C22 (Primary); 82C24; 82C41; 60K35 (Secondary)

Voter models with heterozygosity selection

Authors: Anja Sturm, Jan Swart

Abstract: This paper studies variations of the usual voter model that favor types that are locally less common. Such models are dual to certain systems of branching annihilating random walks that are parity preserving. For both the voter models and their dual branching annihilating systems we determine all homogeneous invariant laws, and we study convergence to these laws started from other initial laws. This paper studies variations of the usual voter model that favor types that are locally less common. Such models are dual to certain systems of branching annihilating random walks that are parity preserving. For both the voter models and their dual branching annihilating systems we determine all homogeneous invariant laws, and we study convergence to these laws started from other initial laws. △ Less

Submitted 25 January, 2008; v1 submitted 19 January, 2007; originally announced January 2007.

Comments: Published in at http://dx.doi.org/10.1214/07-AAP444 the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AAP-AAP444 MSC Class: 82C22 (Primary); 82C41; 60K35 (Secondary)

Journal ref: Annals of Applied Probability 2008, Vol. 18, No. 1, 59-99

The spatial $Λ$-coalescent

Authors: Vlada Limic, Anja Sturm

Abstract: This paper extends the notion of the $\la$-coalescent of Pitman (1999) to the spatial setting. The partition elements of the spatial $Λ$-coalescent migrate in a (finite) geographical space and may only coalesce if located at the same site of the space. We characterize the $Λ$-coalescents that come down from infinity, in an analogous way to Schweinsberg (2000). Surprisingly, all spatial coalescen… ▽ More This paper extends the notion of the $\la$-coalescent of Pitman (1999) to the spatial setting. The partition elements of the spatial $Λ$-coalescent migrate in a (finite) geographical space and may only coalesce if located at the same site of the space. We characterize the $Λ$-coalescents that come down from infinity, in an analogous way to Schweinsberg (2000). Surprisingly, all spatial coalescents that come down from infinity, also come down from infinity in a uniform way. This enables us to study space-time asymptotics of spatial $Λ$-coalescents on large tori in $d\ge 3$ dimensions. Our results generalize and strengthen those of Greven et al. (2005), who studied the spatial Kingman coalescent in this context. △ Less

Submitted 21 November, 2005; originally announced November 2005.

Comments: 29 pages

MSC Class: 60J25; 60K35