Showing 201–250 of 503 results for author: Nowak, M

Conductance through a helical state in an InSb nanowire

Authors: Jakob Kammhuber, Maja C Cassidy, Fei Pei, Michal P Nowak, Adriaan Vuik, Diana Car, Sèbastien R Plissard, Erik P A M Bakkers, Michael Wimmer, Leo P Kouwenhoven

Abstract: The motion of an electron and its spin are generally not coupled. However in a one dimensional material with strong spin-orbit interaction (SOI) a helical state may emerge at finite magnetic fields, where electrons of opposite spin will have opposite momentum. The existence of this helical state has applications for spin filtering and Cooper pair splitter devices and is an essential ingredient for… ▽ More The motion of an electron and its spin are generally not coupled. However in a one dimensional material with strong spin-orbit interaction (SOI) a helical state may emerge at finite magnetic fields, where electrons of opposite spin will have opposite momentum. The existence of this helical state has applications for spin filtering and Cooper pair splitter devices and is an essential ingredient for realizing topologically protected quantum computing using Majorana zero modes. Here we report electrical conductance measurements of a quantum point contact (QPC) formed in an indium antimonide nanowire as a function of magnetic field. At magnetic fields exceeding 3T, the $2e^2/h$ plateau shows a reentrant conductance feature towards $1e^2/h$ which increases linearly in width with magnetic field before enveloping the $1e^2/h$ plateau. Rotating the external magnetic field either parallel or perpendicular to the spin-orbit field allows us to clearly attribute this experimental signature to SOI. We compare our observations with a model of a QPC incorporating SOI and extract a spin-orbit energy of ~6.5meV, which is significantly stronger than the SO energy obtained by other methods. △ Less

Submitted 24 January, 2017; originally announced January 2017.

Journal ref: Nature Communications 8, 478 (2017)

Evolutionary games on cycles with strong selection

Authors: Philipp M. Altrock, Arne Traulsen, Martin A. Nowak

Abstract: Evolutionary games on graphs describe how strategic interactions and population structure determine evolutionary success, quantified by the probability that a single mutant takes over a population. Graph structures, compared to the well-mixed case, can act as amplifiers or suppressors of selection by increasing or decreasing the fixation probability of a beneficial mutant. Properties of the associ… ▽ More Evolutionary games on graphs describe how strategic interactions and population structure determine evolutionary success, quantified by the probability that a single mutant takes over a population. Graph structures, compared to the well-mixed case, can act as amplifiers or suppressors of selection by increasing or decreasing the fixation probability of a beneficial mutant. Properties of the associated mean fixation times can be more intricate, especially when selection is strong. The intuition is that fixation of a beneficial mutant happens fast (in a dominance game), that fixation takes very long (in a coexistence game), and that strong selection eliminates demographic noise. Here we show that these intuitions can be misleading in structured populations. We analyze mean fixation times on the cycle graph under strong frequency-dependent selection for two different microscopic evolutionary update rules (death-birth and birth-death). We establish exact analytical results for fixation times under strong selection, and show that there are coexistence games in which fixation occurs in time polynomial in population size. Depending on the underlying game, we observe inherence of demographic noise even under strong selection, if the process is driven by random death before selection for birth of an offspring (death-birth update). In contrast, if selection for an offspring occurs before random removal (birth-death update), strong selection can remove demographic noise almost entirely. △ Less

Submitted 23 January, 2017; originally announced January 2017.

Comments: 5 Figures. Accepted for publication (Physical Review E)

MSC Class: 60J10; 92D25

Journal ref: Physical Review E 95, 022407 (2017)

Spectra of large time-lagged correlation matrices from Random Matrix Theory

Authors: Maciej A. Nowak, Wojciech Tarnowski

Abstract: We analyze the spectral properties of large, time-lagged correlation matrices using the tools of random matrix theory. We compare predictions of the one-dimensional spectra, based on approaches already proposed in the literature. Employing the methods of free random variables and diagrammatic techniques, we solve a general random matrix problem, namely the spectrum of a matrix… ▽ More We analyze the spectral properties of large, time-lagged correlation matrices using the tools of random matrix theory. We compare predictions of the one-dimensional spectra, based on approaches already proposed in the literature. Employing the methods of free random variables and diagrammatic techniques, we solve a general random matrix problem, namely the spectrum of a matrix $\frac{1}{T}XAX^{\dagger}$, where $X$ is an $N\times T$ Gaussian random matrix and $A$ is \textit{any} $T\times T$, not necessarily symmetric (Hermitian) matrix. As a particular application, we present the spectral features of the large lagged correlation matrices as a function of the depth of the time-lag. We also analyze the properties of left and right eigenvector correlations for the time-lagged matrices. We positively verify our results by the numerical simulations. △ Less

Submitted 30 June, 2017; v1 submitted 20 December, 2016; originally announced December 2016.

Comments: 44 pages, 11 figures; v2 typos corrected, final version

Journal ref: J. Stat. Mech. (2017) 063405

Disorder in the Sachdev-Yee-Kitaev Model

Authors: Yizhuang Liu, Maciej A. Nowak, Ismail Zahed

Abstract: We give qualitative arguments for the mesoscopic nature of the Sachdev-Yee-Kitaev (SYK) model in the holographic regime with $q^2/N\ll 1$ with $N$ Majorana particles coupled by antisymmetric and random interactions of range $q$. Using a stochastic deformation of the SYK model, we show that its characteristic determinant obeys a viscid Burgers equation with a small spectral viscosity in the opposit… ▽ More We give qualitative arguments for the mesoscopic nature of the Sachdev-Yee-Kitaev (SYK) model in the holographic regime with $q^2/N\ll 1$ with $N$ Majorana particles coupled by antisymmetric and random interactions of range $q$. Using a stochastic deformation of the SYK model, we show that its characteristic determinant obeys a viscid Burgers equation with a small spectral viscosity in the opposite regime with $q/N=1/2$, in leading order. The stochastic evolution of the SYK model can be mapped onto that of random matrix theory, with universal Airy oscillations at the edges. A spectral hydrodynamical estimate for the relaxation of the collective modes is made. △ Less

Submitted 27 December, 2016; v1 submitted 15 December, 2016; originally announced December 2016.

Comments: 7 pages, 1 figure

Leveraging High Resolution Spectra to Understand Black Hole Spectra

Authors: Michael A. Nowak

Abstract: For the past 17 years, both XMM-Newton and Chandra have brought the powerful combination of high spatial and spectral resolution to the study of black hole systems. Each of these attributes requires special consideration--- in comparison to lower spatial resolution, CCD quality spectra--- when modeling observations obtained by these spacecraft. A good understanding of the high resolution spectra i… ▽ More For the past 17 years, both XMM-Newton and Chandra have brought the powerful combination of high spatial and spectral resolution to the study of black hole systems. Each of these attributes requires special consideration--- in comparison to lower spatial resolution, CCD quality spectra--- when modeling observations obtained by these spacecraft. A good understanding of the high resolution spectra is in fact required to model properly lower resolution CCD spectra, with the Reflection Grating Spectrometer (RGS) instrument on XMM-Newton maintaining the highest "figure of merit" at soft X-ray energies for all missions flying or currently planned for the next decade. Thanks to its even higher spectral resolution, the use of Chandra-High Energy Transmission Gratings (HETG), albeit with longer integration times, allows for one to bring further clarity to RGS studies. A further promising route for continued studies is the combination of high spectral resolution at soft X-rays, via RGS and/or HETG, with contemporaneous broadband coverage extending to hard X-rays (e.g., NuSTAR or INTEGRAL spectra). Such studies offer special promise for answering fundamental questions about accretion in black hole systems; however, they have received only moderate consideration to date. This may be due in part to the difficulty of analyzing high resolution spectra. In response, we must continue to develop software tools that make the analysis of high resolution X-ray spectra more accessible to the wider astrophysics community. △ Less

Submitted 6 December, 2016; originally announced December 2016.

Comments: 7 pages, 5 figures, to appear in Astronomische Nachrichten, presented at the XMM-Newton: The Next Decade conference, ESAC, Madrid, Spain, 9 - 11 May 2016

Mass-scaling as a method to constrain outflows and particle acceleration from low-luminosity accreting black holes

Authors: R. M. T. Connors, S. Markoff, M. A. Nowak, J. Neilsen, C. Ceccobello, P. Crumley, C. S. Froning, E. Gallo, J. E. Nip

Abstract: The `fundamental plane of black hole accretion' (FP), a relation between the radio luminosities ($L_R$), X-ray luminosities ($L_X$), and masses ($M_{BH}$) of hard/quiescent state black hole binaries and low-luminosity active galactic nuclei, suggests some aspects of black hole accretion may be scale invariant. However, key questions still exist concerning the relationship between the inflow/outflo… ▽ More The `fundamental plane of black hole accretion' (FP), a relation between the radio luminosities ($L_R$), X-ray luminosities ($L_X$), and masses ($M_{BH}$) of hard/quiescent state black hole binaries and low-luminosity active galactic nuclei, suggests some aspects of black hole accretion may be scale invariant. However, key questions still exist concerning the relationship between the inflow/outflow behaviour in the `classic' hard state and quiescence, which may impact this scaling. We show that the broadband spectra of A0620-00 and~\sgra~(the least luminous stellar mass/supermassive black holes on the FP) can be modelled simultaneously with a physically-motivated outflow-dominated model where the jet power and all distances are scaled by the black hole mass. We find we can explain the data of both A0620-00 and~\sgra~(in its non-thermal flaring state) in the context of two outflow-model scenarios: (1) a synchrotron-self-Compton dominated state in which the jet plasma reaches highly sub-equipartition conditions (for the magnetic field with respect to that of the radiating particles), and (2) a synchrotron dominated state in the fast-cooling regime in which particle acceleration occurs within the inner few gravitational radii of the black hole and plasma is close to equipartition. We show that it may be possible to further discriminate between models (1) and (2) through future monitoring of its submm/IR/X-ray emission, in particular via time lags between the variable emission in these bands. △ Less

Submitted 3 December, 2016; originally announced December 2016.

Comments: 19 pages (including appendix), 6 figures, 5 tables. Accepted for publication in MNRAS

Primordial Sex Facilitates the Emergence of Evolution

Authors: Sam Sinai, Jason Olejarz, Iulia A. Neagu, Martin A. Nowak

Abstract: Compartments are ubiquitous throughout biology, yet their importance stretches back to the origin of cells. In the context of origin of life, we assume that a protocell, a compartment enclosing functional components, requires $N$ components to be evolvable. We take interest in the timescale in which a minimal evolvable protocell is produced. We show that when protocells fuse and share information,… ▽ More Compartments are ubiquitous throughout biology, yet their importance stretches back to the origin of cells. In the context of origin of life, we assume that a protocell, a compartment enclosing functional components, requires $N$ components to be evolvable. We take interest in the timescale in which a minimal evolvable protocell is produced. We show that when protocells fuse and share information, the time to produce an evolvable protocell scales algebraically in $N$, in contrast to an exponential scaling in the absence of fusion. We discuss the implications of this result for origins of life, as well as other biological processes. △ Less

