title: Integrability of Black Hole Orbits in Maximal Supergravity

abstract: We study the dynamics of a pair of extremal (half-BPS) black holes in $\mathcal{N}=8$ supergravity, as a potentially solvable model of gravitational dynamics. As a diagnosis of hidden symmetries, we ask whether the perihelion of the orbits precesses over time. We consider black hole charge vectors with arbitrary misalignment. First, we use scattering amplitude methods to compute the leading post-Newtonian correction for general mass ratios. This computation is greatly simplified by introducing a suitable on-shell superspace. Second, we study the probe limit to all orders in velocity and Newton's constant through a ten-dimensional brane setup. In all cases we find no precession. We relate this to the absence of scalar triangle integrals.


title: Dense-gas tracers and carbon isotopes in five 2.5<z<4 lensed dusty star forming galaxies from the SPT SMG sample

abstract: The origin of the high SFR observed in high-z dusty star-forming galaxies is still unknown. Large fractions of dense molecular gas might provide part of the explanation, but there are few observational constraints on the amount of dense gas in high-z systems dominated by star formation. We present the results of our ALMA program targeting dense-gas tracers (HCN(5-4), HCO+(5-4), and HNC(5-4)) in 5 strongly lensed galaxies from the SPT SMG sample. We detected two of these lines (SNR>5) in SPT-125-47 at z=2.51 and tentatively detected all three (SNR~3) in SPT0551-50 at z=3.16. Since a significant fraction of our target lines is not detected, we developed a statistical method to derive unbiased mean properties taking into account both detections and non-detections. On average, the HCN(5-4) and HCO+(5-4) luminosities of our sources are a factor of ~1.7 fainter than expected, based on the local L'HCN(5-4)-LIR relation, but this offset corresponds to only ~2 sigma. We find that both the HCO+/HCN and HNC/HCN flux ratios are compatible with unity. The first ratio is expected for PDRs while the second is consistent with PDRs or XDRs and/or mid-IR pumping of HNC. Our sources are at the high end of the local relation between the star formation efficiency, determined using the LIR/[CI] and LIR/CO ratios, and the dense gas fraction, estimated using the HCN/[CI] and HCN/CO ratios. In SPT0125-47, we found that the velocity profiles of the lines tracing dense (HCN, HCO+) and lower-density (CO, [CI]) gas are similar. In addition to these lines, we obtained one robust and one tentative detection of 13CO(4-3) and found an average I12CO(4-3)/I13CO(4-3) flux ratio of 26.1$_{-3.5}^{+4.5}$, indicating a young but not pristine interstellar medium. We argue that the combination of large and slightly enriched gas reservoirs and high dense-gas fractions could explain the prodigious star formation in these systems.


title: Relations in doubly laced crystal graphs via discrete Morse theory

abstract: We study the combinatorics of crystal graphs given by highest weight representations of types $A_{n}, B_{n}, C_{n}$, and $D_{n}$, uncovering new relations that exist among crystal operators. Much structure in these graphs has been revealed by local relations given by Stembridge and Sternberg. However, there exist relations among crystal operators that are not implied by Stembridge or Sternberg relations. Viewing crystal graphs as edge colored posets, we use poset topology to study them. Using the lexicographic discrete Morse functions of Babson and Hersh, we relate the M\"obius function of a given interval in a crystal poset of simply laced or doubly laced type to the types of relations that can occur among crystal operators within this interval. For a crystal of a highest weight representation of finite classical Cartan type, we show that whenever there exists an interval whose M\"obius function is not equal to -1, 0, or 1, there must be a relation among crystal operators within this interval not implied by Stembridge or Sternberg relations. As an example of an application, this yields relations among crystal operators in type $C_{n}$ that were not previously known. Additionally, by studying the structure of Sternberg relations in the doubly laced case, we prove that crystals of highest weight representations of types $B_{2}$ and $C_{2}$ are not lattices.


title: Prices, Profits, Proxies, and Production

abstract: This paper studies nonparametric identification and counterfactual bounds for heterogeneous firms that can be ranked in terms of productivity. Our approach works when quantities and prices are latent, rendering standard approaches inapplicable. Instead, we require observation of profits or other optimizing-values such as costs or revenues, and either prices or price proxies of flexibly chosen variables. We extend classical duality results for price-taking firms to a setup with discrete heterogeneity, endogeneity, and limited variation in possibly latent prices. Finally, we show that convergence results for nonparametric estimators may be directly converted to convergence results for production sets.


title: On the Balmer spectrum for compact Lie groups

abstract: We study the Balmer spectrum of the category of finite G-spectra for a compact Lie group G, extending the work for finite G by Strickland, Balmer-Sanders, Barthel-Hausmann-Naumann-Nikolaus-Noel-Stapleton and others. We give a description of the underlying set of the spectrum and show that the Balmer topology is completely determined by the inclusions between the prime ideals and the topology on the space of closed subgroups of G. Using this, we obtain a complete description of this topology for all abelian compact Lie groups and consequently a complete classification of thick tensor-ideals. For general compact Lie groups we obtain such a classification away from a finite set of primes p.


title: Blockchain-based Smart Contracts - Applications and Challenges

abstract: A blockchain-based smart contract or a "smart contract" for short, is a computer program intended to digitally facilitate the negotiation or contractual terms directly between users when certain conditions are met. With the advance in blockchain technology, smart contracts are being used to serve a wide range of purposes ranging from self-managed identities on public blockchains to automating business collaboration on permissioned blockchains. In this paper, we present a comprehensive survey of smart contracts with a focus on existing applications and challenges they face.


title: End-to-End Content and Plan Selection for Data-to-Text Generation

abstract: Learning to generate fluent natural language from structured data with neural networks has become an common approach for NLG. This problem can be challenging when the form of the structured data varies between examples. This paper presents a survey of several extensions to sequence-to-sequence models to account for the latent content selection process, particularly variants of copy attention and coverage decoding. We further propose a training method based on diverse ensembling to encourage models to learn distinct sentence templates during training. An empirical evaluation of these techniques shows an increase in the quality of generated text across five automated metrics, as well as human evaluation.


title: Exact results for the infinite supersymmetric extensions of the infinite square well

abstract: One-dimensional potentials defined by $V^{(S)}(x) =S(S+1) \hbar^2 \pi^2 /[2ma^2\sin^2(\pi x/a)]$ (for integer $S$) arise in the repeated supersymmetrization of the infinite square well, here defined over the region $(0,a)$. We review the derivation of this hierarchy of potentials and then use the methods of supersymmetric quantum mechanics, as well as more familiar textbook techniques, to derive compact closed-form expressions for the normalized solutions, $\psi_n^{(S)}(x)$, for all $V^{(S)}(x)$ in terms of well-known special functions in a pedagogically accessible manner. We also note how the solutions can be obtained as a special case of a family of shape-invariant potentials, the trigonometric P\"oschl-Teller potentials, which can be used to confirm our results. We then suggest additional avenues for research questions related to, and pedagogical applications of, these solutions, including the behavior of the corresponding momentum-space wave functions $\phi_n^{(S)}(p)$ for large $|p|$ and general questions about the supersymmetric hierarchies of potentials which include an infinite barrier.


title: Pattern Formation in a Slowly Flattening Spherical Cap: Delayed Bifurcation

abstract: This article describes a reduction of a nonautonomous Brusselator reaction-diffusion system of partial differential equations on a spherical cap with time dependent curvature using the method of centre manifold reduction. Parameter values are chosen such that the change in curvature would cross critical values which would change the stability of the patternless solution in the constant domain case. The evolving domain functions and quasi-patternless solutions are derived as well as a method to obtain this nonautonomous normal form. The coefficients of such a normal form are computed and the reduction solutions are compared to numerical solutions.


title: Deep Inertial Poser: Learning to Reconstruct Human Pose from Sparse Inertial Measurements in Real Time

abstract: We demonstrate a novel deep neural network capable of reconstructing human full body pose in real-time from 6 Inertial Measurement Units (IMUs) worn on the user's body. In doing so, we address several difficult challenges. First, the problem is severely under-constrained as multiple pose parameters produce the same IMU orientations. Second, capturing IMU data in conjunction with ground-truth poses is expensive and difficult to do in many target application scenarios (e.g., outdoors). Third, modeling temporal dependencies through non-linear optimization has proven effective in prior work but makes real-time prediction infeasible. To address this important limitation, we learn the temporal pose priors using deep learning. To learn from sufficient data, we synthesize IMU data from motion capture datasets. A bi-directional RNN architecture leverages past and future information that is available at training time. At test time, we deploy the network in a sliding window fashion, retaining real time capabilities. To evaluate our method, we recorded DIP-IMU, a dataset consisting of $10$ subjects wearing 17 IMUs for validation in $64$ sequences with $330\,000$ time instants; this constitutes the largest IMU dataset publicly available. We quantitatively evaluate our approach on multiple datasets and show results from a real-time implementation. DIP-IMU and the code are available for research purposes.


title: Equivalence of generics

abstract: Given a countable transitive model of set theory and a partial order contained in it, there is a natural countable Borel equivalence relation on generic filters over the model; two are equivalent if they yield the same generic extension. We examine the complexity of this equivalence relation for various partial orders, with particular focus on Cohen and random forcing. We prove, amongst other results, that the former is an increasing union of countably many hyperfinite Borel equivalence relations, while the latter is neither amenable nor treeable.


title: A direct approach to imaging in a waveguide with perturbed geometry

abstract: We introduce a direct, linear sampling approach to imaging in an acoustic waveguide with sound hard walls. The waveguide terminates at one end and has unknown geometry due to compactly supported wall deformations. The goal of imaging is to determine these deformations and to identify localized scatterers in the waveguide, using a remote array of sensors that emits time harmonic probing waves and records the echoes. We present a theoretical analysis of the imaging approach and illustrate its performance with numerical simulations.


title: Low-dissipating push-pull SQUID amplifier for TES detector readout

