-
Evaluating the Role of Security Assurance Cases in Agile Medical Device Development
Authors:
Max Fransson,
Adam Andersson,
Mazen Mohamad,
Jan-Philipp Steghöfer
Abstract:
Cybersecurity issues in medical devices threaten patient safety and can cause harm if exploited. Standards and regulations therefore require vendors of such devices to provide an assessment of the cybersecurity risks as well as a description of their mitigation. Security assurance cases (SACs) capture these elements as a structured argument. Compiling an SAC requires taking domain-specific regulat…
▽ More
Cybersecurity issues in medical devices threaten patient safety and can cause harm if exploited. Standards and regulations therefore require vendors of such devices to provide an assessment of the cybersecurity risks as well as a description of their mitigation. Security assurance cases (SACs) capture these elements as a structured argument. Compiling an SAC requires taking domain-specific regulations and requirements as well as the way of working into account. In this case study, we evaluate CASCADE, an approach for building SAC in the context of a large medical device manufacturer with an established agile development workflow. We investigate the regulatory context as well as the adaptations needed in the development process. Our results show the suitability of SACs in the medical device industry. We identified 17 use cases in which an SAC supports internal and external needs. The connection to safety assurance can be achieved by incorporating information from the risk assessment matrix into the SAC. Integration into the development process can be achieved by introducing a new role and rules for the design review and the release to production as well as additional criteria for the definition of done. We also show that SACs built with CASCADE fulfill the requirements of relevant standards in the medical domain such as ISO 14971.
△ Less
Submitted 10 July, 2024;
originally announced July 2024.
-
Discovery of the optical and radio counterpart to the fast X-ray transient EP240315a
Authors:
J. H. Gillanders,
L. Rhodes,
S. Srivastav,
F. Carotenuto,
J. Bright,
M. E. Huber,
H. F. Stevance,
S. J. Smartt,
K. C. Chambers,
T. -W. Chen,
R. Fender,
A. Andersson,
A. J. Cooper,
P. G. Jonker,
F. J. Cowie,
T. deBoer,
N. Erasmus,
M. D. Fulton,
H. Gao,
J. Herman,
C. -C. Lin,
T. Lowe,
E. A. Magnier,
H. -Y. Miao,
P. Minguez
, et al. (14 additional authors not shown)
Abstract:
Fast X-ray Transients (FXTs) are extragalactic bursts of soft X-rays first identified >10 years ago. Since then, nearly 40 events have been discovered, although almost all of these have been recovered from archival Chandra and XMM-Newton data. To date, optical sky surveys and follow-up searches have not revealed any multi-wavelength counterparts. The Einstein Probe, launched in January 2024, has s…
▽ More
Fast X-ray Transients (FXTs) are extragalactic bursts of soft X-rays first identified >10 years ago. Since then, nearly 40 events have been discovered, although almost all of these have been recovered from archival Chandra and XMM-Newton data. To date, optical sky surveys and follow-up searches have not revealed any multi-wavelength counterparts. The Einstein Probe, launched in January 2024, has started surveying the sky in the soft X-ray regime (0.5-4 keV) and will rapidly increase the sample of FXTs discovered in real time. Here, we report the first discovery of both an optical and radio counterpart to a distant FXT, the fourth source publicly released by the Einstein Probe. We discovered a fast-fading optical transient within the 3 arcmin localisation radius of EP240315a with the all-sky optical survey ATLAS, and our follow-up Gemini spectrum provides a redshift, z=4.859+/-0.002. Furthermore, we uncovered a radio counterpart in the S-band (3.0 GHz) with the MeerKAT radio interferometer. The optical (rest-frame UV) and radio luminosities indicate the FXT most likely originates from either a long gamma-ray burst or a relativistic tidal disruption event. This may be a fortuitous early mission detection by the Einstein Probe or may signpost a mode of discovery for high-redshift, high-energy transients through soft X-ray surveys, combined with locating multi-wavelength counterparts.
△ Less
Submitted 19 June, 2024; v1 submitted 16 April, 2024;
originally announced April 2024.
-
Asymmetric Influence of the Amplitude-Dependent Tune Shift on the Transverse Mode-Coupling Instability
Authors:
Miriam Brosi,
Francis Cullinan,
Åke Andersson,
Jonas Breunlin,
Pedro Fernandes Tavares
Abstract:
At the 3 GeV ring of the MAX IV Laboratory, a fourth generation ring-based synchrotron light source, an asymmetric influence of the sign of the amplitude-dependent tune shift (ADTS) on the transverse mode-coupling instability (TMCI) has been observed. Measurements of the instability, in dedicated single-bunch experiments at low chromaticity, revealed a significant dependence of the dynamics of the…
▽ More
At the 3 GeV ring of the MAX IV Laboratory, a fourth generation ring-based synchrotron light source, an asymmetric influence of the sign of the amplitude-dependent tune shift (ADTS) on the transverse mode-coupling instability (TMCI) has been observed. Measurements of the instability, in dedicated single-bunch experiments at low chromaticity, revealed a significant dependence of the dynamics of the instability above threshold on the sign of the ADTS. While for a negative sign of the ADTS the crossing of the instability does not lead to a loss of beam current, a positive sign results in the loss of 40% or more of the beam current at the threshold. In order to investigate the observed asymmetry, systematic measurement of beam dynamics above the threshold have been conducted in combination with particle tracking simulations with mbtrack2 and theoretical calculations of the Landau damping due to the ADTS. The findings point towards an influence of the Landau damping in combination with the low synchrotron frequency, which indicates that this effect could become relevant in future low-emittance electron storage rings.
△ Less
Submitted 17 April, 2024; v1 submitted 26 January, 2024;
originally announced January 2024.
-
Transcriptome-supervised classification of tissue morphology using deep learning
Authors:
Axel Andersson,
Gabriele Partel,
Leslie Solorzano,
Carolina Wählby
Abstract:
Deep learning has proven to successfully learn variations in tissue and cell morphology. Training of such models typically relies on expensive manual annotations. Here we conjecture that spatially resolved gene expression, e.i., the transcriptome, can be used as an alternative to manual annotations. In particular, we trained five convolutional neural networks with patches of different size extract…
▽ More
Deep learning has proven to successfully learn variations in tissue and cell morphology. Training of such models typically relies on expensive manual annotations. Here we conjecture that spatially resolved gene expression, e.i., the transcriptome, can be used as an alternative to manual annotations. In particular, we trained five convolutional neural networks with patches of different size extracted from locations defined by spatially resolved gene expression. The network is trained to classify tissue morphology related to two different genes, general tissue, as well as background, on an image of fluorescence stained nuclei in a mouse brain coronal section. Performance is evaluated on an independent tissue section from a different mouse brain, reaching an average Dice score of 0.51. Results may indicate that novel techniques for spatially resolved transcriptomics together with deep learning may provide a unique and unbiased way to find genotype-phenotype relationships.
△ Less
Submitted 7 December, 2023;
originally announced December 2023.
-
Cell segmentation of in situ transcriptomics data using signed graph partitioning
Authors:
Axel Andersson,
Andrea Behanova,
Carolina Wählby,
Filip Malmberg
Abstract:
The locations of different mRNA molecules can be revealed by multiplexed in situ RNA detection. By assigning detected mRNA molecules to individual cells, it is possible to identify many different cell types in parallel. This in turn enables investigation of the spatial cellular architecture in tissue, which is crucial for furthering our understanding of biological processes and diseases. However,…
▽ More
The locations of different mRNA molecules can be revealed by multiplexed in situ RNA detection. By assigning detected mRNA molecules to individual cells, it is possible to identify many different cell types in parallel. This in turn enables investigation of the spatial cellular architecture in tissue, which is crucial for furthering our understanding of biological processes and diseases. However, cell typing typically depends on the segmentation of cell nuclei, which is often done based on images of a DNA stain, such as DAPI. Limiting cell definition to a nuclear stain makes it fundamentally difficult to determine accurate cell borders, and thereby also difficult to assign mRNA molecules to the correct cell. As such, we have developed a computational tool that segments cells solely based on the local composition of mRNA molecules. First, a small neural network is trained to compute attractive and repulsive edges between pairs of mRNA molecules. The signed graph is then partitioned by a mutex watershed into components corresponding to different cells. We evaluated our method on two publicly available datasets and compared it against the current state-of-the-art and older baselines. We conclude that combining neural networks with combinatorial optimization is a promising approach for cell segmentation of in situ transcriptomics data.
△ Less
Submitted 7 December, 2023;
originally announced December 2023.
-
Data-driven methods for diffusivity prediction in nuclear fuels
Authors:
Galen T. Craven,
Renai Chen,
Michael W. D. Cooper,
Christopher Matthews,
Jason Rizk,
Walter Malone,
Landon Johnson,
Tammie Gibson,
David A. Andersson
Abstract:
The growth rate of structural defects in nuclear fuels under irradiation is intrinsically related to the diffusion rates of the defects in the fuel lattice. The generation and growth of atomistic structural defects can significantly alter the performance characteristics of the fuel. This alteration of functionality must be accurately captured to qualify a nuclear fuel for use in reactors. Predicti…
▽ More
The growth rate of structural defects in nuclear fuels under irradiation is intrinsically related to the diffusion rates of the defects in the fuel lattice. The generation and growth of atomistic structural defects can significantly alter the performance characteristics of the fuel. This alteration of functionality must be accurately captured to qualify a nuclear fuel for use in reactors. Predicting the diffusion coefficients of defects and how they impact macroscale properties such as swelling, gas release, and creep is therefore of significant importance in both the design of new nuclear fuels and the assessment of current fuel types. In this article, we apply data-driven methods focusing on machine learning (ML) to determine various diffusion properties of two nuclear fuels, uranium oxide and uranium nitride. We show that using ML can increase, often significantly, the accuracy of predicting diffusivity in nuclear fuels in comparison to current analytical models. We also illustrate how ML can be used to quickly develop fuel models with parameter dependencies that are more complex and robust than what is currently available in the literature. These results suggest there is potential for ML to accelerate the design, qualification, and implementation of nuclear fuels.
△ Less
Submitted 7 September, 2023;
originally announced October 2023.
-
GW190425: Pan-STARRS and ATLAS coverage of the skymap and limits on optical emission associated with FRB190425
Authors:
S. J. Smartt,
M. Nicholl,
S. Srivastav,
M. E. Huber,
K. C. Chambers,
K. W. Smith,
D. R. Young,
M. D. Fulton,
J. L. Tonry,
C. W. Stubbs,
L. Denneau,
A. J. Cooper,
A. Aamer,
J. P. Anderson,
A. Andersson,
J. Bulger,
T. -W Chen,
P. Clark,
T. de Boer,
H. Gao,
J. H. Gillanders,
A. Lawrence,
C. C. Lin,
T. B. Lowe,
E. A. Magnier
, et al. (10 additional authors not shown)
Abstract:
GW190425 is the second of only two binary neutron star (BNS) merger events to be significantly detected by the LIGO-Virgo- Kagra gravitational wave detectors. With a detection only in LIGO Livingston, the skymap containing the source was large and no plausible electromagnetic counterpart was found in real time searching in 2019. Here we summarise our ATLAS and Pan-STARRS wide-field optical coverag…
▽ More
GW190425 is the second of only two binary neutron star (BNS) merger events to be significantly detected by the LIGO-Virgo- Kagra gravitational wave detectors. With a detection only in LIGO Livingston, the skymap containing the source was large and no plausible electromagnetic counterpart was found in real time searching in 2019. Here we summarise our ATLAS and Pan-STARRS wide-field optical coverage of the skymap beginning within 1 hour and 3 hours respectively of the GW190425 merger time. More recently, a potential coincidence between GW190425 and a fast radio burst FRB 190425 has been suggested, given their spatial and temporal coincidence. The smaller sky localisation area of FRB 190425 and its dispersion measure have led to the identification of a likely host galaxy, UGC 10667 at a distance of 141 +/- 10 Mpc. Our optical imaging covered the galaxy 6.0 hrs after GW190425 was detected and 3.5 hrs after the FRB 190425. No optical emission was detected and further imaging at +1.2 and +13.2 days also revealed no emission. If the FRB 190425 and GW190425 association were real, we highlight our limits on kilonova emission from a BNS merger in UGC 10667. The model for producing FRB 190425 from a BNS merger involves a supramassive magnetised neutron star spinning down by dipole emission on the timescale of hours. We show that magnetar enhanced kilonova emission is ruled out by optical upper limits. The lack of detected optical emission from a kilonova in UGC 10667 disfavours, but does not disprove, the FRB-GW link for this source.
△ Less
Submitted 20 September, 2023;
originally announced September 2023.
-
Commensal Transient Searches in Eight Short Gamma Ray Burst Fields
Authors:
S. I. Chastain,
A. J. van der Horst,
A. Rowlinson,
L. Rhodes,
A. Andersson,
R. Diretse,
R. P. Fender,
P. A. Woudt
Abstract:
A new generation of radio telescopes with excellent sensitivity, instantaneous {\it uv} coverage, and large fields of view, are providing unprecedented opportunities for performing commensal transient searches. Here we present such a commensal search in deep observations of short gamma-ray burst fields carried out with the MeerKAT radio telescope in South Africa at 1.3 GHz. These four hour observa…
▽ More
A new generation of radio telescopes with excellent sensitivity, instantaneous {\it uv} coverage, and large fields of view, are providing unprecedented opportunities for performing commensal transient searches. Here we present such a commensal search in deep observations of short gamma-ray burst fields carried out with the MeerKAT radio telescope in South Africa at 1.3 GHz. These four hour observations of eight different fields span survey lengths of weeks to months. We also carry out transient searches in time slices of the full observations, at timescales of 15 minutes, and 8 seconds. We find 122 variable sources on the long timescales, of which 52 are likely active galactic nuclei, but there are likely also some radio flaring stars. While the variability is intrinsic in at least two cases, most of it is consistent with interstellar scintillation. In this study, we also place constraints on transient rates based on state-of-the-art transient simulations codes. We place an upper limit of $2\times10^{-4}$ transients per day per square degree for transients with peak flux of 5 mJy, and an upper limit of $2.5\times10^{-2}$ transients per day per square degree for transients with a fluence of 10 Jy ms, the minimum detectable fluence of our survey.
△ Less
Submitted 14 September, 2023;
originally announced September 2023.
-
A Radio Flare in the Long-Lived Afterglow of the Distant Short GRB 210726A: Energy Injection or a Reverse Shock from Shell Collisions?
Authors:
Genevieve Schroeder,
Lauren Rhodes,
Tanmoy Laskar,
Anya Nugent,
Alicia Rouco Escorial,
Jillian C. Rastinejad,
Wen-fai Fong,
Alexander J. van der Horst,
Péter Veres,
Kate D. Alexander,
Alex Andersson,
Edo Berger,
Peter K. Blanchard,
Sarah Chastain,
Lise Christensen,
Rob Fender,
David A. Green,
Paul Groot,
Ian Heywood,
Assaf Horesh,
Luca Izzo,
Charles D. Kilpatrick,
Elmar Körding,
Amy Lien,
Daniele B. Malesani
, et al. (10 additional authors not shown)
Abstract:
We present the discovery of the radio afterglow of the short $γ$-ray burst (GRB) 210726A, localized to a galaxy at a photometric redshift of $z\sim 2.4$. While radio observations commenced $\lesssim 1~$day after the burst, no radio emission was detected until $\sim11$~days. The radio afterglow subsequently brightened by a factor of $\sim 3$ in the span of a week, followed by a rapid decay (a ``rad…
▽ More
We present the discovery of the radio afterglow of the short $γ$-ray burst (GRB) 210726A, localized to a galaxy at a photometric redshift of $z\sim 2.4$. While radio observations commenced $\lesssim 1~$day after the burst, no radio emission was detected until $\sim11$~days. The radio afterglow subsequently brightened by a factor of $\sim 3$ in the span of a week, followed by a rapid decay (a ``radio flare''). We find that a forward shock afterglow model cannot self-consistently describe the multi-wavelength X-ray and radio data, and underpredicts the flux of the radio flare by a factor of $\approx 5$. We find that the addition of substantial energy injection, which increases the isotropic kinetic energy of the burst by a factor of $\approx 4$, or a reverse shock from a shell collision are viable solutions to match the broad-band behavior. At $z\sim 2.4$, GRB\,210726A is among the highest redshift short GRBs discovered to date as well as the most luminous in radio and X-rays. Combining and comparing all previous radio afterglow observations of short GRBs, we find that the majority of published radio searches conclude by $\lesssim 10~$days after the burst, potentially missing these late rising, luminous radio afterglows.
△ Less
Submitted 21 August, 2023;
originally announced August 2023.
-
FRB 20121102A: images of the bursts and the varying radio counterpart
Authors:
L. Rhodes,
M. Caleb,
B. W. Stappers,
A. Andersson,
M. C. Bezuidenhout,
L. N. Driessen,
I. Heywood
Abstract:
As more Fast Radio Bursts (FRBs) are being localised, we are learning that some fraction have persistent radio sources (PRSs). Such a discovery motivates an improvement in our understanding of the nature of those counterparts, the relation to the bursts themselves and why only some FRBs have PRSs. We report on observations made of FRB 20121102A with the MeerKAT radio telescope. Across five epochs,…
▽ More
As more Fast Radio Bursts (FRBs) are being localised, we are learning that some fraction have persistent radio sources (PRSs). Such a discovery motivates an improvement in our understanding of the nature of those counterparts, the relation to the bursts themselves and why only some FRBs have PRSs. We report on observations made of FRB 20121102A with the MeerKAT radio telescope. Across five epochs, we detect the PRS associated with FRB 20121102A. Our observations are split into a cluster of four epochs (MJD 58732 - 58764) and a separate single epoch about 1000days later. The measured flux density is constant across the first four observations but then decays by more than one-third in the final observation. Our observations on MJD 58736 coincided with the detections of 11 bursts from FRB 20121102A by the MeerTRAP backend, seven of which we detected in the image plane. We discuss the importance of image plane detections when considering the commensal transient searches being performed with MeerKAT and other radio facilities. We find that MeerKAT is so sensitive that within a two-second image, we can detect any FRB with a flux density above 2.4mJy at 1.3GHz and so could localise every FRB that has been detected by CHIME to date.
△ Less
Submitted 8 August, 2023;
originally announced August 2023.
-
A new method for short duration transient detection in radio images: Searching for transient sources in MeerKAT data of NGC 5068
Authors:
S. Fijma,
A. Rowlinson,
R. A. M. J. Wijers,
I. de Ruiter,
W. J. G. de Blok,
S. Chastain,
A. J. van der Horst,
Z. S. Meyers,
K. van der Meulen,
R. Fender,
P. A. Woudt,
A. Andersson,
A. Zijlstra,
J. Healy,
F. M. Maccagni
Abstract:
Transient surveys are a vital tool in exploring the dynamic universe, with radio transients acting as beacons for explosive and highly energetic astrophysical phenomena. However, performing commensal transient surveys using radio imaging can require a significant amount of computing power, data storage and time. With the instrumentation available to us, and with new and exciting radio interferomet…
▽ More
Transient surveys are a vital tool in exploring the dynamic universe, with radio transients acting as beacons for explosive and highly energetic astrophysical phenomena. However, performing commensal transient surveys using radio imaging can require a significant amount of computing power, data storage and time. With the instrumentation available to us, and with new and exciting radio interferometers in development, it is essential that we develop efficient methods to probe the radio transient sky. In this paper, we present results from an commensal short duration transient survey, on time scales of 8 seconds, 128 seconds and 1 hour, using data from the MeerKAT radio telescope. The dataset used was obtained as part of a galaxy observing campaign, and we focus on the field of NGC 5068. We present a quick, wide field imaging strategy to enable fast imaging of large datasets, and develop methods to efficiently filter detected transient candidates. No transient candidates were identified on the time scales of 8 seconds, 128 seconds and 1 hour, leading to competitive limits on the transient surface densities of $6.7{\times}10^{-5}$ deg$^{-1}$, $1.1{\times}10^{-3}$ deg$^{-1}$, and $3.2{\times}10^{-2}$ deg$^{-1}$ at sensitivities of 56.4 mJy, 19.2 mJy, and 3.9 mJy for the respective time scales. We find one possible candidate that could be associated with a stellar flare, that was rejected due to strict image quality control. Further short time-scale radio observations of this candidate could give definite results to its origin.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
Bursts from Space: MeerKAT - The first citizen science project dedicated to commensal radio transients
Authors:
Alex Andersson,
Chris Lintott,
Rob Fender,
Joe Bright,
Francesco Carotenuto,
Laura Driessen,
Mathilde Espinasse,
Kelebogile Gaseahalwe,
Ian Heywood,
Alexander J. van der Horst,
Sara Motta,
Lauren Rhodes,
Evangelia Tremou,
David R. A. Williams,
Patrick Woudt,
Xian Zhang,
Steven Bloemen,
Paul Groot,
Paul Vreeswijk,
Stefano Giarratana,
Payaswini Saikia,
Jonas Andersson,
Lizzeth Ruiz Arroyo,
Loïc Baert,
Matthew Baumann
, et al. (18 additional authors not shown)
Abstract:
The newest generation of radio telescopes are able to survey large areas with high sensitivity and cadence, producing data volumes that require new methods to better understand the transient sky. Here we describe the results from the first citizen science project dedicated to commensal radio transients, using data from the MeerKAT telescope with weekly cadence. Bursts from Space: MeerKAT was launc…
▽ More
The newest generation of radio telescopes are able to survey large areas with high sensitivity and cadence, producing data volumes that require new methods to better understand the transient sky. Here we describe the results from the first citizen science project dedicated to commensal radio transients, using data from the MeerKAT telescope with weekly cadence. Bursts from Space: MeerKAT was launched late in 2021 and received ~89000 classifications from over 1000 volunteers in 3 months. Our volunteers discovered 142 new variable sources which, along with the known transients in our fields, allowed us to estimate that at least 2.1 per cent of radio sources are varying at 1.28 GHz at the sampled cadence and sensitivity, in line with previous work. We provide the full catalogue of these sources, the largest of candidate radio variables to date. Transient sources found with archival counterparts include a pulsar (B1845-01) and an OH maser star (OH 30.1-0.7), in addition to the recovery of known stellar flares and X-ray binary jets in our observations. Data from the MeerLICHT optical telescope, along with estimates of long time-scale variability induced by scintillation, imply that the majority of the new variables are active galactic nuclei. This tells us that citizen scientists can discover phenomena varying on time-scales from weeks to several years. The success both in terms of volunteer engagement and scientific merit warrants the continued development of the project, whilst we use the classifications from volunteers to develop machine learning techniques for finding transients.
△ Less
Submitted 27 April, 2023;
originally announced April 2023.
-
Single-photon hot electron ionization of C$_{70}$
Authors:
Åke Andersson,
Luca Schio,
Robert Richter,
Michele Alagia,
Stefano Stranges,
Piero Ferrari,
Klavs Hansen,
Vitali Zhaunerchyk
Abstract:
Gas phase C$_{70}$ molecules have been ionized with single photons of energies between 16 eV and 70 eV and the electron spectra measured with velocity map imaging in coincidence with the ions. The doubly ionized and unfragmented species was present at photon energies of 22 eV and up, and triply charged ions from 55 eV. The low kinetic energy parts of the spectra are explained with thermal emission…
▽ More
Gas phase C$_{70}$ molecules have been ionized with single photons of energies between 16 eV and 70 eV and the electron spectra measured with velocity map imaging in coincidence with the ions. The doubly ionized and unfragmented species was present at photon energies of 22 eV and up, and triply charged ions from 55 eV. The low kinetic energy parts of the spectra are explained with thermal emission of transient hot electrons. Deviations at high photon energies are used to determine a value for the initial electron equilibration time. We propose a generally applicable mechanism, named Resonance Ionization Shadowing, for the creation of hot electrons by absorption of above-threshold energy photons.
△ Less
Submitted 11 October, 2022;
originally announced October 2022.
-
Serendipitous discovery of radio flaring behaviour from a nearby M dwarf with MeerKAT
Authors:
Alex Andersson,
Rob Fender,
Chris Lintott,
David Williams,
Laura Driessen,
Patrick Woudt,
Alexander van der Horst,
David Buckley,
Sara Motta,
Lauren Rhodes,
Nora Eisner,
Rachel Osten,
Paul Vreeswijk,
Steven Bloemen,
Paul Groot
Abstract:
We report on the detection of MKT J174641.0$-$321404, a new radio transient found in untargeted searches of wide-field MeerKAT radio images centred on the black hole X-ray binary H1743$-$322. MKT J174641.0$-$321404 is highly variable at 1.3 GHz and was detected three times during 11 observations of the field in late 2018, reaching a maximum flux density of 590 $\pm$ 60 $μ$Jy. We associate this rad…
▽ More
We report on the detection of MKT J174641.0$-$321404, a new radio transient found in untargeted searches of wide-field MeerKAT radio images centred on the black hole X-ray binary H1743$-$322. MKT J174641.0$-$321404 is highly variable at 1.3 GHz and was detected three times during 11 observations of the field in late 2018, reaching a maximum flux density of 590 $\pm$ 60 $μ$Jy. We associate this radio transient with a high proper motion, M dwarf star SCR~1746$-$3214 12 pc away from the Sun. Multiwavelength observations of this M dwarf indicate flaring activity across the electromagnetic spectrum, consistent with emission expected from dMe stars, and providing upper limits on quiescent brightness in both the radio and X-ray regimes. \textit{TESS} photometry reveals a rotational period for SCR~1746$-$3214 of $0.2292 \pm 0.0025$ days, which at its estimated radius makes the star a rapid rotator, comparable to other low mass systems. Dedicated spectroscopic follow up confirms the star as a mid-late spectral M dwarf with clear magnetic activity indicated by strong H$α$ emission. This transient's serendipitous discovery by MeerKAT, along with multiwavelength characterisation, make it a prime demonstration of both the capabilities of the current generation of radio interferometers and the value of simultaneous observations by optical facilities such as MeerLICHT. Our results build upon the literature of of M dwarfs' flaring behaviour, particularly relevant to the habitability of their planetary systems.
△ Less
Submitted 7 April, 2022;
originally announced April 2022.
-
An energy-based deep splitting method for the nonlinear filtering problem
Authors:
Kasper Bågmark,
Adam Andersson,
Stig Larsson
Abstract:
The purpose of this paper is to explore the use of deep learning for the solution of the nonlinear filtering problem. This is achieved by solving the Zakai equation by a deep splitting method, previously developed for approximate solution of (stochastic) partial differential equations. This is combined with an energy-based model for the approximation of functions by a deep neural network. This res…
▽ More
The purpose of this paper is to explore the use of deep learning for the solution of the nonlinear filtering problem. This is achieved by solving the Zakai equation by a deep splitting method, previously developed for approximate solution of (stochastic) partial differential equations. This is combined with an energy-based model for the approximation of functions by a deep neural network. This results in a computationally fast filter that takes observations as input and that does not require re-training when new observations are received. The method is tested on four examples, two linear in one and twenty dimensions and two nonlinear in one dimension. The method shows promising performance when benchmarked against the Kalman filter and the bootstrap particle filter.
△ Less
Submitted 21 December, 2022; v1 submitted 31 March, 2022;
originally announced March 2022.
-
Search and identification of transient and variable radio sources using MeerKAT observations: a case study on the MAXI J1820+070 field
Authors:
A. Rowlinson,
J. Meijn,
J. Bright,
A. J. van der Horst,
S. Chastain,
S. Fijma,
R. Fender,
I. Heywood,
R. A. M. J. Wijers,
P. A. Woudt,
A. Andersson,
G. R. Sivakoff,
E. Tremou,
L. N. Driessen
Abstract:
Many transient and variable sources detected at multiple wavelengths are also observed to vary at radio frequencies. However, these samples are typically biased towards sources that are initially detected in wide-field optical, X-ray or gamma-ray surveys. Many sources that are insufficiently bright at higher frequencies are therefore missed, leading to potential gaps in our knowledge of these sour…
▽ More
Many transient and variable sources detected at multiple wavelengths are also observed to vary at radio frequencies. However, these samples are typically biased towards sources that are initially detected in wide-field optical, X-ray or gamma-ray surveys. Many sources that are insufficiently bright at higher frequencies are therefore missed, leading to potential gaps in our knowledge of these sources and missing populations that are not detectable in optical, X-rays or gamma-rays. Taking advantage of new state-of-the-art radio facilities that provide high quality wide-field images with fast survey speeds, we can now conduct unbiased surveys for transient and variable sources at radio frequencies. In this paper, we present an unbiased survey using observations obtained by MeerKAT, a mid-frequency ($\sim$1.4 GHz) radio array in South Africa's Karoo Desert. The observations used were obtained as part of a weekly monitoring campaign for X-ray binaries (XRBs) and we focus on the field of MAXI J1820+070. We develop methods to optimally filter transient and variable candidates that can be directly applied to other datasets. In addition to MAXI J1820+070, we identify four likely active galactic nuclei, one source that could be a Galactic source (pulsar or quiescent X-ray binary) or an AGN, and one variable pulsar. No transient sources, defined as being undetected in deep images, were identified leading to a transient surface density of $<3.7\times10^{-2}$ deg$^{-2}$ at a sensitivity of 1 mJy on timescales of one week at 1.4 GHz.
△ Less
Submitted 29 August, 2022; v1 submitted 31 March, 2022;
originally announced March 2022.
-
End-to-end Multiple Instance Learning with Gradient Accumulation
Authors:
Axel Andersson,
Nadezhda Koriakina,
Nataša Sladoje,
Joakim Lindblad
Abstract:
Being able to learn on weakly labeled data, and provide interpretability, are two of the main reasons why attention-based deep multiple instance learning (ABMIL) methods have become particularly popular for classification of histopathological images. Such image data usually come in the form of gigapixel-sized whole-slide-images (WSI) that are cropped into smaller patches (instances). However, the…
▽ More
Being able to learn on weakly labeled data, and provide interpretability, are two of the main reasons why attention-based deep multiple instance learning (ABMIL) methods have become particularly popular for classification of histopathological images. Such image data usually come in the form of gigapixel-sized whole-slide-images (WSI) that are cropped into smaller patches (instances). However, the sheer size of the data makes training of ABMIL models challenging. All the instances from one WSI cannot be processed at once by conventional GPUs. Existing solutions compromise training by relying on pre-trained models, strategic sampling or selection of instances, or self-supervised learning. We propose a training strategy based on gradient accumulation that enables direct end-to-end training of ABMIL models without being limited by GPU memory. We conduct experiments on both QMNIST and Imagenette to investigate the performance and training time, and compare with the conventional memory-expensive baseline and a recent sampled-based approach. This memory-efficient approach, although slower, reaches performance indistinguishable from the memory-expensive baseline.
△ Less
Submitted 8 March, 2022;
originally announced March 2022.
-
Convergence of a robust deep FBSDE method for stochastic control
Authors:
Kristoffer Andersson,
Adam Andersson,
Cornelis W. Oosterlee
Abstract:
In this paper, we propose a deep learning based numerical scheme for strongly coupled FBSDEs, stemming from stochastic control. It is a modification of the deep BSDE method in which the initial value to the backward equation is not a free parameter, and with a new loss function being the weighted sum of the cost of the control problem, and a variance term which coincides with the mean squared erro…
▽ More
In this paper, we propose a deep learning based numerical scheme for strongly coupled FBSDEs, stemming from stochastic control. It is a modification of the deep BSDE method in which the initial value to the backward equation is not a free parameter, and with a new loss function being the weighted sum of the cost of the control problem, and a variance term which coincides with the mean squared error in the terminal condition. We show by a numerical example that a direct extension of the classical deep BSDE method to FBSDEs, fails for a simple linear-quadratic control problem, and motivate why the new method works. Under regularity and boundedness assumptions on the exact controls of time continuous and time discrete control problems, we provide an error analysis for our method. We show empirically that the method converges for three different problems, one being the one that failed for a direct extension of the deep BSDE method.
△ Less
Submitted 9 February, 2023; v1 submitted 18 January, 2022;
originally announced January 2022.
-
Oriented hairy graphs and moduli spaces of curves
Authors:
Assar Andersson,
Thomas Willwacher,
Marko Zivkovic
Abstract:
We discuss a graph complex formed by directed acyclic graphs with external legs. This complex comes in particular with a map to the ribbon graph complex computing the (compactly supported) cohomology of the moduli space of points $\mathcal M_{g,n}$, extending an earlier result of Merkulov-Willwacher. It is furthermore quasi-isomorphic to the hairy graph complex computing the weight 0 part of the c…
▽ More
We discuss a graph complex formed by directed acyclic graphs with external legs. This complex comes in particular with a map to the ribbon graph complex computing the (compactly supported) cohomology of the moduli space of points $\mathcal M_{g,n}$, extending an earlier result of Merkulov-Willwacher. It is furthermore quasi-isomorphic to the hairy graph complex computing the weight 0 part of the compactly supported cohomology of $\mathcal{M}_{g,n}$ according to Chan-Galatius-Payne. Hence we can naturally connect the works Chan-Galatius-Payne and of Merkulov-Willwacher and the ribbon graph complex and obtain a fairly satisfying picture of how all the pieces and various graph complexes fit together, at least in weight zero.
△ Less
Submitted 1 May, 2020;
originally announced May 2020.
-
Hairy graphs to ribbon graphs via a fixed source graph complex
Authors:
Assar Andersson,
Marko Živković
Abstract:
We show that the hairy graph complex $(HGC_{n,n},d)$ appears as an associated graded complex of the oriented graph complex $(OGC_{n+1},d)$, subject to the filtration on the number of targets, or equivalently sources, called the fixed source graph complex. The fixed source graph complex $(OGC_1,d_0)$ maps into the ribbon graph complex $RGC$, which models the moduli space of Riemann surfaces with ma…
▽ More
We show that the hairy graph complex $(HGC_{n,n},d)$ appears as an associated graded complex of the oriented graph complex $(OGC_{n+1},d)$, subject to the filtration on the number of targets, or equivalently sources, called the fixed source graph complex. The fixed source graph complex $(OGC_1,d_0)$ maps into the ribbon graph complex $RGC$, which models the moduli space of Riemann surfaces with marked points. The full differential $d$ on the oriented graph complex $OGC_{n+1}$ corresponds to the deformed differential $d+h$ on the hairy graph complex $HGC_{n,n}$, where $h$ adds a hair. This deformed complex $(HGC_{n,n},d+h)$ is already known to be quasi-isomorphic to standard Kontsevich's graph complex $GC^2_n$. This gives a new connection between the standard and the oriented version of Kontsevich's graph complex.
△ Less
Submitted 16 April, 2020; v1 submitted 19 December, 2019;
originally announced December 2019.
-
From deformation theory of wheeled props to classification of Kontsevich formality maps
Authors:
Assar Andersson,
Sergei Merkulov
Abstract:
We study homotopy theory of the wheeled prop controlling Poisson structures on arbitrary formal graded finite-dimensional manifolds and prove, in particular, that Grothendieck-Teichmueller group acts on that wheeled prop faithfully and homotopy non-trivially. Next we apply this homotopy theory to the study of the deformation complex of an arbitrary Maxim Kontsevich formality map and compute the fu…
▽ More
We study homotopy theory of the wheeled prop controlling Poisson structures on arbitrary formal graded finite-dimensional manifolds and prove, in particular, that Grothendieck-Teichmueller group acts on that wheeled prop faithfully and homotopy non-trivially. Next we apply this homotopy theory to the study of the deformation complex of an arbitrary Maxim Kontsevich formality map and compute the full cohomology group of that deformation complex in terms of the cohomology of a certain graph complex introduced earlier by Maxim Kontsevich in [K1] and studied by Thomas Willwacher in [W1].
△ Less
Submitted 25 November, 2019; v1 submitted 20 November, 2019;
originally announced November 2019.
-
Finite element approximation of Lyapunov equations related to parabolic stochastic PDEs
Authors:
Adam Andersson,
Annika Lang,
Andreas Petersson,
Leander Schroer
Abstract:
A numerical analysis for the fully discrete approximation of an operator Lyapunov equation related to linear SPDEs (stochastic partial differential equations) driven by multiplicative noise is considered. The discretization of the Lyapunov equation in space is given by finite elements and in time by a semiimplicit Euler scheme. The main result is the derivation of the rate of convergence in operat…
▽ More
A numerical analysis for the fully discrete approximation of an operator Lyapunov equation related to linear SPDEs (stochastic partial differential equations) driven by multiplicative noise is considered. The discretization of the Lyapunov equation in space is given by finite elements and in time by a semiimplicit Euler scheme. The main result is the derivation of the rate of convergence in operator norm. Moreover, it is shown that the solution of the equation provides a representation of a quadratic and path dependent functional of the SPDE solution. This fact yields a deterministic numerical method to compute such functionals. As a secondary result, weak error rates are established for a fully discrete finite element approximation of the SPDE with respect to this functional. This is obtained as a consequence of the approximation analysis of the Lyapunov equation. It is the first weak convergence analysis for fully discrete finite element approximations of SPDEs driven by multiplicative noise that obtains double the strong rate of convergence, especially for path dependent functionals and smooth spatial noise. Numerical experiments illustrate the results empirically and it is demonstrated that the deterministic method has advantages over Monte Carlo sampling in a stability context.
△ Less
Submitted 3 May, 2022; v1 submitted 11 October, 2019;
originally announced October 2019.
-
Quantization of Yang--Mills metrics on holomorphic vector bundles
Authors:
Andreas Andersson
Abstract:
We investigate quantization properties of Hermitian metrics on holomorphic vector bundles over homogeneous compact Kähler manifolds. This allows us to study operators on Hilbert function spaces using vector bundles in a new way. We show that Yang--Mills metrics can be quantized in a strong sense and for equivariant vector bundles we deduce a strong stability property which supersedes Gieseker-stab…
▽ More
We investigate quantization properties of Hermitian metrics on holomorphic vector bundles over homogeneous compact Kähler manifolds. This allows us to study operators on Hilbert function spaces using vector bundles in a new way. We show that Yang--Mills metrics can be quantized in a strong sense and for equivariant vector bundles we deduce a strong stability property which supersedes Gieseker-stability. We obtain interesting examples of generalized notions of contractive, isometric, and subnormal operator tuples which have geometric interpretations related to holomorphic vector bundles over coadjoint orbits.
△ Less
Submitted 13 March, 2019;
originally announced March 2019.
-
The Compact Linear Collider (CLIC) - 2018 Summary Report
Authors:
The CLIC,
CLICdp collaborations,
:,
T. K. Charles,
P. J. Giansiracusa,
T. G. Lucas,
R. P. Rassool,
M. Volpi,
C. Balazs,
K. Afanaciev,
V. Makarenko,
A. Patapenka,
I. Zhuk,
C. Collette,
M. J. Boland,
A. C. Abusleme Hoffman,
M. A. Diaz,
F. Garay,
Y. Chi,
X. He,
G. Pei,
S. Pei,
G. Shu,
X. Wang,
J. Zhang
, et al. (671 additional authors not shown)
Abstract:
The Compact Linear Collider (CLIC) is a TeV-scale high-luminosity linear $e^+e^-$ collider under development at CERN. Following the CLIC conceptual design published in 2012, this report provides an overview of the CLIC project, its current status, and future developments. It presents the CLIC physics potential and reports on design, technology, and implementation aspects of the accelerator and the…
▽ More
The Compact Linear Collider (CLIC) is a TeV-scale high-luminosity linear $e^+e^-$ collider under development at CERN. Following the CLIC conceptual design published in 2012, this report provides an overview of the CLIC project, its current status, and future developments. It presents the CLIC physics potential and reports on design, technology, and implementation aspects of the accelerator and the detector. CLIC is foreseen to be built and operated in stages, at centre-of-mass energies of 380 GeV, 1.5 TeV and 3 TeV, respectively. CLIC uses a two-beam acceleration scheme, in which 12 GHz accelerating structures are powered via a high-current drive beam. For the first stage, an alternative with X-band klystron powering is also considered. CLIC accelerator optimisation, technical developments and system tests have resulted in an increased energy efficiency (power around 170 MW) for the 380 GeV stage, together with a reduced cost estimate at the level of 6 billion CHF. The detector concept has been refined using improved software tools. Significant progress has been made on detector technology developments for the tracking and calorimetry systems. A wide range of CLIC physics studies has been conducted, both through full detector simulations and parametric studies, together providing a broad overview of the CLIC physics potential. Each of the three energy stages adds cornerstones of the full CLIC physics programme, such as Higgs width and couplings, top-quark properties, Higgs self-coupling, direct searches, and many precision electroweak measurements. The interpretation of the combined results gives crucial and accurate insight into new physics, largely complementary to LHC and HL-LHC. The construction of the first CLIC energy stage could start by 2026. First beams would be available by 2035, marking the beginning of a broad CLIC physics programme spanning 25-30 years.
△ Less
Submitted 6 May, 2019; v1 submitted 14 December, 2018;
originally announced December 2018.
-
Observatory science with eXTP
Authors:
Jean J. M. in 't Zand,
Enrico Bozzo,
Jinlu Qu,
Xiang-Dong Li,
Lorenzo Amati,
Yang Chen,
Immacolata Donnarumma,
Victor Doroshenko,
Stephen A. Drake,
Margarita Hernanz,
Peter A. Jenke,
Thomas J. Maccarone,
Simin Mahmoodifar,
Domitilla de Martino,
Alessandra De Rosa,
Elena M. Rossi,
Antonia Rowlinson,
Gloria Sala,
Giulia Stratta,
Thomas M. Tauris,
Joern Wilms,
Xuefeng Wu,
Ping Zhou,
Iván Agudo,
Diego Altamirano
, et al. (159 additional authors not shown)
Abstract:
In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to stu…
▽ More
In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.
△ Less
Submitted 10 December, 2018;
originally announced December 2018.
-
Parametrization of LSDA+$U$ for noncollinear magnetic configurations: Multipolar magnetism in UO$_2$
Authors:
S. L. Dudarev,
P. Liu,
D. A. Andersson,
C. R. Stanek,
T. Ozaki,
C. Franchini
Abstract:
To explore the formation of noncollinear magnetic configurations in materials with strongly correlated electrons, we derive a noncollinear LSDA+$U$ model involving only one parameter $U$, as opposed to the difference between the Hubbard and Stoner parameters $U-J$. Computing $U$ in the constrained random phase approximation, we investigate noncollinear magnetism of uranium dioxide UO$_2$ and find…
▽ More
To explore the formation of noncollinear magnetic configurations in materials with strongly correlated electrons, we derive a noncollinear LSDA+$U$ model involving only one parameter $U$, as opposed to the difference between the Hubbard and Stoner parameters $U-J$. Computing $U$ in the constrained random phase approximation, we investigate noncollinear magnetism of uranium dioxide UO$_2$ and find that the spin-orbit coupling (SOC) stabilizes the 3$\textbf{k}$ ordered magnetic ground state. The estimated SOC strength in UO$_2$ is as large as 0.73 eV per uranium atom, making spin and orbital degrees of freedom virtually inseparable. Using a multipolar pseudospin Hamiltonian, we show how octupolar and dipole-dipole exchange coupling help establish the 3$\textbf{k}$ magnetic ground state with canted ordering of uranium $f$-orbitals. The cooperative Jahn-Teller effect does not appear to play a significant part in stabilizing the noncollinear 3$\textbf{k}$ state, which has the lowest energy even in an undistorted lattice. The choice of parameter $U$ in the LSDA+$U$ model has a notable quantitative effect on the predicted properties of UO$_2$, in particular on the magnetic exchange interaction and, perhaps trivially, on the band gap: The value of $U=3.46$ eV computed fully $ab$ $initio$ delivers the band gap of 2.11~eV in good agreement with experiment, and a balanced account of other pertinent energy scales.
△ Less
Submitted 12 August, 2019; v1 submitted 16 November, 2018;
originally announced November 2018.
-
Malliavin regularity and weak approximation of semilinear SPDE with Lévy noise
Authors:
Adam Andersson,
Felix Lindner
Abstract:
We investigate the weak order of convergence for space-time discrete approximations of semilinear parabolic stochastic evolution equations driven by additive square-integrable Lévy noise. To this end, the Malliavin regularity of the solution is analyzed and recent results on refined Malliavin-Sobolev spaces from the Gaussian setting are extended to a Poissonian setting. For a class of path-depende…
▽ More
We investigate the weak order of convergence for space-time discrete approximations of semilinear parabolic stochastic evolution equations driven by additive square-integrable Lévy noise. To this end, the Malliavin regularity of the solution is analyzed and recent results on refined Malliavin-Sobolev spaces from the Gaussian setting are extended to a Poissonian setting. For a class of path-dependent test functions, we obtain that the weak rate of convergence is twice the strong rate.
△ Less
Submitted 26 August, 2018;
originally announced August 2018.
-
Hodge classes of Chern character forms on compact Kähler manifolds
Authors:
Andreas Andersson
Abstract:
In this paper we show that every rational cohomology class of type $(p,p)$ on a compact Kähler manifold can be representated as a differential $(p,p)$-form given by an explicit formula involving a Čech cocycle. First we represent Chern characters of smooth vector bundles by Čech cocycles with values in the sheaf of differential forms. We then consider the behavior of these cocycles with respect to…
▽ More
In this paper we show that every rational cohomology class of type $(p,p)$ on a compact Kähler manifold can be representated as a differential $(p,p)$-form given by an explicit formula involving a Čech cocycle. First we represent Chern characters of smooth vector bundles by Čech cocycles with values in the sheaf of differential forms. We then consider the behavior of these cocycles with respect to the Hodge structure on cohomology when the base manifold is compact Kähler.
△ Less
Submitted 10 August, 2018;
originally announced August 2018.
-
Operational experience of ALBA's Digital LLRF at SOLARIS Light Source
Authors:
P. Borowiec,
L. Dudek,
W. Kitka,
A. Kisiel,
P. Klimczyk,
M. Knafel,
M. Kopec,
A. Wawrzyniak,
F. Perez,
A. Salom,
A. Andersson,
R. Lindvall,
L. Malmgren,
A. Mitrovic
Abstract:
For control of RF cavities installed in Solaris storage ring light source the digital Low Level RF (dLLRF) system was necessary from the beginning of operation. Since there were no expertise at the new constructed facility and no time for development due to funds deadline, almost turn-key dLLRF from Alba has been implemented according to MAXIV selection. Thanks to high flexibility of dLLRF only sm…
▽ More
For control of RF cavities installed in Solaris storage ring light source the digital Low Level RF (dLLRF) system was necessary from the beginning of operation. Since there were no expertise at the new constructed facility and no time for development due to funds deadline, almost turn-key dLLRF from Alba has been implemented according to MAXIV selection. Thanks to high flexibility of dLLRF only small adaptations were needed in terms of interfaces to auxiliary systems and setup of parameters. This paper summarizes operational experience about installation, commissioning, learning-curve from entry-level user, beam operation and future upgrades of this dLLRF.
△ Less
Submitted 23 March, 2018;
originally announced March 2018.
-
Performance Overhead Comparison between Hypervisor and Container based Virtualization
Authors:
Zheng Li,
Maria Kihl,
Qinghua Lu,
Jens A. Andersson
Abstract:
The current virtualization solution in the Cloud widely relies on hypervisor-based technologies. Along with the recent popularity of Docker, the container-based virtualization starts receiving more attention for being a promising alternative. Since both of the virtualization solutions are not resource-free, their performance overheads would lead to negative impacts on the quality of Cloud services…
▽ More
The current virtualization solution in the Cloud widely relies on hypervisor-based technologies. Along with the recent popularity of Docker, the container-based virtualization starts receiving more attention for being a promising alternative. Since both of the virtualization solutions are not resource-free, their performance overheads would lead to negative impacts on the quality of Cloud services. To help fundamentally understand the performance difference between these two types of virtualization solutions, we use a physical machine with "just-enough" resource as a baseline to investigate the performance overhead of a standalone Docker container against a standalone virtual machine (VM). With findings contrary to the related work, our evaluation results show that the virtualization's performance overhead could vary not only on a feature-by-feature basis but also on a job-to-job basis. Although the container-based solution is undoubtedly lightweight, the hypervisor-based technology does not come with higher performance overhead in every case. For example, Docker containers particularly exhibit lower QoS in terms of storage transaction speed.
△ Less
Submitted 4 August, 2017;
originally announced August 2017.
-
The conundrum of relaxation volumes in first-principles calculations of charge defects
Authors:
Anuj Goyal,
Kiran Mathew,
Richard G. Hennig,
Aleksandr Chernatynskiy,
Christopher R. Stank,
Samuel T. Murphy,
David A. Andersson,
Simon R. Phillpot,
Blas P. Uberuaga
Abstract:
The defect relaxation volumes obtained from density-functional theory (DFT) calculations of charged vacancies and interstitials are much larger than their neutral counterparts, seemingly unphysically large. In this work, we investigate the possible reasons for this and revisit the methods that address the calculation of charged defect structures in periodic DFT. We probe the dependence of the prop…
▽ More
The defect relaxation volumes obtained from density-functional theory (DFT) calculations of charged vacancies and interstitials are much larger than their neutral counterparts, seemingly unphysically large. In this work, we investigate the possible reasons for this and revisit the methods that address the calculation of charged defect structures in periodic DFT. We probe the dependence of the proposed energy corrections to charged defect formation energies on relaxation volumes and find that corrections such as the image charge correction and the alignment correction, which can lead to sizable changes in defect formation energies, have an almost negligible effect on the charged defect relaxation volume. We also investigate the volume for the net neutral defect reactions comprised of individual charged defects, and find that the aggregate formation volumes have reasonable magnitudes. This work highlights an important issue that, as for defect formation energies, the defect formation volumes depend on the choice of reservoir. We show that considering the change in volume of the electron reservoir in the formation reaction of the charged defects, analogous to how volumes of atoms are accounted for in defect formation volumes, can renormalize the formation volumes of charged defects such that they are comparable to neutral defects. This approach enables the description of the elastic properties of isolated charged defects within the overall neutral material, beyond the context of the overall defect reactions that produce the charged defect.
△ Less
Submitted 13 April, 2017;
originally announced April 2017.
-
Poisson Malliavin calculus in Hilbert space with an application to SPDE
Authors:
Adam Andersson,
Felix Lindner
Abstract:
In this paper we introduce a Hilbert space-valued Malliavin calculus for Poisson random measures. It is solely based on elementary principles from the theory of point processes and basic moment estimates, and thus allows for a simple treatment of the Malliavin operators. The main part of the theory is developed for general Poisson random measures, defined on a $σ$-finite measure space, with minima…
▽ More
In this paper we introduce a Hilbert space-valued Malliavin calculus for Poisson random measures. It is solely based on elementary principles from the theory of point processes and basic moment estimates, and thus allows for a simple treatment of the Malliavin operators. The main part of the theory is developed for general Poisson random measures, defined on a $σ$-finite measure space, with minimal conditions. The theory is shown to apply to a space-time setting, suitable for studying stochastic partial differential equations. As an application, we analyze the weak order of convergence of space-time approximations for a class of linear equations with $α$-stable noise, $α\in(1,2)$. For a suitable class of test functions, the weak order of convergence is found to be $α$ times the strong order.
△ Less
Submitted 21 March, 2017;
originally announced March 2017.
-
Statistical mechanics for metabolic networks during steady-state growth
Authors:
Daniele De Martino,
Anna MC Andersson,
Tobias Bergmiller,
Călin C Guet,
Gašper Tkačik
Abstract:
Which properties of metabolic networks can be derived solely from stoichiometric information about the network's constituent reactions? Predictive results have been obtained by Flux Balance Analysis (FBA), by postulating that cells set metabolic fluxes within the allowed stoichiometry so as to maximize their growth. Here, we generalize this framework to single cell level using maximum entropy mode…
▽ More
Which properties of metabolic networks can be derived solely from stoichiometric information about the network's constituent reactions? Predictive results have been obtained by Flux Balance Analysis (FBA), by postulating that cells set metabolic fluxes within the allowed stoichiometry so as to maximize their growth. Here, we generalize this framework to single cell level using maximum entropy models from statistical physics. We define and compute, for the core metabolism of Escherichia coli, a joint distribution over all fluxes that yields the experimentally observed growth rate. This solution, containing FBA as a limiting case, provides a better match to the measured fluxes in the wild type and several mutants. We find that E. coli metabolism is close to, but not at, the optimality assumed by FBA. Moreover, our model makes a wide range of predictions: (i) on flux variability, its regulation, and flux correlations across individual cells; (ii) on the relative importance of stoichiometric constraints vs. growth rate optimization; (iii) on quantitative scaling relations for singe-cell growth rate distributions. We validate these scaling predictions using data from individual bacterial cells grown in a microfluidic device at different sub-inhibitory antibiotic concentrations. Under mild dynamical assumptions, fluctuation-response relations further predict the autocorrelation timescale in growth data and growth rate adaptation times following an environmental perturbation.
△ Less
Submitted 6 March, 2017;
originally announced March 2017.
-
Regularity properties for solutions of infinite dimensional Kolmogorov equations in Hilbert spaces
Authors:
Adam Andersson,
Mario Hefter,
Arnulf Jentzen,
Ryan Kurniawan
Abstract:
In this article we establish regularity properties for solutions of infinite dimensional Kolmogorov equations. We prove that if the nonlinear drift coefficients, the nonlinear diffusion coefficients, and the initial conditions of the considered Kolmogorov equations are $n$-times continuously Fréchet differentiable, then so are the generalized solutions at every positive time. In addition, a key co…
▽ More
In this article we establish regularity properties for solutions of infinite dimensional Kolmogorov equations. We prove that if the nonlinear drift coefficients, the nonlinear diffusion coefficients, and the initial conditions of the considered Kolmogorov equations are $n$-times continuously Fréchet differentiable, then so are the generalized solutions at every positive time. In addition, a key contribution of this work is to prove suitable enhanced regularity properties for the derivatives of the generalized solutions of the Kolmogorov equations in the sense that the dominating linear operator in the drift coefficient of the Kolmogorov equation regularizes the higher order derivatives of the solutions. Such enhanced regularity properties are of major importance for establishing weak convergence rates for spatial and temporal numerical approximations of stochastic partial differential equations.
△ Less
Submitted 2 November, 2016;
originally announced November 2016.
-
On the differentiability of solutions of stochastic evolution equations with respect to their initial values
Authors:
Adam Andersson,
Arnulf Jentzen,
Ryan Kurniawan,
Timo Welti
Abstract:
In this article we study the differentiability of solutions of parabolic semilinear stochastic evolution equations (SEEs) with respect to their initial values. We prove that if the nonlinear drift coefficients and the nonlinear diffusion coefficients of the considered SEEs are $n$-times continuously Fréchet differentiable, then the solutions of the considered SEEs are also $n$-times continuously F…
▽ More
In this article we study the differentiability of solutions of parabolic semilinear stochastic evolution equations (SEEs) with respect to their initial values. We prove that if the nonlinear drift coefficients and the nonlinear diffusion coefficients of the considered SEEs are $n$-times continuously Fréchet differentiable, then the solutions of the considered SEEs are also $n$-times continuously Fréchet differentiable with respect to their initial values. Moreover, a key contribution of this work is to establish suitable enhanced regularity properties of the derivative processes of the considered SEE in the sense that the dominating linear operator appearing in the SEE smoothes the higher order derivative processes.
△ Less
Submitted 2 November, 2016;
originally announced November 2016.
-
Characterizing the gamma-ray long-term variability of PKS 2155-304 with H.E.S.S. and Fermi-LAT
Authors:
H. E. S. S. Collaboration,
H. Abdalla,
A. Abramowski,
F. Aharonian,
F. Ait Benkhali,
A. G. Akhperjanian T. Andersson,
M. Arrieta,
P. Aubert,
M. Backes,
A. Balzer,
M. Barnard,
Y. Becherini,
J. Becker Tjus,
D. Berge S. Bernhard,
K. Bernlöhr,
R. Blackwell,
M. Böttcher,
C. Boisson,
J. Bolmont,
P. Bordas,
J. Bregeon,
F. Brun,
P. Brun,
M. Bryan,
T. Bulik
, et al. (222 additional authors not shown)
Abstract:
Studying the temporal variability of BL Lac objects at the highest energies provides unique insights into the extreme physical processes occurring in relativistic jets and in the vicinity of super-massive black holes. To this end, the long-term variability of the BL Lac object PKS 2155-304 is analyzed in the high (HE, 100 MeV < E < 300 GeV) and very high energy (VHE, E > 200 GeV) gamma-ray domain.…
▽ More
Studying the temporal variability of BL Lac objects at the highest energies provides unique insights into the extreme physical processes occurring in relativistic jets and in the vicinity of super-massive black holes. To this end, the long-term variability of the BL Lac object PKS 2155-304 is analyzed in the high (HE, 100 MeV < E < 300 GeV) and very high energy (VHE, E > 200 GeV) gamma-ray domain. Over the course of ~9 yr of H.E.S.S observations the VHE light curve in the quiescent state is consistent with a log-normal behavior. The VHE variability in this state is well described by flicker noise (power-spectral-density index ß_VHE = 1.10 +0.10 -0.13) on time scales larger than one day. An analysis of 5.5 yr of HE Fermi LAT data gives consistent results (ß_HE = 1.20 +0.21 -0.23, on time scales larger than 10 days) compatible with the VHE findings. The HE and VHE power spectral densities show a scale invariance across the probed time ranges. A direct linear correlation between the VHE and HE fluxes could neither be excluded nor firmly established. These long-term-variability properties are discussed and compared to the red noise behavior (ß ~ 2) seen on shorter time scales during VHE-flaring states. The difference in power spectral noise behavior at VHE energies during quiescent and flaring states provides evidence that these states are influenced by different physical processes, while the compatibility of the HE and VHE long-term results is suggestive of a common physical link as it might be introduced by an underlying jet-disk connection.
△ Less
Submitted 11 October, 2016;
originally announced October 2016.
-
Updated baseline for a staged Compact Linear Collider
Authors:
The CLIC,
CLICdp collaborations,
:,
M. J. Boland,
U. Felzmann,
P. J. Giansiracusa,
T. G. Lucas,
R. P. Rassool,
C. Balazs,
T. K. Charles,
K. Afanaciev,
I. Emeliantchik,
A. Ignatenko,
V. Makarenko,
N. Shumeiko,
A. Patapenka,
I. Zhuk,
A. C. Abusleme Hoffman,
M. A. Diaz Gutierrez,
M. Vogel Gonzalez,
Y. Chi,
X. He,
G. Pei,
S. Pei,
G. Shu
, et al. (493 additional authors not shown)
Abstract:
The Compact Linear Collider (CLIC) is a multi-TeV high-luminosity linear e+e- collider under development. For an optimal exploitation of its physics potential, CLIC is foreseen to be built and operated in a staged approach with three centre-of-mass energy stages ranging from a few hundred GeV up to 3 TeV. The first stage will focus on precision Standard Model physics, in particular Higgs and top-q…
▽ More
The Compact Linear Collider (CLIC) is a multi-TeV high-luminosity linear e+e- collider under development. For an optimal exploitation of its physics potential, CLIC is foreseen to be built and operated in a staged approach with three centre-of-mass energy stages ranging from a few hundred GeV up to 3 TeV. The first stage will focus on precision Standard Model physics, in particular Higgs and top-quark measurements. Subsequent stages will focus on measurements of rare Higgs processes, as well as searches for new physics processes and precision measurements of new states, e.g. states previously discovered at LHC or at CLIC itself. In the 2012 CLIC Conceptual Design Report, a fully optimised 3 TeV collider was presented, while the proposed lower energy stages were not studied to the same level of detail. This report presents an updated baseline staging scenario for CLIC. The scenario is the result of a comprehensive study addressing the performance, cost and power of the CLIC accelerator complex as a function of centre-of-mass energy and it targets optimal physics output based on the current physics landscape. The optimised staging scenario foresees three main centre-of-mass energy stages at 380 GeV, 1.5 TeV and 3 TeV for a full CLIC programme spanning 22 years. For the first stage, an alternative to the CLIC drive beam scheme is presented in which the main linac power is produced using X-band klystrons.
△ Less
Submitted 27 March, 2017; v1 submitted 26 August, 2016;
originally announced August 2016.
-
Existence, uniqueness, and regularity for stochastic evolution equations with irregular initial values
Authors:
Adam Andersson,
Arnulf Jentzen,
Ryan Kurniawan
Abstract:
In this article we develop a framework for studying parabolic semilinear stochastic evolution equations (SEEs) with singularities in the initial condition and singularities at the initial time of the time-dependent coefficients of the considered SEE. We use this framework to establish existence, uniqueness, and regularity results for mild solutions of parabolic semilinear SEEs with singularities a…
▽ More
In this article we develop a framework for studying parabolic semilinear stochastic evolution equations (SEEs) with singularities in the initial condition and singularities at the initial time of the time-dependent coefficients of the considered SEE. We use this framework to establish existence, uniqueness, and regularity results for mild solutions of parabolic semilinear SEEs with singularities at the initial time. We also provide several counterexample SEEs that illustrate the optimality of our results.
△ Less
Submitted 6 October, 2016; v1 submitted 21 December, 2015;
originally announced December 2015.
-
Electromagnetism in terms of quantum measurements
Authors:
Andreas Andersson
Abstract:
We consider the question whether electromagnetism can be derived from quantum physics of measurements. It turns out that this is possible, both for quantum and classical electromagnetism, if we use more recent innovations such as smearing of observables and simultaneous measurability. In this way we justify the use of von Neumann-type measurement models for physical processes.
We apply operation…
▽ More
We consider the question whether electromagnetism can be derived from quantum physics of measurements. It turns out that this is possible, both for quantum and classical electromagnetism, if we use more recent innovations such as smearing of observables and simultaneous measurability. In this way we justify the use of von Neumann-type measurement models for physical processes.
We apply operational quantum measurement theory to gain insight in fundamental aspects of quantum physics. Interactions of von Neumann type make the Heisenberg evolution of observables describable using explicit operator deformations. In this way one can obtain quantized electromagnetism as a measurement of a system by another. The relevant deformations (Rieffel deformations) have a mathematically well-defined "classical" limit which is indeed classical electromagnetism for our choice of interaction.
△ Less
Submitted 16 September, 2015;
originally announced September 2015.
-
Mean-square convergence of the BDF2-Maruyama and backward Euler schemes for SDE satisfying a global monotonicity condition
Authors:
Adam Andersson,
Raphael Kruse
Abstract:
In this paper the numerical approximation of stochastic differential equations satisfying a global monotonicity condition is studied. The strong rate of convergence with respect to the mean square norm is determined to be $\frac{1}{2}$ for the two-step BDF-Maruyama scheme and for the backward Euler-Maruyama method. In particular, this is the first paper which proves a strong convergence rate for a…
▽ More
In this paper the numerical approximation of stochastic differential equations satisfying a global monotonicity condition is studied. The strong rate of convergence with respect to the mean square norm is determined to be $\frac{1}{2}$ for the two-step BDF-Maruyama scheme and for the backward Euler-Maruyama method. In particular, this is the first paper which proves a strong convergence rate for a multi-step method applied to equations with possibly superlinearly growing drift and diffusion coefficient functions. We also present numerical experiments for the $\tfrac32$-volatility model from finance, which verify our results in practice and indicate that the BDF2-Maruyama method offers advantages over Euler-type methods if the stochastic differential equation is stiff or driven by a noise with small intensity.
△ Less
Submitted 21 June, 2016; v1 submitted 2 September, 2015;
originally announced September 2015.
-
Berezin quantization of noncommutative projective varieties
Authors:
Andreas Andersson
Abstract:
We use operator algebras and operator theory to obtain new result concerning Berezin quantization of compact Kähler manifolds. Our main tool is the notion of subproduct systems of finite-dimensional Hilbert spaces, which enables all involved objects, such as the Toeplitz operators, to be very conveniently expressed in terms of shift operators compressed to a subspace of full Fock space. This subsp…
▽ More
We use operator algebras and operator theory to obtain new result concerning Berezin quantization of compact Kähler manifolds. Our main tool is the notion of subproduct systems of finite-dimensional Hilbert spaces, which enables all involved objects, such as the Toeplitz operators, to be very conveniently expressed in terms of shift operators compressed to a subspace of full Fock space. This subspace is not required to be contained in the symmetric Fock space, so from finite-dimensional matrix algebras we can construct noncommutative manifolds with extra structure generalizing that of a projective variety endowed with a positive Hermitian line bundle and a canonical Kähler metric in the class of the line bundle. Even in the commutative setting these constructions are very fruitful. Firstly, we show that the algebra of smooth functions on any smooth projective variety can be quantized in a strong sense of inductive limits, as was previously only accomplished for homogeneous manifolds. In this way the Kähler manifold is recovered exactly from quantization and not just approximately. Secondly, we obtain a strict quantization also for singular varieties. Thirdly, we show that the Arveson conjecture is true in full generality for shift operators compressed to the subspace of symmetric Fock space associated with any homogeneous ideal. For noncommutative examples we consider homogeneous spaces for compact matrix quantum groups which generalize $q$-deformed projective spaces, and we show that these can be obtained as the cores of Cuntz--Pimsner algebras constructed solely from the representation theory of the quantum group. We also discuss interesting connections with noncommutative random walks.
△ Less
Submitted 4 February, 2018; v1 submitted 3 June, 2015;
originally announced June 2015.
-
Dequantization via quantum channels
Authors:
Andreas Andersson
Abstract:
For a unital completely positive map $Φ$ ("quantum channel") governing the time propagation of a quantum system, the Stinespring representation gives an enlarged system evolving unitarily. We argue that the Stinespring representations of each power $Φ^m$ of the single map together encode the structure of the original quantum channel and provides an interaction-dependent model for the bath. The sam…
▽ More
For a unital completely positive map $Φ$ ("quantum channel") governing the time propagation of a quantum system, the Stinespring representation gives an enlarged system evolving unitarily. We argue that the Stinespring representations of each power $Φ^m$ of the single map together encode the structure of the original quantum channel and provides an interaction-dependent model for the bath. The same bath model gives a "classical limit" at infinite time $m\to\infty$ in the form of a noncommutative "manifold" determined by the channel. In this way a simplified analysis of the system can be performed by making the large-$m$ approximation. These constructions are based on a noncommutative generalization of Berezin quantization. The latter is shown to involve very fundamental aspects of quantum-information theory, which are thereby put in a completely new light.
△ Less
Submitted 3 June, 2015;
originally announced June 2015.
-
Detailed balance as a quantum-group symmetry of Kraus operators
Authors:
Andreas Andersson
Abstract:
A unital completely positive map governing the time evolution of a quantum system is usually called a quantum channel, and it can be represented by a tuple of operators which are then referred to as the Kraus operators of the channel. We look at states of the system whose correlations with respect to the channel have a certain symmetry. Then we show that detailed balance holds if the Kraus operato…
▽ More
A unital completely positive map governing the time evolution of a quantum system is usually called a quantum channel, and it can be represented by a tuple of operators which are then referred to as the Kraus operators of the channel. We look at states of the system whose correlations with respect to the channel have a certain symmetry. Then we show that detailed balance holds if the Kraus operators satisfy a very interesting algebraic relation which plays an important role in the representation theory of any compact quantum group.
△ Less
Submitted 20 July, 2015; v1 submitted 1 June, 2015;
originally announced June 2015.
-
Weak error analysis for semilinear stochastic Volterra equations with additive noise
Authors:
Adam Andersson,
Mihály Kovács,
Stig Larsson
Abstract:
We prove a weak error estimate for the approximation in space and time of a semilinear stochastic Volterra integro-differential equation driven by additive space-time Gaussian noise. We treat this equation in an abstract framework, in which parabolic stochastic partial differential equations are also included as a special case. The approximation in space is performed by a standard finite element m…
▽ More
We prove a weak error estimate for the approximation in space and time of a semilinear stochastic Volterra integro-differential equation driven by additive space-time Gaussian noise. We treat this equation in an abstract framework, in which parabolic stochastic partial differential equations are also included as a special case. The approximation in space is performed by a standard finite element method and in time by an implicit Euler method combined with a convolution quadrature. The weak rate of convergence is proved to be twice the strong rate, as expected. Our convergence result concerns not only functionals of the solution at a fixed time but also more complicated functionals of the entire path and includes convergence of covariances and higher order statistics. The proof does not rely on a Kolmogorov equation. Instead it is based on a duality argument from Malliavin calculus.
△ Less
Submitted 14 March, 2016; v1 submitted 24 November, 2014;
originally announced November 2014.
-
Modeling and simulations of quantum phase slips in ultrathin superconducting wires
Authors:
Andreas Andersson,
Jack Lidmar
Abstract:
We study quantum phase slips (QPS) in ultrathin superconducting wires. Starting from an effective one-dimensional microscopic model, which includes electromagnetic fluctuations, we map the problem to a (1+1)-dimensional gas of interacting instantons. We introduce a method to calculate the tunneling amplitude of quantum phase slips directly from Monte Carlo simulations. This allows us to go beyond…
▽ More
We study quantum phase slips (QPS) in ultrathin superconducting wires. Starting from an effective one-dimensional microscopic model, which includes electromagnetic fluctuations, we map the problem to a (1+1)-dimensional gas of interacting instantons. We introduce a method to calculate the tunneling amplitude of quantum phase slips directly from Monte Carlo simulations. This allows us to go beyond the dilute instanton gas approximation and study the problem without any limitations of the density of QPS. We find that the tunneling amplitude shows a characteristic scaling behavior near the superconductor-insulator transition. We also calculate the voltage-charge relation of the insulating state, which is the dual of the Josephson current-phase relation in ordinary superconducting weak links. This evolves from a sinusoidal form in the regime of dilute QPS to more exotic shapes for higher QPS densities, where interactions are important.
△ Less
Submitted 16 April, 2015; v1 submitted 14 November, 2014;
originally announced November 2014.
-
Index pairings for $\mathbb{R}^n$-actions and Rieffel deformations
Authors:
Andreas Andersson
Abstract:
With an action $α$ of $\mathbb{R}^n$ on a $C^*$-algebra $A$ and a skew-symmetric $n\times n$ matrix $Θ$ one can consider the Rieffel deformation $A_Θ$ of $A$, which is a $C^*$-algebra generated by the $α$-smooth elements of $A$ with a new multiplication. The purpose of this paper is to obtain explicit formulas for $K$-theoretical quantities defined by elements of $A_Θ$. We assume that there is a d…
▽ More
With an action $α$ of $\mathbb{R}^n$ on a $C^*$-algebra $A$ and a skew-symmetric $n\times n$ matrix $Θ$ one can consider the Rieffel deformation $A_Θ$ of $A$, which is a $C^*$-algebra generated by the $α$-smooth elements of $A$ with a new multiplication. The purpose of this paper is to obtain explicit formulas for $K$-theoretical quantities defined by elements of $A_Θ$. We assume that there is a densely defined trace on $A$, invariant under the action. We give an explicit realization of Thom class in $KK$ in any dimension $n$, and use it in the index pairings. When $n$ is odd, for example, we give a formula for the index of operators of the form $Pπ^Θ(u)P$, where $π^Θ(u)$ is the operator of left Rieffel multiplication by an invertible element $u$ over the unitization of $A$, and $P$ is projection onto the nonnegative eigenspace of a Dirac operator constructed from the action $α$. The results are new also for the undeformed case $Θ=0$. The construction relies on two approaches to Rieffel deformations in addition to Rieffel's original one: "Kasprzak deformation" and "warped convolution". We end by outlining potential applications in mathematical physics.
△ Less
Submitted 23 November, 2017; v1 submitted 16 June, 2014;
originally announced June 2014.
-
Duality in refined Sobolev-Malliavin spaces and weak approximations of SPDE
Authors:
Adam Andersson,
Raphael Kruse,
Stig Larsson
Abstract:
We introduce a new family of refined Sobolev-Malliavin spaces that capture the integrability in time of the Malliavin derivative. We consider duality in these spaces and derive a Burkholder type inequality in a dual norm. The theory we develop allows us to prove weak convergence with essentially optimal rate for numerical approximations in space and time of semilinear parabolic stochastic evolutio…
▽ More
We introduce a new family of refined Sobolev-Malliavin spaces that capture the integrability in time of the Malliavin derivative. We consider duality in these spaces and derive a Burkholder type inequality in a dual norm. The theory we develop allows us to prove weak convergence with essentially optimal rate for numerical approximations in space and time of semilinear parabolic stochastic evolution equations driven by Gaussian additive noise. In particular, we combine a standard Galerkin finite element method with backward Euler timestepping. The method of proof does not rely on the use of the Kolmogorov equation or the Itō formula and is therefore non-Markovian in nature. Test functions satisfying polynomial growth and mild smoothness assumptions are allowed, meaning in particular that we prove convergence of arbitrary moments with essentially optimal rate.
△ Less
Submitted 1 June, 2015; v1 submitted 20 December, 2013;
originally announced December 2013.
-
CONCOCT: Clustering cONtigs on COverage and ComposiTion
Authors:
Johannes Alneberg,
Brynjar Smari Bjarnason,
Ino de Bruijn,
Melanie Schirmer,
Joshua Quick,
Umer Z. Ijaz,
Nicholas J. Loman,
Anders F. Andersson,
Christopher Quince
Abstract:
Metagenomics enables the reconstruction of microbial genomes in complex microbial communities without the need for culturing. Since assembly typically results in fragmented genomes the grouping of genome fragments (contigs) belonging to the same genome, a process referred to as binning, remains a major informatics challenge. Here we present CONCOCT, a computer program that combines three types of…
▽ More
Metagenomics enables the reconstruction of microbial genomes in complex microbial communities without the need for culturing. Since assembly typically results in fragmented genomes the grouping of genome fragments (contigs) belonging to the same genome, a process referred to as binning, remains a major informatics challenge. Here we present CONCOCT, a computer program that combines three types of information - sequence composition, coverage across multiple sample, and read-pair linkage - to automatically bin contigs into genomes. We demonstrate high recall and precision rates of the program on artificial as well as real human gut metagenome datasets.
△ Less
Submitted 14 December, 2013;
originally announced December 2013.
-
Fourier Methods for Harmonic Scalar Waves in General Waveguides
Authors:
Anders Andersson,
Borje Nilsson,
Thomas Biro
Abstract:
A set of semi-analytical techniques based on Fourier analysis is used to solve wave scattering problems in variously shaped waveguides with varying normal admittance boundary conditions. Key components are newly developed conformal mapping methods, wave splitting, Fourier series expansions in eigen-functions to non-normal operators, the building block method or the cascade technique, Dirichlet-to-…
▽ More
A set of semi-analytical techniques based on Fourier analysis is used to solve wave scattering problems in variously shaped waveguides with varying normal admittance boundary conditions. Key components are newly developed conformal mapping methods, wave splitting, Fourier series expansions in eigen-functions to non-normal operators, the building block method or the cascade technique, Dirichlet-to-Neumann operators, and reformulation in terms of stable differential equations for reflection and transmission matrices. For an example the results show good correspondence with a finite element method solution to the same problem in the low and medium frequency domain. The Fourier method complements finite element analysis as a waveguide simulation tool. For inverse engineering involving tuning of straight waveguide parts joining complicated waveguide elements, the Fourier method is an attractive alternative including time aspects. The prime motivation for the Fourier method is its added physical understanding primarily at low frequencies.
△ Less
Submitted 1 June, 2015; v1 submitted 30 July, 2013;
originally announced July 2013.
-
Operator Deformations in Quantum Measurement Theory
Authors:
Andreas Andersson
Abstract:
We describe rigorous quantum measurement theory in the Heisenberg picture by applying operator deformation techniques previously used in noncommutative quantum field theory. This enables the conventional observables (represented by unbounded operators) to play a role also in the more general setting.
We describe rigorous quantum measurement theory in the Heisenberg picture by applying operator deformation techniques previously used in noncommutative quantum field theory. This enables the conventional observables (represented by unbounded operators) to play a role also in the more general setting.
△ Less
Submitted 21 March, 2014; v1 submitted 9 April, 2013;
originally announced April 2013.