Search

We fabricate paper-supported semiconducting devices by rubbing a layered molybdenum disulfide (MoS2) crystal onto a piece of paper, similarly to the action of drawing/writing with a pencil on paper. We show that the abrasion between the MoS2 crystal and the paper substrate efficiently exfoliates the crystals, breaking the weak van der Waals interlayer bonds and leading to the deposition of a film of interconnected MoS2 platelets. Employing this simple method, that can be easily extended to other 2D materials, we fabricate MoS2-on-paper broadband photodectectors with spectral sensitivity from the ultraviolet (UV) to the near-infrared (NIR). We also used these paper-based photodetectors to acquire pictures of objects by mounting the photodetectors in a homebuilt single-pixel camera setup., Comment: 6 main text figures + 4 Supp Info figures

Assessments are usually thought of as ways for instructors to get information from students. In this work, we flip this perspective and explore how assessments communicate information to students. Specifically, we consider how assessments may provide information about what faculty and/or researchers think it means to know and do physics, i.e., their epistemologies. Using data from students completing assessment questions during one-on-one think aloud interviews, we explore how assessment features did (or did not) impact student engagement with the assessment problems. We analyze video recordings and transcripts to infer the epistemological framings and resources students use while completing introductory-level physics problems. Students’ framings tended to be fairly stable, but when shifts occurred, they were triggered by a shift in epistemological resource, which can be activated by assessment feature. This work extends existing work on epistemological framing into the realm of assessment and allows us to consider the effects of assessments on our students’ understanding of physics teaching and learning.

The paper alluded to in the title contains the following striking result: Let $I$ be the unit interval and $\Delta$ the Cantor set. If $X$ is a quasi Banach space containing no copy of $c_0$ which is isomorphic to a closed subspace of a space with a basis and $C(I, X)$ is linearly homeomorphic to $C(\Delta, X)$, then $X$ is locally convex, i.e., a Banach space. It is shown that Kalton result is sharp by exhibiting non locally convex quasi Banach spaces X with a basis for which $C(I, X)$ and $C(\Delta, X)$ are isomorphic. Our examples are rather specific and actually in all cases X is isomorphic to $C(\Delta, X)$ if $K$ is a metric compactum of finite covering dimension., Comment: 4 pages

This letter solves an open question of origami paper spring risen by Yoneda et al.(Phys. Rev. E 2019). By using both dimensional analysis and data fitting, a universal scaling law of a paper spring is formulated. The scaling law shows that origami spring force obeys power square law of spring extension, however strong nonlinear to the total twist angle. The study has also successfully generalized the scaling law from the Poisson ratio 0.3 to an arbitrary Poisson's ratio with the help of dimensional analysis., Comment: 3 pages, 3 figues and 1 table

We investigate a modified spatial stochastic Lotka-Volterra formulation of the rock-paper-scissors model using off-lattice stochastic simulations. In this model one of the species moves preferentially in a specific direction -- the level of preference being controlled by a noise strength parameter $\eta \in [0, 1]$ ($\eta = 0$ and $\eta = 1$ corresponding to total preference and no preference, respectively) -- while the other two species have no referred direction of motion. We study the behaviour of the system starting from random initial conditions, showing that the species with asymmetric mobility has always an advantage over its predator. We also determine the optimal value of the noise strength parameter which gives the maximum advantage to that species. Finally, we find that the critical number of individuals, below which the probability of extinction becomes significant, decreases as the noise level increases, thus showing that the addition of a preferred mobility direction studied in the present paper does not favour coexistence., Comment: 5 pages, 7 figures

We investigate the problem of the predominance and survival of "weak" species in the context of the simplest generalization of the spatial stochastic rock-paper-scissors model to four species by considering models in which one, two, or three species have a reduced predation probability. We show, using lattice based spatial stochastic simulations with random initial conditions, that if only one of the four species has its probability reduced then the most abundant species is the prey of the "weakest" (assuming that the simulations are large enough for coexistence to prevail). Also, among the remaining cases, we present examples in which "weak" and "strong" species have similar average abundances and others in which either of them dominates -- the most abundant species being always a prey of a "weak" species with which it maintains a unidirectional predator-prey interaction. However, in contrast to the three-species model, we find no systematic difference in the global performance of "weak" and "strong" species, and we conjecture that the same result will hold if the number of species is further increased. We also determine the probability of single species survival and coexistence as a function of the lattice size, discussing its dependence on initial conditions and on the change to the dynamics of the model which results from the extinction of one of the species., Comment: 7 pages, 6 figures. arXiv admin note: text overlap with arXiv:1907.11280

We present a sensitive search with the Karl G. Jansky Very Large Array (VLA) for the radio counterpart of the gravitational wave candidate S191216ap, classified as a binary black hole merger, and suggested to be a possible multi-messenger event, based on the detection of a high energy neutrino and a TeV photon. We carried out a blind search at C band (4--8 GHz) over 0.3 deg$^2$ of the gamma-ray counterpart of S191216ap reported by the High-Altitude Water Cherenkov Observatory (HAWC). Our search, spanning three epochs over 130 days post-merger and having mean source-detection threshold of 75$\mu$Jy/beam (4$\sigma$), yielded 5 variable sources associated with AGN activity and no definitive counterpart of S191216ap. We find $, Comment: 9 pages, 5 figures, 2 tables, submitted to ApJ

One of the simplest examples of a robust heteroclinic cycle involves three saddle equilibria: each one is unstable to the next in turn, and connections from one to the next occur within invariant subspaces. Such a situation can be described by a third-order ordinary differential equation (ODE), and typical trajectories approach each equilibrium point in turn, spending progressively longer to cycle around the three points but never stopping. This cycle has been invoked as a model of cyclic competition between populations adopting three strategies, characterised as Rock, Paper and Scissors. When spatial distribution and mobility of the populations is taken into account, waves of Rock can invade regions of Scissors, only to be invaded by Paper in turn. The dynamics is described by a set of partial differential equations (PDEs) that has travelling wave (in one dimension) and spiral (in two dimensions) solutions. In this paper, we explore how the robust heteroclinic cycle in the ODE manifests itself in the PDEs. Taking the wavespeed as a parameter, and moving into a travelling frame, the PDEs reduce to a sixth-order set of ODEs, in which travelling waves are created in a Hopf bifurcation and are destroyed in three different heteroclinic bifurcations, depending on parameters, as the travelling wave approaches the heteroclinic cycle. We explore the three different heteroclinic bifurcations, none of which have been observed in the context of robust heteroclinic cycles previously. These results are an important step towards a full understanding of the spiral patterns found in two dimensions, with possible application to travelling waves and spirals in other population dynamics models., Comment: 36 pages, 8 figures