Submitted 7 December, 2017; v1 submitted 2 December, 2016; originally announced December 2016.

Comments: Second version with minor reference updates

Journal ref: Journal of The Royal Society Interface 15.139 (2018): 20180003

Reaching sub-milimag photometric precision on Beta Pictoris with a nanosat: the PicSat mission

Authors: M. Nowak, S. Lacour, V. Lapeyrère, L. David, A. Crouzier, C. Dufoing, H. Faiz, T. Lemoult, P. Trebuchet

Abstract: PicSat is a nanosatellite currently being developed to observe the transit of the giant planet \b{eta} Pictoris, expected some time between July 2017 and June 2018. The mission is based on a Cubesat architecture, with a small but ambitious 2 kg opto-mechanical payload specifically designed for high precision photometry. The satellite will be launched in early 2017, probably on a 600 km Sun synchro… ▽ More PicSat is a nanosatellite currently being developed to observe the transit of the giant planet \b{eta} Pictoris, expected some time between July 2017 and June 2018. The mission is based on a Cubesat architecture, with a small but ambitious 2 kg opto-mechanical payload specifically designed for high precision photometry. The satellite will be launched in early 2017, probably on a 600 km Sun synchronous orbit. The main objective of the mission is the constant monitoring of the brightness of Pic at an unprecedented combination of reliability and precision (200 ppm per hour, with interruptions of at most 30 minutes) to finely characterize the transiting exoplanet and detect exocomets in the Pictoris system. To achieve this difficult objective, the payload is designed with a 3.5 cm effective aperture telescope which injects the light in a single-mode optical fiber linked to an avalanche photodioode. A two-axis piezoelectric actuation system, driven by a tailor-made feedback loop control algorithm, is used to lock the fiber on the center of the star in the focal plane. These actuators complement the attitude determination and control system of the satellite to maintain the sub-arcsecond pointing accuracy required to reach the excellent level of photometric precision. Overall, the mission raises multiple very difficult challenges: high temperature stability of the avalanche detector (achieved with a thermoelectric colling device), high pointing accuracy and stability, and short timeframe for the development. △ Less

Submitted 25 October, 2016; originally announced October 2016.

Comments: 7 pages, 6 figures

Journal ref: Proceedings of the SPIE, Volume 9904, id. 99044L 7 pp. (2016)

Squared eigenvalue condition numbers and eigenvector correlations from the single ring theorem

Authors: Serban Belinschi, Maciej A. Nowak, Roland Speicher, Wojciech Tarnowski

Abstract: We extend the so-called "single ring theorem"[1], also known as the Haagerup-Larsen theorem[2], by showing that in the limit when the size of the matrix goes to infinity a particular correlator between left and right eigenvectors of the relevant non-hermitian matrix $X$, being the spectral density weighted by the squared eigenvalue condition number, is given by a simple formula involving only the… ▽ More We extend the so-called "single ring theorem"[1], also known as the Haagerup-Larsen theorem[2], by showing that in the limit when the size of the matrix goes to infinity a particular correlator between left and right eigenvectors of the relevant non-hermitian matrix $X$, being the spectral density weighted by the squared eigenvalue condition number, is given by a simple formula involving only the radial spectral cumulative distribution function of $X$. We show that this object allows to calculate the conditional expectation of the squared eigenvalue condition number. We give examples and we provide cross-check of the analytic prediction by the large scale numerics. △ Less

Submitted 22 December, 2016; v1 submitted 17 August, 2016; originally announced August 2016.

Comments: 9 pages, 2 figures

Journal ref: J. Phys. A: Math. Th. 50 (10) 105204 (2017)

Transparent Semiconductor-Superconductor Interface and Induced Gap in an Epitaxial Heterostructure Josephson Junction

Authors: M. Kjaergaard, H. J. Suominen, M. P. Nowak, A. R. Akhmerov, J. Shabani, C. J. Palmstrøm, F. Nichele, C. M. Marcus

Abstract: Measurement of multiple Andreev reflection (MAR) in a Josephson junction made from an InAs heterostructure with epitaxial aluminum is used to quantify the highly transparent semiconductor-superconductor interface, indicating near-unity transmission. The observed temperature dependence of MAR does not follow a conventional BCS form, but instead agrees with a model in which the density of states in… ▽ More Measurement of multiple Andreev reflection (MAR) in a Josephson junction made from an InAs heterostructure with epitaxial aluminum is used to quantify the highly transparent semiconductor-superconductor interface, indicating near-unity transmission. The observed temperature dependence of MAR does not follow a conventional BCS form, but instead agrees with a model in which the density of states in the quantum well acquires an effective induced gap, in our case 180 μeV, close to that of the epitaxial superconductor. Carrier density dependence of MAR is investigated using a depletion gate, revealing the subband structure of the semiconductor quantum well, consistent with magnetotransport experiment of the bare InAs performed on the same wafer. △ Less

Submitted 4 August, 2016; v1 submitted 14 July, 2016; originally announced July 2016.

Comments: Includes supplementary material

Report number: NBI QDEV 2016

Journal ref: Phys. Rev. Applied 7, 034029 (2017)

Games of multicellularity

Authors: Kamran Kaveh, Carl Veller, Martin A. Nowak

Abstract: Evolutionary game dynamics are often studied in the context of different population structures. Here we propose a new population structure that is inspired by simple multicellular life forms. In our model, cells reproduce but can stay together after reproduction. They reach complexes of a certain size, n, before producing single cells again. The cells within a complex derive payoff from an evoluti… ▽ More Evolutionary game dynamics are often studied in the context of different population structures. Here we propose a new population structure that is inspired by simple multicellular life forms. In our model, cells reproduce but can stay together after reproduction. They reach complexes of a certain size, n, before producing single cells again. The cells within a complex derive payoff from an evolutionary game by interacting with each other. The reproductive rate of cells is proportional to their payoff. We consider all two-strategy games. We study deterministic evolutionary dynamics with mutations, and derive exact conditions for selection to favor one strategy over another. Our main result has the same symmetry as the well-known sigma condition, which has been proven for stochastic game dynamics and weak selection. For a maximum complex size of n=2 our result holds for any intensity of selection. For n > 2 it holds for weak selection. As specific examples we study the prisoner's dilemma and hawk-dove games. Our model advances theoretical work on multicellularity by allowing for frequency-dependent interactions within groups. △ Less

Submitted 24 May, 2016; originally announced May 2016.

Comments: 44 pages, 11 figues

Journal ref: Journal of Theoretical Biology (2016), pp. 143-158

Evolutionary dynamics on any population structure

Authors: Benjamin Allen, Gabor Lippner, Yu-Ting Chen, Babak Fotouhi, Naghmeh Momeni, Martin A. Nowak, Shing-Tung Yau

Abstract: Evolution occurs in populations of reproducing individuals. The structure of a biological population affects which traits evolve. Understanding evolutionary game dynamics in structured populations is difficult. Precise results have been absent for a long time, but have recently emerged for special structures where all individuals have the same number of neighbors. But the problem of determining wh… ▽ More Evolution occurs in populations of reproducing individuals. The structure of a biological population affects which traits evolve. Understanding evolutionary game dynamics in structured populations is difficult. Precise results have been absent for a long time, but have recently emerged for special structures where all individuals have the same number of neighbors. But the problem of determining which trait is favored by selection in the natural case where the number of neighbors can vary, has remained open. For arbitrary selection intensity, the problem is in a computational complexity class which suggests there is no efficient algorithm. Whether there exists a simple solution for weak selection was unanswered. Here we provide, surprisingly, a general formula for weak selection that applies to any graph or social network. Our method uses coalescent theory and relies on calculating the meeting times of random walks. We can now evaluate large numbers of diverse and heterogeneous population structures for their propensity to favor cooperation. We can also study how small changes in population structure---graph surgery---affect evolutionary outcomes. We find that cooperation flourishes most in societies that are based on strong pairwise ties. △ Less

Submitted 22 December, 2016; v1 submitted 20 May, 2016; originally announced May 2016.

Comments: 68 pages, 10 figures

MSC Class: 92D15; 05C81

Journal ref: Nature 544 (2017) 227-230

The Goodness of Simultaneous Fits in ISIS

Authors: Matthias Kühnel, Sebastian Falkner, Christoph Grossberger, Ralf Ballhausen, Thomas Dauser, Fritz-Walter Schwarm, Ingo Kreykenbohm, Michael A. Nowak, Katja Pottschmidt, Carlo Ferrigno, Richard E. Rothschild, Silvia Martínez-Núñez, José Miguel Torrejón, Felix Fürst, Dmitry Klochkov, Rüdiger Staubert, Peter Kretschmar, Jörn Wilms

Abstract: In a previous work, we introduced a tool for analyzing multiple datasets simultaneously, which has been implemented into ISIS. This tool was used to fit many spectra of X-ray binaries. However, the large number of degrees of freedom and individual datasets raise an issue about a good measure for a simultaneous fit quality. We present three ways to check the goodness of these fits: we investigate t… ▽ More In a previous work, we introduced a tool for analyzing multiple datasets simultaneously, which has been implemented into ISIS. This tool was used to fit many spectra of X-ray binaries. However, the large number of degrees of freedom and individual datasets raise an issue about a good measure for a simultaneous fit quality. We present three ways to check the goodness of these fits: we investigate the goodness of each fit in all datasets, we define a combined goodness exploiting the logical structure of a simultaneous fit, and we stack the fit residuals of all datasets to detect weak features. These tools are applied to all RXTE-spectra from GRO 1008-57, revealing calibration features that are not detected significantly in any single spectrum. Stacking the residuals from the best-fit model for the Vela X-1 and XTE J1859+083 data evidences fluorescent emission lines that would have gone undetected otherwise. △ Less

Submitted 6 May, 2016; originally announced May 2016.

Comments: 6 pages, 6 figures

Journal ref: Acta Polytechnica 56(1):41-46, 2016

Fixation probabilities for any configuration of two strategies on regular graphs

Authors: Yu-Ting Chen, Alex McAvoy, Martin A. Nowak

