title: On the persistence of lytic bacteriophages in vivo and its consequences for bacteriophage therapy

abstract: Bacteriophages are viruses infecting bacteria and archaea. Many phage species cause infections which lead to the certain death of the infected prokaryotic host cell and the release of a large batch of phage progeny, yet they have been able to stably coexist with their bacterial hosts over the eons in nature, as well as in the majority of laboratory experiments reported in the literature. This possibility of a stable coexistence between populations of bacteria and bacteriophages can be suspected to critically reduce the chances of a successful therapeutic application of bacteriophages for combatting pathogenic bacterial infections in vivo. Here, we extend an established differential equation model describing the interaction of bacteria, bacteriophages, and the host immune system, modelling different degrees of spatial heterogeneity of the host organism by introducing a scaling parameter which alters the encounter rates of the different cell populations. We demonstrate by rigorous mathematical analysis that, depending on the degree of spatial heterogeneity, the system will either converge to the desired state of bacterial and phage extinction, to a state of persisting bacterial infection with a completely eliminated bacteriophage population or to an equally undesirable state of a stable long-term bacteria-bacteriophage coexistence. We additionally provide numerical solutions of the model to illustrate the emerging dynamics and discuss some implications for bacteriophage therapy.


title: Catalan words avoiding pairs of length three patterns

abstract: Catalan words are particular growth-restricted words counted by the eponymous integer sequence. In this article we consider Catalan words avoiding a pair of patterns of length 3, pursuing the recent initiating work of the first and last authors and of S. Kirgizov where (among other things) the enumeration of Catalan words avoiding a patterns of length 3 is completed. More precisely, we explore systematically the structural properties of the sets of words under consideration and give enumerating results by means of recursive decomposition, constructive bijections or bivariate generating functions with respect to the length and descent number. Some of the obtained enumerating sequences are known, and thus the corresponding results establish new combinatorial interpretations for them.


title: Three-dimensional and four-dimensional scalar, vector, tensor cosmological fluctuations and the cosmological decomposition theorem

abstract: In cosmological perturbation theory it is convenient to use the scalar, vector, tensor (SVT) basis as defined according to how these components transform under 3-dimensional rotations. In attempting to solve the fluctuation equations that are automatically written in terms of gauge-invariant combinations of these components, the equations are taken to break up into separate SVT sectors, the decomposition theorem. Here, without needing to specify a gauge, we solve the fluctuation equations exactly for some standard cosmologies, to show that in general the various gauge-invariant combinations only separate at a higher-derivative level. To achieve separation at the level of the fluctuation equations themselves one has to assume boundary conditions for the higher-derivative equations. While asymptotic conditions suffice for fluctuations around a dS background or a $k=0$ RW background, for fluctuations around a $k\neq 0$ RW background one additionally has to require that the fluctuations be well-behaved at the origin. We show that in certain cases the gauge-invariant combinations themselves involve both scalars and vectors. For such cases there is no decomposition theorem for the individual SVT components themselves, but for the gauge-invariant combinations there still can be. Given the lack of manifest covariance in defining a basis with respect to 3-dimensional rotations, we introduce an alternate SVT basis whose components are defined according to how they transform under 4-dimensional general coordinate transformations. With this basis the fluctuation equations greatly simplify, and while one can again break them up into separate gauge-invariant sectors at the higher-derivative level, in general we find that even with boundary conditions we do not obtain a decomposition theorem in which the fluctuations separate at the level of the fluctuation equations themselves.


title: Artifacts and Errors in Cross-Spectrum Phase Noise Measurements

abstract: This article deals with the erratic and inconsistent phase-noise spectra often seen in low-noise oscillators, whose floor is of the order of $-180$ dBc/Hz or less. Such oscillators are generally measured with two-channel instruments based on averaging two simultaneous and statistically independent measures. Our new method consists of inserting a dissipative attenuator between the oscillator under test and the phase-noise analyzer. The thermal noise of the attenuator introduces a controlled amount of phase noise. We compare the phase noise floor to the theoretical expectation with different values of the attenuation in small steps. The analysis reveals a negative bias (underestimation of phase noise) due to the thermal energy of the internal power splitter at the instrument input, and an uncertainty due to crosstalk between the two channels. In not-so-rare unfortunate cases, the bias results in a negative phase-noise spectrum, which is an obvious nonsense. Similar results are observed separately in three labs with instruments from the two major brands. We give experimental evidence, full theory, and suggestions to mitigate the problem. Our multiple-attenuators method provides quantitative information about the correlation phenomena inside the instrument.


title: Formation and critical dynamics of topological defects in Lifshitz holography

abstract: We examine the formation and critical dynamics of topological defects via Kibble-Zurek mechanism in a (2+1)-dimensional quantum critical point, which is conjectured to dual to a Lifshitz geometry. Quantized magnetic fluxoids are spontaneously generated and trapped in the cores of order parameter vortices, which is a feature of type II superconductors. Time evolution of the average condensate is found to lag behind the instantaneous equilibrium value of the order parameter, a typical phenomenon in nonequilibrium dynamics. Scalings of vortex number density and the "freeze-out" time match the predictions from Kibble-Zurek mechanism. From these scalings, the dynamic and static critical exponents in the boundary field theory are found, at least at finite temperature, to be irrespective of the Lifshitz exponent in the bulk.


title: Long-time behavior of a reaction-diffusion model with strong Allee effect and free boundary: effect of a protection zone

abstract: This paper concerns the effect of the (separated/connected) protection zone for the evolution of an endangered species on the reaction-diffusion equation with strong Allee effect and free boundary. We give a description of the long-time dynamical behavior of the problem of two types protection zones with the same length. Furthermore, the asymptotic profiles of solutions and the asymptotic spreading speed are estimated when spreading happens. Our results, together with those in previous papers [8, 12] on two other closely related models, show that the protection zone and the free boundary play an important role in the evolution of the endangered species.


title: Tune-out and magic wavelengths for ground-state $^{23}$Na$^{40}$K molecules

abstract: We demonstrate a versatile, rotational-state dependent trapping scheme for the ground and first excited rotational states of $^{23}$Na$^{40}$K molecules. Close to the rotational manifold of a narrow electronic transition, we determine tune-out frequencies where the polarizability of one state vanishes while the other remains finite, and a magic frequency where both states experience equal polarizability. The proximity of these frequencies of only 10 GHz allows for dynamic switching between different trap configurations in a single experiment, while still maintaining sufficiently low scattering rates.


title: Effects of momentum-dependent quasiparticle renormalization on the gap structure of iron-based superconductors

abstract: We discuss the influence of momentum-dependent correlations on the superconducting gap structure in iron-based superconductors. Within the weak coupling approach including self-energy effects at the one-loop spin-fluctuation level, we construct a dimensionless pairing strength functional which includes the effects of quasiparticle renormalization. The stationary solution of this equation determines the gap function at $T_c$. The resulting equations represent the simplest generalization of spin fluctuation pairing theory to include the effects of an anisotropic quasiparticle weight. We obtain good agreement with experimentally observed anisotropic gap structures in LiFeAs, indicating that the inclusion of quasiparticle renormalization effects in the existing weak-coupling theories can account for the observed anomalies in the gap structure of Fe-based superconductors.


title: On the Initialization of Long Short-Term Memory Networks

abstract: Weight initialization is important for faster convergence and stability of deep neural networks training. In this paper, a robust initialization method is developed to address the training instability in long short-term memory (LSTM) networks. It is based on a normalized random initialization of the network weights that aims at preserving the variance of the network input and output in the same range. The method is applied to standard LSTMs for univariate time series regression and to LSTMs robust to missing values for multivariate disease progression modeling. The results show that in all cases, the proposed initialization method outperforms the state-of-the-art initialization techniques in terms of training convergence and generalization performance of the obtained solution.


title: Dissipative hydrodynamics of relativistic shock waves in a Quark Gluon Plasma: comparing and benchmarking alternate numerical methods

abstract: This paper presents numerical cross-comparisons and benchmark results for two different kinetic numerical methods, capable of describing relativistic dissipative fluid dynamics in a wide range of kinematic regimes, typical of relevant physics applications, such as transport phenomena in quark-gluon plasmas. We refer to relativistic lattice Boltzmann versus Montecarlo Test-Particle methods. Lacking any realistic option for accurate validation vis-a-vis experimental data, we check the consistency of our results against established simulation packages available in the literature. We successfully cross-compare the results of the two aforementioned numerical approaches for momentum integrated quantities like the hydrostatic and dynamical pressure profiles, the collective flow and the heat flux. These results corroborate the confidence on the robustness and correctness of these computational methods and on the accurate calibration of their numerical parameters with respect to the physical transport coefficients. Our numerical results are made available as supplemental material, with the aim of establishing a reference benchmark for other numerical approaches.


title: LB-1 is inconsistent with the X-ray source population and pulsar-black hole binary searches in the Milky Way

abstract: If confirmed, a wide binary system of 70 $M_{\odot}$ black hole (BH) and an 8 $M_{\odot}$ main sequence star (LB-1) is observed to reside in the Milky Way (MW). We show that long term evolution of an 8 $M_{\odot}$ star around a BH with mass between 5-70 $M_{\odot}$ makes them visible as ultra-luminous X-ray (ULX) sources in the sky. Given the expected ULX phase lifetime ($\approx0.1$ Myr) and their lack of detection in the MW, we conclude that the frequency of an 8-20 $M_{\odot}$ star to be in binary around a stellar mass BH should be less than ($f<2\times10^{-3}$). This is in tension with Liu et al. (2019) claimed detection frequency of LB-1 like system around 8-20$M_{\odot}$ stars ($f\approx3\times10^{-2}$). Moreover, the 8 $M_{\odot}$ star is likely to end as a neutron star (NS) born with a very small kick from an electron capture supernova (ECSN), leaving behind a wide NS-BH binary. So far less than 1\% of all the detectable pulsars in the MW are mapped and there has been no detection of any pulsars in binary systems around BHs which sets an upper bound of about 100 possible pulsar-BH systems in the MW. We show if the NS is born from ECSN, a frequency upper limit of ($f=5\times10^{-4}$) for stars with masses $\approx 8-20~M_{\odot}$ in the MW to have a BH companion. The rate discrepancy will further increase as more pulsars are mapped in the MW, yet these searches would not be able to rule out the Liu et al. detection frequency if NSs are instead born in core collapse SNe with the commonly inferred high kick velocities.