We present the completed KMOS$^\mathrm{3D}$ survey $-$ an integral field spectroscopic survey of 739, $\log(M_{\star}/M_{\odot})>9$, galaxies at $0.610.5$). We confirm our first year results of a high rotation dominated fraction (monotonic velocity gradient and $v_\mathrm{rot}$/$\sigma_0 > \sqrt{3.36}$) of 77% for the full KMOS$^\mathrm{3D}$ H$\alpha$sample. The rotation-dominated fraction is a function of both stellar mass and redshift with the strongest evolution measured over the redshift range of the survey for galaxies with $\log(M_{\star}/M_{\odot}), Comment: 26 pages, 18 figures, 8 tables; re-submitted after minor revisions to ApJ; associated data release at: http://www.mpe.mpg.de/ir/KMOS3D

While 3-SAT is NP-hard, 2-SAT is solvable in polynomial time. Austrin, Guruswami, and H\r{a}stad roved a result known as "$(2+\varepsilon)$-SAT is NP-hard" [FOCS'14/SICOMP'17]. They showed that the problem of distinguishing k-CNF formulas that are g-satisfiable (i.e. some assignment satisfies at least g literals in every clause) from those that are not even 1-satisfiable is NP-hard if $\frac{g}{k} < \frac{1}{2}$ and is in P otherwise. We study a generalisation of SAT on arbitrary finite domains, with clauses that are disjunctions of unary constraints, and establish analogous behaviour. Thus we give a dichotomy for a natural fragment of promise constraint satisfaction problems (PCSPs) on arbitrary finite domains. The hardness side is proved using the algebraic approach, via a new general NP-hardness criterion on polymorphisms of the problem, based on a gap version of the Layered Label Cover problem. We show that previously used criteria are insufficient -- the problem hence gives an interesting benchmark of algebraic techniques for proving hardness of approximation problems such as PCSPs., Comment: Full version of an ICALP 2020 paper

We present a simple resolution proof system for higher-order constrained Horn clauses (HoCHC)-a system of higher-order logic modulo theories-and prove its soundness and refutational completeness w.r.t. both standard and Henkin semantics. As corollaries, we obtain the compactness theorem and semi-decidability of HoCHC for semi-decidable background theories, and we prove that HoCHC satisfies a canonical model property. Moreover a variant of the well-known translation from higher-order to 1st-order logic is shown to be sound and complete for HoCHC in both semantics. We illustrate how to transfer decidability results for (fragments of) 1st-order logic modulo theories to our higher-order setting, using as example the Bernays-Schonflukel-Ramsey fragment of HoCHC modulo a restricted form of Linear Integer Arithmetic.

We describe a method to construct smooth and compactly supported solutions of 3D incompressible Euler equations and related models. The method is based on localizable Grad–Shafranov equations and is inspired by the recent result (Gavrilov in A steady Euler flow with compact support. Geom Funct Anal 29(1):90–197, [Gav19]).

Dirty paper codes are a powerful tool for combating known interference. However, there is a significant difference between knowing the transmitted interference sequence and knowing the received interference sequence, especially when the channel modifying the interference is uncertain. We present an upper bound on the capacity of a compound vector dirty paper channel where although an additive Gaussian sequence is known to the transmitter, the channel matrix between the interferer and receiver is uncertain but known to lie within a bounded set. Our bound is tighter than previous bounds in the low-SIR regime for the scalar version of the compound dirty paper channel and employs a construction that focuses on the relationship between the dimension of the message-bearing signal and the dimension of the additive state sequence. Additionally, a bound on the high-SNR behavior of the system is established., Comment: to be presented at ISIT 2013

We characterize functions of finite energy in the plane in terms of their traces on the lines that make up "graph paper" with squares of side length $mn$ for all $n$, and certain $\12-$order Sobolev norms on the graph paper lines. We also obtain analogous results for functions of finite energy on two classical fractals: the Sierpinski gasket and the Sierpinski carpet.

We develop a theoretical framework for rigid origami, and show how this framework can be used to connect rigid origami and results from cognate areas, such as the rigidity theory, graph theory, linkage folding and computer science. First, we give definitions on important concepts in rigid origami, then focus on how to describe the configuration space of a creased paper. The shape and 0-connectedness of the configuration space are analyzed using algebraic, geometric and numeric methods, where the key results from each method are gathered and reviewed., Comment: 29 pages, 7 figures, submitted to the Proceedings of the Royal Society A

The "Carbon Copy onto Dirty Paper" (CCDP) channel is the compound "writing on dirty paper" channel in which the channel output is obtained as the sum of the channel input, white Gaussian noise and a Gaussian state sequence randomly selected among a set possible realizations. The transmitter has non-causal knowledge of the set of possible state sequences but does not know which sequence is selected to produce the channel output. We study the capacity of the CCDP channel for two scenarios: (i) the state sequences are independent and identically distributed, and (ii) the state sequences are scaled versions of the same sequence. In the first scenario, we show that a combination of superposition coding, time-sharing and Gel'fand-Pinsker binning is sufficient to approach the capacity to within three bits per channel use for any number of possible state realizations. In the second scenario, we derive capacity to within four bits-per-channel-use for the case of two possible state sequences. This result is extended to the CCDP channel with any number of possible state sequences under certain conditions on the scaling parameters which we denote as "strong fading" regime. We conclude by providing some remarks on the capacity of the CCDP channel in which the state sequences have any jointly Gaussian distribution.

Automatic summarisation is a popular approach to reduce a document to its main arguments. Recent research in the area has focused on neural approaches to summarisation, which can be very data-hungry. However, few large datasets exist and none for the traditionally popular domain of scientific publications, which opens up challenging research avenues centered on encoding large, complex documents. In this paper, we introduce a new dataset for summarisation of computer science publications by exploiting a large resource of author provided summaries and show straightforward ways of extending it further. We develop models on the dataset making use of both neural sentence encoding and traditionally used summarisation features and show that models which encode sentences as well as their local and global context perform best, significantly outperforming well-established baseline methods., Comment: 11 pages, 6 figures