Abstract: Population structure and spatial heterogeneity are integral components of evolutionary dynamics, in general, and of evolution of cooperation, in particular. Structure can promote the emergence of cooperation in some populations and suppress it in others. Here, we provide results for weak selection to favor cooperation on regular graphs for any configuration, meaning any arrangement of cooperators… ▽ More Population structure and spatial heterogeneity are integral components of evolutionary dynamics, in general, and of evolution of cooperation, in particular. Structure can promote the emergence of cooperation in some populations and suppress it in others. Here, we provide results for weak selection to favor cooperation on regular graphs for any configuration, meaning any arrangement of cooperators and defectors. Our results extend previous work on fixation probabilities of single, randomly placed mutants. We find that for any configuration cooperation is never favored for birth-death (BD) updating. In contrast, for death-birth (DB) updating, we derive a simple, computationally tractable formula for weak selection to favor cooperation when starting from any configuration containing any number of cooperators and defectors. This formula elucidates two important features: (i) the takeover of cooperation can be enhanced by the strategic placement of cooperators and (ii) adding more cooperators to a configuration can sometimes suppress the evolution of cooperation. These findings give a formal account for how selection acts on all transient states that appear in evolutionary trajectories. They also inform the strategic design of initial states in social networks to maximally promote cooperation. We also derive general results that characterize the interaction of any two strategies, not only cooperation and defection. △ Less

Submitted 5 December, 2016; v1 submitted 4 May, 2016; originally announced May 2016.

Comments: 28 pages; final version

Journal ref: Scientific Reports 6, 39181 (2016)

An Ultra-Fast X-ray Disk Wind in the Neutron Star Binary GX 340+0

Authors: J. M. Miller, J. Raymond, E. Cackett, V. Grinberg, M. Nowak

Abstract: We present a spectral analysis of a brief Chandra/HETG observation of the neutron star low-mass X-ray binary GX~340+0. The high-resolution spectrum reveals evidence of ionized absorption in the Fe K band. The strongest feature, an absorption line at approximately 6.9 keV, is required at the 5 sigma level of confidence via an F-test. Photoionization modeling with XSTAR grids suggests that the line… ▽ More We present a spectral analysis of a brief Chandra/HETG observation of the neutron star low-mass X-ray binary GX~340+0. The high-resolution spectrum reveals evidence of ionized absorption in the Fe K band. The strongest feature, an absorption line at approximately 6.9 keV, is required at the 5 sigma level of confidence via an F-test. Photoionization modeling with XSTAR grids suggests that the line is the most prominent part of a disk wind with an apparent outflow speed of v = 0.04c. This interpretation is preferred at the 4 sigma level over a scenario in which the line is H-like Fe XXVI at a modest red-shift. The wind may achieve this speed owing to its relatively low ionization, enabling driving by radiation pressure on lines; in this sense, the wind in GX 340+0 may be the stellar-mass equivalent of the flows in broad absorption line quasars (BALQSOs). If the gas has a unity volume filling factor, the mass ouflow rate in the wind is over 10^-5 Msun/year, and the kinetic power is nearly 10^39 erg/s (or, 5-6 times the radiative Eddington limit for a neutron star). However, geometrical considerations - including a small volume filling factor and low covering factor - likely greatly reduce these values. △ Less

Submitted 12 April, 2016; originally announced April 2016.

Comments: Accepted for publication in ApJ Letters

Chandra X-Ray Spectroscopy of the Focused Wind in the Cygnus X-1 System. II. The Nondip Spectrum in the Low/Hard State - Modulations with Orbital Phase

Authors: Ivica Miškovičová, Natalie Hell, Manfred Hanke, Michael A. Nowak, Katja Pottschmidt, Norbert S. Schulz, Victoria Grinberg, Refiz Duro, Oliwia K. Madej, Anne M. Lohfink, Jérôme Rodriguez, Marion Cadolle Bel, Arash Bodaghee, John A. Tomsick, Julia C. Lee, Gregory V. Brown, Jörn Wilms

Abstract: The accretion onto the black hole in the system HDE 226868/Cygnus X-1 is powered by the strong line driven stellar wind of the O-type donor star. We study the X-ray properties of the stellar wind in the hard state of Cyg X-1 as determined with data from the Chandra High Energy Transmission Gratings. Large density and temperature inhomogeneities are present in the wind, with a fraction of the wind… ▽ More The accretion onto the black hole in the system HDE 226868/Cygnus X-1 is powered by the strong line driven stellar wind of the O-type donor star. We study the X-ray properties of the stellar wind in the hard state of Cyg X-1 as determined with data from the Chandra High Energy Transmission Gratings. Large density and temperature inhomogeneities are present in the wind, with a fraction of the wind consisting of clumps of matter with higher density and lower temperature embedded in a photoionized gas. Absorption dips observed in the light curve are believed to be caused by these clumps. This work concentrates on the non-dip spectra as a function of orbital phase. The spectra show lines of H-like and He-like ions of S, Si, Na, Mg, Al and highly ionized Fe (Fe xvii-Fe xxiv). We measure velocity shifts, column densities, and thermal broadening of the line series. The excellent quality of these five observations allows us to investigate the orbital phase dependence of these parameters. We show that the absorber is located close to the black hole. Doppler shifted lines point at a complex wind structure in this region, while emission lines seen in some observations are from a denser medium than the absorber. The observed line profiles are phase dependent. Their shapes vary from pure, symmetric absorption at the superior conjunction to P Cygni profiles at the inferior conjunction of the black hole. △ Less

Submitted 1 April, 2016; originally announced April 2016.

Comments: A&A, accepted (26 pages, 18 figures)

Journal ref: A&A 590, A114 (2016)

Evolution of worker policing

Authors: Jason W. Olejarz, Benjamin Allen, Carl Veller, Raghavendra Gadagkar, Martin A. Nowak

Abstract: Workers in insect societies are sometimes observed to kill male eggs of other workers, a phenomenon known as worker policing. We perform a mathematical analysis of the evolutionary dynamics of policing. We investigate the selective forces behind policing for both dominant and recessive mutations for different numbers of matings of the queen. The traditional, relatedness-based argument suggests tha… ▽ More Workers in insect societies are sometimes observed to kill male eggs of other workers, a phenomenon known as worker policing. We perform a mathematical analysis of the evolutionary dynamics of policing. We investigate the selective forces behind policing for both dominant and recessive mutations for different numbers of matings of the queen. The traditional, relatedness-based argument suggests that policing evolves if the queen mates with more than two males, but does not evolve if the queen mates with a single male. We derive precise conditions for the invasion and stability of policing alleles. We find that the relatedness-based argument is not robust with respect to small changes in colony efficiency caused by policing. We also calculate evolutionarily singular strategies and determine when they are evolutionarily stable. We use a population genetics approach that applies to dominant or recessive mutations of any effect size. △ Less

Submitted 10 May, 2016; v1 submitted 16 March, 2016; originally announced March 2016.

Quantized conductance doubling and hard gap in a two-dimensional semiconductor-superconductor heterostructure

Authors: M. Kjaergaard, F. Nichele, H. J. Suominen, M. P. Nowak, M. Wimmer, A. R. Akhmerov, J. A. Folk, K. Flensberg, J. Shabani, C. J. Palmstrom, C. M. Marcus

Abstract: The prospect of coupling a two-dimensional (2D) semiconductor heterostructure to a superconductor opens new research and technology opportunities, including fundamental problems in mesoscopic superconductivity, scalable superconducting electronics, and new topological states of matter. For instance, one route toward realizing topological matter is by coupling a 2D electron gas (2DEG) with strong s… ▽ More The prospect of coupling a two-dimensional (2D) semiconductor heterostructure to a superconductor opens new research and technology opportunities, including fundamental problems in mesoscopic superconductivity, scalable superconducting electronics, and new topological states of matter. For instance, one route toward realizing topological matter is by coupling a 2D electron gas (2DEG) with strong spin-orbit interaction to an s-wave superconductor. Previous efforts along these lines have been hindered by interface disorder and unstable gating. Here, we report measurements on a gateable InGaAs/InAs 2DEG with patterned epitaxial Al, yielding multilayer devices with atomically pristine interfaces between semiconductor and superconductor. Using surface gates to form a quantum point contact (QPC), we find a hard superconducting gap in the tunneling regime, overcoming the soft-gap problem in 2D superconductor-semiconductor hybrid systems. With the QPC in the open regime, we observe a first conductance plateau at 4e^2/h, as expected theoretically for a normal-QPC-superconductor structure. The realization of a hard-gap semiconductor-superconductor system that is amenable to top-down processing provides a means of fabricating scalable multicomponent hybrid systems for applications in low-dissipation electronics and topological quantum information. △ Less

Submitted 29 November, 2016; v1 submitted 6 March, 2016; originally announced March 2016.

Comments: includes main text, supplementary information and code for simulations. Published version

Report number: NBI QDEV CMT 2016

Journal ref: Nat. Commun. 7, 12841 (2016)

Revealing the broad iron Kalpha line in Cygnus X-1 through simultaneous XMM-Newton, RXTE, and INTEGRAL observations

Authors: Refiz Duro, Thomas Dauser, Victoria Grinberg, Ivica Miškovičová, Jérôme Rodriguez, John Tomsick, Manfred Hanke, Katja Pottschmidt, Michael A. Nowak, Sonja Kreykenbohm, Marion Cadolle Bel, Arash Bodaghee, Anne Lohfink, Christopher S. Reynolds, Eckhard Kendziorra, Marcus G. F. Kirsch, Rüdiger Staubert, Jörn Wilms

Abstract: We report on the analysis of the broad Fe Kalpha line feature of Cygnus X-1 in the spectra of four simultaneous hard intermediate state observations made with the X-ray Multiple Mirror mission (XMM-Newton), the Rossi X-ray Timing Explorer (RXTE), and the International Gamma-Ray Astrophysics Laboratory (INTEGRAL). The high quality of the XMM-Newton data taken in the Modified Timing Mode of the EPIC… ▽ More We report on the analysis of the broad Fe Kalpha line feature of Cygnus X-1 in the spectra of four simultaneous hard intermediate state observations made with the X-ray Multiple Mirror mission (XMM-Newton), the Rossi X-ray Timing Explorer (RXTE), and the International Gamma-Ray Astrophysics Laboratory (INTEGRAL). The high quality of the XMM-Newton data taken in the Modified Timing Mode of the EPIC-pn camera provides a great opportunity to investigate the broadened Fe Kalpha reflection line at 6.4keV with a very high signal to noise ratio. The 4-500keV energy range is used to constrain the underlying continuum and the reflection at higher energies. We first investigate the data by applying a phenomenological model that consists of the sum of an exponentially cutoff power law and relativistically smeared reflection. Additionally, we apply a more physical approach and model the irradiation of the accretion disk directly from the lamp post geometry. All four observations show consistent values for the black hole parameters with a spin of $a\sim 0.9$, in agreement with recent measurements from reflection and disk continuum fitting. The inclination is found to be $i\sim30^\circ$, consistent with the orbital inclination and different from inclination measurements made during the soft state, which show a higher inclination. We speculate that the difference between the inclination measurements is due to changes in the inner region of the accretion disk. △ Less

Submitted 28 February, 2016; originally announced February 2016.