title: The Need to Reconcile Concepts that Characterize Systems Withstanding Threats

abstract: Desirable system performance in the face of threats and disruptions has been characterized by various management concepts. Through semi-structured interviews with editors of journals in the fields of emergency response and systems management, a literature review, and professional judgment, we identified nine related and often interchangeably-used system performance concepts: adaptability, agility, reliability, resilience, resistance, robustness, safety, security, and sustainability. We analysed expert responses and reviewed the linguistic definitions and mathematical framing of these concepts to understand their applications. We found a lack of consensus on their usage between interview subjects, but using a mathematical framing enriched the linguistic definitions and enabled formulating comparative visualizations and system-specific definitions for the concepts. We propose a conceptual framing to relate the concepts for management purposes. A better understanding of these concepts will allow system planners to pursue management strategies best suited to their unique system dynamics and specific objectives of 'goodness' that all these concepts bring.


title: Emotion Recognition from Speech

abstract: In this work, we conduct an extensive comparison of various approaches to speech based emotion recognition systems. The analyses were carried out on audio recordings from Ryerson Audio-Visual Database of Emotional Speech and Song (RAVDESS). After pre-processing the raw audio files, features such as Log-Mel Spectrogram, Mel-Frequency Cepstral Coefficients (MFCCs), pitch and energy were considered. The significance of these features for emotion classification was compared by applying methods such as Long Short Term Memory (LSTM), Convolutional Neural Networks (CNNs), Hidden Markov Models (HMMs) and Deep Neural Networks (DNNs). On the 14-class (2 genders x 7 emotions) classification task, an accuracy of 68% was achieved with a 4-layer 2 dimensional CNN using the Log-Mel Spectrogram features. We also observe that, in emotion recognition, the choice of audio features impacts the results much more than the model complexity.


title: A Novel Hybrid Opportunistic Scalable Energy Efficient Routing Design For Low Power, Lossy Wireless Sensor Networks

abstract: Opportunistic Routing (OR) scheme increases the transmission reliability despite the lossy wireless radio links by exploiting the broadcast nature of the wireless medium. However, OR schemes in low power Wireless Sensor Network (WSN) leads to energy drain in constrained sensor nodes due to constant overhearing, periodic beaconing for Neighbourhood Management (NM) and increase in packet header length to append priority wise sorted Forwarding Candidates Set (FCS) prior to data transmission. The focus of this work is to propose a hybrid opportunistic energy efficient routing design for large scale, low power and lossy WSN. This design avoids periodic 'hello' beacons for NM, limits constant overhearing and increase in packet header length. There are two modes of operation i) opportunistic ii) unicast mode. The sender node adopts opportunistic forwarding for its initial data packet transmission and instead of pre-computing the FCS, it is dynamically computed in a completely distributed manner. The nodes part of FCS based on cross layered multi-metrics and fuzzy decision logic determines its priority level to compute Dynamic Holding Delay (DHD) for effective timer coordination. The differentiated back off implementation along with DHD enables the higher priority candidate that had received data packet to forward the packet first and facilitates others to cancel its timer upon overhearing. The sender node switches to unicast mode for successive transmissions to achieve energy efficiency till link-level transmission error occur. Opportunistic mode is used as the first and last resort to increase the transmission reliability. Simulation results in NS2 show significant increase in Packet Delivery Ratio (PDR), decrease in both average energy consumption per node and Normalized Energy Consumption (NEC) per packet in comparison with existing protocols.


title: Quantum gravitational wave function for the interior of a black hole and the generalized uncertainty principle

abstract: We investigate the internal structures of a Schwarzschild black hole by solving the Wheeler-DeWitt equation. The generic bounded wave function has a bouncing point around $r \simeq M$, where $M$ is the black hole mass. Due to this quantum bouncing, there appears an ambiguity to define the arrow of time. If we introduce two arrows of time, one can then interpret that two classical spacetime is annihilated around the bouncing point. Finally, we provide a conceptual explanation based on the generalized uncertainty principle.


title: Dark origin of quark flavor hierarchy and mixing

abstract: We provide a dynamical mechanism for the generation of the Cabibbo-Kobayashi-Maskawa matrix in the context of a recently proposed model of flavor. The framework, based on the paradigm that Yukawa couplings are effective low energy couplings generated radiatively by dark sector interactions, is here extended to include a new scalar field which plays the role of a {\it dark flavon}. Being singlet under the Standard Model gauge group but charged under the $U(1)$ symmetry of the dark sector, this particle sources new two- and three-loop diagrams that result in off-diagonal Yukawa interactions, thereby providing a simple explanation for the observed structure of the CKM matrix. By using an effective parametrization of these new loop contributions, we show that the CKM matrix elements can be correctly reproduced for perturbative values of the involved couplings. We outline the rich phenomenology predicted by the scenario and detail its implications for the LHC and future collider experiments.


title: Lattice points in spherical segments

abstract: We study lattice points in d-dimensional spheres, and count their number in thin spherical segments. We found an upper bound depending only on the radius of the sphere and opening angle of the segment. To obtain this bound we slice the segment by hyperplanes of rational direction, and then cover an arbitrary segment with one having rational direction. Diophantine approximation can be used to obtain the best rational direction possible.


title: Stochastic Recursive Optimal Control Problem with Mixed Delay under Viscosity Solution's Framework

abstract: This paper is concerned with the stochastic recursive optimal control problem with mixed delay. The connection between Pontryagin's maximum principle and Bellman's dynamic programming principle is discussed. Without containing any derivatives of the value function, relations among the adjoint processes and the value function are investigated by employing the notions of super- and sub-jets introduced in defining the viscosity solutions. Stochastic verification theorem is also given to verify whether a given admissible control is really optimal.


title: Storage Control for Carbon Emission Reduction: Opportunities and Challenges

abstract: Storage is vital to power systems as it provides the urgently needed flexibility to the system. Meanwhile, it can contribute more than flexibility. In this paper, we study the possibility of utilizing storage system for carbon emission reduction. The opportunity arises due to the pending implementation of carbon tax throughout the world. Without the right incentive, most system operators have to dispatch the generators according to the merit order of the fuel costs, without any control for carbon emissions. However, we submit that storage may provide necessary flexibility in carbon emission reduction even without carbon tax. We identify the non-convex structure to conduct storage control for this task and propose an easy to implement dynamic programming algorithm to investigate the value of storage in carbon emission reduction.


title: Topological full groups of ultragraph groupoids as an isomorphism invariant

abstract: We prove two isomorphism-invariance theorems for groupoids associated with ultragraphs. These theorems characterize ultragraphs for which the topological full group of an associated groupoid is an isomorphism invariant. These results extend those of graph groupoids to ultragraph groupoids while providing another concrete example where the topological full group of a groupoid is a complete isomorphism invariant.


title: Broadband optical parametric amplification by two-dimensional semiconductors

abstract: Optical parametric amplification is a second-order nonlinear process whereby an optical signal is amplified by a pump via the generation of an idler field. It is the key ingredient of tunable sources of radiation that play an important role in several photonic applications. This mechanism is inherently related to spontaneous parametric down-conversion that currently constitutes the building block for entangled photon pair generation, which has been exploited in modern quantum technologies ranging from computing to communications and cryptography. Here we demonstrate single-pass optical parametric amplification at the ultimate thickness limit; using semiconducting transition-metal dichalcogenides, we show that amplification can be attained over a propagation through a single atomic layer. Such a second-order nonlinear interaction at the 2D limit bypasses phase-matching requirements and achieves ultrabroad amplification bandwidths. The amplification process is independent on the in-plane polarization of the impinging signal and pump fields. First-principle calculations confirm the observed polarization invariance and linear relationship between idler and pump powers. Our results pave the way for the development of atom-sized tunable sources of radiation with applications in nanophotonics and quantum information technology.


title: Two novel results on the existence of $3$-kernels in digraphs

abstract: Let $D$ be a digraph. We call a subset $N$ of $V(D)$ $k$-independent if for every pair of vertices $u,v \in N$, $d(u,v) \geq k$; and we call it $\ell$-absorbent if for every vertex $u \in V(D) \setminus N$, there exists $v \in N$ such that $d(u,v) \leq \ell$. A $(k,\ell)$-kernel of $D$ is a subset of vertices which is $k$-independent and $\ell$-absorbent. A $k$-kernel is a $(k,k-1)$-kernel. In this report, we present the main results from our master's research regarding kernel theory. We prove that if a digraph $D$ is strongly connected and every cycle $C$ of $D$ satisfies: $(i)$ if $C \equiv 0 \pmod 3$, then $C$ has a short chord and $(ii)$ if $C \not \equiv 0 \pmod 3$, then $C$ has three short chords: two consecutive and a third crossing one of the former, then $D$ has a $3$-kernel. Moreover, we introduce a modification of the substitution method, proposed by Meyniel and Duchet in 1983, for $3$-kernels and use it to prove that a quasi-$3$-kernel-perfect digraph $D$ is $3$-kernel-perfect if every circuit of length not dividable by three has four short chords.


title: A biology-inspired approach to the positive integral control of positive systems -- the antithetic, exponential, and logistic integral controllers

abstract: The integral control of positive systems using nonnegative control input is an important problem arising, among others, in biochemistry, epidemiology and ecology. An immediate solution is to use an ON-OFF nonlinearity between the controller and the system. However, this solution is only available when controllers are implemented in computer systems. When this is not the case, like in biology, alternative approaches need to be explored. Based on recent research in the control of biological systems, we propose to develop a theory for the integral control of positive systems using nonnegative controls based on the so-called \emph{antithetic integral controller} and two \emph{positively regularized integral controllers}, the so-called \emph{exponential integral controller} and \emph{logistic integral controller}. For all these controllers, we establish several qualitative results, which we connect to standard results on integral control. We also obtain additional results which are specific to the type of controllers. For instance, we show an interesting result stipulating that if the gain of the antithetic integral controller is suitably chosen, then the local stability of the equilibrium point of the closed-loop system does not depend on the choice for the coupling parameter, an additional parameter specific to this controller. Conversely, we also show that if the coupling parameter is suitably chosen, then the equilibrium point of the closed-loop system is locally stable regardless the value of the gain. For the exponential integral controller, we can show that the local stability of the equilibrium point of the closed-loop system is independent of the gain of the controller and the gain of the system. The stability only depends on the exponential rate of the controller, again a parameter that is specific to this type of controllers. Several examples are given for illustration.