The rock-paper-scissors game is a paradigmatic model for biodiversity, with applications ranging from microbial populations to human societies. Research has shown, however, that mobility jeopardizes biodiversity by promoting the formation of spiral waves, especially if there is no conservation law in place for the total number of competing players. Firstly, we show that even if such a conservation law applies, mobility still jeopardizes biodiversity in the spatial rock-paper-scissors game if only a small fraction of links of the square lattice is randomly rewired. Secondly, we show that zealots are very effective in taming the amplitude of oscillations that emerge due to mobility and/or interaction randomness, and this regardless of whether the later is quenched or annealed. While even a tiny fraction of zealots brings significant benefits, at 5\% occupancy zealots practically destroy all oscillations regardless of the intensity of mobility, and regardless of the type and strength of randomness in the interaction structure. Interestingly, by annealed randomness the impact of zealots is qualitatively the same as by mobility, which highlights that fast diffusion does not necessarily destroy the coexistence of species, and that zealotry thus helps to recover the stable mean-field solution. Our results strengthen the important role of zealots in models of cyclic dominance, and they reveal fascinating evolutionary outcomes in structured populations that are a unique consequence of such uncompromising behavior., Comment: 6 two-column pages, 4 figures; accepted for publication in Physical Review E

Answer Set Programming (ASP) is a popular logic programming paradigm that has been applied for solving a variety of complex problems. Among the most challenging real-world applications of ASP are two industrial problems defined by Siemens: the Partner Units Problem (PUP) and the Combined Configuration Problem (CCP). The hardest instances of PUP and CCP are out of reach for state-of-the-art ASP solvers. Experiments show that the performance of ASP solvers could be significantly improved by embedding domain-specific heuristics, but a proper effective integration of such criteria in off-the-shelf ASP implementations is not obvious. In this paper the combination of ASP and domain-specific heuristics is studied with the goal of effectively solving real-world problem instances of PUP and CCP. As a byproduct of this activity, the ASP solver WASP was extended with an interface that eases embedding new external heuristics in the solver. The evaluation shows that our domain-heuristic-driven ASP solver finds solutions for all the real-world instances of PUP and CCP ever provided by Siemens. This paper is under consideration for acceptance in TPLP., Comment: Paper presented at the 32nd International Conference on Logic Programming (ICLP 2016), New York City, USA, 16-21 October 2016, 15 pages, LaTeX, 3 PDF figures

We present the first major data release of the largest single key-project in area carried out in open time with the Herschel Space Observatory. The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) is a survey of 600 deg^2 in five photometric bands - 100, 160, 250, 350 and 500 um - with the PACS and SPIRE cameras. In this paper and a companion paper (Bourne et al. 2016) we present the survey of three fields on the celestial equator, covering a total area of 161.6 deg^2 and previously observed in the Galaxy and Mass Assembly (GAMA) spectroscopic survey. This paper describes the Herschel images and catalogues of the sources detected on the SPIRE 250 um images. The 1-sigma noise for source detection, including both confusion and instrumental noise, is 7.4, 9.4 and 10.2 mJy at 250, 350 and 500 um. Our catalogue includes 120230 sources in total, with 113995, 46209 and 11011 sources detected at >4-sigma at 250, 350 and 500 um. The catalogue contains detections at >3-sigma at 100 and 160 um for 4650 and 5685 sources, and the typical noise at these wavelengths is 44 and 49 mJy. We include estimates of the completeness of the survey and of the effects of flux bias and also describe a novel method for determining the true source counts. The H-ATLAS source counts are very similar to the source counts from the deeper HerMES survey at 250 and 350 um, with a small difference at 500 um. Appendix A provides a quick start in using the released datasets, including instructions and cautions on how to use them., Accepted by MNRAS; data available at www.h-atlas.org

We consider (N,r) games of predation with N species and r < N prey and predators, acting in a cyclic way. Further basic reactions include reproduction, decay and diffusion over a regular grid, without a hard constraint on the occupation number per site, so in a "bosonic" implementation. For special combinations of N and r and appropriate parameter choices we observe games within games, that is, different coexisting games, depending on the spatial resolution. As concrete and simplest example we analyze the (6,3) game. Once the players segregate from a random initial distribution, domains emerge, which effectively play a (2,1)-game on the coarse scale of domain diameters, while agents inside the domains play (3,1) (rock-paper-scissors), leading to spiral formation with species chasing each other. As the (2,1)-game has a winner in the end, the coexistence of domains is transient, while agents inside the remaining domain coexist, until demographic fluctuations lead to extinction of all but one species in the very end. This means that we observe a dynamical generation of multiple space and time scales with emerging re-organization of players upon segregation, starting from a simple set of rules on the smallest scale and changed rules from the coarser perspective. These observations are based on Gillespie simulations. We discuss the deterministic limit derived from a van Kampen expansion. The linear stability analysis predicts the number of forming domains, their composition in terms of species; it also explains the instability of interfaces between domains, which drives their extinction; spiral patterns are identified as motion along heteroclinic cycles. The numerical solutions reproduce the observed patterns of the Gillespie simulations including even extinction events, so that the mean-field analysis here is very conclusive, which is due to the specific implementation of rules., 21 pages main text, 16 pages supplementary material (merged)