abstract: We suggest SQUID amplifier configuration intended for millikelvin refrigerators where cooling power is limited and hence high power efficiency is desirable. The configuration operates the SQUIDs in class-AB rather than the traditional class-A. A proof-of principle demonstration at T = 4.2 K is presented.


title: An Inductive Logic Programming Approach to Validate Hexose Binding Biochemical Knowledge

abstract: Hexoses are simple sugars that play a key role in many cellular pathways, and in the regulation of development and disease mechanisms. Current protein-sugar computational models are based, at least partially, on prior biochemical findings and knowledge. They incorporate different parts of these findings in predictive black-box models. We investigate the empirical support for biochemical findings by comparing Inductive Logic Programming (ILP) induced rules to actual biochemical results. We mine the Protein Data Bank for a representative data set of hexose binding sites, non-hexose binding sites and surface grooves. We build an ILP model of hexose-binding sites and evaluate our results against several baseline machine learning classifiers. Our method achieves an accuracy similar to that of other black-box classifiers while providing insight into the discriminating process. In addition, it confirms wet-lab findings and reveals a previously unreported Trp-Glu amino acids dependency.


title: First-principles investigation of the effect of substitution and surface adsorption on magnetostrictive properties of Fe-Ga alloys

abstract: Materials with large magnetostriction are widely used in sensors, actuators, micro electromechanical systems, and energy-harvesters. Binary Fe-Ga alloys (Galfenol) are the most promising rare-earth-free candidates combining numerous advantages such as low saturation magnetic field (~200 Oe), excellent ductility and low cost, while further improving their performance is imperative for practical applications. Using density functional theory calculation, we report results of the effect of substituting small amount of additional elements X (eg. X = Ag, Pd and Cu) on magnetostriction of Fe-Ga alloys, and find that it may double the magnetostriction with a substitutional percentage of only 1.6%. Moreover, adsorbents with high chemical activity (eg. O or Os atoms) may affect the surface energy of different face-orientations of Fe-Ga alloys, indicating proper surface treatments are necessary to tune the alignment of Fe-Ga grains to achieve better performance. These results may be helpful to further optimize the magnetostrictive properties of Fe-Ga alloys for device applications.


title: Velocity Dependence from Resonant Self-Interacting Dark Matter

abstract: The dark matter density distribution in small-scale astrophysical objects may indicate that dark matter is self-interacting, while observations from clusters of galaxies suggest that the corresponding cross section depends on the velocity. Using a model-independent approach, we show that resonant self-interacting dark matter (RSIDM) can naturally explain such a behavior. In contrast to what is often assumed, this does not require a light mediator. We present explicit realizations of this mechanism and discuss the corresponding astrophysical constraints.


title: Ramanujan complexes and Golden Gates in PU(3)

abstract: In a seminal series of papers from the 80's, Lubotzky, Phillips and Sarnak applied the Ramanujan-Petersson Conjecture for $GL_{2}$ (Deligne's theorem), to a special family of arithmetic lattices, which act simply-transitively on the Bruhat-Tits trees associated with $SL_{2}(\mathbb{Q}_{p})$. As a result, they obtained explicit Ramanujan Cayley graphs from $PSL_{2}\left(\mathbb{F}_{p}\right)$, as well as optimal topological generators ("Golden Gates") for the compact Lie group $PU(2)$. In higher dimension, the naive generalization of the Ramanujan Conjecture fails, due to the phenomenon of endoscopic lifts. In this paper we overcome this problem for $PU_{3}$ by constructing a family of arithmetic lattices which act simply-transitively on the Bruhat-Tits buildings associated with $SL_{3}(\mathbb{Q}_{p})$ and $SU_{3}(\mathbb{Q}_{p})$, while at the same time do not admit any representation which violates the Ramanujan Conjecture. This gives us Ramanujan complexes from $PSL_{3}(\mathbb{F}_{p})$ and $PSU_{3}(\mathbb{F}_{p})$, as well as golden gates for $PU(3)$.


title: The role of nonlinearities in topological protection: using magnetically coupled fidget spinners

abstract: We investigate and experimentally observe the existence of topologically protected interface modes in a one-dimensional mechanical lattice, and we report on the effect of nonlinearities on topological protection. The lattice consists of a one-dimensional array of spinners with nearest neighbor coupling resulting from magnetic interactions. The distance between the spinners is spatially modulated to obtain a diatomic configuration, and to produce a non-trivial interface by breaking spatial inversion symmetry. For small amplitudes of motion, the interactions are approximately linear, and the system supports topologically protected interface modes at frequencies inside the bulk bandgaps of the lattice. Nonlinearities induced by increasing amplitude of motion cause the interface modes to shift outside the bandgaps and merge with the bulk bands. The resulting edge-to-bulk transition causes the extinction of the topologically protected interface mode and extends it to the entire length of the chain. Such transition is predicted by analytical calculations and verified by experimental observations. The paper thus investigates the existence of topologically protected interface modes obtained through broken spatial inversion symmetry, and documents their lack of robustness in the presence of nonlinearities.


title: Estimation of cosmic-muon flux attenuation by Monserrate Hill in Bogota

abstract: When designing an experimental setup for measuring the flux of cosmic ray muons through a geological structure, it is crucial to make detailed estimates about the expected muon flux attenuation by the structure studied. In this way, the coverage area of the detectors, the data acquisition system and the time length of data taking can be planned accordingly to the aimed accuracy of the measurements. We present here a computational study of the cosmic muon flux at the location of Monserrate Hill in Bogota, Colombia, using the CORSIKA cosmic ray simulator. We use ArcGIS package to estimate the mean depth muons have to travel through the mountain to reach the detector. We consider two different uniform densities of rock for the mountain to make estimates of the muon flux. Through the interpolation of stopping power available data, one can then determine the attenuation as a function of the muon incident energy. Our final result provides an estimation for the muon flux attenuation as a function of the location of the mountain where the muons pass through. We also supply estimated values for the time needed to reach an accuracy of about 3$\%$ in the measurement of the muon flux through Monserrate Hill. Our strategy can be easily adapted for muon flux studies at any other location in the world.


title: Fingerprint of the first stars: multi-enriched extremely metal-poor stars in the TOPoS survey

abstract: Extremely metal poor (EMP) stars in the Milky Way inherited the chemical composition of the gas out of which they formed. They therefore carry the chemical fingerprint of the first stars in their spectral lines. It is commonly assumed that EMP stars form from gas that was enriched by only one progenitor supernova ('mono-enriched'). However, recent numerical simulations show that the first stars form in small clusters. Consequently, we expect several supernovae to contribute to the abundances of an EMP star ('multi-enriched'). We analyse seven recently observed EMP stars from the TOPoS survey by applying the divergence of the chemical displacement and find that J1035+0641 is mono-enriched ($p_{mono}=53\%$) and J1507+0051 is multi-enriched ($p_{mono}=4\%$). For the remaining five stars we can not make a distinct prediction ($p_{mono} \lesssim 50\%$) due to theoretical and observational uncertainties. Further observations in the near-UV will help to improve our diagnostic and therefore contribute to constrain the nature of the first stars.


title: Training Generative Adversarial Networks with Binary Neurons by End-to-end Backpropagation

abstract: We propose the BinaryGAN, a novel generative adversarial network (GAN) that uses binary neurons at the output layer of the generator. We employ the sigmoid-adjusted straight-through estimators to estimate the gradients for the binary neurons and train the whole network by end-to-end backpropogation. The proposed model is able to directly generate binary-valued predictions at test time. We implement such a model to generate binarized MNIST digits and experimentally compare the performance for different types of binary neurons, GAN objectives and network architectures. Although the results are still preliminary, we show that it is possible to train a GAN that has binary neurons and that the use of gradient estimators can be a promising direction for modeling discrete distributions with GANs. For reproducibility, the source code is available at https://github.com/salu133445/binarygan .


title: Homogeneous and Isotropic Spacetime in Conformal Scalar-Tensor Gravity