Comments: (25 pages, 16 figures, Astronomy and Astrophysics, accepted)

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

Chiral Random Matrix Model at Finite Chemical Potential: Characteristic Determinant and Edge Universality

Authors: Yizhuang Liu, Maciej A. Nowak, Ismail Zahed

Abstract: We derive an exact formula for the stochastic evolution of the characteristic determinant of a class of deformed Wishart matrices following from a chiral random matrix model of QCD at finite chemical potential. In the WKB approximation, the characteristic determinant describes a sharp droplet of eigenvalues that deforms and expands at large stochastic times. Beyond the WKB limit, the edges of the… ▽ More We derive an exact formula for the stochastic evolution of the characteristic determinant of a class of deformed Wishart matrices following from a chiral random matrix model of QCD at finite chemical potential. In the WKB approximation, the characteristic determinant describes a sharp droplet of eigenvalues that deforms and expands at large stochastic times. Beyond the WKB limit, the edges of the droplet are fuzzy and described by universal edge functions. At the chiral point, the characteristic determinant in the microscopic limit is universal. Remarkably, the physical chiral condensate at finite chemical potential may be extracted from current and quenched lattice Dirac spectra using the universal edge scaling laws, without having to solve the QCD sign problem. △ Less

Submitted 2 February, 2016; originally announced February 2016.

Comments: 16 pages, 4 figures

Comparing reactive and memory-one strategies of direct reciprocity

Authors: Seung Ki Baek, Hyeong-Chai Jeong, Christian Hilbe, Martin A. Nowak

Abstract: Direct reciprocity is a mechanism for the evolution of cooperation based on repeated interactions. When individuals meet repeatedly, they can use conditional strategies to enforce cooperative outcomes that would not be feasible in one-shot social dilemmas. Direct reciprocity requires that individuals keep track of their past interactions and find the right response. However, there are natural boun… ▽ More Direct reciprocity is a mechanism for the evolution of cooperation based on repeated interactions. When individuals meet repeatedly, they can use conditional strategies to enforce cooperative outcomes that would not be feasible in one-shot social dilemmas. Direct reciprocity requires that individuals keep track of their past interactions and find the right response. However, there are natural bounds on strategic complexity: Humans find it difficult to remember past interactions accurately, especially over long timespans. Given these limitations, it is natural to ask how complex strategies need to be for cooperation to evolve. Here, we study stochastic evolutionary game dynamics in finite populations to systematically compare the evolutionary performance of reactive strategies, which only respond to the co-player's previous move, and memory-one strategies, which take into account the own and the co-player's previous move. In both cases, we compare deterministic strategy and stochastic strategy spaces. For reactive strategies and small costs, we find that stochasticity benefits cooperation, because it allows for generous-tit-for-tat. For memory one strategies and small costs, we find that stochasticity does not increase the propensity for cooperation, because the deterministic rule of win-stay, lose-shift works best. For memory one strategies and large costs, however, stochasticity can augment cooperation. △ Less

Submitted 23 May, 2016; v1 submitted 28 January, 2016; originally announced January 2016.

Comments: 18 pages, 7 figures

Journal ref: Sci. Rep. 6, 25676 (2016)

Nulling Data Reduction and On-Sky Performance of the Large Binocular Telescope Interferometer

Authors: D. Defrère, P. M. Hinz, B. Mennesson, W. F. Hoffmann, R. Millan-Gabet, A. J. Skemer, V. Bailey, W. C. Danchi, E. C. Downey, O. Durney, P. Grenz, J. M. Hill, T. J. McMahon, M. Montoya, E. Spalding, A. Vaz, O. Absil, P. Arbo, H. Bailey, G. Brusa, G. Bryden, S. Esposito, A. Gaspar, C. A. Haniff, G. M. Kennedy , et al. (14 additional authors not shown)

Abstract: The Large Binocular Telescope Interferometer (LBTI) is a versatile instrument designed for high-angular resolution and high-contrast infrared imaging (1.5-13 microns). In this paper, we focus on the mid-infrared (8-13 microns) nulling mode and present its theory of operation, data reduction, and on-sky performance as of the end of the commissioning phase in March 2015. With an interferometric base… ▽ More The Large Binocular Telescope Interferometer (LBTI) is a versatile instrument designed for high-angular resolution and high-contrast infrared imaging (1.5-13 microns). In this paper, we focus on the mid-infrared (8-13 microns) nulling mode and present its theory of operation, data reduction, and on-sky performance as of the end of the commissioning phase in March 2015. With an interferometric baseline of 14.4 meters, the LBTI nuller is specifically tuned to resolve the habitable zone of nearby main-sequence stars, where warm exozodiacal dust emission peaks. Measuring the exozodi luminosity function of nearby main-sequence stars is a key milestone to prepare for future exoEarth direct imaging instruments. Thanks to recent progress in wavefront control and phase stabilization, as well as in data reduction techniques, the LBTI demonstrated in February 2015 a calibrated null accuracy of 0.05% over a three-hour long observing sequence on the bright nearby A3V star beta Leo. This is equivalent to an exozodiacal disk density of 15 to 30 zodi for a Sun-like star located at 10pc, depending on the adopted disk model. This result sets a new record for high-contrast mid-infrared interferometric imaging and opens a new window on the study of planetary systems. △ Less

Submitted 25 January, 2016; originally announced January 2016.