title: Classification of Vertex-Transitive Zonotopes

abstract: We give a full classification of vertex-transitive zonotopes. We prove that a vertex-transitive zonotope is a $\Gamma$-permutahedron for some finite reflection group $\Gamma\subset\mathrm{O}(\mathbb R^d)$. The same holds true for zonotopes in which all vertices are on a common sphere, and all edges are of the same length (which we call homogeneous zonotopes). The classification of these then follows from the classification of finite reflection groups. We proof that root systems can be characterized as those centrally symmetric sets of vectors, for which all intersections with half-spaces, that contain exactly half the vectors, are congruent. We provide a further sufficient condition for a centrally symmetric set being a root system.


title: On the topology of real Lagrangians in toric symplectic manifolds

abstract: We explore the topology of real Lagrangian submanifolds in a toric symplectic manifold which come from involutive symmetries on its moment polytope. We establish a real analog of the Delzant construction for those real Lagrangians, which says that their diffeomorphism type is determined by combinatorial data. As an application, we realize all possible diffeomorphism types of connected real Lagrangians in toric symplectic del Pezzo surfaces.


title: Evidence of a dynamically evolving Galactic warp

abstract: In a cosmological setting, the disc of a galaxy is expected to continuously experience gravitational torques and perturbations from a variety of sources, which can cause the disc to wobble, flare and warp. Specifically, the study of galactic warps and their dynamical nature can potentially reveal key information on the formation history of galaxies and the mass distribution of their halos. Our Milky Way presents a unique case study for galactic warps, thanks to the detailed knowledge of its stellar distribution and kinematics. Using a simple model of how the warp's orientation is changing with time, we here measure the precession rate of the Milky Way's warp using 12 million giant stars from Gaia Data Release 2, finding that it is precessing at $10.86 \pm 0.03_{stat} \pm 3.20_{sys}$ km/s/kpc in the direction of Galactic rotation, about one third the angular velocity at the Sun's position in the Galaxy. The direction and magnitude of the warp's precession rate favour the scenario that the warp is the results of a recent or ongoing encounter with a satellite galaxy, rather than the relic of the ancient assembly history of the Galaxy.


title: Testing the equality of multivariate means when $p>n$ by combining the Hoteling and Simes tests

abstract: We propose a method of testing the shift between mean vectors of two multivariate Gaussian random variables in a high-dimensional setting incorporating the possible dependency and allowing $p > n$. This method is a combination of two well-known tests: the Hotelling test and the Simes test. The tests are integrated by sampling several dimensions at each iteration, testing each using the Hotelling test, and combining their results using the Simes test. We prove that this procedure is valid asymptotically. This procedure can be extended to handle non-equal covariance matrices by plugging in the appropriate extension of the Hotelling test. Using a simulation study, we show that the proposed test is advantageous over state-of-the-art tests in many scenarios and robust to violation of the Gaussian assumption.


title: Sharp asymptotics in a fractional Sturm-Liouville problem

abstract: The current research of fractional Sturm-Liouville boundary value problems focuses on the qualitative theory and numerical methods, and much progress has been recently achieved in both directions. The objective of this paper is to explore a different route, namely, construction of explicit asymptotic approximations for the solutions. As a study case, we consider a problem with left and right Riemann-Liouville derivatives, for which our analysis yields asymptotically sharp estimates for the sequence of eigenvalues and eigenfunctions.


title: Fluctuation theory for spectrally positive additive L\'evy fields

abstract: A spectrally positive additive L\'evy field is a multidimensional field obtained as the sum $\mathbf{X}_{\rm t}={\rm X}^{(1)}_{t_1}+{\rm X}^{(2)}_{t_2}+\dots+{\rm X}^{(d)}_{t_d}$, ${\rm t}=(t_1,\dots,t_d)\in\mathbb{R}_+^d$, where ${\rm X}^{(j)}={}^t (X^{1,j},\dots,X^{d,j})$, $j=1,\dots,d$, are $d$ independent $\mathbb{R}^d$-valued L\'evy processes issued from 0, such that $X^{i,j}$ is non decreasing for $i\neq j$ and $X^{j,j}$ is spectrally positive. It can also be expressed as $\mathbf{X}_{\rm t}=\mathbb{X}_{\rm t}{\bf 1}$, where ${\bf 1}={}^t(1,1,\dots,1)$ and $\mathbb{X}_{\rm t}=(X^{i,j}_{t_j})_{1\leq i,j\leq d}$. The main interest of spaLf's lies in the Lamperti representation of multitype continuous state branching processes. In this work, we study the law of the first passage times $\mathbf{T}_{\rm r}$ of such fields at levels $-{\rm r}$, where ${\rm r}\in\mathbb{R}_+^d$. We prove that the field $\{(\mathbf{T}_{\rm r},\mathbb{X}_{\mathbf{T}_{\rm r}}),{\rm r}\in\mathbb{R}_+^d\}$ has stationary and independent increments and we describe its law in terms of this of the spaLf $\mathbf{X}$. In particular, the Laplace exponent of $(\mathbf{T}_{\rm r},\mathbb{X}_{\mathbf{T}_{\rm r}})$ solves a functional equation leaded by the Laplace exponent of $\mathbf{X}$. This equation extends in higher dimension a classical fluctuation identity satisfied by the Laplace exponents of the ladder processes. Then we give an expression of the distribution of $\{(\mathbf{T}_{\rm r},\mathbb{X}_{\mathbf{T}_{\rm r}}),{\rm r}\in\mathbb{R}_+^d\}$ in terms of the distribution of $\{\mathbb{X}_{\rm t},{\rm t}\in\mathbb{R}_+^d\}$ by the means of a Kemperman-type formula, well-known for spectrally positive L\'evy processes.


title: Multi-type Dirac fermions protected by orthogonal glide symmetries in a noncentrosymmetric system

abstract: Compared to inversion-symmetric systems, emerging bulk Dirac point (DP) in noncentrosymmetric systems is much harder and its symmetry protection mechanism is poorly understood. In this work, we propose that orthogonal glide symmetries can protect two distinct types of bulk anisotropic DPs. One is the ordinary anisotropic DP that is doubly degenerate only along one invariant axis, which is resulted from symmetry-protected accidental band crossing. The other one is a symmetry-enforced DP fixed at a non-time-reversal-invariant point, which is doubly degenerate at three orthogonal directions. This unique topological phase is exemplified by KSnSe$_{2}$ in space group 108. Our work not only unveils an unique symmetry protection mechanism but also provides the first material candidate for exploring multi-type Dirac fermions in noncentrosymmetric systems.


title: Cloud model inversions of strong chromospheric absorption lines using principal component analysis

abstract: High-resolution spectroscopy of strong chromospheric absorption lines delivers nowadays several millions of spectra per observing day, when using fast scanning devices to cover large regions on the solar surface. Therefore, fast and robust inversion schemes are needed to explore the large data volume. Cloud Model (CM) inversions of the chromospheric H$\alpha$ line are commonly employed to investigate various solar features including filaments, prominences, surges, jets, mottles, and (macro-)spicules. The choice of the CM was governed by its intuitive description of complex chromospheric structures as clouds suspended above the solar surface by magnetic fields. This study is based on observations of active region NOAA 11126 in H$\alpha$, which were obtained 2010 November 18-23 with the echelle spectrograph of the Vacuum Tower Telescope (VTT) at the Observatorio del Teide, Spain. Principal Component Analysis (PCA) reduces the dimensionality of spectra and conditions noise-stripped spectra for CM inversions. Modeled H$\alpha$ intensity and contrast profiles as well as CM parameters are collected in a database, which facilitates efficient processing of the observed spectra. Physical maps are computed representing the line-core and continuum intensity, absolute contrast, equivalent width, and Doppler velocities, among others. Noise-free spectra expedite the analysis of bisectors. The data processing is evaluated in the context of "big data", in particular with respect to automatic classification of spectra.


title: Electroweak Phase Transition in Non-Minimal Higgs Sectors

abstract: Higgs sector extensions beyond the Standard Model (BSM) provide additional sources of CP violation and further scalar states that help to trigger a strong first order electroweak phase transition (SFOEWPT) required to generate the observed baryon asymmetry of the Universe through electroweak baryogenesis. We investigate the CP-violating 2-Higgs-Doublet Model (C2HDM) and the Next-to-Minimal 2-Higgs-Doublet Model (N2HDM) with respect to their potential to generate an SFOEWPT while being compatible with all relevant and recent theoretical and experimental constraints. The implications of an SFOEWPT on the collider phenomenology of the two models are analysed in detail in particular with respect to Higgs pair production. We provide benchmark points for parameter points that are compatible with an SFOEWPT and provide distinct di-Higgs signatures.


title: Domain structured chaos for discrete random processes

abstract: We introduce the notion of domain-structured chaos and apply it to establish a connection between stochastic dynamics and deterministic chaos.


title: Facial Synthesis from Visual Attributes via Sketch using Multi-Scale Generators

abstract: Automatic synthesis of faces from visual attributes is an important problem in computer vision and has wide applications in law enforcement and entertainment. With the advent of deep generative convolutional neural networks (CNNs), attempts have been made to synthesize face images from attributes and text descriptions. In this paper, we take a different approach, where we formulate the original problem as a stage-wise learning problem. We first synthesize the facial sketch corresponding to the visual attributes and then we generate the face image based on the synthesized sketch. The proposed framework, is based on a combination of two different Generative Adversarial Networks (GANs) - (1) a sketch generator network which synthesizes realistic sketch from the input attributes, and (2) a face generator network which synthesizes facial images from the synthesized sketch images with the help of facial attributes. Extensive experiments and comparison with recent methods are performed to verify the effectiveness of the proposed attribute-based two-stage face synthesis method.