We have commenced a multi-year program, the Caltech-NRAO Stripe 82 Survey (CNSS), to search for radio transients with the Jansky VLA in the SDSS Stripe 82 region. The CNSS will deliver five epochs over the entire $\sim$270 deg$^2$ of Stripe 82, an eventual deep combined map with a rms noise of $\sim$40 $\mu$Jy and catalogs at a frequency of 3 GHz, and having a spatial resolution of 3". This first paper presents the results from an initial pilot survey of a 50 deg$^2$ region of Stripe 82, involving four epochs spanning logarithmic timescales between one week and 1.5 years, with the combined map having a median rms noise of 35 $\mu$Jy. This pilot survey enabled the development of the hardware and software for rapid data processing, as well as transient detection and follow-up, necessary for the full 270 deg$^2$ survey. Classification of variable and transient sources relied heavily on the wealth of multi-wavelength data in the Stripe 82 region, supplemented by repeated mapping of the region by the Palomar Transient Factory. $3.9^{+0.5}_{-0.9}$% of the detected point sources were found to vary by greater than 30%, consistent with similar studies at 1.4 GHz and 5 GHz. Multi-wavelength photometric data and light curves suggest that the variability is mostly due to shock-induced flaring in the jets of AGN. Although this was only a pilot survey, we detected two bona fide transients, associated with an RS CVn binary and a dKe star. Comparison with existing radio survey data revealed additional highly variable and transient sources on timescales between 5-20 years, largely associated with renewed AGN activity. The rates of such AGN possibly imply episodes of enhanced accretion and jet activity occurring once every $\sim$40,000 years in these galaxies. We compile the revised radio transient rates and make recommendations for future transient surveys and joint radio-optical experiments. (Abridged), Comment: 26 pages, 22 figures, 4 tables. Accepted for publication in the Astrophysical Journal. Data products (images, catalogs, tables, and light curves) available at http://tauceti.caltech.edu/stripe82 . A regularly-updated compilation of radio transient surveys is available at http://www.tauceti.caltech.edu/kunal/radio-transient-surveys/index.html