abstract: The background field equations for homogeneous and isotropic spacetime are derived in conformal scalar-tensor gravity. The background temporal evolution is entirely driven by the dynamical evolution of the scalar field, i.e. particle masses, and satisfies an equation which is identical in form to the Friedmann equation of the standard cosmological model in general relativity. In a static background spacetime the scalar field (logarithmic) time-derivative replaces the `Hubble function'. It is also shown that linear perturbations are governed by equations which are identical to those obtained in general relativity, but with their evolution stemming from the scalar field dynamics.


title: The geometro-hydrodynamical formalism of quantum spinning particle

abstract: We present the development of the realistic geometro-hydrodynamical formalism of quantum mechanics for the spinning particle, that involves the vortical flows and is based on the idea, that the spinor wave represents a new type of physical field, propagating in space-time and influencing the corpuscule embedded in the wave. We assume that this new field created by the intrinsic spin of the particle and can produce new intrinsic forces and torques affecting the particle with spin. We derive and discuss physical effects of this field.


title: Detecting Core-Periphery Structures by Surprise

abstract: Detecting the presence of mesoscale structures in complex networks is of primary importance. This is especially true for financial networks, whose structural organization deeply affects their resilience to events like default cascades, shocks propagation, etc. Several methods have been proposed, so far, to detect communities, i.e. groups of nodes whose connectivity is significantly large. Communities, however do not represent the only kind of mesoscale structures characterizing real-world networks: other examples are provided by bow-tie structures, core-periphery structures and bipartite structures. Here we propose a novel method to detect statistically-signifcant bimodular structures, i.e. either bipartite or core-periphery ones. It is based on a modification of the surprise, recently proposed for detecting communities. Our variant allows for bimodular nodes partitions to be revealed, by letting links to be placed either 1) within the core part and between the core and the periphery parts or 2) just between the (empty) layers of a bipartite network. From a technical point of view, this is achieved by employing a multinomial hypergeometric distribution instead of the traditional (binomial) hypergeometric one; as in the latter case, this allows a p-value to be assigned to any given (bi)partition of the nodes. To illustrate the performance of our method, we report the results of its application to several real-world networks, including social, economic and financial ones.


title: A Stochastic Approximation Approach for Foresighted Task Scheduling in Cloud Computing

abstract: With the increasing and elastic demand for cloud resources, finding an optimal task scheduling mechanism become a challenge for cloud service providers. Due to the time-varying nature of resource demands in length and processing over time and dynamics and heterogeneity of cloud resources, existing myopic task scheduling solutions intended to maximize the performance of task scheduling are inefficient and sacrifice the long-time system performance in terms of resource utilization and response time. In this paper, we propose an optimal solution for performing foresighted task scheduling in a cloud environment. Since a-priori knowledge from the dynamics in queue length of virtual machines is not known in run time, an online reinforcement learning approach is proposed for foresighted task allocation. The evaluation results show that our method not only reduce the response time and makespan of submitted tasks, but also increase the resource efficiency. So in this thesis a scheduling method based on reinforcement learning is proposed. Adopting with environment conditions and responding to unsteady requests, reinforcement learning can cause a long-term increase in system's performance. The results show that this proposed method can not only reduce the response time and makespan but also increase resource efficiency as a minor goal.


title: Fully Supervised Speaker Diarization

abstract: In this paper, we propose a fully supervised speaker diarization approach, named unbounded interleaved-state recurrent neural networks (UIS-RNN). Given extracted speaker-discriminative embeddings (a.k.a. d-vectors) from input utterances, each individual speaker is modeled by a parameter-sharing RNN, while the RNN states for different speakers interleave in the time domain. This RNN is naturally integrated with a distance-dependent Chinese restaurant process (ddCRP) to accommodate an unknown number of speakers. Our system is fully supervised and is able to learn from examples where time-stamped speaker labels are annotated. We achieved a 7.6% diarization error rate on NIST SRE 2000 CALLHOME, which is better than the state-of-the-art method using spectral clustering. Moreover, our method decodes in an online fashion while most state-of-the-art systems rely on offline clustering.


title: Toward Resolved Simulations of Burning Fronts in Thermonuclear X-ray Bursts

abstract: We discuss the challenges of modeling X-ray bursts in multi-dimensions, review the different calculations done to date, and discuss our new set of ongoing simulations. We also describe algorithmic improvements that may help in the future to offset some of the expense of these simulations, and describe what may be possible with exascale computing.


title: Measurements of jets in heavy ion collisions

abstract: The Quark-Gluon Plasma (QGP) is created in high energy heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). This medium is transparent to electromagnetic probes but nearly opaque to colored probes. Hard partons produced early in the collision fragment and hadronize into a collimated spray of particles called a jet. The partons lose energy as they traverse the medium, a process called jet quenching. Most of the lost energy is still correlated with the parent parton, contributing to particle production at larger angles and lower momenta relative to the parent parton than in proton-proton collisions. This partonic energy loss can be measured through several observables, each of which give different insights into the degree and mechanism of energy loss. The measurements to date are summarized and the path forward is discussed.


title: A van Benthem Theorem for Quantitative Probabilistic Modal Logic

abstract: In probabilistic transition systems, behavioural metrics provide a more fine-grained and stable measure of system equivalence than crisp notions of bisimilarity. They correlate strongly to quantitative probabilistic logics, and in fact the distance induced by a probabilistic modal logic taking values in the real unit interval has been shown to coincide with behavioural distance. For probabilistic systems, probabilistic modal logic thus plays an analogous role to that of Hennessy-Milner logic on classical labelled transition systems. In the quantitative setting, invariance of modal logic under bisimilarity becomes non-expansivity of formula evaluation w.r.t. behavioural distance. In the present paper, we provide a characterization of the expressive power of probabilistic modal logic based on this observation: We prove a probabilistic analogue of the classical van Benthem theorem, which states that modal logic is precisely the bisimulation-invariant fragment of first-order logic. Specifically, we show that quantitative probabilistic modal logic lies dense in the bisimulation-invariant fragment, in the indicated sense of non-expansive formula evaluation, of quantitative probabilistic first-order logic; more precisely, bisimulation-invariant first-order formulas are approximable by modal formulas of bounded rank. For a description logic perspective on the same result, see arXiv:1906.00784.


title: Tight Dimension Independent Lower Bound on the Expected Convergence Rate for Diminishing Step Sizes in SGD

abstract: We study the convergence of Stochastic Gradient Descent (SGD) for strongly convex objective functions. We prove for all $t$ a lower bound on the expected convergence rate after the $t$-th SGD iteration; the lower bound is over all possible sequences of diminishing step sizes. It implies that recently proposed sequences of step sizes at ICML 2018 and ICML 2019 are {\em universally} close to optimal in that the expected convergence rate after {\em each} iteration is within a factor $32$ of our lower bound. This factor is independent of dimension $d$. We offer a framework for comparing with lower bounds in state-of-the-art literature and when applied to SGD for strongly convex objective functions our lower bound is a significant factor $775\cdot d$ larger compared to existing work.


title: Search for a Stable Six-Quark State at BABAR

abstract: Recent investigations have suggested that the six-quark combination uuddss could be a deeply bound state (S) that has eluded detection so far, and a potential dark matter candidate. We report the first search for a stable, doubly strange six-quark state in Upsilon -> S anti-Lambda anti-Lambda decays based on a sample of 90 million Upsilon(2S) and 110 million Upsilon(3S) decays collected by the BABAR experiment. No signal is observed, and 90% confidence level limits on the combined Upsilon(2S,3S) -> S anti-Lambda anti-Lambda branching fraction in the range (1.2-1.4)x10^-7 are derived for m_S < 2.05 GeV. These bounds set stringent limits on the existence of such exotic particles.


title: Inference for Volatility Functionals of Multivariate It\^o Semimartingales Observed with Jump and Noise

abstract: This paper presents the nonparametric inference for nonlinear volatility functionals of general multivariate It\^o semimartingales, in high-frequency and noisy setting. Pre-averaging and truncation enable simultaneous handling of noise and jumps. Second-order expansion reveals explicit biases and a pathway to bias correction. Estimators based on this framework achieve the optimal convergence rate. A class of stable central limit theorems are attained with estimable asymptotic covariance matrices. This paper form a basis for infill asymptotic results of, for example, the realized Laplace transform, the realized principal component analysis, the continuous-time linear regression, and the generalized method of integrated moments, hence helps to extend the application scopes to more frequently sampled noisy data.


title: Alternative stable states in a model of microbial community limited by multiple essential nutrients

abstract: Microbial communities routinely have several alternative stable states observed for the same environmental parameters. Sudden and irreversible transitions between these states make external manipulation of these systems more complicated. To better understand the mechanisms and origins of multistability in microbial communities, we introduce and study a model of a microbial ecosystem colonized by multiple specialist species selected from a fixed pool. Growth of each species can be limited by essential nutrients of two types, e.g. carbon and nitrogen, each represented in the environment by multiple metabolites. We demonstrate that our model has an exponentially large number of potential stable states realized for different environmental parameters. Using game theoretical methods adapted from the stable marriage problem we predict all of these states based only on ranked lists of competitive abilities of species for each of the nutrients. We show that for every set of nutrient influxes, several mutually uninvadable stable states are generally feasible and we distinguish them based upon their dynamic stability. We further explore an intricate network of discontinuous transitions (regime shifts) between these alternative states both in the course of community assembly, or upon changes of nutrient influxes.


title: The genomic Schur function is fundamental-positive

abstract: In work with A. Yong, the author introduced genomic tableaux to prove the first positive combinatorial rule for the Littlewood-Richardson coefficients in torus-equivariant $K$-theory of Grassmannians. We then studied the genomic Schur function $U_\lambda$, a generating function for such tableaux, showing that it is non-trivially a symmetric function, although generally not Schur-positive. Here we show that $U_\lambda$ is, however, positive in the basis of fundamental quasisymmetric functions. We give a positive combinatorial formula for this expansion in terms of gapless increasing tableaux; this is, moreover, the first finite expression for $U_\lambda$. Combined with work of A. Garsia and J. Remmel, this yields a compact combinatorial (but necessarily non-positive) formula for the Schur expansion of $U_\lambda$.


title: Systematic Quantum Cluster Typical Medium Method For the Study of Localization in Strongly Disordered Electronic Systems

abstract: Great progress has been made in the last several years towards understanding the properties of disordered electronic systems. In part, this is made possible by recent advances in quantum effective medium methods which enable the study of disorder and electron-electronic interactions on equal footing. They include dynamical mean field theory and the coherent potential approximation, and their cluster extension, the dynamical cluster approximation. Despite their successes, these methods do not enable the first-principles study of the strongly disordered regime, including the effects of electronic localization. The main focus of this review is the recently developed typical medium dynamical cluster approximation for disordered electronic systems. This method has been constructed to capture disorder-induced localization, and is based on a mapping of a lattice onto a quantum cluster embedded in an effective typical medium, which is determined self-consistently. Here we provide an overview of various recent applications of the typical medium dynamical cluster approximation to a variety of models and systems, including single and multi-band Anderson model, and models with local and off-diagonal disorder. We then present the application of the method to realistic systems in the framework of the density functional theory.


title: Buckling of geometrically confined shells

abstract: We study the periodic buckling patterns that emerge when elastic shells are subjected to geometric confinement. Residual swelling provides access to range of shapes (saddles, rolled sheets, cylinders, and spherical sections) which vary in their extrinsic and intrinsic curvatures. Our experimental and numerical data show that when these structures are radially confined, a single geometric parameter -- the ratio of the total shell radius to the amount of unconstrained material -- predicts the number of lobes formed. We then generalize our model to account for a relaxation of the in-plane constraint by interpreting this parameter as an effective foundation stiffness. Experimentally, we show that reducing the transverse confinement of saddles causes the lobe number to decrease according to our effective stiffness model. Hence, one geometric parameter captures the wave number through a wide range of radial and transverse confinement, connecting the shell shape to the shape of the boundary that confines it. We expect these results to be relevant for an expanse of shell shapes, and thus apply to the design of shape--shifting materials and the swelling and growth of soft structures.


title: Leveraging local network communities to predict academic performance

abstract: For more than 20 years, social network analysis of student collaboration networks has focused on a student's centrality to predict academic performance. And even though a growing amount of sociological literature has supported that academic success is contagious, identifying central students in the network alone does not capture how peer interactions facilitate the spread of academic success throughout the network. Consequently, we propose novel predictors that treat academic success as a contagion by identifying a student's learning community, consisting of the peers that are most likely to influence a student's performance in a course. We evaluate the importance of these learning communities by predicting academic outcomes in an introductory college statistics course with 103 students. In particular, we observe that by including these learning community predictors, the resulting model is 68 times more likely to be the correct model than the current state-of-the-art centrality network models in the literature.


title: Confirmation of the radial velocity super-Earth K2-18c with HARPS and CARMENES

abstract: In an earlier campaign to characterize the mass of the transiting temperate super-Earth K2-18b with HARPS, a second, non-transiting planet was posited to exist in the system at $\sim 9$ days. Further radial velocity follow-up with the CARMENES spectrograph visible channel revealed a much weaker signal at 9 days which also appeared to vary chromatically and temporally leading to the conclusion that the origin of the 9 day signal was more likely to be related to stellar activity than to being planetary. Here we conduct a detailed re-analysis of all available RV time-series, including a set of 31 previously unpublished HARPS measurements, to investigate the effects of time-sampling and of simultaneous modelling of planetary + activity signals on the existence and origin of the curious 9 day signal. We conclude that the 9 day signal is real and was initially seen to be suppressed in the CARMENES data due to a small number of anomalous measurements, although the exact cause of these anomalies remains unknown. Investigation of the signal's evolution in time, with wavelength, and detailed model comparison reveals that the 9 day signal is most likely planetary in nature. By this analysis, we reconcile the conflicting HARPS and CARMENES results and measure precise and self-consistent planet masses of $m_{p,b} = 8.63 \pm 1.35$ and $m_{p,c}\sin{i_c}=5.62 \pm 0.84$ M$_{\oplus}$. This work, along with the previously published RV papers on the K2-18 planetary system, highlight the importance of understanding one's time-sampling and of simultaneous planet + stochastic activity modelling, particularly when searching for sub-Neptune-sized planets with radial velocities.


title: Averages and higher moments for the $\ell$-torsion in class groups

abstract: We prove upper bounds for the average size of the $\ell$-torsion $\text{Cl}_K[\ell]$ of the class group of $K$, as $K$ runs through certain natural families of number fields and $\ell$ is a positive integer. We refine a key argument, used in almost all results of this type, which links upper bounds for $\text{Cl}_K[\ell]$ to the existence of many primes splitting completely in $K$ that are small compared to the discriminant of $K$. Our improvements are achieved through the introduction of a new family of specialised invariants of number fields to replace the discriminant in this argument, in conjunction with new counting results for these invariants. This leads to significantly improved upper bounds for the average and sometimes even higher moments of $\text{Cl}_K[\ell]$ for many families of number fields $K$ considered in the literature, for example, for the families of all degree-$d$-fields for $d\in\{2,3,4,5\}$ (and non-$D_4$ if $d=4$). As an application of the case $d=2$ we obtain the best upper bounds for the number of $D_p$-fields of bounded discriminant, for primes $p>3$.


title: The fluid mechanics of bubbly drinks

abstract: Bubbly drinks are surprisingly attractive. There is something about the nature of the these beverages that make them preferable among other choices. In this article we explore the physics involved in this particular kind of two-phase, mass-transfer-driven flows.


title: Model Selection of Nested and Non-Nested Item Response Models using Vuong Tests

abstract: In this paper, we apply Vuong's (1989) general approach of model selection to the comparison of nested and non-nested unidimensional and multidimensional item response theory (IRT) models. Vuong's approach of model selection is useful because it allows for formal statistical tests of both nested and non-nested models. However, only the test of non-nested models has been applied in the context of IRT models to date. After summarizing the statistical theory underlying the tests, we investigate the performance of all three distinct Vuong tests in the context of IRT models using simulation studies and real data. In the non-nested case we observed that the tests can reliably distinguish between the graded response model and the generalized partial credit model. In the nested case, we observed that the tests typically perform as well as or sometimes better than the traditional likelihood ratio test. Based on these results, we argue that Vuong's approach provides a useful set of tools for researchers and practitioners to effectively compare competing nested and non-nested IRT models.


title: Towards the Targeted Environment-Specific Evolution of Robot Components

abstract: This research considers the task of evolving the physical structure of a robot to enhance its performance in various environments, which is a significant problem in the field of Evolutionary Robotics. Inspired by the fields of evolutionary art and sculpture, we evolve only targeted parts of a robot, which simplifies the optimisation problem compared to traditional approaches that must simultaneously evolve both (actuated) body and brain. Exploration fidelity is emphasised in areas of the robot most likely to benefit from shape optimisation, whilst exploiting existing robot structure and control. Our approach uses a Genetic Algorithm to optimise collections of Bezier splines that together define the shape of a legged robot's tibia, and leg performance is evaluated in parallel in a high-fidelity simulator. The leg is represented in the simulator as 3D-printable file, and as such can be readily instantiated in reality. Provisional experiments in three distinct environments show the evolution of environment-specific leg structures that are both high-performing and notably different to those evolved in the other environments. This proof-of-concept represents an important step towards the environment-dependent optimisation of performance-critical components for a range of ubiquitous, standard, and already-capable robots that can carry out a wide variety of tasks.


title: Strangeness Enhancement due to String Fluctuations

abstract: We study string fragmentation in high multiplicity proton-proton collisions in a model where the string tension fluctuates. These fluctuations produce exponential pion spectra which are fitted to the transverse momentum distributions of charged particles for different multiplicities. For each multiplicity the so obtained hadronic slope parameter defines the magnitude of the string fluctuations which in turn determines the produced ratio of strange to light quarks. PYTHIA string decay simulations are used to convert each ratio of strange to light quarks to the appropriate ratio of strange hadrons to pions.


title: Compactification of the moduli space of instanton sheaves on the Fano threefold $V_5$

abstract: We study semistable sheaves of rank $2$ with Chern classes $c_1=0$, $c_2=2$ and $c_3=0$ on the Fano 3-fold $V_5$ of Picard number $1$, degree $5$ and index $2$. We show that the moduli space of such sheaves has a component that is isomorphic to $\mathbb{P}^5$ by identifying it with the moduli space of semistable quiver representations. This provides a natural smooth compactification of the moduli space of minimal instantons, as well as Ulrich bundles of rank $2$ on $V_5$.


title: Probabilistic Clustering Using Maximal Matrix Norm Couplings

abstract: In this paper, we present a local information theoretic approach to explicitly learn probabilistic clustering of a discrete random variable. Our formulation yields a convex maximization problem for which it is NP-hard to find the global optimum. In order to algorithmically solve this optimization problem, we propose two relaxations that are solved via gradient ascent and alternating maximization. Experiments on the MSR Sentence Completion Challenge, MovieLens 100K, and Reuters21578 datasets demonstrate that our approach is competitive with existing techniques and worthy of further investigation.


title: Compactification of the moduli space of minimal instantons on the Fano threefold $V_4$

abstract: We study semi-stable sheaves of rank $2$ with Chern class $c_1=0$, $c_2=2$ and $c_3=0$ on the Fano 3-folds $V_4$ of Picard number $1$, degree $4$ and index $2$. We show the moduli space of such sheaves is isomorphic to the moduli space of semi-stable rank $2$ even degree vector bundles on a genus $2$ curve. This provides a natural smooth compatification of the moduli space of Ulrich bundles of rank $2$ on $V_4$.


title: T-dual solutions of the Hull-Strominger system on non-K\"ahler threefolds

abstract: We construct new examples of solutions of the Hull-Strominger system on non-K\"ahler torus bundles over K3 surfaces, with the property that the connection $\nabla$ on the tangent bundle is Hermite-Yang-Mills. With this ansatz for the connection $\nabla$, we show that the existence of solutions reduces to known results about moduli spaces of slope-stable sheaves on a K3 surface, combined with elementary analytical methods. We apply our construction to find the first examples of T-dual solutions of the Hull-Strominger system on compact non-K\"ahler manifolds with different topology.


title: Mirach's Goblin: Discovery of a dwarf spheroidal galaxy behind the Andromeda galaxy

abstract: It is of broad interest for galaxy formation theory to carry out a full inventory of the numbers and properties of dwarf galaxies in the Local Volume, both satellites and isolated ones. Ultra-deep imaging in wide areas of the sky with small amateur telescopes can help to complete the census of these hitherto unknown low surface brightness galaxies, which cannot be detected by the current resolved stellar population and HI surveys. We report the discovery of Donatiello I, a dwarf spheroidal galaxy located one degree from the star Mirach (Beta And) in a deep image taken with an amateur telescope. The color--magnitude diagram obtained from follow-up observations obtained with the Gran Telescopio Canarias (La Palma, Spain) reveals that this system is beyond the Local Group and is mainly composed of old stars. The absence of young stars and HI emission in the ALFALFA survey are typical of quenched dwarf galaxies. Our photometry suggests a distance modulus for this galaxy of (m-M)=27.6 +/- 0.2 (3.3 Mpc), although this distance cannot yet be established securely owing to the crowding effects in our color--magnitude diagram. At this distance, Donatiello I's absolute magnitude (M_v =-8.3), surface brightness (mu_v=26.5 mag arcsec^-2) and stellar content are similar to the "classical" Milky Way companions Draco or Ursa Minor. The projected position and distance of Donatiello I are consistent with being a dwarf satellite of the closest S0-type galaxy NGC 404 ("Mirach's Ghost"). Alternatively, it could be one of the most isolated quenched dwarf galaxies reported so far behind the Andromeda galaxy.


title: Jordan totient quotients

abstract: The Jordan totient $J_k(n)$ can be defined by $J_k(n)=n^k\prod_{p\mid n}(1-p^{-k})$. In this paper, we study the average behavior of fractions $P/Q$ of two products $P$ and $Q$ of Jordan totients, which we call Jordan totient quotients. To this end, we describe two general and ready-to-use methods that allow one to deal with a larger class of totient functions. The first one is elementary and the second one uses an advanced method due to Balakrishnan and P\'etermann. As an application, we determine the average behavior of the Jordan totient quotient, the $k^{th}$ normalized derivative of the $n^{th}$ cyclotomic polynomial $\Phi_n(z)$ at $z=1$, the second normalized derivative of the $n^{th}$ cyclotomic polynomial $\Phi_n(z)$ at $z=-1$, and the average order of the Schwarzian derivative of $\Phi_n(z)$ at $z=1$.


title: Affine-null metric formulation of General Relativity at two intersecting null hypersurfaces

abstract: We revisit Winicour's affine-null metric initial value formulation of General Relativity, where the characteristic initial value formulation is set up with a null metric having two affine parameters. In comparison to past work, where the application of the formulation was aimed for the timelike-null initial value problem, we consider here a boundary surface that is a null hypersurface. All of the initial data are either metric functions or first derivatives of the metric. Given such a set of initial data, Einstein equations can be integrated in a hierarchical manner, where first a set of equations is solved hierarchically on the null hypersurface serving as a boundary. Second, with the obtained boundary values, a set of differential equations, similar to the equations of the Bondi-Sachs formalism, comprising of hypersurface and evolution equations is solved hierarchically to find the entire space-time metric. An example is shown how the double null Israel black hole solution arises after specification to spherical symmetry and vacuum. This black hole solution is then discussed to with respect to Penrose conformal compactification of spacetime.


title: Generalized Ziggurat Algorithm for Unimodal and Unbounded Probability Density Functions with Zest

abstract: We present a modified Ziggurat algorithm that could generate a random number from all unimodal and unbounded PDFs. For PDFs that have unbounded density and/or unbounded support we use a combination of nonlinear mapping function and rejection sampling to generate a random number from the peak and/or the tail distribution. A family of mapping functions and their corresponding acceptance probability functions are presented (along with the criteria for their use and their efficiency) that could be used to generate random numbers from infinite tails and unbounded densities. The Zest library which is a C++ implementation of this algorithm is also presented. Zest can efficiently generate normal, exponential, cauchy, gamma, Weibull, log-normal, chi-squared, student's t and Fisher's f variates. The user can also define their custom PDF as a class and supply it as a template argument to our library's class without modifying any part of the library. Performance of Zest is compared against performance of random modules of (GCC's implementation of) Standard Template Library (STL) and Boost. The presented results show that Zest is faster than both in most cases, sometimes by a factor of more than 10. We also present a C++ implementation of a uniform floating-point random number generator (RNG) which is capable of producing all representable floating-point numbers in $[0,1)$ which will be used in the Ziggurat algorithm near unbounded peaks and tails. The common method of dividing a random integer by the range of the RNG can not produce random floating-point numbers with fully random fraction bits and very small random numbers. The presented uniform floating-point RNG is very efficient and in the case of producing double precision floating-point numbers it's even faster than simply multiplying a 64-bit integer by $2^{-64}$.


title: Thickness-Dependent Differential Reflectance Spectra of Monolayer and Few-Layer MoS2, MoSe2, WS2 and WSe2

abstract: The research field of two dimensional (2D) materials strongly relies on optical microscopy characterization tools to identify atomically thin materials and to determine their number of layers. Moreover, optical microscopy-based techniques opened the door to study the optical properties of these nanomaterials. We presented a comprehensive study of the differential reflectance spectra of 2D semiconducting transition metal dichalcogenides (TMDCs), MoS2, MoSe2, WS2, and WSe2, with thickness ranging from one layer up to six layers. We analyzed the thickness-dependent energy of the different excitonic features, indicating the change in the band structure of the different TMDC materials with the number of layers. Our work provided a route to employ differential reflectance spectroscopy for determining the number of layers of MoS2, MoSe2, WS2, and WSe2.


title: Stability and Erd\H{o}s--Stone type results for $F$-free graphs with a fixed number of edges

abstract: A fundamental problem of extremal graph theory is to ask, 'What is the maximum number of edges in an $F$-free graph on $n$ vertices?' Recently Alon and Shikhelman proposed a more general, subgraph counting, version of this question. They considered the question of determining the maximum number of copies of a fixed graph $T$ in an $F$-free graph on $n$ vertices. In this more general context, where we are no longer counting edges, it is also natural to ask what is the maximum number of copies of $T$ in an $F$-free graph with $m$ edges and no restriction on the number of vertices. Frohmader, in a different context, determined the answer when $T$ and $F$ are both complete graphs. We prove results for this problem analogous to the Erd\H{o}s--Stone theorem, the Erd\H{o}s--Simonovits theorem, and the stability theorem of Erd\H{o}s--Simonovits.


title: Flavored Dark Sectors running to low energy

abstract: We consider the effective field theory generated by a heavy mediator that connects Standard Model particles to a Dark Sector, considering explicitly the flavor structure of the operators. In particular, we study the model independent running and mixing between operators, as well as their matching at the electroweak scale. In addition to the explicit expression of the Renormalization Group Equations, we show the numerical solutions as well as some approximate analytical expressions that help understanding these solutions. At low energy, our results are particularly important in the case of light dark sectors communicating to the b quark, and can be immediately applied to flavored Dark Matter.


title: Empirical Bayes to assess ecological diversity and similarity with overdispersion in multivariate counts

abstract: The assessment of diversity and similarity is relevant in monitoring the status of ecosystems. The respective indicators are based on the taxonomic composition of biological communities of interest, currently estimated through the proportions computed from sampling multivariate counts. In this work we present a novel method able to work with only one sample to estimate the taxonomic composition when the data are affected by overdispersion. The presence of overdispersion in taxonomic counts may be the result of significant environmental factors which are often unobservable but influence communities. Following the empirical Bayes approach, we combine a Bayesian model with the marginal likelihood method to jointly estimate the taxonomic proportions and the level of overdispersion from one sample of multivariate counts. Our proposal is compared to the classical maximum likelihood method in an extensive simulation study with different realistic scenarios. An application to real data from aquatic biomonitoring is also presented. In both the simulation study and the real data application, we consider communities characterized by a large number of taxonomic categories, such as aquatic macroinvertebrates or bacteria which are often overdispersed. The applicative results demonstrate an overall superiority of the empirical Bayes method in almost all examined cases, for both assessments of diversity and similarity. We would recommend practitioners in biomonitoring to use the proposed approach in addition to the traditional procedures. The empirical Bayes estimation allows to better control the error propagation due to the presence of overdispersion in biological data, with a more efficient managerial decision making.


title: Exploring Potential Environmental Benefits of Asteroid Mining

abstract: Asteroid mining has been proposed as an approach to complement Earth-based supplies of rare earth metals and supplying resources in space, such as water. Existing research on asteroid mining has mainly looked into its economic viability, technological feasibility, cartography of asteroids, and legal aspects. More recently, potential environmental benefits for asteroid mining have been considered. However, no quantitative estimate of these benefits has been given. This paper attempts to determine if and under which conditions asteroid mining would have environmental benefits, compared to either Earth-based mining or launching equipment and resources into space. We focus on two cases: Water supply to cis-lunar orbit and platinum mining. First, we conduct a state-of-the-art of current environmental life cycle assessment for the space domain and platinum mining. Second, a first order environmental life cycle assessment is conducted, including goal and scope definition, inventory analysis, and impact assessment. We compare water supply to cis-lunar orbit with and without asteroid mining and go on to compare terrestrial with space-based platinum mining. The results indicate that asteroid water mining would have environmental benefits, as soon as the amount of water supplied via mining is larger than the mass of the spacecraft used for mining. For platinum mining, we find that by comparing the operations phase of terrestrial and space mining, space mining would have a lower environmental impact, if the spacecraft is able to return between 0.3 to 7% of its mass in platinum to Earth, assuming 100% primary platinum or 100% secondary platinum, respectively. For future work, we propose a more detailed analysis, based on a more precise inventory and a larger system boundary, including the production of the launcher and spacecraft.


title: Non-reciprocity and zero reflection in nonlinear cavities with tailored loss

abstract: We demonstrate how to tailor the losses of nonlinear cavities in order to suppress their reflection and enhance their non-reciprocal transmission. We derive analytical expressions predicting the existence of zero-reflection channels in single and coupled nonlinear cavities, depending on the driving frequency and loss rates. While suppressing the reflection from a single cavity imposes a stringent condition on the input-output leakage rates, we demonstrate that this condition can be significantly relaxed in systems of coupled cavities. In particular, zero-reflection and non-reciprocity can be achieved across a range of driving frequencies in coupled cavities by tuning the output leakage rate alone. Numerical calculations based on the driven-dissipative Gross-Pitaevksii equation, usually employed to describe microcavity polaritons, reveal the spatial phenomenology associated with zero-reflection states and provide design guidelines for the construction of nonlinear optical isolators.


title: Non-local interactions in the $(d,p)$ surrogate method for $(n,\gamma)$ reactions

abstract: Single-neutron transfer reactions populating states in the continuum are interesting both for structure and astrophysics. In their description often global optical potentials are used for the nucleon-target interactions, and these interactions are typically local. In our work, we study the effects of nonlocality in $(d,p)$ reactions populating continuum states. This work is similar to that of [1] but now for transfer to the continuum. A theory for computing cross sections for inclusive processes $A (d,p) X$ was explored in [2]. Therein, local optical potentials were used to describe the nucleon-target effective interaction. The goal of the present work is to extend the theory developed in [2] to investigate the effects of including nonlocality in the effective interaction on the relevant reaction observables. We implement the R--matrix method to solve the non--local equations both for the nucleon wavefunctions and the propagator. We then apply the method to systematically study the inclusive process of $(d,p)$ on $^{16}$O, $^{40}$Ca, $^{48}$Ca and $^{208}$Pb at 10, 20 and 50 MeV. We compare the results obtained when non--local interactions are used with those obtained when local equivalent interactions are included. We find that nonlocality affects different pieces of the model in complex ways. Depending on the beam energy and the target, the non-elastic breakup can either increase or decrease. While the non-elastic transfer cross section for each final spin state can change considerably, the main prediction of the model [2], namely the shape of the spin distributions, remains largely unaltered by nonlocality. [1] L. J. Titus and F. M. Nunes. Testing the Perey effect. Phys. Rev. C, 89:034609, Mar 2014. [2] G. Potel, F. M. Nunes, and I. J. Thompson. Establishing a theory for deuteron-induced surrogate reactions. Phys. Rev. C, 92:034611, Sep 2015.


title: Deep Recurrent Level Set for Segmenting Brain Tumors

abstract: Variational Level Set (VLS) has been a widely used method in medical segmentation. However, segmentation accuracy in the VLS method dramatically decreases when dealing with intervening factors such as lighting, shadows, colors, etc. Additionally, results are quite sensitive to initial settings and are highly dependent on the number of iterations. In order to address these limitations, the proposed method incorporates VLS into deep learning by defining a novel end-to-end trainable model called as Deep Recurrent Level Set (DRLS). The proposed DRLS consists of three layers, i.e, Convolutional layers, Deconvolutional layers with skip connections and LevelSet layers. Brain tumor segmentation is taken as an instant to illustrate the performance of the proposed DRLS. Convolutional layer learns visual representation of brain tumor at different scales. Since brain tumors occupy a small portion of the image, deconvolutional layers are designed with skip connections to obtain a high quality feature map. Level-Set Layer drives the contour towards the brain tumor. In each step, the Convolutional Layer is fed with the LevelSet map to obtain a brain tumor feature map. This in turn serves as input for the LevelSet layer in the next step. The experimental results have been obtained on BRATS2013, BRATS2015 and BRATS2017 datasets. The proposed DRLS model improves both computational time and segmentation accuracy when compared to the the classic VLS-based method. Additionally, a fully end-to-end system DRLS achieves state-of-the-art segmentation on brain tumors.


title: Energy returns in global AdS_4

abstract: Recent studies of the weakly nonlinear dynamics of probe fields in global AdS$_4$ (and of the nonrelativistic limit of AdS resulting in the Gross-Pitaevskii equation) have revealed a number of cases with exact dynamical returns for two-mode initial data. In this paper, we address the question whether similar exact returns are present in the weakly nonlinear dynamics of gravitationally backreacting perturbations in global AdS$_4$. In the literature, approximate returns were first pointed out numerically and with limited precision. We first provide a thorough numerical analysis and discover returns that are so accurate that it would be tantalizing to sign off the small imperfections as an artifact of numerics. To clarify the situation, we introduce a systematic analytic approach by focusing on solutions with spectra localized around one of the two lowest modes. This allows us to demonstrate that in the gravitational case the returns are not exact. Furthermore, our analysis predicts and explains specific integer numbers of direct-reverse cascade sequences that result in particularly accurate energy returns (elaborate hierarchies of more and less precise returns arise if one waits for appropriate longer multiple periods in this manner). In addition, we explain, at least in this regime, the ubiquitous appearance of direct-reverse cascades in the weakly nonlinear dynamics of AdS-like systems.


title: Efficient Tensor Decomposition with Boolean Factors

abstract: Tensor decomposition has been extensively used as a tool for exploratory analysis. Motivated by neuroscience applications, we study tensor decomposition with Boolean factors. The resulting optimization problem is challenging due to the non-convex objective and the combinatorial constraints. We propose Binary Matching Pursuit (BMP), a novel generalization of the matching pursuit strategy to decompose the tensor efficiently. BMP iteratively searches for atoms in a greedy fashion. The greedy atom search step is solved efficiently via a MAXCUT-like boolean quadratic program. We prove that BMP is guaranteed to converge sublinearly to the optimal solution and recover the factors under mild identifiability conditions. Experiments demonstrate the superior performance of our method over baselines on synthetic and real datasets. We also showcase the application of BMP in quantifying neural interactions underlying high-resolution spatiotemporal ECoG recordings.


title: Leveraging Textual Specifications for Grammar-based Fuzzing of Network Protocols

abstract: Grammar-based fuzzing is a technique used to find software vulnerabilities by injecting well-formed inputs generated following rules that encode application semantics. Most grammar-based fuzzers for network protocols rely on human experts to manually specify these rules. In this work we study automated learning of protocol rules from textual specifications (i.e. RFCs). We evaluate the automatically extracted protocol rules by applying them to a state-of-the-art fuzzer for transport protocols and show that it leads to a smaller number of test cases while finding the same attacks as the system that uses manually specified rules.


title: Stochastic Synthesis for Stochastic Computing

abstract: Stochastic computing (SC) is an emerging computing technique which offers higher computational density, and lower power over binary-encoded (BE) computation. Unlike BE computation, SC encodes values as probabilistic bitstreams which makes designing new circuits unintuitive. Existing techniques for synthesizing SC circuits are limited to specific classes of functions such as polynomial evaluation or constant scaling. In this paper, we propose using stochastic synthesis, which is originally a program synthesis technique, to automate the task of synthesizing new SC circuits. Our results show stochastic synthesis is more general than past techniques and can synthesize manually designed SC circuits as well as new ones such as an approximate square root unit.


title: Accretion in Common Envelope Evolution

abstract: Common envelope evolution (CEE) occurs in some binary systems involving asymptotic giant branch (AGB) or red giant branch (RGB) stars, and understanding this process is crucial for understanding the origins of various transient phenomena. CEE has been shown to be highly asymmetrical and global 3D simulations are needed to help understand the dynamics. We perform and analyze hydrodynamic CEE simulations with the adaptive mesh refinement (AMR) code AstroBEAR, and focus on the role of accretion onto the companion star. We bracket the range of accretion rates by comparing a model that removes mass and pressure using a subgrid accretion prescription with one that does not. Provided a pressure-release valve, such as a bipolar jet, is available, super-Eddington accretion could be common. Finally, we summarize new results pertaining to the energy budget, and discuss the overall implications relating to the feasibility of unbinding the envelope in CEE simulations.


title: Hybrid KNN-Join: Parallel Nearest Neighbor Searches Exploiting CPU and GPU Architectural Features

abstract: K Nearest Neighbor (KNN) joins are used in scientific domains for data analysis, and are building blocks of several well-known algorithms. KNN-joins find the KNN of all points in a dataset. This paper focuses on a hybrid CPU/GPU approach for low-dimensional KNN-joins, where the GPU may not yield substantial performance gains over parallel CPU algorithms. We utilize a work queue that prioritizes computing data points in high density regions on the GPU, and low density regions on the CPU, thereby taking advantage of each architecture's relative strengths. Our approach, HybridKNN-Join, effectively augments a state-of-the-art multi-core CPU algorithm. We propose optimizations that $(i)$ maximize GPU query throughput by assigning the GPU large batches of work; $(ii)$ increase workload granularity to optimize GPU utilization; and, $(iii)$ limit load imbalance between CPU and GPU architectures. We compare HybridKNN-Join to one GPU and two parallel CPU reference implementations. Compared to the reference implementations, we find that the hybrid algorithm performs best on larger workloads (dataset size and K). The methods employed in this paper show promise for the general division of work in other hybrid algorithms.


title: Integrating electricity markets: Impacts of increasing trade on prices and emissions in the western United States

abstract: This paper presents empirically-estimated average hourly relationships between regional electricity trade in the United States and prices, emissions, and generation from 2015 through 2018. Consistent with economic theory, the analysis finds a negative relationship between electricity prices in California and regional trade, conditional on local demand. Each 1 gigawatt-hour increase in California electricity imports is associated with an average $0.15 per megawatt-hour decrease in the California Independent System Operator's wholesale electricity price. There is a net-negative short term relationship between carbon dioxide emissions in California and electricity imports that is partially offset by positive emissions from exporting neighbors. Specifically, each 1 GWh increase in regional trade is associated with a net 70-ton average decrease in CO2 emissions across the western U.S., conditional on demand levels. The results provide evidence that electricity imports mostly displace natural gas generation on the margin in the California electricity market. A small positive relationship is observed between short-run SO2 and NOx emissions in neighboring regions and California electricity imports. The magnitude of the SO2 and NOx results suggest an average increase of 0.1 MWh from neighboring coal plants is associated with a 1 MWh increase in imports to California.


title: Towards Differentially Private Truth Discovery for Crowd Sensing Systems

abstract: Nowadays, crowd sensing becomes increasingly more popular due to the ubiquitous usage of mobile devices. However, the quality of such human-generated sensory data varies significantly among different users. To better utilize sensory data, the problem of truth discovery, whose goal is to estimate user quality and infer reliable aggregated results through quality-aware data aggregation, has emerged as a hot topic. Although the existing truth discovery approaches can provide reliable aggregated results, they fail to protect the private information of individual users. Moreover, crowd sensing systems typically involve a large number of participants, making encryption or secure multi-party computation based solutions difficult to deploy. To address these challenges, in this paper, we propose an efficient privacy-preserving truth discovery mechanism with theoretical guarantees of both utility and privacy. The key idea of the proposed mechanism is to perturb data from each user independently and then conduct weighted aggregation among users' perturbed data. The proposed approach is able to assign user weights based on information quality, and thus the aggregated results will not deviate much from the true results even when large noise is added. We adapt local differential privacy definition to this privacy-preserving task and demonstrate the proposed mechanism can satisfy local differential privacy while preserving high aggregation accuracy. We formally quantify utility and privacy trade-off and further verify the claim by experiments on both synthetic data and a real-world crowd sensing system.


title: O(D,D) gauge fields in the T-dual string Lagrangian

abstract: We present the string Lagrangian with manifest T-duality. Not only zero-modes but also all string modes are doubled. The gravitational field is an O(D,D) gauge field. We give a Lagrangian version of the section condition for the gauge invariance which compensates the O(D,D) transformation from the gravitational field and the GL(2D) coordinate transformation. We also show the gauge invariance of the line element of the manifest T-duality space and the O(D,D) condition on the background. Different sections describe dual spaces.


title: The Connection between Starspots and Flares on Main-Sequence Kepler Stars

abstract: Starspots and flares are indicators of stellar magnetic activity that can both be studied in greater detail by utilizing the long-term, space-based archive of the Kepler satellite. Here, we aim to investigate a subset of the Kepler archive to reveal a connection between the starspots and the stellar flares, in order to provide insight into the overall stellar magnetic field. We use the flare-finding algorithm FLATW'RM in conjunction with a new suite of algorithms that aim to locate the local minima caused by starspot groups. We compare the phase difference between the time of maximum flux of a flare and the time of minimum stellar flux due to a starspot group. The strongest flares do not appear to be correlated to the largest starspot group present, but are also not uniformly distributed in phase with respect to the starspot group. The weaker flares, however, do show an increased occurrence close to the starspot groups.


title: Recursion Schemes, the MSO Logic, and the U quantifier

abstract: We study the model-checking problem for recursion schemes: does the tree generated by a given higher-order recursion scheme satisfy a given logical sentence. The problem is known to be decidable for sentences of the MSO logic. We prove decidability for an extension of MSO in which we additionally have an unbounding quantifier U, saying that a subformula is true for arbitrarily large finite sets. This quantifier can be used only for subformulae in which all free variables represent finite sets (while an unrestricted use of the quantifier leads to undecidability). We also show that the logic has the properties of reflection and effective selection for trees generated by recursion schemes.


title: Support theorems for degenerate stochastic differential equations with jumps and applications

abstract: In the paper, we are concerned with degenerate stochastic differential equations with jumps. Firstly, we establish two support theorems for the solutions of the degenerate stochastic equations, under different (sufficient) conditions. Secondly, we apply one of our support theorems to a class of degenerate stochastic evolution equations (i.e., infinite-dimensional stochastic differential equations) with jumps to get a characterisation of path-independence for the densities of their Girsanov transformations.


title: Frank-Wolfe Method is Automatically Adaptive to Error Bound Condition

abstract: Error bound condition has recently gained revived interest in optimization. It has been leveraged to derive faster convergence for many popular algorithms, including subgradient methods, proximal gradient method and accelerated proximal gradient method. However, it is still unclear whether the Frank-Wolfe (FW) method can enjoy faster convergence under error bound condition. In this short note, we give an affirmative answer to this question. We show that the FW method (with a line search for the step size) for optimization over a strongly convex set is automatically adaptive to the error bound condition of the problem. In particular, the iteration complexity of FW can be characterized by $O(\max(1/\epsilon^{1-\theta}, \log(1/\epsilon)))$ where $\theta\in[0,1]$ is a constant that characterizes the error bound condition. Our results imply that if the constrained set is characterized by a strongly convex function and the objective function can achieve a smaller value outside the considered domain, then the FW method enjoys a fast rate of $O(1/t^2)$.


title: An edge-based pressure stabilisation technique for finite elements on arbitrarily anisotropic meshes

abstract: In this article, we analyse a stabilised equal-order finite element approximation for the Stokes equations on anisotropic meshes. In particular, we allow arbitrary anisotropies in a sub-domain, for example along the boundary of the domain, with the only condition that a maximum angle is fulfilled in each element.This discretisation is motivated by applications on moving domains as arising e.g. in fluid-structure interaction or multiphase-flow problems. To deal with the anisotropies, we define a modification of the original Continuous Interior Penalty stabilisation approach. We show analytically the discrete stability of the method and convergence of order ${\cal O}(h^{3/2})$ in the energy norm and ${\cal O}(h^{5/2})$ in the $L^2$-norm of the velocities. We present numerical examples for a linear Stokes problem and for a non-linear fluid-structure interaction problem, that substantiate the analytical results and show the capabilities of the approach.


title: The RHIC Beam Energy Scan Phase II: Physics and Upgrades

abstract: The exploration of the QCD phase diagram has been one of the main drivers of contemporary nuclear physics. The Relativistic Heavy Ion Collider (RHIC) at BNL is uniquely suited for this task through its Beam Energy Scan (BES) program which allowed for a large range in baryon chemical potential $\mu_B$ as was successfully demonstrated after the completion of Phase 1 in 2014. Phase 2 of the BES at RHIC is scheduled to start in 2019 and will explore with precision measurements the intermediate-to-high $\mu_B$ region of the QCD phase diagram, five energies $\sqrt{s_{NN}}$ from 7.7 to 19.6 GeV in collider mode and eight energies $\sqrt{s{_{NN}}}$ from 3.0 to 7.7 GeV in fixed-target mode. Some of the key measurements are: the net-protons kurtosis that could pinpoint the position of a critical point, the directed flow that might prove a softening of the EOS, and the chiral restoration in the dielectron channel. These measurements will be possible with an order of magnitude better statistics provided by the electron cooling upgrade of RHIC and with the detector upgrades planned to improve STAR's acceptance. These proceedings review the BES Phase-2 program and the physics opportunities enabled by these upgrades.


title: Reconstruction of Function Fields from their pro-l abelian divisorial Inertia

abstract: Let $\Pi^c_K\to\Pi_K$ be the maximal pro-$\ell$ abelian-by-central, respectively abelian, Galois groups of a function field $K|k$ with $k$ algebraically closed and ${\rm char}\neq\ell$. We show that $K|k$ can be functorially reconstructed by group theoretical recipes from $\Pi^c_K$ endowed with the set of divisorial inertia ${\rm Inrdiv}(K)\subset\Pi_K$. As applications, one has: (i) A group theoretical recipe to reconstruct $K|k$ from $\Pi^c_K$, provided either ${\rm Tr.deg}(K|k)>{\rm dim}(k)+1$ or ${\rm tr.deg}(K|k) >{\rm dim}(k)>1$, where ${\rm dim}(k)$ is the Kronecker dimension; (ii) An application to the pro-$\ell\!$ abelian-by-central I/OM (Ihara's question / Oda-Matsumoto conjecture), which in the cases considered here implies the classical I/OM.


title: Localization in Khovanov homology

abstract: We construct equivariant Khovanov spectra for periodic links, using the Burnside functor construction introduced by Lawson, Lipshitz, and Sarkar. By identifying the fixed-point sets, we obtain rank inequalities for odd and even Khovanov homologies, and their annular filtrations, for prime-periodic links in $S^3$.


title: Smoothly embedding Seifert fibered spaces in $S^4$

abstract: Using an obstruction based on Donaldson's theorem, we derive strong restrictions on when a Seifert fibered space $Y = F(e; \frac{p_1}{q_1}, \ldots, \frac{p_k}{q_k})$ over an orientable base surface $F$ can smoothly embed in $S^4$. This allows us to classify precisely when $Y$ smoothly embeds provided $e > k/2$, where $e$ is the normalized central weight and $k$ is the number of singular fibers. Based on these results and an analysis of the Neumann-Siebenmann invariant $\overline{\mu}$, we make some conjectures concerning Seifert fibered spaces which embed in $S^4$. Finally, we also provide some applications to doubly slice Montesinos links, including a classification of the smoothly doubly slice odd pretzel knots up to mutation.


title: The mod-$p$ homology of the classifying spaces of certain gauge groups

abstract: Let $G$ be a simply-connected, simple compact Lie group of type $\{n_{1},\ldots,n_{\ell}\}$, where $n_{1}\le\cdots \le n_{\ell}$. Let $\mathcal{G}_k$ be the gauge group of the principal $G$-bundle (namedright{P}{}{S^{4}}) whose isomorphism class is determined by the the second Chern class having value $k\in\mathbb{Z}$. We calculate the mod-$p$ homology of the classifying space $B\mathcal{G}_k$ provided that $n_{\ell}<p-1$ and $p$ does not divide $k$.


title: On the cover time of dense graphs

abstract: We consider arbitrary graphs $G$ with $n$ vertices and minimum degree at least $\delta n$ where $\delta>0$ is constant. If the conductance of $G$ is sufficiently large then we obtain an asymptotic expression for the cover time $C_G$ of $G$ as the solution to an explicit transcendental equation. Failing this, if the mixing time of a random walk on $G$ is of a lesser magnitude than the cover time, then we can obtain an asymptotic deterministic estimate via a decomposition into a bounded number of dense sub-graphs with high conductance. Failing this we give a deterministic asymptotic (2+o(1))-approximation of $C_G$.


title: The IFF Foundation for Ontological Knowledge Organization

abstract: This paper discusses an axiomatic approach for the integration of ontologies, an approach that extends to first order logic a previous approach (Kent 2000) based on information flow. This axiomatic approach is represented in the Information Flow Framework (IFF), a metalevel framework for organizing the information that appears in digital libraries, distributed databases and ontologies (Kent 2001). The paper argues that the integration of ontologies is the two-step process of alignment and unification. Ontological alignment consists of the sharing of common terminology and semantics through a mediating ontology. Ontological unification, concentrated in a virtual ontology of community connections, is fusion of the alignment diagram of participant community ontologies - the quotient of the sum of the participant portals modulo the ontological alignment structure.


title: Distributed Conceptual Structures

abstract: The theory of distributed conceptual structures, as outlined in this paper, is concerned with the distribution and conception of knowledge. It rests upon two related theories, Information Flow and Formal Concept Analysis, which it seeks to unify. Information Flow (IF) is concerned with the distribution of knowledge. The foundations of Information Flow is explicitly based upon a mathematical theory known as the Chu Construction in *-autonomous categories and implicitly based upon the mathematics of closed categories. Formal Concept Analysis (FCA) is concerned with the conception and analysis of knowledge. In this paper we connect these two studies by extending the basic theorem of Formal Concept Analysis to the distributed realm of Information Flow. The main results are the categorical equivalence between classifications and concept lattices at the level of functions, and the categorical equivalence between bonds and complete adjoints at the level of relations. With this we hope to accomplish a rapprochement between Information Flow and Formal Concept Analysis.


title: Straintronics beyond homogeneous deformation

abstract: We present a continuum theory of graphene treating on an equal footing both homogeneous Cauchy-Born (CB) deformation, as well as the microscopic degrees of freedom associated with the two sublattices. While our theory recovers all extant results from homogeneous continuum theory, the Dirac-Weyl equation is found to be augmented by new pseudo-gauge and chiral fields fundamentally different from those that result from homogeneous deformation. We elucidate three striking electronic consequences: (i) non-CB deformations allow for the transport of valley polarized charge over arbitrarily long distances e.g. along a designed ridge; (ii) the triaxial deformations required to generate an approximately uniform magnetic field are unnecessary with non-CB deformation; and finally (iii) the vanishing of the effects of a one dimensional corrugation seen in \emph{ab-initio} calculation upon lattice relaxation are explained as a compensation of CB and non-CB deformation.


title: Probabilistic Safety Analysis using Traffic Microscopic Simulation

abstract: Traffic microscopic simulation applications are a common tool in road transportation analysis and several attempts to perform road safety assessments have recently been carried out. However, these approaches often ignore causal relationships between different levels of vehicle interactions and/or accident types and they lack a physical representation of the accident phenomena itself. In this paper, a new generic probabilistic safety assessment framework for traffic microscopic simulation tools is proposed. The probability of a specific accident occurring is estimated by an accident propensity function that consists of a deterministic safety score component and a random component. The formulation of the safety score depends on the type of occurrence, on detailed vehicle interactions and maneuvers and on its representation in a simulation environment. This generic model is applied to the case of an urban motorway and specified to four types of outcomes: non-accident events and three types of accidents in a nested structure: rear-end, lane-changing, and run-off-road accidents. The model was estimated and validated using simulated microscopic data. To obtained the consistent simulated data, a two-step simulation calibration procedure was adopted: (1) using real trajectories collected on site for detailed behavior representation; and (2) using aggregate data from each event used in safety model estimation. The final estimated safety model is able to identify and interpret several simulated vehicle interactions. The fact that these outcomes were extracted from simulated analysis shows the real potential of calibrated traffic microscopic simulation for detailed safety assessments.


title: Rao-Blackwellized Stochastic Gradients for Discrete Distributions

abstract: We wish to compute the gradient of an expectation over a finite or countably infinite sample space having $K \leq \infty$ categories. When $K$ is indeed infinite, or finite but very large, the relevant summation is intractable. Accordingly, various stochastic gradient estimators have been proposed. In this paper, we describe a technique that can be applied to reduce the variance of any such estimator, without changing its bias---in particular, unbiasedness is retained. We show that our technique is an instance of Rao-Blackwellization, and we demonstrate the improvement it yields on a semi-supervised classification problem and a pixel attention task.


title: Offline Multi-Action Policy Learning: Generalization and Optimization

abstract: In many settings, a decision-maker wishes to learn a rule, or policy, that maps from observable characteristics of an individual to an action. Examples include selecting offers, prices, advertisements, or emails to send to consumers, as well as the problem of determining which medication to prescribe to a patient. While there is a growing body of literature devoted to this problem, most existing results are focused on the case where data comes from a randomized experiment, and further, there are only two possible actions, such as giving a drug to a patient or not. In this paper, we study the offline multi-action policy learning problem with observational data and where the policy may need to respect budget constraints or belong to a restricted policy class such as decision trees. We build on the theory of efficient semi-parametric inference in order to propose and implement a policy learning algorithm that achieves asymptotically minimax-optimal regret. To the best of our knowledge, this is the first result of this type in the multi-action setup, and it provides a substantial performance improvement over the existing learning algorithms. We then consider additional computational challenges that arise in implementing our method for the case where the policy is restricted to take the form of a decision tree. We propose two different approaches, one using a mixed integer program formulation and the other using a tree-search based algorithm.


title: Redirect2Own: Protecting the Intellectual Property of User-uploaded Content through Off-site Indirect Access

abstract: Social networking services have attracted millions of users, including individuals, professionals, and companies, that upload massive amounts of content, such as text, pictures, and video, every day. Content creators retain the intellectual property (IP) rights on the content they share with these networks, however, very frequently they implicitly grant them, a sometimes, overly broad license to use that content, which enables the services to use it in possibly undesirable ways. For instance, Facebook claims a transferable, sub-licensable, royalty-free, worldwide license on all user-provided content. Professional content creators, like photographers, are particularly affected. In this paper we propose a design for decoupling user data from social networking services without any loss of functionality for the users. Our design suggests that user data are kept off the social networking service, in third parties that enable the hosting of user-generated content under terms of service and overall environment (e.g., a different location) that better suit the user's needs and wishes. At the same time, indirection schemata are seamlessly integrated in the social networking service, without any cooperation from the server side necessary, so that users can transparently access the off-site data just as they would if hosted in-site. We have implemented our design as an extension for the Chrome Web browser, called Redirect2Own, and show that it incurs negligible overhead on accessing 'redirected' content. We offer the extension as free software and its code as an open-source project.


title: Common pulse retrieval algorithm: a fast and universal method to retrieve ultrashort pulses

abstract: We present a common pulse retrieval algorithm (COPRA) that can be used for a broad category of ultrashort laser pulse measurement schemes including frequency-resolved optical gating (FROG), interferometric FROG, dispersion scan, time domain ptychography, and pulse shaper assisted techniques such as multiphoton intrapulse interference phase scan (MIIPS). We demonstrate its properties in comprehensive numerical tests and show that it is fast, reliable and accurate in the presence of Gaussian noise. For FROG it outperforms retrieval algorithms based on generalized projections and ptychography. Furthermore, we discuss the pulse retrieval problem as a nonlinear least-squares problem and demonstrate the importance of obtaining a least-squares solution for noisy data. These results improve and extend the possibilities of numerical pulse retrieval. COPRA is faster and provides more accurate results in comparison to existing retrieval algorithms. Furthermore, it enables full pulse retrieval from measurements for which no retrieval algorithm was known before, e.g., MIIPS measurements.


title: Ultra-Narrow Linewidth Brillouin Lasers with Nanokelvin Thermometry

abstract: Ultrastable lasers serve as the backbone for some of the most advanced scientific experiments today and enable the ability to perform atomic spectroscopy and laser interferometry at the highest levels of precision possible. With the recent and increasing interest in applying these systems outside of the laboratory, it remains an open question as how to realize a laser source that can reach the extraordinary levels of narrow linewidth required and yet still remain sufficiently compact and portable for field use. Critical to the development of this ideal laser source is the necessity for the laser to be insensitive to both short- and long-term fluctuations in temperature, which ultimately broaden the laser linewidth and cause drift in the laser's center frequency. We show here that the use of a large mode-volume optical resonator, which acts to suppress the resonator's fast thermal fluctuations, together with the stimulated Brillouin scattering (SBS) optical nonlinearity presents a powerful combination that enables the ability to lase with an ultra-narrow linewidth of 20 Hz. To address the laser's long-term temperature drift, we apply the narrow Brillouin line as a metrological tool that precisely senses a minute change in the resonator's temperature at the level of 85 nK. The precision afforded by this temperature measurement enables new possibilities for the stabilization of resonators against environmental perturbation.


title: Maximum Likelihood Analysis of Reaction Coordinates during Solidification in Ni

abstract: Understanding the underlying mechanism of crystal nucleation during solidification is a fundamental aspect in the prediction and control of materials properties. Classical nucleation theory (CNT) assumes that homogeneous nucleation occurs via random fluctuations within the supercooled liquid, that the structure of the growing clusters resembles the most stable bulk phase, and that the nucleus size is the sole reaction coordinate (RC) of the process. Many materials are, however, known to exhibit multiple steps during crystallization, forming different polymorphs. As a consequence, more complex RCs are often required to capture all relevant information about the process. In this work, we employ transition path sampling together with a maximum likelihood analysis of candidate order parameters to identify suitable reaction coordinates for the nucleation mechanism during solidification in Ni. In contrast to CNT, the analysis of the reweighted path ensemble shows that a pre-structured liquid region that surrounds the crystal cluster is a relevant order parameter that enhances the RC and therefore plays a key role in the description of the growing nucleus and the interfacial free energy. We demonstrate that pre-structured liquid clusters that emerge within the liquid act as precursors of the crystallization in a non-classical two-step mechanism which predetermines the coordination of the polymorphs that are being selected.


title: Stability, convergence, and limit cycles in some human physiological processes

abstract: Mathematical models for physiological processes aid qualitative understanding of the impact of various parameters on the underlying process. We analyse two such models for human physiological processes: the Mackey-Glass and the Lasota equations, which model the change in the concentration of blood cells in the human body. We first study the local stability of these models, and derive bounds on various model parameters and the feedback delay for the concentration to equilibrate. We then deduce conditions for non-oscillatory convergence of the solutions, which could ensure that the blood cell concentration does not oscillate. Further, we define the convergence characteristics of the solutions which govern the rate at which the concentration equilibrates when the system is stable. Owing to the possibility that physiological parameters can seldom be estimated precisely, we also derive bounds for robust stability\textemdash which enable one to ensure that the blood cell concentration equilibrates despite parametric uncertainty. We also highlight that when the necessary and sufficient condition for local stability is violated, the system transits into instability via a Hopf bifurcation, leading to limit cycles in the blood cell concentration. We then outline a framework to characterise the type of the Hopf bifurcation and determine the asymptotic orbital stability of limit cycles. The analysis is complemented with numerical examples, stability charts and bifurcation diagrams. The insights into the dynamical properties of the mathematical models may serve to guide the study of dynamical diseases.