title: Parseval Proximal Neural Networks

abstract: The aim of this paper is twofold. First, we show that a certain concatenation of a proximity operator with an affine operator is again a proximity operator on a suitable Hilbert space. Second, we use our findings to establish so-called proximal neural networks (PNNs) and stable tight frame proximal neural networks. Let $\mathcal H$ and $\mathcal K$ be real Hilbert spaces, $b\in\mathcal K$ and $T\in\mathcal{B}(\mathcal H,\mathcal K)$ have closed range and Moore-Penrose inverse $T^\dagger$. Based on the well-known characterization of proximity operators by Moreau, we prove that for any proximity operator $\text{Prox}\colon\mathcal K\to\mathcal K$ the operator $T^\dagger\,\text{Prox} (T\cdot +b)$ is a proximity operator on $\mathcal H$ equipped with a suitable norm. In particular, it follows for the frequently applied soft shrinkage operator $\text{Prox} = S_{\lambda}\colon\ell_2 \rightarrow\ell_2$ and any frame analysis operator $T\colon\mathcal H\to\ell_2$ that the frame shrinkage operator $T^\dagger\, S_\lambda\,T$ is a proximity operator on a suitable Hilbert space. The concatenation of proximity operators on $\mathbb R^d$ equipped with different norms establishes a PNN. If the network arises from tight frame analysis or synthesis operators, then it forms an averaged operator. Hence, it has Lipschitz constant 1 and belongs to the class of so-called Lipschitz networks, which were recently applied to defend against adversarial attacks. Moreover, due to its averaging property, PNNs can be used within so-called Plug-and-Play algorithms with convergence guarantee. In case of Parseval frames, we call the networks Parseval proximal neural networks (PPNNs). Then, the involved linear operators are in a Stiefel manifold and corresponding minimization methods can be applied for training. Finally, some proof-of-the concept examples demonstrate the performance of PPNNs.


title: A Systematic Comparison of Bayesian Deep Learning Robustness in Diabetic Retinopathy Tasks