Comments: 17 pages, 18 figures (resubmitted to ApJ with referee's comments)

Ornstein-Uhlenbeck diffusion of hermitian and non-hermitian matrices - unexpected links

Authors: Jean-Paul Blaizot, Jacek Grela, Maciej A. Nowak, Wojciech Tarnowski, Piotr Warchoł

Abstract: We compare the Ornstein-Uhlenbeck process for the Gaussian Unitary Ensemble to its non-hermitian counterpart - for the complex Ginibre ensemble. We exploit the mathematical framework based on the generalized Green's functions, which involves a new, hidden complex variable, in comparison to the standard treatment of the resolvents. This new variable turns out to be crucial to understand the pattern… ▽ More We compare the Ornstein-Uhlenbeck process for the Gaussian Unitary Ensemble to its non-hermitian counterpart - for the complex Ginibre ensemble. We exploit the mathematical framework based on the generalized Green's functions, which involves a new, hidden complex variable, in comparison to the standard treatment of the resolvents. This new variable turns out to be crucial to understand the pattern of the evolution of non-hermitian systems. The new feature is the emergence of the coupling between the flow of eigenvalues and that of left/right eigenvectors. We analyze local and global equilibria for both systems. Finally, we highlight some unexpected links between both ensembles. △ Less

Submitted 4 January, 2016; v1 submitted 21 December, 2015; originally announced December 2015.

Comments: 21 pages, 4 figures, submitted to J. Stat. Phys

K-Stacker, a new way of detecting and characterizing exoplanets with high contrast imaging instruments

Authors: H. Le Coroller, M. Nowak, L. Arnold, K. Dohlen, T. Fusco, J. F. Sauvage, A. Vigan

Abstract: This year, a second generation of coronagraphs dedicated to high-contrast direct imaging of exoplanets is starting operations. Among them, SPHERE, installed at the focus of the UT3 Very Large Telescope, reaches unprecedented contrast ratios up to $10^{-6}$ -$ 10^{-7}$, using eXtreme Adaptive Optics and the Angular Differential Imaging (ADI) techniques. In this paper, we present a new method call… ▽ More This year, a second generation of coronagraphs dedicated to high-contrast direct imaging of exoplanets is starting operations. Among them, SPHERE, installed at the focus of the UT3 Very Large Telescope, reaches unprecedented contrast ratios up to $10^{-6}$ -$ 10^{-7}$, using eXtreme Adaptive Optics and the Angular Differential Imaging (ADI) techniques. In this paper, we present a new method called Keplerian-Stacker that improves the detection limit of high contrast instruments like SPHERE, by up to a factor of 10. It consists of observing a star on a long enough period to let a hypothetical planet around that star move along its orbit. Even if in each individual observation taken during one night, we do not detect anything, we show that it is possible, using an optimization algorithm, to re-center the images according to keplerian motions (ex: 10-100 images taken over a long period of typically 1-10 years) and detect planets otherwise unreachable. This method can be used in combination with the ADI technics (or possibly any other high contrast data reduction method) to improve the Signal to Noise Ratio in each individual image, and to further improve the global detection limit. It also directly provides orbital parameters of the detected planets, as a by-product of the optimization algorithm. △ Less

Submitted 21 October, 2015; originally announced October 2015.

Comments: Twenty years of giant exoplanets - Proceedings of the Haute Provence Observatory Colloquium, 5-9 October 2015 Edited by I. Boisse, O. Demangeon, F. Bouchy & L. Arnold

Counterexamples in theory of fractal dimension for fractal structures

Authors: Magdalena Nowak, Manuel Fernández-Martínez, Miguel Angel Sánchez-Granero

Abstract: Fractal dimension constitutes the main tool to test for fractal patterns in Euclidean contexts. For this purpose, it is always used the box dimension, since it is easy to calculate, though the Hausdorff dimension, which is the oldest and also the most accurate fractal dimension, presents the best analytical properties. Additionally, fractal structures provide an appropriate topological context whe… ▽ More Fractal dimension constitutes the main tool to test for fractal patterns in Euclidean contexts. For this purpose, it is always used the box dimension, since it is easy to calculate, though the Hausdorff dimension, which is the oldest and also the most accurate fractal dimension, presents the best analytical properties. Additionally, fractal structures provide an appropriate topological context where new models of fractal dimension for a fractal structure could be developed in order to generalize the classical models of fractal dimension. In this paper, we provide some counterexamples regarding these new models of fractal dimension in order to show the reader how they behave mathematically with respect to the classical models, and also to point out which features of such models can be exploited to powerful effect in applications. △ Less

Submitted 12 October, 2015; originally announced October 2015.

Counterexamples for IFS-attractors

Authors: Magdalena Nowak, Manuel Fernandez-Martinez

Abstract: In this paper, we deal with the part of Fractal Theory related to finite families of (weak) contractions, called iterated function systems (IFS, herein). An attractor is a compact set which remains invariant for such a family. Thus, we consider spaces homeomorphic to attractors of either IFS or weak IFS, as well, which we will refer to as Banach and topological fractals, respectively. We present a… ▽ More In this paper, we deal with the part of Fractal Theory related to finite families of (weak) contractions, called iterated function systems (IFS, herein). An attractor is a compact set which remains invariant for such a family. Thus, we consider spaces homeomorphic to attractors of either IFS or weak IFS, as well, which we will refer to as Banach and topological fractals, respectively. We present a collection of counterexamples in order to show that all the presented definitions are essential, though they are not equivalent in general. △ Less

Submitted 12 October, 2015; originally announced October 2015.

As above, so below: exploiting mass scaling in black hole accretion to break degeneracies in spectral interpretation

Authors: Sera Markoff, Michael Nowak, Elena Gallo, Robert Hynes, Jörn Wilms, Richard M. Plotkin, Dipankar Maitra, Catia V. Silva, Samia Drappeau

Abstract: Over the last decade, the evidence is mounting that several aspects of black hole accretion physics proceed in a mass-invariant way. One of the best examples of this scaling is the empirical "Fundamental Plane of Black Hole Accretion" relation linking mass, radio and X-ray luminosity over eight orders of magnitude in black hole mass. The currently favored theoretical interpretation of this relatio… ▽ More Over the last decade, the evidence is mounting that several aspects of black hole accretion physics proceed in a mass-invariant way. One of the best examples of this scaling is the empirical "Fundamental Plane of Black Hole Accretion" relation linking mass, radio and X-ray luminosity over eight orders of magnitude in black hole mass. The currently favored theoretical interpretation of this relation is that the physics governing power output in weakly accreting black holes depends more on relative accretion rate than on mass. In order to test this theory, we explore whether a mass-invariant approach can simultaneously explain the broadband spectral energy distributions from two black holes at opposite ends of the mass scale but at similar Eddington accretion fractions. We find that the same model, with the same value of several fitted physical parameters expressed in mass-scaling units to enforce self-similarity, can provide a good description of two datasets from V404 Cyg and M81*, a stellar and supermassive black hole, respectively. Furthermore, only one of several potential emission scenarios for the X-ray band is successful, suggesting it is the dominant process driving the Fundamental Plane relation at this accretion rate. This approach thus holds promise for breaking current degeneracies in the interpretation of black hole high-energy spectra, and for constructing better prescriptions of black hole accretion for use in various local and cosmological feedback applications. △ Less

Submitted 8 October, 2015; originally announced October 2015.

Comments: 6 pages, 5 Figures, 1 Table, accepted for publication in ApJL

Seeking for a fingerprint: analysis of point processes in actigraphy recording

Authors: E. Gudowska-Nowak, J. K. Ochab, K. Oles, E. Beldzik, D. R. Chialvo, A. Domagalik, M. Fafrowicz, T. Marek, M. A. Nowak, H. Oginska, J. Szwed, J. Tyburczyk

Abstract: Motor activity of humans displays complex temporal fluctuations which can be characterized by scale-invariant statistics, thus documenting that structure and fluctuations of such kinetics remain similar over a broad range of time scales. Former studies on humans regularly deprived of sleep or suffering from sleep disorders predicted change in the invariant scale parameters with respect to those re… ▽ More Motor activity of humans displays complex temporal fluctuations which can be characterized by scale-invariant statistics, thus documenting that structure and fluctuations of such kinetics remain similar over a broad range of time scales. Former studies on humans regularly deprived of sleep or suffering from sleep disorders predicted change in the invariant scale parameters with respect to those representative for healthy subjects. In this study we investigate the signal patterns from actigraphy recordings by means of characteristic measures of fractional point processes. We analyse spontaneous locomotor activity of healthy individuals recorded during a week of regular sleep and a week of chronic partial sleep deprivation. Behavioural symptoms of lack of sleep can be evaluated by analysing statistics of duration times during active and resting states, and alteration of behavioural organization can be assessed by analysis of power laws detected in the event count distribution, distribution of waiting times between consecutive movements and detrended fluctuation analysis of recorded time series. We claim that among different measures characterizing complexity of the actigraphy recordings and their variations implied by chronic sleep distress, the exponents characterizing slopes of survival functions in resting states are the most effective biomarkers distinguishing between healthy and sleep-deprived groups. △ Less

Submitted 26 February, 2016; v1 submitted 31 July, 2015; originally announced July 2015.

Comments: Communicated at UPON2015, 14-17 July 2015, Barcelona. 21 pages, 11 figures; updated: figures 4-7, text revised, expanded Sec. 1,3,5

How we move is universal: scaling in the average shape of human activity

Authors: Dante R. Chialvo, Ana Maria Gonzalez Torrado, Ewa Gudowska-Nowak, Jeremi K. Ochab, Pedro Montoya, Maciej A. Nowak, Enzo Tagliazucchi

Abstract: Human motor activity is constrained by the rhythmicity of the 24 hours circadian cycle, including the usual 12-15 hours sleep-wake cycle. However, activity fluctuations also appear over a wide range of temporal scales, from days to a few seconds, resulting from the concatenation of a myriad of individual smaller motor events. Furthermore, individuals present different propensity to wakefulness and… ▽ More Human motor activity is constrained by the rhythmicity of the 24 hours circadian cycle, including the usual 12-15 hours sleep-wake cycle. However, activity fluctuations also appear over a wide range of temporal scales, from days to a few seconds, resulting from the concatenation of a myriad of individual smaller motor events. Furthermore, individuals present different propensity to wakefulness and thus to motor activity throughout the circadian cycle. Are activity fluctuations across temporal scales intrinsically different, or is there a universal description encompassing them? Is this description also universal across individuals, considering the aforementioned variability? Here we establish the presence of universality in motor activity fluctuations based on the empirical study of a month of continuous wristwatch accelerometer recordings. We study the scaling of average fluctuations across temporal scales and determine a universal law characterized by critical exponents $α$, $τ$ and $1/μ$. Results are highly reminiscent of the universality described for the average shape of avalanches in systems exhibiting crackling noise. Beyond its theoretical relevance, the present results can be important for developing objective markers of healthy as well as pathological human motor behavior. △ Less

Submitted 22 November, 2015; v1 submitted 22 June, 2015; originally announced June 2015.

Comments: Communicated to the Granada Seminar, "Physics Meets the Social Sciences: Emergent cooperative phenomena, from bacterial to human group behavior". June 14-19, 2015. La Herradura, Spain

Journal ref: Papers in Physics, vol. 7, art. 070017 (2015)

Lord of the Rings: A Kinematic Distance to Circinus X-1 from a Giant X-Ray Light Echo

Authors: Sebastian Heinz, Michael G. Burton, Catherine Braiding, William N. Brandt, Peter G. Jonker, Paul Sell, Robert P. Fender, Michael A. Nowak, Norbert S. Schulz

Abstract: Circinus X-1 exhibited a bright X-ray flare in late 2013. Follow-up observations with Chandra and XMM-Newton from 40 to 80 days after the flare reveal a bright X-ray light echo in the form of four well-defined rings with radii from 5 to 13 arcminutes, growing in radius with time. The large fluence of the flare and the large column density of interstellar dust towards Circinus X-1 make this the lar… ▽ More Circinus X-1 exhibited a bright X-ray flare in late 2013. Follow-up observations with Chandra and XMM-Newton from 40 to 80 days after the flare reveal a bright X-ray light echo in the form of four well-defined rings with radii from 5 to 13 arcminutes, growing in radius with time. The large fluence of the flare and the large column density of interstellar dust towards Circinus X-1 make this the largest and brightest set of rings from an X-ray light echo observed to date. By deconvolving the radial intensity profile of the echo with the MAXI X-ray lightcurve of the flare we reconstruct the dust distribution towards Circinus X-1 into four distinct dust concentrations. By comparing the peak in scattering intensity with the peak intensity in CO maps of molecular clouds from the Mopra Southern Galactic Plane CO Survey we identify the two innermost rings with clouds at radial velocity ~ -74 km/s and ~ -81 km/s, respectively. We identify a prominent band of foreground photoelectric absorption with a lane of CO gas at ~ -32 km/s. From the association of the rings with individual CO clouds we determine the kinematic distance to Circinus X-1 to be $D_{Cir X-1} = 9.4^{+0.8}_{-1.0}$ kpc. This distance rules out earlier claims of a distance around 4 kpc, implies that Circinus X-1 is a frequent super-Eddington source, and places a lower limit of $Γ\gtrsim 22$ on the Lorentz factor and an upper limit of $θ_{jet} \lesssim 3^{\circ}$ on the jet viewing angle. △ Less

Submitted 19 June, 2015; originally announced June 2015.

Comments: 20 pages, 21 figures, Astrophysical Journal, in print

The complex accretion geometry of GX 339-4 as seen by NuSTAR and Swift

Authors: F. Fuerst, M. A. Nowak, J. A. Tomsick, J. M. Miller, S. Corbel, M. Bachetti, S. E. Boggs, F. E. Christensen, W. W. Craig, A. C. Fabian, P. Gandhi, V. Grinberg, C. J. Hailey, F. A. Harrison, E. Kara, J. A. Kennea, K. K. Madsen, K. Pottschmidt, D. Stern, D. J. Walton, J. Wilms, W. W. Zhang

Abstract: We present spectral analysis of five NuSTAR and Swift observations of GX 339-4 taken during a failed outburst in summer 2013. These observations cover Eddington luminosity fractions in the range ~0.9-6%. Throughout this outburst, GX 339-4 stayed in the hard state, and all five observations show similar X-ray spectra with a hard power-law with a photon index near 1.6 and significant contribution fr… ▽ More We present spectral analysis of five NuSTAR and Swift observations of GX 339-4 taken during a failed outburst in summer 2013. These observations cover Eddington luminosity fractions in the range ~0.9-6%. Throughout this outburst, GX 339-4 stayed in the hard state, and all five observations show similar X-ray spectra with a hard power-law with a photon index near 1.6 and significant contribution from reflection. Using simple reflection models we find unrealistically high iron abundances. Allowing for different photon indices for the continuum incident on the reflector relative to the underlying observed continuum results in a statistically better fit and reduced iron abundances. With a photon index around 1.3, the input power-law on the reflector is significantly harder than that which is directly observed. We study the influence of different emissivity profiles and geometries and consistently find an improvement when using separate photon indices. The inferred inner accretion disk radius is strongly model dependent, but we do not find evidence for a truncation radius larger than 100 r_g in any model. The data do not allow independent spin constraints but the results are consistent with the literature (i.e., a>0). Our best-fit models indicate an inclination angle in the range 40-60 degrees, consistent with limits on the orbital inclination but higher than reported in the literature using standard reflection models. The iron line around 6.4 keV is clearly broadened, and we detect a superimposed narrow core as well. This core originates from a fluorescence region outside the influence of the strong gravity of the black hole and we discuss possible geometries. △ Less

Submitted 3 June, 2015; originally announced June 2015.

Comments: 11 pages, 6 figures, 6 tables, plus 9 tables in the appendix. Submitted to ApJ

Investigating the nature of IGR J17454-2919 using X-ray and Near-Infrared observations

Authors: A. Paizis, M. A. Nowak, J. Rodriguez, A. Segreto, S. Chaty, A. Rau, J. Chenevez, M. Del Santo, J. Greiner, S. Schmidl

Abstract: IGR J17454-2919 is a hard X-ray transient discovered by INTEGRAL on 2014 September 27. We report on our 20ks Chandra observation of the source, performed about five weeks after the discovery, as well as on INTEGRAL and Swift monitoring long-term observations. X-ray broad-band spectra of the source are compatible with an absorbed power-law, $Γ\sim$1.6-1.8, ${\rm N_H}\sim$(10-12)… ▽ More IGR J17454-2919 is a hard X-ray transient discovered by INTEGRAL on 2014 September 27. We report on our 20ks Chandra observation of the source, performed about five weeks after the discovery, as well as on INTEGRAL and Swift monitoring long-term observations. X-ray broad-band spectra of the source are compatible with an absorbed power-law, $Γ\sim$1.6-1.8, ${\rm N_H}\sim$(10-12)$\times 10^{22}\,{\rm cm}^{-2}$, with no trace of a cut-off in the data up to about 100keV, and with an average absorbed 0.5-100keV flux of about (7.1-9.7)${\times 10^{-10}~erg~cm^{-2}~s^{-1}}$. With Chandra, we determine the most accurate X-ray position of IGR J17454-2919, $α_{J2000}$=17$^{h}$ 45$^{m}$ 27$^{s}$.69, $δ_{J2000}$= $-$29$^{\circ}$ 19$^{\prime}$ 53$^{\prime \prime}$.8 (90% uncertainty of 0$^{\prime\prime}$.6), consistent with the NIR source 2MASS J17452768-2919534. We also include NIR investigations from our observations of the source field on 2014 October 6 with GROND. With the multi-wavelength information at hand, we discuss the possible nature of IGR J17454-2919. △ Less

Submitted 3 June, 2015; originally announced June 2015.

Comments: 9 pages, 4 figures, accepted for publication in ApJ

NuSTAR and Suzaku observations of the hard state in Cygnus X-1: locating the inner accretion disk

Authors: M. L. Parker, J. A. Tomsick, J. M. Miller, K. Yamaoka, A. Lohfink, M. Nowak, A. C. Fabian, W. N. Alston, S. E. Boggs, F. E. Christensen, W. W. Craig, F. Fuerst, P. Gandhi, B. W. Grefenstette, V. Grinberg, C. J. Hailey, F. A. Harrison, E. Kara, A. L. King, D. Stern, D. J. Walton, J. Wilms, W. W. Zhang

Abstract: We present simultaneous Nuclear Spectroscopic Telescope Array (NuSTAR ) and Suzaku observations of the X-ray binary Cygnus X-1 in the hard state. This is the first time this state has been observed in Cyg X-1 with NuSTAR, which enables us to study the reflection and broad-band spectra in unprecedented detail. We confirm that the iron line cannot be fit with a combination of narrow lines and absorp… ▽ More We present simultaneous Nuclear Spectroscopic Telescope Array (NuSTAR ) and Suzaku observations of the X-ray binary Cygnus X-1 in the hard state. This is the first time this state has been observed in Cyg X-1 with NuSTAR, which enables us to study the reflection and broad-band spectra in unprecedented detail. We confirm that the iron line cannot be fit with a combination of narrow lines and absorption features, and instead requires a relativistically blurred profile in combination with a narrow line and absorption from the companion wind. We use the reflection models of Garcia et al. (2014) to simultaneously measure the black hole spin, disk inner radius, and coronal height in a self-consistent manner. Detailed fits to the iron line profile indicate a high level of relativistic blurring, indicative of reflection from the inner accretion disk. We find a high spin, a small inner disk radius, and a low source height, and rule out truncation to greater than three gravitational radii at the 3σ confidence level. In addition, we find that the line profile has not changed greatly in the switch from soft to hard states, and that the differences are consistent with changes in the underlying reflection spectrum rather than the relativistic blurring. We find that the blurring parameters are consistent when fitting either just the iron line or the entire broad-band spectrum, which is well modelled with a Comptonized continuum plus reflection model. △ Less

Submitted 29 May, 2015; originally announced June 2015.

Comments: 12 pages, 7 figures, accepted for publication in ApJ

Spatial model predicts dispersal and cell turnover cause reduced intra-tumor heterogeneity

Authors: Bartlomiej Waclaw, Ivana Bozic, Meredith E. Pittman, Ralph H. Hruban, Bert Vogelstein, Martin A. Nowak

Abstract: Most cancers in humans are large, measuring centimeters in diameter, composed of many billions of cells. An equivalent mass of normal cells would be highly heterogeneous as a result of the mutations that occur during each cell division. What is remarkable about cancers is their homogeneity - virtually every neoplastic cell within a large cancer contains the same core set of genetic alterations, wi… ▽ More Most cancers in humans are large, measuring centimeters in diameter, composed of many billions of cells. An equivalent mass of normal cells would be highly heterogeneous as a result of the mutations that occur during each cell division. What is remarkable about cancers is their homogeneity - virtually every neoplastic cell within a large cancer contains the same core set of genetic alterations, with heterogeneity confined to mutations that have emerged after the last clonal expansions. How such clones expand within the spatially-constrained three dimensional architecture of a tumor, and come to dominate a large, pre-existing lesion, has never been explained. We here describe a model for tumor evolution that shows how short-range migration and cell turnover can account for rapid cell mixing inside the tumor. With it, we show that even a small selective advantage of a single cell within a large tumor allows the descendants of that cell to replace the precursor mass in a clinically relevant time frame. We also demonstrate that the same mechanisms can be responsible for the rapid onset of resistance to chemotherapy. Our model not only provides novel insights into spatial and temporal aspects of tumor growth but also suggests that targeting short range cellular migratory activity could have dramatic effects on tumor growth rates. △ Less

Submitted 24 March, 2015; originally announced March 2015.

Comments: 37 pages, 14 figures

Unveiling the significance of eigenvectors in diffusing non-hermitian matrices by identifying the underlying Burgers dynamics

Authors: Zdzislaw Burda, Jacek Grela, Maciej A. Nowak, Wojciech Tarnowski, Piotr Warchoł

Abstract: Following our recent letter, we study in detail an entry-wise diffusion of non-hermitian complex matrices. We obtain an exact partial differential equation (valid for any matrix size $N$ and arbitrary initial conditions) for evolution of the averaged extended characteristic polynomial. The logarithm of this polynomial has an interpretation of a potential which generates a Burgers dynamics in quate… ▽ More Following our recent letter, we study in detail an entry-wise diffusion of non-hermitian complex matrices. We obtain an exact partial differential equation (valid for any matrix size $N$ and arbitrary initial conditions) for evolution of the averaged extended characteristic polynomial. The logarithm of this polynomial has an interpretation of a potential which generates a Burgers dynamics in quaternionic space. The dynamics of the ensemble in the large $N$ is completely determined by the coevolution of the spectral density and a certain eigenvector correlation function. This coevolution is best visible in an electrostatic potential of a quaternionic argument built of two complex variables, the first of which governs standard spectral properties while the second unravels the hidden dynamics of eigenvector correlation function. We obtain general large $N$ formulas for both spectral density and 1-point eigenvector correlation function valid for any initial conditions. We exemplify our studies by solving three examples, and we verify the analytic form of our solutions with numerical simulations. △ Less

Submitted 23 March, 2015; originally announced March 2015.

Comments: 24 pages, 11 figures

Journal ref: Nucl. Phys. B 897, 421 (2015)

Detecting topological and Banach fractals among zero-dimensional spaces

Authors: Taras Banakh, Magdalena Nowak, Filip Strobin

Abstract: A topological space $X$ is called a topological fractal if $X=\bigcup_{f\in\mathcal F}f(X)$ for a finite system $\mathcal F$ of continuous self-maps of $X$, which is topologically contracting in the sense that for every open cover $\mathcal U$ of $X$ there is a number $n\in\mathbb N$ such that for any functions $f_1,\dots,f_n\in \mathcal F$, the set $f_1\circ\dots\circ f_n(X)$ is contained in some… ▽ More A topological space $X$ is called a topological fractal if $X=\bigcup_{f\in\mathcal F}f(X)$ for a finite system $\mathcal F$ of continuous self-maps of $X$, which is topologically contracting in the sense that for every open cover $\mathcal U$ of $X$ there is a number $n\in\mathbb N$ such that for any functions $f_1,\dots,f_n\in \mathcal F$, the set $f_1\circ\dots\circ f_n(X)$ is contained in some set $U\in\mathcal U$. If, in addition, all functions $f\in\mathcal F$ have Lipschitz constant $<1$ with respect to some metric generating the topology of $X$, then the space $X$ is called a Banach fractal. It is known that each topological fractal is compact and metrizable. We prove that a zero-dimensional compact metrizable space $X$ is a topological fractal if and only if $X$ is a Banach fractal if and only if $X$ is either uncountable or $X$ is countable and its scattered height $\hbar(X)$ is a successor ordinal. For countable compact spaces this classification was recently proved by M.Nowak. △ Less

Submitted 31 March, 2015; v1 submitted 22 March, 2015; originally announced March 2015.

Comments: 7 pages

MSC Class: 28A80; 37C25; 37C70

Journal ref: Topology Appl. 196:A (2015), 22-30

Long term variability of Cygnus X-1: VII. Orbital variability of the focussed wind in Cyg X-1 / HDE 226868 system

Authors: V. Grinberg, M. A. Leutenegger, N. Hell, K. Pottschmidt, M. Böck, J. A. García, M. Hanke, M. A. Nowak, J. O. Sundqvist, R. H. D. Townsend, J. Wilms

Abstract: Binary systems with an accreting compact object are a unique chance to investigate the strong, clumpy, line-driven winds of early type supergiants by using the compact object's X-rays to probe the wind structure. We analyze the two-component wind of HDE 226868, the O9.7Iab giant companion of the black hole Cyg X-1 using 4.77 Ms of RXTE observations of the system taken over the course of 16 years.… ▽ More Binary systems with an accreting compact object are a unique chance to investigate the strong, clumpy, line-driven winds of early type supergiants by using the compact object's X-rays to probe the wind structure. We analyze the two-component wind of HDE 226868, the O9.7Iab giant companion of the black hole Cyg X-1 using 4.77 Ms of RXTE observations of the system taken over the course of 16 years. Absorption changes strongly over the 5.6 d binary orbit, but also shows a large scatter at a given orbital phase, especially at superior conjunction. The orbital variability is most prominent when the black hole is in the hard X-ray state. Our data are poorer for the intermediate and soft state, but show signs for orbital variability of the absorption column in the intermediate state. We quantitatively compare the data in the hard state to a toy model of a focussed Castor-Abbott-Klein-wind: as it does not incorporate clumping, the model does not describe the observations well. A qualitative comparison to a simplified simulation of clumpy winds with spherical clumps shows good agreement in the distribution of the equivalent hydrogen column density for models with a porosity length on the order of the stellar radius at inferior conjunction; we conjecture that the deviations between data and model at superior conjunction could be either due to lack of a focussed wind component in the model or a more complicated clump structure. △ Less

Submitted 25 February, 2015; originally announced February 2015.

Comments: proposed for acceptance in A&A, 11 pages, 11 figures (two in appendix)

Journal ref: A&A 576, A117 (2015)

Signal from noise retrieval from one and two-point Green's function - comparison

Authors: Zbigniew Drogosz, Jerzy Jurkiewicz, Grzegorz Łukaszewski, Maciej A. Nowak

Abstract: We compare two methods of eigen-inference from large sets of data, based on the analysis of one-point and two-point Green's functions, respectively. Our analysis points at the superiority of eigen-inference based on one-point Green's function. First, the applied by us method based on Pad?e approximants is orders of magnitude faster comparing to the eigen-inference based on uctuations (two-point Gr… ▽ More We compare two methods of eigen-inference from large sets of data, based on the analysis of one-point and two-point Green's functions, respectively. Our analysis points at the superiority of eigen-inference based on one-point Green's function. First, the applied by us method based on Pad?e approximants is orders of magnitude faster comparing to the eigen-inference based on uctuations (two-point Green's functions). Second, we have identified the source of potential instability of the two-point Green's function method, as arising from the spurious zero and negative modes of the estimator for a variance operator of the certain multidimensional Gaussian distribution, inherent for the two-point Green's function eigen-inference method. Third, we have presented the cases of eigen-inference based on negative spectral moments, for strictly positive spectra. Finally, we have compared the cases of eigen-inference of real-valued and complex-valued correlated Wishart distributions, reinforcing our conclusions on an advantage of the one-point Green's function method. △ Less

Submitted 9 January, 2015; originally announced January 2015.

Comments: 14 pages, 8 figures, 3 tables

Journal ref: Phys. Rev. E 92, 022111 (2015)

The X-ray Flux Distribution of Sagittarius A* as Seen by Chandra

Authors: J. Neilsen, S. Markoff, M. A. Nowak, J. Dexter, G. Witzel, N. Barrière, Y. Li, F. K. Baganoff, N. Degenaar, P. C. Fragile, C. Gammie, A. Goldwurm, N. Grosso, D. Haggard

Abstract: We present a statistical analysis of the X-ray flux distribution of Sgr A* from the Chandra X-ray Observatory's 3 Ms Sgr A* X-ray Visionary Project (XVP) in 2012. Our analysis indicates that the observed X-ray flux distribution can be decomposed into a steady quiescent component, represented by a Poisson process with rate $Q=(5.24\pm0.08)\times10^{-3}$ cts s$^{-1},$ and a variable component, repre… ▽ More We present a statistical analysis of the X-ray flux distribution of Sgr A* from the Chandra X-ray Observatory's 3 Ms Sgr A* X-ray Visionary Project (XVP) in 2012. Our analysis indicates that the observed X-ray flux distribution can be decomposed into a steady quiescent component, represented by a Poisson process with rate $Q=(5.24\pm0.08)\times10^{-3}$ cts s$^{-1},$ and a variable component, represented by a power law process ($dN/dF\propto F^{-ξ},$ $ξ=1.92_{-0.02}^{+0.03}$). This slope matches our recently-reported distribution of flare luminosities. The variability may also be described by a log-normal process with a median unabsorbed 2-8 keV flux of $1.8^{+0.9}_{-0.6}\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$ and a shape parameter $σ=2.4\pm0.2,$ but the power law provides a superior description of the data. In this decomposition of the flux distribution, all of the intrinsic X-ray variability of Sgr A* (spanning at least three orders of magnitude in flux) can be attributed to flaring activity, likely in the inner accretion flow. We confirm that at the faint end, the variable component contributes ~10% of the apparent quiescent flux, as previously indicated by our statistical analysis of X-ray flares in these Chandra observations. Our flux distribution provides a new and important observational constraint on theoretical models of Sgr A*, and we use simple radiation models to explore the extent to which a statistical comparison of the X-ray and infrared can provide insights into the physics of the X-ray emission mechanism. △ Less

Submitted 9 December, 2014; originally announced December 2014.

Comments: 13 pages, 7 figures, accepted for publication in ApJ. Comments welcome

Spatial heterogeneity in drug concentrations can facilitate the emergence of resistance to cancer therapy

Authors: Feng Fu, Martin A. Nowak, Sebastian Bonhoeffer

Abstract: Acquired resistance is one of the major barriers to successful cancer therapy. The development of resistance is commonly attributed to genetic heterogeneity. However, heterogeneity of drug penetration of the tumor microenvironment both on the microscopic level within solid tumors as well as on the macroscopic level across metastases may also contribute to acquired drug resistance. Here we use math… ▽ More Acquired resistance is one of the major barriers to successful cancer therapy. The development of resistance is commonly attributed to genetic heterogeneity. However, heterogeneity of drug penetration of the tumor microenvironment both on the microscopic level within solid tumors as well as on the macroscopic level across metastases may also contribute to acquired drug resistance. Here we use mathematical models to investigate the effect of drug heterogeneity on the probability of escape from treatment and time to resistance. Specifically we address scenarios with sufficiently efficient therapies that suppress growth of all preexisting genetic variants in the compartment with highest drug concentration. To study the joint effect of drug heterogeneity, growth rate, and evolution of resistance we analyze a multitype stochastic branching process describing growth of cancer cells in two compartments with different drug concentration and limited migration between compartments. We show that resistance is more likely to arise first in the low drug compartment and from there populate the high drug compartment. Moreover, we show that only below a threshold rate of cell migration does spatial heterogeneity accelerate resistance evolution, otherwise deterring drug resistance with excessively high migration rates. Our results provide new insights into understanding why cancers tend to quickly become resistant, and that cell migration and the presence of sanctuary sites with little drug exposure are essential to this end. △ Less

Submitted 24 November, 2014; originally announced November 2014.

Comments: Collaborations on further followup work are extremely welcome. Please see contact details at http://www.tb.ethz.ch/people/person-detail.html?persid=189998

Single-electron shell occupation and effective $g$-factor in few-electron nanowire quantum dots

Authors: M. P. Nowak, B. Szafran

Abstract: Nanowire double quantum dots occupied by an even number of electrons are investigated in the context of energy level structure revealed by electric dipole spin resonance measurements. We use numerically exact configuration interaction approach up to 6 electrons for systems tuned to Pauli spin blockade regime. We point out the differences between the spectra of systems with two and a greater number… ▽ More Nanowire double quantum dots occupied by an even number of electrons are investigated in the context of energy level structure revealed by electric dipole spin resonance measurements. We use numerically exact configuration interaction approach up to 6 electrons for systems tuned to Pauli spin blockade regime. We point out the differences between the spectra of systems with two and a greater number of electrons. For two-electrons the unequal length of the dots results in a different effective $g$-factors in the dots as observed by the recent experiments. For an increased number of electrons the $g$-factor difference between the dots appears already for symmetric systems and it is greatly amplified when the dots are of unequal length. We find that in the lowest energy states all except of two electrons form a frozen core of singlet shells that can be separated away for description of the properties of the system. △ Less

Submitted 27 January, 2015; v1 submitted 23 November, 2014; originally announced November 2014.

Journal ref: Phys. Rev. B 91, 085102 (2015)

Assembly of repetitive regions using next-generation sequencing data

Authors: Robert M. Nowak

Abstract: High read depth can be used to assemble short sequence repeats. The existing genome assemblers fail in repetitive regions of longer than average read. I propose a new algorithm for a DNA assembly which uses the relative frequency of reads to properly reconstruct repetitive sequences. The mathematical model shows the upper limits of accuracy of the results as a function of read coverage. For high… ▽ More High read depth can be used to assemble short sequence repeats. The existing genome assemblers fail in repetitive regions of longer than average read. I propose a new algorithm for a DNA assembly which uses the relative frequency of reads to properly reconstruct repetitive sequences. The mathematical model shows the upper limits of accuracy of the results as a function of read coverage. For high coverage, the estimation error depends linearly on repetitive sequence length and inversely proportional to the sequencing coverage. The algorithm requires high read depth, provided by the next-generation sequencers and could use the existing data. The tests on errorless reads, generated in silico from several model genomes, pointed the properly reconstructed repetitive sequences, where existing assemblers fail. △ Less

Submitted 3 November, 2014; originally announced November 2014.

Comments: 11 pages, 5 figures, 6 tables. The C++ sources, the Python scripts and the additional data are available at http://dnaasm.sourceforge.org

Shell-Shocked: The Interstellar Medium Near Cygnus X-1

Authors: P. H. Sell, S. Heinz, E. Richards, T. J. Maccarone, D. M. Russell, E. Gallo, R. Fender, S. Markoff, M. Nowak

Abstract: We conduct a detailed case-study of the interstellar shell near the high-mass X-ray binary, Cygnus X-1. We present new WIYN optical spectroscopic and Chandra X-ray observations of this region, which we compare with detailed MAPPINGS III shock models, to investigate the outflow powering the shell. Our analysis places improved, physically motivated constraints on the nature of the shockwave and the… ▽ More We conduct a detailed case-study of the interstellar shell near the high-mass X-ray binary, Cygnus X-1. We present new WIYN optical spectroscopic and Chandra X-ray observations of this region, which we compare with detailed MAPPINGS III shock models, to investigate the outflow powering the shell. Our analysis places improved, physically motivated constraints on the nature of the shockwave and the interstellar medium (ISM) it is plowing through. We find that the shock is traveling at less than a few hundred km/s through a low-density ISM (< 5 cm^-3). We calculate a robust, 3 sigma upper limit to the total, time-averaged power needed to drive the shockwave and inflate the bubble, < 2 x 10^38 erg/s. We then review possible origins of the shockwave. We find that a supernova origin to the shockwave is unlikely and that the black hole jet and/or O-star wind can both be central drivers of the shockwave. We conclude that the source of the Cygnus X-1 shockwave is far from solved. △ Less

Submitted 30 October, 2014; originally announced October 2014.

Comments: 16 pages, 8 figures, accepted for publication in MNRAS

Spitzer/infrared spectrograph investigation of MIPSGAL 24 μm compact bubbles : Low resolution observations

Authors: M. Nowak, N. Flagey, A. Noriega-Crespo, N. Billot, S. J. Carey, R. Paladini, S. D. Van Dyk

Abstract: We present Spitzer/IRS low resolution observations of 11 compact circumstellar bubbles from the MIPSGAL 24 μm Galactic Plane Survey. We find that this set of MIPSGAL bubbles (MBs) is divided into two categories, and that this distinction correlates with the morphologies of the MBs in the mid- IR. The four MBs with central sources in the mid-IR exhibit dust-rich, low excitation spectra, and their 2… ▽ More We present Spitzer/IRS low resolution observations of 11 compact circumstellar bubbles from the MIPSGAL 24 μm Galactic Plane Survey. We find that this set of MIPSGAL bubbles (MBs) is divided into two categories, and that this distinction correlates with the morphologies of the MBs in the mid- IR. The four MBs with central sources in the mid-IR exhibit dust-rich, low excitation spectra, and their 24 μm emission is accounted for by the dust continuum. The seven MBs without central sources in the mid-IR have spectra dominated by high excitation gas lines (e.g., [O IV] 26.0 μm, [Ne V] 14.3 and 24.3 μm, [Ne III] 15.5 μm), and the [O IV] line accounts for 50 to almost 100% of the 24 μm emission in five of them. In the dust-poor MBs, the [Ne V] and [Ne III] line ratios correspond to high excitation conditions. Based on comparisons with published IRS spectra, we suggest that the dust-poor MBs are highly excited planetary nebulae with peculiar white dwarfs (e.g., [WR], novae) at their centers. The central stars of the four dust-rich MBs are all massive star candidates. Dust temperatures range from 40 to 100 K in the outer shells. We constrain the extinction along the lines of sight from the IRS spectra. We then derive distance, dust masses, and dust production rate estimates for these objects. These estimates are all consistent with the nature of the central stars. We summarize the identifications of MBs made to date and discuss the correlation between their mid-IR morphologies and natures. Candidate Be/B[e]/LBV and WR stars are mainly "rings" with mid-IR central sources, whereas PNe are mostly "disks" without mid-IR central sources. Therefore we expect that most of the 300 remaining unidentified MBs will be classified as PNe. △ Less

Submitted 22 October, 2014; originally announced October 2014.

Comments: 22 pages, 10 figures, accepted for publication in ApJ

The molecular clock of neutral evolution can be accelerated or slowed by asymmetric spatial structure

Authors: Benjamin Allen, Christine Sample, Yulia A. Dementieva, Ruben C. Medeiros, Christopher Paoletti, Martin A. Nowak

Abstract: Over time, a population acquires neutral genetic substitutions as a consequence of random drift. A famous result in population genetics asserts that the rate, $K$, at which these substitutions accumulate in the population coincides with the mutation rate, $u$, at which they arise in individuals: $K=u$. This identity enables genetic sequence data to be used as a "molecular clock" to estimate the ti… ▽ More Over time, a population acquires neutral genetic substitutions as a consequence of random drift. A famous result in population genetics asserts that the rate, $K$, at which these substitutions accumulate in the population coincides with the mutation rate, $u$, at which they arise in individuals: $K=u$. This identity enables genetic sequence data to be used as a "molecular clock" to estimate the timing of evolutionary events. While the molecular clock is known to be perturbed by selection, it is thought that $K=u$ holds very generally for neutral evolution. Here we show that asymmetric spatial population structure can alter the molecular clock rate for neutral mutations, leading to either $K<u$ or $K>u$. Deviations from $K=u$ occur because mutations arise unequally at different sites and have different probabilities of fixation depending on where they arise. If birth rates are uniform across sites, then $K \leq u$. In general, $K$ can take any value between 0 and $Nu$. Our model can be applied to a variety of population structures. In one example, we investigate the accumulation of genetic mutations in the small intestine. In another application, we analyze over 900 Twitter networks to study the effect of network topology on the fixation of neutral innovations in social evolution. △ Less

Submitted 18 September, 2014; originally announced September 2014.

Comments: 31 pages, 9 figures

MSC Class: 92D15

Journal ref: PLoS Comput Biol (2015) 11(2): e1004108

The Large Area Detector of LOFT: the Large Observatory for X-ray Timing

Authors: S. Zane, D. Walton, T. Kennedy, M. Feroci, J. -W. Den Herder, M. Ahangarianabhari, A. Argan, P. Azzarello, G. Baldazzi, M. Barbera, D. Barret, G. Bertuccio, P. Bodin, E. Bozzo, L. Bradley, F. Cadoux, P. Cais, R. Campana, J. Coker, A. Cros, E. Del Monte, A. De Rosa, S. Di Cosimo, I. Donnarumma, Y. Evangelista , et al. (34 additional authors not shown)

Abstract: LOFT (Large Observatory for X-ray Timing) is one of the five candidates that were considered by ESA as an M3 mission (with launch in 2022-2024) and has been studied during an extensive assessment phase. It is specifically designed to perform fast X-ray timing and probe the status of the matter near black holes and neutron stars. Its pointed instrument is the Large Area Detector (LAD), a 10 m 2 -cl… ▽ More LOFT (Large Observatory for X-ray Timing) is one of the five candidates that were considered by ESA as an M3 mission (with launch in 2022-2024) and has been studied during an extensive assessment phase. It is specifically designed to perform fast X-ray timing and probe the status of the matter near black holes and neutron stars. Its pointed instrument is the Large Area Detector (LAD), a 10 m 2 -class instrument operating in the 2-30keV range, which holds the capability to revolutionise studies of variability from X-ray sources on the millisecond time scales. The LAD instrument has now completed the assessment phase but was not down-selected for launch. However, during the assessment, most of the trade-offs have been closed leading to a robust and well documented design that will be re- proposed in future ESA calls. In this talk, we will summarize the characteristics of the LAD design and give an overview of the expectations for the instrument capabilities. △ Less

Submitted 27 August, 2014; originally announced August 2014.

Comments: Proc. SPIE 9144, Space Telescopes and Instrumentation 2014: Ultraviolet to Gamma Ray, 91442W

The Large Observatory For x-ray Timing

Authors: M. Feroci, J. W. den Herder, E. Bozzo, D. Barret, S. Brandt, M. Hernanz, M. van der Klis, M. Pohl, A. Santangelo, L. Stella, A. Watts, J. Wilms, S. Zane, M. Ahangarianabhari, C. Albertus, M. Alford, A. Alpar, D. Altamirano, L. Alvarez, L. Amati, C. Amoros, N. Andersson, A. Antonelli, A. Argan, R. Artigue , et al. (320 additional authors not shown)

Abstract: The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost… ▽ More The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study. △ Less

Submitted 29 August, 2014; v1 submitted 27 August, 2014; originally announced August 2014.

Comments: Proc. SPIE 9144, Space Telescopes and Instrumentation 2014: Ultraviolet to Gamma Ray, 91442T

Universality of fixation probabilities in randomly structured populations

Authors: Ben Adlam, Martin A. Nowak

Abstract: The stage of evolution is the population of reproducing individuals. The structure of the population is know to affect the dynamics and outcome of evolutionary processes, but analytical results for generic random structures have been lacking. The most general result so far, the isothermal theorem, assumes the propensity for change in each position is exactly the same, but realistic biological stru… ▽ More The stage of evolution is the population of reproducing individuals. The structure of the population is know to affect the dynamics and outcome of evolutionary processes, but analytical results for generic random structures have been lacking. The most general result so far, the isothermal theorem, assumes the propensity for change in each position is exactly the same, but realistic biological structures are always subject to variation and noise. We consider a population of finite size $n$ under constant selection whose structure is given by a wide variety of weighted, directed, random graphs; vertices represent individuals and edges interactions between individuals. By establishing a robustness result for the isothermal theorem and using large deviation estimates to understand the typical structure of random graphs, we prove that for a generalization of the Erdős-Rényi model the fixation probability of an invading mutant is approximately the same as that of a mutant of equal fitness in a well-mixed population with high probability. Simulations of perturbed lattices, small-world networks, and scale-free networks behave similarly. We conjecture that the fixation probability in a well-mixed population, $(1-r^{-1})/(1-r^{-n})$, is universal: for many random graph models, the fixation probability approaches the above function uniformly as the graphs become large. △ Less

Submitted 9 July, 2014; originally announced July 2014.

Spectral density of generalized Wishart matrices and free multiplicative convolution

Authors: Wojciech Mlotkowski, Maciej A. Nowak, Karol A. Penson, Karol Zyczkowski

Abstract: We investigate the level density for several ensembles of positive random matrices of a Wishart--like structure, $W=XX^{\dagger}$, where $X$ stands for a nonhermitian random matrix. In particular, making use of the Cauchy transform, we study free multiplicative powers of the Marchenko-Pastur (MP) distribution, ${\rm MP}^{\boxtimes s}$, which for an integer $s$ yield Fuss-Catalan distributions corr… ▽ More We investigate the level density for several ensembles of positive random matrices of a Wishart--like structure, $W=XX^{\dagger}$, where $X$ stands for a nonhermitian random matrix. In particular, making use of the Cauchy transform, we study free multiplicative powers of the Marchenko-Pastur (MP) distribution, ${\rm MP}^{\boxtimes s}$, which for an integer $s$ yield Fuss-Catalan distributions corresponding to a product of $s$ independent square random matrices, $X=X_1\cdots X_s$. New formulae for the level densities are derived for $s=3$ and $s=1/3$. Moreover, the level density corresponding to the generalized Bures distribution, given by the free convolution of arcsine and MP distributions is obtained. We also explain the reason of such a curious convolution. The technique proposed here allows for the derivation of the level densities for several other cases. △ Less

Submitted 24 June, 2015; v1 submitted 4 July, 2014; originally announced July 2014.

Comments: 10 latex pages including 4 figures, Ver 4, minor improvements and references update

Journal ref: Phys. Rev. E 92, 012121 (2015)

Risk management of solitary and eusocial reproduction

Authors: Feng Fu, Sarah D. Kocher, Martin A. Nowak

Abstract: Social insect colonies can be seen as a distinct form of biological organization because they function as superorganisms. Understanding how natural selection acts on the emergence and maintenance of these colonies remains a major question in evolutionary biology and ecology. Here, we explore this by using multi-type branching processes to calculate the basic reproductive ratios and the extinction… ▽ More Social insect colonies can be seen as a distinct form of biological organization because they function as superorganisms. Understanding how natural selection acts on the emergence and maintenance of these colonies remains a major question in evolutionary biology and ecology. Here, we explore this by using multi-type branching processes to calculate the basic reproductive ratios and the extinction probabilities for solitary versus eusocial reproductive strategies. In order to derive precise mathematical results, we use a simple haploid, asexual model. In general, we show that eusocial reproductive strategies are unlikely to materialize unless large fitness advantages are gained by the production of only a few workers. These fitness advantages are maximized through obligate rather than facultative eusocial strategies. Furthermore, we find that solitary reproduction is `unbeatable' as long as the solitary reproductive ratio exceeds a critical value. In these cases, no eusocial parameters exist that would reduce their probability of extinction. Our results help to explain why the number of solitary species exceeds that of eusocial ones: eusociality is a high risk, high reward strategy, while solitary reproduction is better at risk management. △ Less

Submitted 1 July, 2014; originally announced July 2014.

Comments: Comments are highly appreciated, and collaborations on further followup work are also extremely welcome. For contact details see http://www.tb.ethz.ch/people/fuf/index