The origin of the low luminosity radio emission in radio-quiet AGN, is unknown. The detection of a positive correlation between the radio and X-ray emission would imply a jet-like origin, similar to that seen in `hard state' X-ray binary systems. In our previous work, we found no believable radio variability in the well known X-ray bright Seyfert 1 galaxy NGC 4051, despite large amplitude X-ray variability. In this study we have carefully re-analysed radio and X-ray observations using the same methods as our previous work, we again find no evidence for core radio variability. In direct contrast to our findings, another study claim significant radio variability and a distinctive anti-correlation between radio and X-ray data for the same source. The other study report only integral flux values and do not consider the effect of the changing array on the synthesised beam. In both our studies of NGC 4051 we have taken great care to account for the effect that the changing beam size has on the measured radio flux and as a result we are confident that our method gives more accurate values for the intrinsic core radio flux. However, the lack of radio variability we find is hard to reconcile because radio images of NGC 4051 do show jet-like structure. We suggest that the radio structures observed are likely the result of a previous period of higher radio activity and that the current level of radio emission from a compact nuclear jet is low., Comment: 13 pages, 8 Figures

We revisit two old and apparently little known papers by Basuev [2] [3] and show that the results contained there yield strong improvements on current lower bounds of the convergence radius of the Mayer series for continuous particle systems interacting via a very large class of stable and tempered potentials which includes the Lennard-Jones type potentials. In particular we analyze the case of the classical Lennard-Jones gas under the light of the Basuev scheme and, using also some new results [33] on this model recently obtained by one of us, we provide a new lower bound for the Mayer series convergence radius of the classical Lennard-Jones gas which improves by a factor of the order $10^5$ on the current best lower bound recently obtained in [17]., Comment: Final version as will appear in Comm. Math. Phys

We present observations from the Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa, observed in May and September 2010. Using two nights of drift scanning with PAPER's 60\arcdeg\ FWHM beam we have made a map covering the entire sky below +10 degrees declination with an effective center frequency of 145 MHz, a 70-MHz bandwidth, and a resolution of 26\arcmin. A 4800 square-degree region of this large map with the lowest Galactic emission reaches an RMS of 0.7 Jy. We establish an absolute flux scale using sources from the 160-MHz Culgoora catalog. Using the 408-MHz Molonglo Reference Catalog (MRC) as a finding survey, we identify counterparts to 480 sources in our maps, and compare our fluxes to the MRC and to 332 sources in the Culgoora catalog. For both catalogs, the ratio of PAPER to catalog flux averages to 1, with a standard deviation of 50%. This measured variation is consistent with comparisons between independent catalogs observed at different bands. The PAPER data represent new 145-MHz flux measurements for a large number of sources in the band expected to encompass cosmic reionization, and represents a significant step toward establishing a model for removing foregrounds to the reionization signal., Comment: 3 figures, 2 tables, machine readable table. Accepted by ApJL

An important goal for digital libraries is to enable researchers to more easily explore related work. While citation data is often used as an indicator of relatedness, in this paper we demonstrate that digital access records (e.g. http-server logs) can be used as indicators as well. In particular, we show that measures based on co-access provide better coverage than co-citation, that they are available much sooner, and that they are more accurate for recent papers., Comment: 2 pages, 3 postscript figures, to appear in proceedings of JCDL'07, additional and more detailed results can be found in arXiv:0704.2963v1

We present spectral fits to five epochs of the typical Type Ib supernova 1999dn using the generalized, non-LTE, stellar atmospheres code PHOENIX. Our goal is threefold: to determine basic physical properties of the supernova ejecta, such as velocity, temperature, and density gradients; to reproduce He I absorption lines by invoking non-thermal excitation; and, to investigate possible spectral signatures of hydrogen, especially a feature around 6200 Angstrom, which has been attributed to high velocity $H_\alpha$. Our models assume an atmosphere with uniform composition devoid of any hydrogen. Our model spectra fit the observed spectra well, successfully reproducing most of the features, including the prominent He I absorptions. The most plausible alternative to $H_\alpha$ as the source of the 6200 Angstrom feature is a blend of Fe II and Si II lines, which can be made stronger to fit the observed feature better by increasing the metallicity of the ejecta. High-metallicity models fit well at early epochs, but not as well as solar-metallicity models after maximum light. While this blend of metal lines is a reasonable explanation of the source of the 6200 Angstrom feature, it is still important to investigate hydrogen as the source; therefore, a second paper will present models that include a thin shell of hydrogen around the main composition structure., Comment: 24 pages, 11 figures, 2 tables, submitted to ApJ

Biomedical photoacoustic tomography, which can provide high resolution 3D soft tissue images based on the optical absorption, has advanced to the stage at which translation from the laboratory to clinical settings is becoming possible. The need for rapid image formation and the practical restrictions on data acquisition that arise from the constraints of a clinical workflow are presenting new image reconstruction challenges. There are many classical approaches to image reconstruction, but ameliorating the effects of incomplete or imperfect data through the incorporation of accurate priors is challenging and leads to slow algorithms. Recently, the application of Deep Learning, or deep neural networks, to this problem has received a great deal of attention. This paper reviews the literature on learned image reconstruction, summarising the current trends, and explains how these new approaches fit within, and to some extent have arisen from, a framework that encompasses classical reconstruction methods. In particular, it shows how these new techniques can be understood from a Bayesian perspective, providing useful insights. The paper also provides a concise tutorial demonstration of three prototypical approaches to learned image reconstruction. The code and data sets for these demonstrations are available to researchers. It is anticipated that it is in in vivo applications - where data may be sparse, fast imaging critical and priors difficult to construct by hand - that Deep Learning will have the most impact. With this in mind, the paper concludes with some indications of possible future research directions.

Observational longitudinal data on treatments and covariates are increasingly used to investigate treatment effects, but are often subject to time‐dependent confounding. Marginal structural models (MSMs), estimated using inverse probability of treatment weighting or the g‐formula, are popular for handling this problem. With increasing development of advanced causal inference methods, it is important to be able to assess their performance in different scenarios to guide their application. Simulation studies are a key tool for this, but their use to evaluate causal inference methods has been limited. This paper focuses on the use of simulations for evaluations involving MSMs in studies with a time‐to‐event outcome. In a simulation, it is important to be able to generate the data in such a way that the correct forms of any models to be fitted to those data are known. However, this is not straightforward in the longitudinal setting because it is natural for data to be generated in a sequential conditional manner, whereas MSMs involve fitting marginal rather than conditional hazard models. We provide general results that enable the form of the correctly specified MSM to be derived based on a conditional data generating procedure, and show how the results can be applied when the conditional hazard model is an Aalen additive hazard or Cox model. Using conditional additive hazard models is advantageous because they imply additive MSMs that can be fitted using standard software. We describe and illustrate a simulation algorithm. Our results will help researchers to effectively evaluate causal inference methods via simulation.

An Omnimagnet is an electromagnetic device that enables remote magnetic manipulation of devices such as medical implants and microrobots. It is comprised of three orthogonal nested solenoids with a ferromagnetic core at the center. Electrical current within the solenoids leads to Joule heating, resulting in undesired temperature increase within the Omnimagnet. If the temperature exceeds the melting point of the wire insulation, device failure will occur. Thus, a study of heat transfer within an Omnimagnet is a necessity, particularly to maximize the performance of the device. For the first time, a transient heat transfer model, that incorporates all three heat transfer modes, is proposed and validated with experimental data for an Omnimagnet with maximum root mean square error equal to 8${\%}$ (4$^{\circ}$C). This transient model is not computationally expensive. It is relatively easy to apply to Omnimagnets with different structures. The accuracy of this model depends on the accuracy of the input data. The code is applied to calculate the maximum safe operational time at a fixed input current or the maximum safe input current for a fixed time interval. The maximum safe operational time and maximum safe input current depend on size and structure of the Omnimagnet and the lowest melting point of all the Omnimagnet materials. A parametric study shows that increasing convective heat transfer during cooling, and during heating with low input currents, is an effective method to increase the maximum operational time of the Omnimagnet.The thermal model is also presented in a state-space equation format that can be used in a real-time Kalman filter current controller to avoid device failure due to excessive heating., Comment: This paper is under review in Asme journal of thermal science and engineering applications

The identification of the right order of the equation in applied fractional modeling plays an important role. In this paper we consider an inverse problem for determining the order of time fractional derivative in a subdiffusion equation with an arbitrary second order elliptic differential operator. We prove that the additional information about the solution at a fixed time instant at a monitoring location, as "the observation data", identifies uniquely the order of the fractional derivative., Comment: Version 1

Significance: Histopathological analysis of tissues is an essential tool for grading, staging, diagnosing and resecting cancers and other malignancies. Current histopathological techniques require substantial sample processing prior to staining with hematoxylin and eosin (H&E) dyes, to highlight nuclear and cellular morphology. Sample preparation and staining is resource-intensive and introduces potential for variability during sample preparation. Aim: We present a novel method for direct label-free histopathological assessment of formalin fixed paraffin embedded tissue blocks and thin tissue sections using a dual contrast photoacoustic remote sensing (PARS) microscopy system. Approach: To emulate the nuclear and cellular contrast of H&E staining, we leverage unique properties of the PARS system. Here the ultraviolet excitation of the PAARS microscope takes advantage of DNA's unique optical absorption to provide nuclear contrast analogous to hematoxylin staining of cell nucelli. Concurrently, the optical scattering contrast of the PARS detection system is leveraged to provide bulk tissue contrast analogous to eosin staining of cell membranes. Results: We demonstrate the efficacy of this technique by imaging human breast tissue and human skin tissues in formalin fixed paraffin embedded tissue blocks and frozen sections respectively. Salient nuclear and extra-nuclear features such as cancerous cells, glands and ducts, adipocytes, and stromal structures such as collagen. Conclusions. The presented dual contrast PARS microscope enables label-free visualizations of tissue with contrast and quality analogous to the current gold standard for histopathological analysis. Thus, the proposed system is well positioned to augment existing histopathological workflows, providing histological imaging directly on unstained tissue blocks and sections., 13 Pages, 3 figures

Reliable methods for the classification and quantification of quantum entanglement are fundamental to understanding its exploitation in quantum technologies. One such method, known as Separable Neural Network Quantum States (SNNS), employs a neural network inspired parameterisation of quantum states whose entanglement properties are explicitly programmable. Combined with generative machine learning methods, this ansatz allows for the study of very specific forms of entanglement which can be used to infer/measure entanglement properties of target quantum states. In this work, we extend the use of SNNS to mixed, multipartite states, providing a versatile and efficient tool for the investigation of intricately entangled quantum systems. We illustrate the effectiveness of our method through a number of examples, such as the computation of novel tripartite entanglement measures, and the approximation of ultimate upper bounds for qudit channel capacities., Comment: 14 pages, 7 figures

We show that bounded self-contracted curves are rectifiable in metric spaces with weak lower curvature bound in a sense we introduce in this article. This class of spaces is wide and includes, for example, finite-dimensional Alexandrov spaces of curvature bounded below and Berwald spaces of nonnegative flag curvature. (To be more precise, our condition is regarded as a strengthened doubling condition and holds also for a certain class of metric spaces with upper curvature bound.) We also provide the non-embeddability of large snowflakes into (balls in) metric spaces in the same class. We follow the strategy of the last author's previous paper based on the small rough angle condition, where spaces with upper curvature bound are considered. The results in this article show that such a strategy applies to spaces with lower curvature bound as well., Comment: A goofy mistake in formulations of Theorem 2 and 4 is fixed

Merging photonic structures and optoelectronic sensors into a single chip may yield a sensor‐on‐chip spectroscopic device that can measure the spectrum of matter. In this work, an on‐chip concurrent multiwavelength infrared (IR) sensor, which consists of a set of narrowband wavelength‐selective plasmonic perfect absorbers combined with pyroelectric sensors, where the response of each pyroelectric sensor is boosted only at the resonance of the nanostructured absorber, is proposed and realized. The proposed absorber, which is based on Wood's anomaly absorption from a 2D plasmonic square lattice, shows a narrowband polarization‐independent resonance (quality factor – Q of 73) with a nearly perfect absorptivity as high as 0.99 at normal incidence. The fabricated quad‐wavelength IR sensors exhibit four different narrowband spectral responses at normal incidence following the predesigned resonances in the mid‐wavelength infrared region that corresponds to the atmospheric window. The device can be applied for practical spectroscopic applications such as nondispersive IR sensors, IR chemical imaging devices, pyrometers, and spectroscopic thermography imaging., An approach to realize on‐chip multiwavelength infrared (IR) sensors is reported. As a proof of concept, the fabricated quad‐wavelength IR sensors exhibit four different narrowband spectral responses at normal incidence. The multiwavelength IR sensors can be applied for portable spectroscopic devices such as nondispersive IR sensors, IR chemical imaging, pyrometers, and thermography cameras.

2D materials with nontrivial energy bands are highly desirable for exploring various topological phases of matter, as low dimensionality opens unprecedented opportunities for manipulating the quantum states. Here, it is reported that monolayer (ML) dialkali-metal monoxides, in the well-known 2H-MoS$_2$ type lattice, host multiple symmetry-protected topological phases with emergent fermions, which can be effectively tuned by strain engineering. Based on first-principles calculations, it is found that in the equilibrium state, ML Na$_2$O is a 2D double Weyl semimetal, while ML K$_2$O is a 2D pseudospin-1 metal. These exotic topological states exhibit a range of fascinating effects, including universal optical absorbance, super Klein tunneling, and super collimation effect. By introducing biaxial or uniaxial strain, a series of quantum phase transitions between 2D double Weyl semimetal, 2D Dirac semimetal, 2D pseudospin-1 metal, and semiconductor phases can be realized. The results suggest monolayer dialkali-metal monoxides as a promising platform to explore fascinating physical phenomena associated with novel 2D emergent fermions., Comment: 7 pages, 4 figures

Funder: John Simon Guggenheim Memorial Foundation; doi: https://doi.org/10.13039/100005851, In this work we study the global solvability of the primitive equations for the atmosphere coupled to moisture dynamics with phase changes for warm clouds, where water is present in the form of water vapor and in the liquid state as cloud water and rain water. This moisture model contains closures for the phase changes condensation and evaporation, as well as the processes of auto conversion of cloud water into rainwater and the collection of cloud water by the falling rain droplets. It has been used by Klein and Majda in [17] and corresponds to a basic form of the bulk microphysics closure in the spirit of Kessler [16] and Grabowski and Smolarkiewicz [12]. The moisture balances are strongly coupled to the thermodynamic equation via the latent heat associated to the phase changes. In [14] we assumed the velocity field to be given and proved rigorously the global existence and uniqueness of uniformly bounded solutions of the moisture balances coupled to the thermodynamic equation. In this paper we present the solvability of a full moist atmospheric flow model, where the moisture model is coupled to the primitive equations of atmospherical dynamics governing the velocity field. For the derivation of a priori estimates for the velocity field we thereby use the ideas of Cao and Titi [6], who succeeded in proving the global solvability of the primitive equations.

UID/AMB/04085/2019 UID/FIS/00068/2019 PTDC/FIS-NAN/0909/2014 SFRH/BD/114674/2016 CTQ2017-90659-REDT (MCIUN, Spain). FEDER Andalucia 2014–2020” – UMA18-FEDERJA-279. Degradation technologies applied to emerging organic contaminants from human activities are one of the major water challenges in the contamination legacy. Triclosan is an emerging contaminant, commonly used as antibacterial agent in personal care products. Triclosan is stable, lipophilic and it is proved to have ecotoxicology effects in organics. This induces great concern since its elimination in wastewater treatment plants is not efficient and its by-products (e.g. methyl-triclosan, 2,4-dichlorophenol or 2,4,6-trichlorophenol) are even more hazardous to several environmental compartments. This work provides understanding of two different electrochemical reactors for the degradation of triclosan and its derivative by-products in effluent. A batch reactor and a flow reactor (mimicking a secondary settling tank in a wastewater treatment plant) were tested with two different working anodes: Ti/MMO and Nb/BDD. The degradation efficiency and kinetics were evaluated to find the best combination of current density, electrodes and set-up design. For both reactors the best electrode combination was achieved with Ti/MMO as anode. The batch reactor at 7 mA/cm2 during 4 h attained degradation rates below the detection limit for triclosan and 2,4,6-trichlorophenol and, 94% and 43% for 2,4-dichlorophenol and methyl triclosan, respectively. The flow reactor obtained, in approximately 1 h, degradation efficiencies between 41% and 87% for the four contaminants. This study suggests an alternative technology for emerging organic contaminants degradation, since the combination of a low current density with the flow and matrix induced disturbance increases and speeds up the compounds’ elimination in a real environmental matrix. proof published

Significance: Spatial frequency-domain imaging (SFDI) is a powerful technique for mapping tissue oxygen saturation over a wide field of view. However, current SFDI methods either require a sequence of several images with different illumination patterns or, in the case of single-snapshot optical properties (SSOP), introduce artifacts and sacrifice accuracy. Aim: We introduce OxyGAN, a data-driven, content-aware method to estimate tissue oxygenation directly from single structured-light images. Approach: OxyGAN is an end-to-end approach that uses supervised generative adversarial networks. Conventional SFDI is used to obtain ground truth tissue oxygenation maps for ex vivo human esophagi, in vivo hands and feet, and an in vivo pig colon sample under 659- and 851-nm sinusoidal illumination. We benchmark OxyGAN by comparing it with SSOP and a two-step hybrid technique that uses a previously developed deep learning model to predict optical properties followed by a physical model to calculate tissue oxygenation. Results: When tested on human feet, cross-validated OxyGAN maps tissue oxygenation with an accuracy of 96.5%. When applied to sample types not included in the training set, such as human hands and pig colon, OxyGAN achieves a 93% accuracy, demonstrating robustness to various tissue types. On average, OxyGAN outperforms SSOP and a hybrid model in estimating tissue oxygenation by 24.9% and 24.7%, respectively. Finally, we optimize OxyGAN inference so that oxygenation maps are computed ∼10 times faster than previous work, enabling video-rate, 25-Hz imaging. Conclusions: Due to its rapid acquisition and processing speed, OxyGAN has the potential to enable real-time, high-fidelity tissue oxygenation mapping that may be useful for many clinical applications.

As the COVID-19 outbreak is developing the two most frequently reported statistics seem to be the raw confirmed case and case fatalities counts. Focusing on Italy, one of the hardest hit countries, we look at how these two values could be put in perspective to reflect the dynamics of the virus spread. In particular, we find that merely considering the confirmed case counts would be very misleading. The number of daily tests grows, while the daily fraction of confirmed cases to total tests has a change point. It (depending on region) generally increases with strong fluctuations till (around, depending on region) 15th-22nd March and then decreases linearly after. Combined with the increasing trend of daily performed tests, the raw confirmed case counts are not representative of the situation and are confounded with the sampling effort. This we observe when regressing on time the logged fraction of positive tests and for comparison the logged raw confirmed count. Hence, calibrating model parameters for this virus's dynamics should not be done based only on confirmed case counts (without rescaling by the number of tests), but take also fatalities and hospitalization count under consideration as variables not prone to be distorted by testing efforts. Furthermore, reporting statistics on the national level does not say much about the dynamics of the disease, which are taking place at the regional level. These findings are based on the official data of total death counts up to 15th April 2020 released by ISTAT and up to 10th May 2020 for the number of cases. In this work we do not fit models but we rather investigate whether this task is possible at all. This work also informs about a new tool to collect and harmonize official statistics coming from different sources in the form of a package for the R statistical environment and presents the COVID-19 Data Hub., Comment: updated references

Indirect imaging problems in biomedical optics generally require repeated evaluation of forward models of radiative transport, for which Monte Carlo is accurate yet computationally costly. We develop a novel approach to reduce this bottleneck which has significant implications for quantitative tomographic imaging in a variety of medical and industrial applications. Using Monte Carlo we compute a fully stochastic gradient of an objective function for a given imaging problem. Leveraging techniques from the machine learning community we then adaptively control the accuracy of this gradient throughout the iterative inversion scheme, in order to substantially reduce computational resources at each step. For example problems of Quantitative Photoacoustic Tomography and Ultrasound Modulated Optical Tomography, we demonstrate that solutions are attainable using a total computational expense that is comparable to (or less than) that which is required for a single high accuracy forward run of the same Monte Carlo model. This approach demonstrates significant computational savings when approaching the full non-linear inverse problem of optical property estimation using stochastic methods., Comment: 24 Pages, 11 figures

According to the current perception, symptomatic, presymptomatic, and asymptomatic infectious persons can infect the healthy population susceptible to the SARS-Cov-2. More importantly, various reports indicate that the number of asymptomatic cases can be several-fold higher than the reported symptomatic cases. In this article, we take the reported cases in India and various states within the country till September 1, as the specimen to understand the progression of the COVID-19. Employing a modified SEIRD model, we predict the spread of COVID-19 by the symptomatic as well as asymptomatic infectious population. Considering reported infection primarily due to symptomatic we compare the model predicted results with the available data to estimate the dynamics of the asymptomatically infected population. Our data indicate that in the absence of the asymptomatic infectious population, the number of symptomatic cases would have been much less. Therefore, the current progress of the symptomatic infection can be reduced by quarantining the asymptomatically infectious population via extensive or random testing. This study is motivated strictly towards academic pursuit; this theoretical investigation is not meant for influencing policy decisions or public health practices.

We present the CHAOS-7 model of the time-dependent near-Earth geomagnetic field between 1999 and 2020 based on magnetic field observations collected by the low-Earth orbit satellites Swarm, CryoSat-2, CHAMP, SAC-C and Ørsted, and on annual differences of monthly means of ground observatory measurements. The CHAOS-7 model consists of a time-dependent internal field up to spherical harmonic degree 20, a static internal field which merges to the LCS-1 lithospheric field model above degree 25, a model of the magnetospheric field and its induced counterpart, estimates of Euler angles describing the alignment of satellite vector magnetometers, and magnetometer calibration parameters for CryoSat-2. Only data from dark regions satisfying strict geomagnetic quiet-time criteria (including conditions on IMF $$B_z$$ B z and $$B_y$$ B y at all latitudes) were used in the field estimation. Model parameters were estimated using an iteratively reweighted regularized least-squares procedure; regularization of the time-dependent internal field was relaxed at high spherical harmonic degree compared with previous versions of the CHAOS model. We use CHAOS-7 to investigate recent changes in the geomagnetic field, studying the evolution of the South Atlantic weak field anomaly and rapid field changes in the Pacific region since 2014. At Earth’s surface a secondary minimum of the South Atlantic Anomaly is now evident to the south west of Africa. Green’s functions relating the core–mantle boundary radial field to the surface intensity show this feature is connected with the movement and evolution of a reversed flux feature under South Africa. The continuing growth in size and weakening of the main anomaly is linked to the westward motion and gathering of reversed flux under South America. In the Pacific region at Earth’s surface between 2015 and 2018 a sign change has occurred in the second time derivative (acceleration) of the radial component of the field. This acceleration change took the form of a localized, east–west oriented, dipole. It was clearly recorded on ground, for example at the magnetic observatory at Honolulu, and was seen in Swarm observations over an extended region in the central and western Pacific. Downward continuing to the core–mantle boundary, we find this event originated in field acceleration changes at low latitudes beneath the central and western Pacific in 2017.

Motivation The location, timing and abundance of gene expression (both mRNA and proteins) within a tissue define the molecular mechanisms of cell functions. Recent technology breakthroughs in spatial molecular profiling, including imaging-based technologies and sequencing-based technologies, have enabled the comprehensive molecular characterization of single cells while preserving their spatial and morphological contexts. This new bioinformatics scenario calls for effective and robust computational methods to identify genes with spatial patterns. Results We represent a novel Bayesian hierarchical model to analyze spatial transcriptomics data, with several unique characteristics. It models the zero-inflated and over-dispersed counts by deploying a zero-inflated negative binomial model that greatly increases model stability and robustness. Besides, the Bayesian inference framework allows us to borrow strength in parameter estimation in a de novo fashion. As a result, the proposed model shows competitive performances in accuracy and robustness over existing methods in both simulation studies and two real data applications. Availability and implementation The related R/C++ source code is available at https://github.com/Minzhe/BOOST-GP. Supplementary information Supplementary data are available at Bioinformatics online.

The application of network analysis has found great success in a wide variety of disciplines; however, the popularity of these approaches has revealed the difficulty in handling networks whose complexity scales rapidly. One of the main interests in network analysis is the online detection of anomalous behaviour. To overcome the curse of dimensionality, we introduce a network surveillance method bringing together network modelling and statistical process control. Our approach is to apply multivariate control charts based on exponential smoothing and cumulative sums in order to monitor networks determined by temporal exponential random graph models (TERGM). This allows us to account for temporal dependence, while simultaneously reducing the number of parameters to be monitored. The performance of the proposed charts is evaluated by calculating the average run length for both simulated and real data. To prove the appropriateness of the TERGM to describe network data, some measures of goodness of fit are inspected. We demonstrate the effectiveness of the proposed approach by an empirical application, monitoring daily flights in the United States to detect anomalous patterns.

Many attempts have been made to provide Quantum Field Theory with conceptually clear and mathematically rigorous foundations; remarkable examples are the Bohmian and the algebraic perspectives respectively. In this essay we introduce the dissipative approach to QFT, a new alternative formulation of the theory explaining the phenomena of particle creation and annihilation starting from nonequilibrium thermodynamics. It is shown that DQFT presents a rigorous mathematical structure, and a clear particle ontology, taking the best from the mentioned perspectives. Finally, after the discussion of its principal implications and consequences, we compare it with the main Bohmian QFTs implementing a particle ontology., European Journal for Philosophy of Science, 11 (1)

The sample paths of white noise are proved to be elements of certain Besov spaces with dominating mixed smoothness. Unlike in isotropic spaces, here the regularity does not get worse with increasing space dimension. Consequently, white noise is actually much smoother than the known sharp regularity results in isotropic spaces suggest. An application of our techniques yields new results for the regularity of solutions of Poisson and heat equation on the half space with boundary noise. The main novelty is the flexible treatment of the interplay between the singularity at the boundary and the smoothness in tangential, normal and time direction., Comment: 29 pages

A large scientific community depends on the precise modelling of complex processes in particle cascades in various types of matter. These models are used most prevalently in cosmic-ray physics, astrophysical-neutrino physics, and gamma-ray astronomy. In this white paper, we summarize the necessary steps to ensure the evolution and future availability of optimal simulation tools. The purpose of this document is not to act as a strict blueprint for next-generation software, but to provide guidance for the vital aspects of its design. The topics considered here are driven by physics and scientific applications. Furthermore, the main consequences of implementation decisions on performance are outlined. We highlight the computational performance as an important aspect guiding the design since future scientific applications will heavily depend on an efficient use of computational resources., Comment: 17 pages, 2 figures, content matches published version

In this note, we consider the antibandwidth problem, also known as dual bandwidth problem, separation problem and maximum differential coloring problem. Given a labeled graph (i.e., a numbering of the vertices of a graph), the antibandwidth of a node is defined as the minimum absolute difference of its labeling to the labeling of all its adjacent vertices. The goal in the antibandwidth problem is to find a labeling maximizing the antibandwidth. The problem is NP-hard in general graphs and has applications in diverse areas like scheduling, radio frequency assignment, obnoxious facility location and map-coloring. There has been much work on deriving theoretical bounds for the problem and also in the design of metaheuristics in recent years. However, the optimality gaps between the best known solution values and reported upper bounds for the HarwellBoeing Matrix-instances, which are the commonly used benchmark instances for this problem, are often very large (e.g., up to 577%). Moreover, only for three of these 24 instances, the optimal solution is known, leading the authors of a state-of-the-art heuristic to conclude “HarwellBoeing instances are actually a challenge for modern heuristic methods”. The upper bounds reported in literature are based on the theoretical bounds involving simple graph characteristics, i.e., size, order and degree, and a mixed-integer programming (MIP) model. We present new MIP models for the problem, together with valid inequalities, and design a branch-and-cut algorithm and an iterative solution algorithm based on them. These algorithms also include two starting heuristics and a primal heuristic. We also present a constraint programming approach, and calculate upper bounds based on the stability number and chromatic number. Our computational study shows that the developed approaches allow to find the proven optimal solution for eight instances from literature, where the optimal solution was unknown and also provide reduced gaps for eleven additional instances, including improved solution values for seven instances, the largest optimality gap is now 46%.

Compound refractive lenses (CRL) offer outstanding potential for probe formation in X-ray micro- and nano-analysis as well as for imaging applications, but effective aberration-free focusing of X-rays is a key consideration when high spatial resolution and intensities are required. Available refractive optics have non-negligible aberrations which degrade their final performance. We describe a framework based on physical-optics for simulating the effect of imperfect CRL upon an X-ray beam, taking into account measured phase errors obtained from at-wavelength metrology. We model, with increasing complexity, a CRL stack as a single thin phase element, then as a more realistic compound-element including absorption and thickness effects and finally, we add realistic optical imperfections to the CRL. Coherent and partially-coherent simulations using the "Synchrotron Radiation Workshop" (SRW) are used to evaluate the different models, the effects of the phase errors and to check the validity of the design equations and suitability of the figures of merit., Comment: 14 pages, 10 figures, 5 tables