abstract: Evaluation of Bayesian deep learning (BDL) methods is challenging. We often seek to evaluate the methods' robustness and scalability, assessing whether new tools give `better' uncertainty estimates than old ones. These evaluations are paramount for practitioners when choosing BDL tools on-top of which they build their applications. Current popular evaluations of BDL methods, such as the UCI experiments, are lacking: Methods that excel with these experiments often fail when used in application such as medical or automotive, suggesting a pertinent need for new benchmarks in the field. We propose a new BDL benchmark with a diverse set of tasks, inspired by a real-world medical imaging application on \emph{diabetic retinopathy diagnosis}. Visual inputs (512x512 RGB images of retinas) are considered, where model uncertainty is used for medical pre-screening---i.e. to refer patients to an expert when model diagnosis is uncertain. Methods are then ranked according to metrics derived from expert-domain to reflect real-world use of model uncertainty in automated diagnosis. We develop multiple tasks that fall under this application, including out-of-distribution detection and robustness to distribution shift. We then perform a systematic comparison of well-tuned BDL techniques on the various tasks. From our comparison we conclude that some current techniques which solve benchmarks such as UCI `overfit' their uncertainty to the dataset---when evaluated on our benchmark these underperform in comparison to simpler baselines. The code for the benchmark, its baselines, and a simple API for evaluating new BDL tools are made available at https://github.com/oatml/bdl-benchmarks.


title: Excellent electronic transport in heterostructures of graphene and monoisotopic boron-nitride grown at atmospheric pressure

abstract: Hexagonal boron nitride (BN), one of the very few layered insulators, plays a crucial role in 2D materials research. In particular, BN grown with a high pressure technique has proven to be an excellent substrate material for graphene and related 2D materials, but at the same time very hard to replace. Here we report on a method of growth at atmospheric pressure as a true alternative for producing BN for high quality graphene/BN heterostructures. The process is not only more scalable, but also allows to grow isotopically purified BN crystals. We employ Raman spectroscopy, cathodoluminescence, and electronic transport measurements to show the high-quality of such monoisotopic BN and its potential for graphene-based heterostructures. The excellent electronic performance of our heterostructures is demonstrated by well developed fractional quantum Hall states, ballistic transport over distances around $10\,\mathrm{\mu m}$ at low temperatures and electron-phonon scattering limited transport at room temperature.


title: Vector Representations of Graphs and Distinguishing Quantum Product States with One-way LOCC

abstract: Distinguishing sets of quantum states shared by two parties using only local operations and classical communication measurements is a fundamental topic in quantum communication and quantum information theory. We introduce a graph-theoretic approach, based on the theory of vector representations of graphs, to the core problem of distinguishing product states with one-way LOCC. We establish a number of results that show how distinguishing such states can be framed in terms of properties of the underlying graphs associated with a set of vector product states. We also present a number of illustrative examples.


title: Stability for inverse source problems by Carleman estimates

abstract: In this article, we provide a modified argument for proving conditional stability for inverse problems of determining spatially varying functions in evolution equations by Carleman estimates. Our method needs not any cut-off procedures and can simplify the existing proofs. We establish the conditional stability for inverse source problems for a hyperbolic equation and a parabolic equation, and our method is widely applicable to various evolution equations.


title: Energy-Aware Multi-Server Mobile Edge Computing: A Deep Reinforcement Learning Approach

abstract: We investigate the problem of computation offloading in a mobile edge computing architecture, where multiple energy-constrained users compete to offload their computational tasks to multiple servers through a shared wireless medium. We propose a multi-agent deep reinforcement learning algorithm, where each server is equipped with an agent, observing the status of its associated users and selecting the best user for offloading at each step. We consider computation time (i.e., task completion time) and system lifetime as two key performance indicators, and we numerically demonstrate that our approach outperforms baseline algorithms in terms of the trade-off between computation time and system lifetime.


title: A Polynomial Time Algorithm for the $k$-Disjoint Shortest Paths Problem

abstract: The disjoint paths problem is a fundamental problem in algorithmic graph theory and combinatorial optimization. For a given graph $G$ and a set of $k$ pairs of terminals in $G$, it asks for the existence of $k$ vertex-disjoint paths connecting each pair of terminals. The proof of Robertson and Seymour [JCTB 1995] of the existence of an $n^3$ algorithm for any fixed $k$ is one of the highlights of their Graph Minors project. In this paper, we focus on the version of the problem where all the paths are required to be shortest paths. This problem, called the disjoint shortest paths problem, was introduced by Eilam-Tzoreff [DAM 1998] where she proved that the case $k = 2$ admits a polynomial time algorithm. This problem has received some attention lately, especially since the proof of the existence of a polynomial time algorithm in the directed case when $k = 2$ by B\'erczi and Kobayashi [ESA 2017]. However, the existence of a polynomial algorithm when $k = 3$ in the undirected version remained open since 1998. In this paper we show that for any fixed $k$, the disjoint shortest paths problem admits a polynomial time algorithm. In fact for any fixed $C$, the algorithm can be extended to treat the case where each path connecting the pair $(s,t)$ has length at most $d(s,t) + C$.


title: Exploring Dark Matter, Neutrino mass and flavour anomalies in L_{\mu}-L_{\tau} model

abstract: We investigate Majorana dark matter in a new variant of $U(1)_{L_{\mu}-L_{\tau}}$ gauge extension of Standard Model, containing three additional neutral fermions $N_{e}, N_{\mu}, N_{\tau}$, along with a $(\bar{3},1,1/3)$ scalar Leptoquark (SLQ) and an inert doublet, to study the phenomenology of dark matter, neutrino mass generation and flavour anomalies on a single platform. The lightest mass eigenstate of the $N_{\mu}, N_{\tau}$ neutral fermions plays the role of dark matter. We compute the WIMP-nucleon cross section in leptoquark portal and the relic density mediated by inert doublet components, leptoquark and the new $Z^{\prime}$ boson. We constrain the parameter space consistent with Planck limit on relic density, PICO-60 and LUX bounds on spin-dependent direct detection cross section. We also discuss about the neutrino mass generation at one-loop level and the viable parameter region to explain current neutrino oscillation data. The $Z^\prime$ gauge boson of extended $U(1)$ symmetry and the SLQ play an important role in settling the known issues of flavor sector.


title: Efficient and Convergent Sequential Pseudo-Likelihood Estimation of Dynamic Discrete Games

abstract: We propose a new sequential Efficient Pseudo-Likelihood (k-EPL) estimator for dynamic discrete choice games of incomplete information. k-EPL considers the joint behavior of multiple players simultaneously, as opposed to individual responses to other agents' equilibrium play. This, in addition to reframing the problem from conditional choice probability (CCP) space to value function space, yields a computationally tractable, stable, and efficient estimator. We show that each iteration in the k-EPL sequence is consistent and asymptotically efficient, so the first-order asymptotic properties do not vary across iterations. Furthermore, we show the sequence achieves higher-order equivalence to the finite-sample maximum likelihood estimator with iteration and that the sequence of estimators converges almost surely to the maximum likelihood estimator at a nearly-superlinear rate when the data are generated by any regular Markov perfect equilibrium, including equilibria that lead to inconsistency of other sequential estimators. When utility is linear in parameters, k-EPL iterations are computationally simple, only requiring that the researcher solve linear systems of equations to generate pseudo-regressors which are used in a static logit/probit regression. Monte Carlo simulations demonstrate the theoretical results and show k-EPL's good performance in finite samples in both small- and large-scale games, even when the game admits spurious equilibria in addition to one that generated the data. We apply the estimator to study the role of competition in the U.S. wholesale club industry.


title: Recurrent Feedback Improves Feedforward Representations in Deep Neural Networks

abstract: The abundant recurrent horizontal and feedback connections in the primate visual cortex are thought to play an important role in bringing global and semantic contextual information to early visual areas during perceptual inference, helping to resolve local ambiguity and fill in missing details. In this study, we find that introducing feedback loops and horizontal recurrent connections to a deep convolution neural network (VGG16) allows the network to become more robust against noise and occlusion during inference, even in the initial feedforward pass. This suggests that recurrent feedback and contextual modulation transform the feedforward representations of the network in a meaningful and interesting way. We study the population codes of neurons in the network, before and after learning with feedback, and find that learning with feedback yielded an increase in discriminability (measured by d-prime) between the different object classes in the population codes of the neurons in the feedforward path, even at the earliest layer that receives feedback. We find that recurrent feedback, by injecting top-down semantic meaning to the population activities, helps the network learn better feedforward paths to robustly map noisy image patches to the latent representations corresponding to important visual concepts of each object class, resulting in greater robustness of the network against noises and occlusion as well as better fine-grained recognition.


title: Learning Improved Representations by Transferring Incomplete Evidence Across Heterogeneous Tasks

abstract: Acquiring ground truth labels for unlabelled data can be a costly procedure, since it often requires manual labour that is error-prone. Consequently, the available amount of labelled data is increasingly reduced due to the limitations of manual data labelling. It is possible to increase the amount of labelled data samples by performing automated labelling or crowd-sourcing the annotation procedure. However, they often introduce noise or uncertainty in the labelset, that leads to decreased performance of supervised deep learning methods. On the other hand, weak supervision methods remain robust during noisy labelsets or can be effective even with low amounts of labelled data. In this paper we evaluate the effectiveness of a representation learning method that uses external categorical evidence called "Evidence Transfer", against low amount of corresponding evidence termed as incomplete evidence. Evidence transfer is a robust solution against external unknown categorical evidence that can introduce noise or uncertainty. In our experimental evaluation, evidence transfer proves to be effective and robust against different levels of incompleteness, for two types of incomplete evidence.


title: Building general Langevin models from discrete data sets

abstract: Many living and complex systems exhibit second order emergent dynamics. Limited experimental access to the configurational degrees of freedom results in data that appears to be generated by a non-Markovian process. This poses a challenge in the quantitative reconstruction of the model from experimental data, even in the simple case of equilibrium Langevin dynamics of Hamiltonian systems. We develop a novel Bayesian inference approach to learn the parameters of such stochastic effective models from discrete finite length trajectories. We first discuss the failure of naive inference approaches based on the estimation of derivatives through finite differences, regardless of the time resolution and the length of the sampled trajectories. We then derive, adopting higher order discretization schemes, maximum likelihood estimators for the model parameters that provide excellent results even with moderately long trajectories. We apply our method to second order models of collective motion and show that our results also hold in the presence of interactions.


title: Pooled variable scaling for cluster analysis

abstract: We propose a new approach for scaling prior to cluster analysis based on the concept of pooled variance. Unlike available scaling procedures such as the standard deviation and the range, our proposed scale avoids dampening the beneficial effect of informative clustering variables. We confirm through an extensive simulation study and applications to well known real data examples that the proposed scaling method is safe and generally useful. Finally, we use our approach to cluster a high dimensional genomic dataset consisting of gene expression data for several specimens of breast cancer cells tissue.


title: Active Learning for Segmentation Based on Bayesian Sample Queries

abstract: Segmentation of anatomical structures is a fundamental image analysis task for many applications in the medical field. Deep learning methods have been shown to perform well, but for this purpose large numbers of manual annotations are needed in the first place, which necessitate prohibitive levels of resources that are often unavailable. In an active learning framework of selecting informed samples for manual labeling, expert clinician time for manual annotation can be optimally utilized, enabling the establishment of large labeled datasets for machine learning. In this paper, we propose a novel method that combines representativeness with uncertainty in order to estimate ideal samples to be annotated, iteratively from a given dataset. Our novel representativeness metric is based on Bayesian sampling, by using information-maximizing autoencoders. We conduct experiments on a shoulder magnetic resonance imaging (MRI) dataset for the segmentation of four musculoskeletal tissue classes. Quantitative results show that the annotation of representative samples selected by our proposed querying method yields an improved segmentation performance at each active learning iteration, compared to a baseline method that also employs uncertainty and representativeness metrics. For instance, with only 10% of the dataset annotated, our method reaches within 5% of Dice score expected from the upper bound scenario of all the dataset given as annotated (an impractical scenario due to resource constraints), and this gap drops down to a mere 2% when less than a fifth of the dataset samples are annotated. Such active learning approach to selecting samples to annotate enables an optimal use of the expert clinician time, being often the bottleneck in realizing machine learning solutions in medicine.


title: A unified framework for X-shaped radio galaxies

abstract: We propose a radically different picture for X-shaped radio galaxies compared to existing models as transition objects between cold mode accreting low spinning retrograde black holes and low spinning prograde black holes. The model explains their smaller average black hole masses, their general aversion for cluster compared to isolated environments, the negligible difference in spectral index between primary and secondary jets despite a time difference in their formation, their absence among the most powerful radio quasars and radio galaxies, and their connection to the elusive FRI quasar class, among others. The key to their formation is cold gas accretion onto Schwarzschild black holes.


title: Improved success probability with greater circuit depth for the quantum approximate optimization algorithm

abstract: Present-day, noisy, small or intermediate-scale quantum processors---although far from fault-tolerant---support the execution of heuristic quantum algorithms, which might enable a quantum advantage, for example, when applied to combinatorial optimization problems. On small-scale quantum processors, validations of such algorithms serve as important technology demonstrators. We implement the quantum approximate optimization algorithm (QAOA) on our hardware platform, consisting of two superconducting transmon qubits and one parametrically modulated coupler. We solve small instances of the NP-complete exact-cover problem, with 96.6% success probability, by iterating the algorithm up to level two.


title: Discussion of "Unbiased Markov chain Monte Carlo with couplings" by Pierre E. Jacob, John O'Leary and Yves F. Atchad\'e

abstract: This is a contribution for the discussion on "Unbiased Markov chain Monte Carlo with couplings" by Pierre E. Jacob, John O'Leary and Yves F. Atchad\'e to appear in the Journal of the Royal Statistical Society Series B.


title: Approximate Maximum Matching in Random Streams

abstract: In this paper, we study the problem of finding a maximum matching in the semi-streaming model when edges arrive in a random order. In the semi-streaming model, an algorithm receives a stream of edges and it is allowed to have a memory of $\tilde{O}(n)$ where $n$ is the number of vertices in the graph. A recent inspiring work by Assadi et al. shows that there exists a streaming algorithm with the approximation ratio of $\frac{2}{3}$ that uses $\tilde{O}(n^{1.5})$ memory. However, the memory of their algorithm is much larger than the memory constraint of the semi-streaming algorithms. In this work, we further investigate this problem in the semi-streaming model, and we present simple algorithms for approximating maximum matching in the semi-streaming model. Our main results are as follows. We show that there exists a single-pass deterministic semi-streaming algorithm that finds a $\frac{3}{5} (= 0.6)$ approximation of the maximum matching in bipartite graphs using $\tilde{O}(n)$ memory. This result significantly outperforms the state-of-the-art result of Konrad that finds a $0.539$ approximation of the maximum matching using $\tilde{O}(n)$ memory. By giving a black-box reduction from finding a matching in general graphs to finding a matching in bipartite graphs, we show there exists a single-pass deterministic semi-streaming algorithm that finds a $\frac{6}{11} (\approx 0.545)$ approximation of the maximum matching in general graphs, improving upon the state-of-art result $0.506$ approximation by Gamlath et al.


title: The Statistics of Extended Debris Disks Measured with Gaia and Planck

abstract: Thermal emission from debris disks around stars has been measured using targeted and resolved observations. We present an alternative, likelihood-based approach in which temperature maps from the Planck CMB survey at 857 and 545 GHz are analyzed in conjunction with stellar positions from Gaia to estimate the fraction of stars hosting disks and the thermal emission from the disks. The debris disks are not resolved (or even necessarily detected individually) but their statistical properties and the correlations with stellar properties are measured for several thousand stars. We compare our findings with higher sensitivity surveys of smaller samples of stars. For dimmer stars, in particular K and M-dwarfs, we find about 10 percent of stars within 80 pc have emission consistent with debris disks. We also report on 80 candidate disks, the majority of which are not previously identified. We have previously constrained the properties of Exo-Oort clouds using Planck data -- with future CMB surveys both components can be measured for different stellar types, providing a new avenue to study the outer parts of planetary systems.

categories: astro-ph.EP astro-ph.CO astro-ph.GA astro-ph.SR


title: Applying the Quantum Approximate Optimization Algorithm to the Tail Assignment Problem

abstract: Airlines today are faced with a number of large scale scheduling problems. One such problem is the tail assignment problem, which is the task of assigning individual aircraft to a given set of flights, minimizing the overall cost. Each aircraft is identified by the registration number on its tail fin. In this article, we simulate the Quantum Approximate Optimization Algorithm (QAOA) applied to instances of this problem derived from real world data. The QAOA is a variational hybrid quantum-classical algorithm recently introduced and likely to run on near-term quantum devices. The instances are reduced to fit on quantum devices with 8, 15 and 25 qubits. The reduction procedure leaves only one feasible solution per instance, which allows us to map the tail assignment problem onto the Exact Cover problem. We find that repeated runs of the QAOA identify the feasible solution with close to unit probability for all instances. Furthermore, we observe patterns in the variational parameters such that an interpolation strategy can be employed which significantly simplifies the classical optimization part of the QAOA. Finally, we empirically find a relation between the connectivity of the problem graph and the single-shot success probability of the algorithm.


title: End point of the first-order phase transition of QCD in the heavy quark region by reweighting from quenched QCD

abstract: We study the end point of the first-order deconfinement phase transition in two and 2+1 flavor QCD in the heavy quark region of the quark mass parameter space. We determine the location of critical point at which the first-order deconfinement phase transition changes to crossover, and calculate the pseudo-scalar meson mass at the critical point. Performing quenched QCD simulations on lattices with the temporal extents Nt=6 and 8, the effects of heavy quarks are determined using the reweighting method. We adopt the hopping parameter expansion to evaluate the quark determinants in the reweighting factor. We estimate the truncation error of the hopping parameter expansion by comparing the results of leading and next-to-leading order calculations, and study the lattice spacing dependence as well as the spatial volume dependence of the result for the critical point. The overlap problem of the reweighting method is also examined. Our results for Nt=4 and 6 suggest that the critical quark mass decreases as the lattice spacing decreases and increases as the spatial volume increases.


title: Axial anomaly generation by domain wall motion in Weyl semimetals

abstract: A space-time dependent node separation in Weyl semimetals acts as an axial vector field. Coupled with domain wall motion in magnetic Weyl semimetals, this induces axial electric and magnetic fields localized at the domain wall. We show how these fields can activate the axial (chiral) anomaly and provide a direct experimental signature of it. Specifically, a domain wall provides a spatially dependent Weyl node separation and an axial magnetic field $\textbf{B}_5$, and domain wall movement, driven by an external magnetic field, gives the Weyl node separation a time dependence, inducing an axial electric field $\textbf{E}_5$. At magnetic fields beyond the Walker breakdown, $\textbf{E}_5\cdot\textbf{B}_5$ becomes nonzero and activates the axial anomaly that induces a finite axial charge density -- imbalance in the number of left- and right-handed fermions -- moving with the domain wall. This axial density, in turn, produces, via the chiral magnetic effect, an oscillating current flowing along the domain wall plane, resulting in a characteristic radiation of electromagnetic waves emanating from the domain wall. A detection of this radiation would constitute a direct measurement of the axial anomaly induced by axial electromagnetic fields.


title: Broad spectral lines in AGNs and supermassive black hole mass measurements

abstract: The mass measurement of supermassive black holes (SMBHs) is a very complex task. Between several methods for SMBH mass measurements, some of them use the spectral lines, which indicate the motion of the emitting/absorbing material around an SMBH. Mostly there is assumption of virialization of line emitting gas in the region close to SMBHs that is used for their mass measurements. In this paper we will give an overview of methods for the SMBH mass measurements using broad emission spectral lines observed in Type 1 AGNs. First we give the basic idea to use the parameters of broad lines to SMBH mass measurements. After that we give an overview of broad lines from X-ray (Fe k{\alpha}) to the IR (Pashen and Brecket lines) which have been used for SMBH mass estimates. Additionally, we describe and discuss a new method for SMBH mass measurements using the polarization in the broad lines emitted from Type 1 AGNs.


title: Rapid Whole-Heart CMR with Single Volume Super-resolution

abstract: Background: Three-dimensional, whole heart, balanced steady state free precession (WH-bSSFP) sequences provide delineation of intra-cardiac and vascular anatomy. However, they have long acquisition times. Here, we propose significant speed ups using a deep learning single volume super resolution reconstruction, to recover high resolution features from rapidly acquired low resolution WH-bSSFP images. Methods: A 3D residual U-Net was trained using synthetic data, created from a library of high-resolution WH-bSSFP images by simulating 0.5 slice resolution and 0.5 phase resolution. The trained network was validated with synthetic test data, as well as prospective low-resolution data. Results: Synthetic low-resolution data had significantly better image quality after super-resolution reconstruction. Qualitative image scores showed super-resolved images had better edge sharpness, fewer residual artefacts and less image distortion than low-resolution images, with similar scores to high-resolution data. Quantitative image scores showed super-resolved images had significantly better edge sharpness than low-resolution or high-resolution images, with significantly better signal-to-noise ratio than high-resolution data. Vessel diameters measurements showed over-estimation in the low-resolution measurements, compared to the high-resolution data. No significant differences and no bias was found in the super-resolution measurements. Conclusion: This paper demonstrates the potential of using a residual U-Net for super-resolution reconstruction of rapidly acquired low-resolution whole heart bSSFP data within a clinical setting. The resulting network can be applied very quickly, making these techniques particularly appealing within busy clinical workflow. Thus, we believe that this technique may help speed up whole heart CMR in clinical practice.


title: Independent Study and Spectral Classification of a Sample of Poorly Studied High Proper Motion M-dwarf Candidate Stars

abstract: We report an independent spectral classification of a sample of poorly studied M-dwarf candidate stars observed with the OSIRIS instrument at GTC. Our project was carried out as an independent test of the spectral classification. It is crucial for the studies of extrasolar planets orbiting M-dwarfs, since properties of the host star are directly related to understanding the planet properties and possible habitability. Understanding of the statistical properties of the dwarf stars is also crucial for the Simple Stellar Population models that play a major role in the modern astrophysics. Halpha emission was detected in 33% of the sample with evidence of Halpha variability in one object.


title: Sheared Amorphous Packings Display Two Separate Particle Transport Mechanisms

abstract: Shearing granular materials induces non-affine displacements. Such non-affine displacements have been studied extensively, and are known to correlate with plasticity and other mechanical features of amorphous packings. A well known example is shear transformation zones as captured by the local deviation from affine deformation, $D^2_{min}$, and their relevance to failure and stress fluctuations. We analyze sheared frictional athermal disk packings and show that there exists at least one additional mesoscopic transport mechanism that superimposes itself on top of local diffusive motion. We evidence this second transport mechanism in a homogeneous system via a diffusion tensor analysis and show that the trace of the diffusion tensor equals the classic $D^2_{min}$ when this second mesoscopic transport is corrected for. The new transport mechanism is consistently observed over a wide range of volume fractions and even for particles with different friction coefficients and is consistently observed also upon shear reversal, hinting at its relevance for memory effects.


title: Ab initio modeling of the energy landscape for screw dislocations in body-centered cubic high-entropy alloys

abstract: In traditional body-centered cubic (bcc) metals, the core properties of screw dislocations play a critical role in plastic deformation at low temperatures. Recently, much attention has been focused on refractory high-entropy alloys (RHEAs), which also possess bcc crystal structures. However, unlike face-centered cubic high-entropy alloys (HEAs), there have been far fewer investigations on bcc HEAs, specifically on the possible effects of chemical short-range order (SRO) in these multiple principal element alloys on dislocation mobility. Here, using density functional theory, we investigate the distribution of dislocation core properties in MoNbTaW RHEAs alloys, and how they are influenced by SRO. The average values of the core energies in the RHEA are found to be larger than those in the corresponding pure constituent bcc metals, and are relatively insensitive to the degree of SRO. However, the presence of SRO is shown to have a large effect on narrowing the distribution of dislocation core energies and decreasing the spatial heterogeneity of dislocation core energies in the RHEA. It is argued that the consequences for the mechanical behavior of HEAs is a change in the energy landscape of the dislocations which would likely heterogeneously inhibit their motion.


title: Resolvent conditions and growth of powers of operators

abstract: Following Berm\'udez et al. (ArXiv: 1706.03638v1), we study the rate of growth of the norms of the powers of a linear operator, under various resolvent conditions or Ces\`aro boundedness assumptions. We show that $T$ is power-bounded if (and only if) both $T$ and $T^*$ are absolutely Ces\`aro bounded. In Hilbert spaces, we prove that if $T$ satisfies the Kreiss condition, $\|T^n\|=O(n/\sqrt {\log n})$; if $T$ is absolutely Ces\`aro bounded, $\|T^n\|=O(n^{1/2 -\varepsilon})$ for some $\varepsilon >0$ (which depends on $T$); if $T$ is strongly Kreiss bounded, then $\|T^n\|=O((\log n)^\kappa)$ for some $\kappa >0$. We show that a Kreiss bounded operator on a reflexive space is Abel ergodic, and its Ces\`aro means of order $\alpha$ converge strongly when $\alpha >1$.


title: Direct and Indirect Effects -- An Information Theoretic Perspective

abstract: Information theoretic (IT) approaches to quantifying causal influences have experienced some popularity in the literature, in both theoretical and applied (e.g. neuroscience and climate science) domains. While these causal measures are desirable in that they are model agnostic and can capture non-linear interactions, they are fundamentally different from common statistical notions of causal influence in that they (1) compare distributions over the effect rather than values of the effect and (2) are defined with respect to random variables representing a cause rather than specific values of a cause. We here present IT measures of direct, indirect, and total causal effects. The proposed measures are unlike existing IT techniques in that they enable measuring causal effects that are defined with respect to specific values of a cause while still offering the flexibility and general applicability of IT techniques. We provide an identifiability result and demonstrate application of the proposed measures in estimating the causal effect of the El Ni\~no-Southern Oscillation on temperature anomalies in the North American Pacific Northwest.


title: AdS$_3$ solutions with $\mathcal{N}=(3,0)$ from S$^3\times$S$^3$ fibrations

abstract: We study warped AdS$_3$ solutions in massive IIA supergravity preserving $\mathcal{N}=(3,0)$ and $\mathcal{N}=(1,0)$ supersymmetry. We consider solutions whose internal spaces decompose as an S$^3\times$S$^3$ fibration and a interval over which the rest of the solution is foliated. We present necessary and sufficient conditions for these solutions to exist, in terms of systems of ordinary differential equations and find several new analytic and numerical examples with internal spaces bounded between D-brane and O-plane behaviors.


title: A Survey of Nuclear Pasta in the Intermediate Density Regime: Structure Functions for Neutrino Scattering

abstract: Background: Nuclear pasta matter, emerging due to the competition between the long-range Coulomb force and the short-range strong force, is believed to be present in astrophysical scenarios, such as neutron stars and core-collapse supernovae. Its structure can have a high impact e.g. on neutrino transport or the tidal deformability of neutron stars. Purpose: We investigate the impact of nuclear pasta on neutrino interactions and compare the results to uniform matter. Method: We calculate the elastic and inelastic static structure factors for nuclear pasta matter using density functional theory (DFT), which contain the main nuclear input for neutrino scattering. Results: Each pasta structure leaves a unique imprint in the elastic structure factor and it is largely enhanced. The inelastic structure factors are very similar for all configurations. Conclusion: Nuclear pasta has a noticeable impact on neutrino neutral-current scattering opacities. While for inelastic reactions the cross section is reduced, the elastic coherent scattering increases dramatically. The effect can be of importance for the cooling of neutron stars as well as for core-collapse supernova models.


title: A structure theorem for fundamental solutions of analytic multipliers in $\mathbb{R}^n$

abstract: Using a version of Hironaka's resolution of singularities for real-analytic functions, any elliptic multiplier $\mathrm{Op}(p)$ of order $d>0$, real-analytic near $p^{-1}(0)$, has a fundamental solution $\mu_0$. We give an integral representation of $\mu_0$ in terms of the resolutions supplied by Hironaka's theorem. This $\mu_0$ is weakly approximated in $H^t_{\mathrm{loc}}(\mathbb{R}^n)$ for $t<d-\frac{n}{2}$ by a sequence from a Paley-Wiener space. In special cases of global symmetry, the obtained integral representation can be made fully explicit, and we use this to compute fundamental solutions for two non-polynomial symbols.


title: Symmetric cluster expansions with tensor networks

abstract: Cluster expansions for the exponential of local operators are constructed using tensor networks. In contrast to other approaches, the cluster expansion does not break any spatial or internal symmetries and exhibits a very favourable prefactor to the error scaling versus bond dimension. This is illustrated by time evolving a matrix product state using very large time steps, and by constructing a novel robust algorithm for finding ground states of 2-dimensional Hamiltonians using projected entangled pair states as fixed points of 2-dimensional transfer matrices.


title: Observation of magnetic solitons in two-component Bose-Einstein condensates

abstract: We experimentally investigate the dynamics of spin solitary waves (magnetic solitons) in a harmonically trapped, binary superfluid mixture. We measure the in-situ density of each pseudospin component and their relative local phase via an interferometric technique we developed, and as such, fully characterise the magnetic solitons while they undergo oscillatory motion in the trap. Magnetic solitons exhibit non-dispersive, dissipationless long-time dynamics. By imprinting multiple magnetic solitons in our ultracold gas sample, we engineer binary collisions between solitons of either same or opposite magnetisation and map out their trajectories.


title: Tag-less Back-Translation

abstract: An effective method to generate a large number of parallel sentences for training improved neural machine translation (NMT) systems is the use of the back-translations of the target-side monolingual data. The standard back-translation method has been shown to be unable to efficiently utilize the available huge amount of existing monolingual data because of the inability of translation models to differentiate between the authentic and synthetic parallel data during training. Tagging, or using gates, has been used to enable translation models to distinguish between synthetic and authentic data, improving standard back-translation and also enabling the use of iterative back-translation on language pairs that underperformed using standard back-translation. In this work, we approach back-translation as a domain adaptation problem, eliminating the need for explicit tagging. In the approach -- \emph{tag-less back-translation} -- the synthetic and authentic parallel data are treated as out-of-domain and in-domain data respectively and, through pre-training and fine-tuning, the translation model is shown to be able to learn more efficiently from them during training. Experimental results have shown that the approach outperforms the standard and tagged back-translation approaches on low resource English-Vietnamese and English-German neural machine translation.


title: Bringing Belief Base Change into Dynamic Epistemic Logic

abstract: AGM's belief revision is one of the main paradigms in the study of belief change operations. In this context, belief bases (prioritised bases) have been primarily used to specify the agent's belief state. While the connection of iterated AGM-like operations and their encoding in dynamic epistemic logics have been studied before, few works considered how well-known postulates from iterated belief revision theory can be characterised by means of belief bases and their counterpart in dynamic epistemic logic. Particularly, it has been shown that some postulates can be characterised through transformations in priority graphs, while others may not be represented that way. This work investigates changes in the semantics of Dynamic Preference Logic that give rise to an appropriate syntactic representation for its models that allow us to represent and reason about iterated belief base change in this logic.


title: Universal $\mathbb{F}_q-$Family of $v-$adic Multiple Zeta Values over Function Fields

abstract: This paper aims to study the $\mathbb{F}_q-$linear relations between interpolated $v-$adic multiple zeta values over function fields. We proved a universal family of linear relations of interpolated $v-$adic MZVs, which is conjectured to generate all linear relations over $\mathbb{F}_q$. At the end of the paper, we also show that these relations can be generalized to all multiple harmonic type sums.


title: Memgames

abstract: In this article, we study the structure, and in particular the Grundy values, of a family of games known as memgames.


title: Theta divisors whose Gauss map has a fiber of positive dimension

abstract: We construct families of principally polarized abelian varieties whose theta divisor is irreducible and contains an abelian subvariety. These families are used to construct examples when the Gauss map of the theta divisor is only generically finite and not finite. That is, the Gauss map in these cases has at least one positive-dimensional fiber. We also obtain lower-bounds on the dimension of Andreotti-Mayer loci.


title: Non-Orthogonal eMBB-URLLC Radio Access for Cloud Radio Access Networks with Analog Fronthauling

abstract: This paper considers the coexistence of Ultra Reliable Low Latency Communications (URLLC) and enhanced Mobile BroadBand (eMBB) services in the uplink of Cloud Radio Access Network (C-RAN) architecture based on the relaying of radio signals over analog fronthaul links. While Orthogonal Multiple Access (OMA) to the radio resources enables the isolation and the separate design of different 5G services, Non-Orthogonal Multiple Access (NOMA) can enhance the system performance by sharing wireless and fronthaul resources. This paper provides an information-theoretic perspective in the performance of URLLC and eMBB traffic under both OMA and NOMA. The analysis focuses on standard cellular models with additive Gaussian noise links and a finite inter-cell interference span, and it accounts for different decoding strategies such as puncturing, Treating Interference as Noise (TIN) and Successive Interference Cancellation (SIC). Numerical results demonstrate that, for the considered analog fronthauling C-RAN architecture, NOMA achieves higher eMBB rates with respect to OMA, while guaranteeing reliable low-rate URLLC communication with minimal access latency. Moreover, NOMA under SIC is seen to achieve the best performance, while, unlike the case with digital capacity-constrained fronthaul links, TIN always outperforms puncturing.


title: Lieb-Schultz-Mattis type constraints on Fractonic Matter

abstract: The Lieb-Schultz-Mattis (LSM) theorem and its descendants impose strong constraints on the low-energy behavior of interacting quantum systems. In this paper, we formulate LSM-type constraints for lattice translation invariant systems with generalized U(1) symmetries which have recently appeared in the context of fracton phases: U(1) polynomial shift and subsystem symmetries. Starting with a generic interacting system with conserved dipole moment, we examine the conditions under which it supports a symmetric, gapped, and non-degenerate ground state, which we find requires that both the filling fraction and the bulk charge polarization take integer-values. Similar constraints are derived for systems with higher moment conservation laws or subsystem symmetries, in addition to lower bounds on the ground state degeneracy when certain conditions are violated. Finally, we discuss the mapping between LSM-type constraints for subsystem symmetries and the anomalous symmetry action at boundaries of subsystem symmetric topological (SSPT) states.


title: Viewing Computer Science through Citation Analysis; Salton and Bergmark Redux

abstract: Computer science has experienced dramatic growth and diversification over the last twenty years. Towards a current understanding of the structure of this discipline, we analyze a cohort of the computer science literature using the DBLP database. For insight on the features of this cohort and the relationship within its components, we constructed article level clusters based on either direct citations or co-citations, and reconciled them to major and minor subject categories in the Scopus All Science Journal Classification (ASJC). We described complementary insights from clustering by direct citation and co-citation, and both point to the increase in computer science publications and their scope. Our analysis shows cross-category clusters, some that interact with external fields, such as the biological sciences, while others remain inward looking.


title: Acoustic versus electromagnetic field theory: scalar, vector, spinor representations and the emergence of acoustic spin

abstract: We construct a novel Lagrangian representation of acoustic field theory that describes the local vector properties of longitudinal (curl-free) acoustic fields. In particular, this approach accounts for the recently-discovered nonzero spin angular momentum density in inhomogeneous sound fields in fluids or gases. The traditional acoustic Lagrangian representation with a ${\it scalar}$ potential is unable to describe such vector properties of acoustic fields adequately, which are however observable via local radiation forces and torques on small probe particles. By introducing a displacement ${\it vector}$ potential analogous to the electromagnetic vector potential, we derive the appropriate canonical momentum and spin densities as conserved Noether currents. The results are consistent with recent theoretical analyses and experiments. Furthermore, by an analogy with dual-symmetric electromagnetic field theory that combines electric- and magnetic-potential representations, we put forward an acoustic ${\it spinor}$ representation combining the scalar and vector representations. This approach also includes naturally coupling to sources. The strong analogies between electromagnetism and acoustics suggest further productive inquiry, particularly regarding the nature of the apparent spacetime symmetries inherent to acoustic fields.


title: Modeling Hessian-vector products in nonlinear optimization: New Hessian-free methods

abstract: In this paper, we suggest two ways of calculating interpolation models for unconstrained smooth nonlinear optimization when Hessian-vector products are available. The main idea is to interpolate the objective function using a quadratic on a set of points around the current one and concurrently using the curvature information from products of the Hessian times appropriate vectors, possibly defined by the interpolating points. These enriched interpolating conditions form then an affine space of model Hessians or model Newton directions, from which a particular one can be computed once an equilibrium or least secant principle is defined. A first approach consists of recovering the Hessian matrix satisfying the enriched interpolating conditions, from which then a Newton direction model can be computed. In a second approach we pose the recovery problem directly in the Newton direction. These techniques can lead to a significant reduction in the overall number of Hessian-vector products when compared to the inexact or truncated Newton method, although simple implementations may pay a cost in linear algebra or number of function evaluations.


title: BNS Invariants and Algebraic Fibrations of Group Extensions

abstract: Let $G$ be a finitely generated group that can be written as an extension \[ 1 \longrightarrow K \stackrel{i}{\longrightarrow} G \stackrel{f}{\longrightarrow} \Gamma \longrightarrow 1 \] where $K$ is a finitely generated group. By a study of the BNS invariants we prove that if $b_1(G) > b_1(\Gamma) > 0$, then $G$ algebraically fibers, i.e. admits an epimorphism to $\Bbb{Z}$ with finitely generated kernel. An interesting case of this occurrence is when $G$ is the fundamental group of a surface bundle over a surface $F \hookrightarrow X \rightarrow B$ with Albanese dimension $a(X) = 2$. As an application, we show that if $X$ has virtual Albanese dimension $va(X) = 2$ and base and fiber have genus greater that $1$, $G$ is noncoherent. This answers for a broad class of bundles a question of J. Hillman.


title: Re-Identification and Growth Detection of Pulmonary Nodules without Image Registration Using 3D Siamese Neural Networks

abstract: Lung cancer follow-up is a complex, error prone, and time consuming task for clinical radiologists. Several lung CT scan images taken at different time points of a given patient need to be individually inspected, looking for possible cancerogenous nodules. Radiologists mainly focus their attention in nodule size, density, and growth to assess the existence of malignancy. In this study, we present a novel method based on a 3D siamese neural network, for the re-identification of nodules in a pair of CT scans of the same patient without the need for image registration. The network was integrated into a two-stage automatic pipeline to detect, match, and predict nodule growth given pairs of CT scans. Results on an independent test set reported a nodule detection sensitivity of 94.7%, an accuracy for temporal nodule matching of 88.8%, and a sensitivity of 92.0% with a precision of 88.4% for nodule growth detection.


title: Building and Testing Yield Curve Generators for P&C Insurance

abstract: Interest-rate risk is a key factor for property-casualty insurer capital. P&C companies tend to be highly leveraged, with bond holdings much greater than capital. For GAAP capital, bonds are marked to market but liabilities are not, so shifts in the yield curve can have a significant impact on capital. Yield-curve scenario generators are one approach to quantifying this risk. They produce many future simulated evolutions of the yield curve, which can be used to quantify the probabilities of bond-value changes that would result from various maturity-mix strategies. Some of these generators are provided as black-box models where the user gets only the projected scenarios. One focus of this paper is to provide methods for testing generated scenarios from such models by comparing to known distributional properties of yield curves. P&C insurers hold bonds to maturity and manage cash-flow risk by matching asset and liability flows. Derivative pricing and stochastic volatility are of little concern over the relevant time frames. This requires different models and model testing than what is common in the broader financial markets. To complicate things further, interest rates for the last decade have not been following the patterns established in the sixty years following WWII. We are now coming out of the period of very low rates, yet are still not returning to what had been thought of as normal before that. Modeling and model testing are in an evolving state while new patterns emerge. Our analysis starts with a review of the literature on interest-rate model testing, with a P&C focus, and an update of the tests for current market behavior. We then discuss models, and use them to illustrate the fitting and testing methods. The testing discussion does not require the model-building section.


title: Ramsey theory and topological dynamics for first order theories

abstract: We investigate interactions between Ramsey theory, topological dynamics, and model theory. We introduce various Ramsey-like properties for first order theories and characterize them in terms of the appropriate dynamical properties of the theories in question (such as [extreme] amenability of a theory or some properties of the associated Ellis semigroups). Then we relate them to profiniteness and triviality of the Ellis groups of first order theories. In particular, we find various criteria for [pro]finiteness and for triviality of the Ellis group of a given theory from which we obtain wide classes of examples of theories with [pro]finite or trivial Ellis groups. As an initial motivation, we note that profiniteness of the Ellis group of a theory implies that the Kim-Pillay Galois group of this theory is also profinite, which in turn is equivalent to the equality of the Shelah and Kim-Pillay strong types. We also find several concrete examples illustrating the lack of implications between some fundamental properties. In the appendix, we give a full computation of the Ellis group of the theory of the random hypergraph with one binary and one 4-ary relation. This example shows that the assumption of NIP in the version of Newelski's conjecture for amenable theories (proved in [15]) cannot be dropped.


title: Emotion Recognition Using Wearables: A Systematic Literature Review Work in progress

abstract: Wearables like smartwatches or wrist bands equipped with pervasive sensors enable us to monitor our physiological signals. In this study, we address the question whether they can help us to recognize our emotions in our everyday life for ubiquitous computing. Using the systematic literature review, we identified crucial research steps and discussed the main limitations and problems in the domain.


title: Improved Central Limit Theorem and bootstrap approximations in high dimensions

abstract: This paper deals with the Gaussian and bootstrap approximations to the distribution of the max statistic in high dimensions. This statistic takes the form of the maximum over components of the sum of independent random vectors and its distribution plays a key role in many high-dimensional econometric problems. Using a novel iterative randomized Lindeberg method, the paper derives new bounds for the distributional approximation errors. These new bounds substantially improve upon existing ones and simultaneously allow for a larger class of bootstrap methods.


title: Analysis of a Sinclair-type domain decomposition solver for atomistic/continuum coupling

abstract: The "flexible boundary condition" method, introduced by Sinclair and coworkers in the 1970s, remains among the most popular methods for simulating isolated two-dimensional crystalline defects, embedded in an effectively infinite atomistic domain. In essence, the method can be characterized as a domain decomposition method which iterates between a local anharmonic and a global harmonic problem, where the latter is solved by means of the lattice Green function of the ideal crystal. This local/global splitting gives rise to tremendously improved convergence rates over related alternating Schwarz methods. In a previous publication (Hodapp et al., 2019, Comput. Methods in Appl. Mech. Eng. 348), we have shown that this method also applies to large-scale three-dimensional problems, possibly involving hundreds of thousands of atoms, using fast summation techniques exploiting the low-rank nature of the asymptotic lattice Green function. Here, we generalize the Sinclair method to bounded domains and develop an implementation using a discrete boundary element method to correct the infinite solution with respect to a prescribed far-field condition, thus preserving the advantage of the original method of not requiring a global spatial discretization. Moreover, we present a detailed convergence analysis and show for a one-dimensional problem that the method is unconditionally stable under physically motivated assumptions. To further improve the convergence behavior, we develop an acceleration technique based on a relaxation of the transmission conditions between the two subproblems. Numerical examples for linear and nonlinear problems are presented to validate the proposed methodology.


title: The CoCo-Beholder: Enabling Comprehensive Evaluation of Congestion Control Algorithms

abstract: The recent endeavors of the research community to unite efforts on the design and evaluation of congestion control algorithms have created a growing collection of congestion control schemes called Pantheon. However, the virtual network emulator that comes with the collection has very limited capabilities: it can run flows of only one scheme at once, and the flows cannot have individual network settings, as the topology is point-to-point. This thesis addresses those limitations and presents CoCo-Beholder, a human-friendly emulator providing the popular dumbbell topology of any size, each link of which may have individual rate, delay, and queue size. The central link of the topology may also have a variable delay with optional jitter. Flows of different schemes may run between the halves of the topology at once, and for each flow, the direction and starting time can be chosen. CoCo-Beholder's reliability is ensured by testing schemes in the dumbbell topology of size one and comparing the results against those of Pantheon emulator. With CoCo-Beholder, the thesis successfully reproduces experiments from a recent paper that evaluated the fairness and RTT-fairness of schemes using a real hardware dumbbell testbed. Finally, the thesis explores the behavior of schemes under the square-wave delay using CoCo-Beholder's variable delay feature.


title: Banks-Zaks Cosmology, Inflation, and the Big Bang Singularity

abstract: We consider the thermodynamical behavior of Banks-Zaks theory close to the conformal point in a cosmological setting. Due to the anomalous dimension, the resulting pressure and energy density deviate from that of radiation and result in various interesting cosmological scenarios. Specifically, for a given range of parameters one avoids the cosmological singularity. We provide a full "phase diagram" of possible Universe evolution for the given parameters. For a certain range of parameters, the thermal averaged Banks-Zaks theory alone results in an exponentially contracting universe followed by a non-singular bounce and an exponentially expanding universe, i.e. \textit{Inflation without a Big Bang singularity}, or shortly termed "dS Bounce". The temperature of such a universe is bounded from above and below. The result is a theory violating the classical Null Energy Condition (NEC). Considering the Banks-Zaks theory with an additional perfect fluid, yields an even richer phase diagram that includes the standard Big Bang model, stable single "normal" bounce, dS Bounce and stable cyclic solutions. The bouncing and cyclic solutions are with no singularities, and the violation of the NEC happens only near the bounce. We also provide simple analytical conditions for the existence of these non-singular solutions. Hence, within effective field theory, we have a new alternative non-singular cosmology based on the anomalous dimension of Bank-Zaks theory that may include inflation and without resorting to scalar fields.


title: Inferring Parameters of GW170502: The Loudest Intermediate-mass Black Hole Trigger in LIGO's O1/O2 data

abstract: Gravitational wave (GW) measurements provide the most robust constraints of the mass of astrophysical black holes. Using state-of-the-art GW signal models and a unique parameter estimation technique, we infer the source parameters of the loudest marginal trigger, GW170502, found by LIGO from 2015 to 2017. If this trigger is assumed to be a binary black hole merger, we find it corresponds to a total mass in the source frame of $157^{+55}_{-41}~\rm{M}_\odot$ at redshift $z=1.37^{+0.93}_{-0.64}$. The primary and secondary black hole masses are constrained to $94^{+44}_{-28}~\rm{M}_{\odot}$ and $62^{+30}_{-25}~\rm{M}_{\odot}$ respectively, with 90\% confidence. Across all signal models, we find $\gtrsim 70\%$ probability for the effective spin parameter $\chi_\mathrm{eff}>0.1$. Furthermore, we find that the inclusion of higher-order modes in the analysis narrows the confidence region for the primary black hole mass by 10\%, however, the evidence for these modes in the data remains negligible. The techniques outlined in this study could lead to robust inference of the physical parameters for all intermediate-mass black hole binary candidates $(\gtrsim100~\mathrm{M}_\odot)$ in the current GW network.


title: Designing Radiation Transport Tests: Simulation-Driven Uncertainty-Quantification of the COAX Temperature Diagnostic

abstract: One of the difficulties in developing accurate numerical models of radiation flow in a coupled radiation-hydrodynamics setting is accurately modeling the transmission across a boundary layer. The COAX experiment is a platform design to test this transmission including standard radiograph and flux diagnostics as well as a temperature diagnostic measuring the population of excitation levels and ionization states of a dopant embedded within the target material. Using a broad range of simulations, we study the experimental errors in this temperature diagnostic. We conclude with proposed physics experiments that show features that are much stronger than the experimental errors and provide the means to study transport models.