title: Axiogenesis from $SU(2)_R$ phase transition

abstract: The baryon asymmetry of the universe may be explained by rotations of the QCD axion in field space and baryon number violating processes. We consider the minimal extension of the Standard Model by a non-Abelian gauge interaction, $SU(2)_R$, whose sphaleron process violates baryon number. Assuming that axion dark matter is also created from the axion rotation by the kinetic misalignment mechanism, the mass scale of the $SU(2)_R$ gauge boson is fixed as a function of the QCD axion decay constant, and vise versa. Significant portion of the parameter space has already been excluded by new gauge boson searches, and the high-luminocity LHC will further probe the viable parameter space.


title: The Fornax3D project: Planetary nebulae catalogue and independent distance measurements to Fornax cluster galaxies

abstract: Extragalactic planetary nebulae (PNe) offer a way to determine the distance to their host galaxies thanks to the nearly universal shape of the planetary nebulae luminosity function (PNLF). Accurate PNe distance measurements rely on obtaining well-sampled PNLFs and the number of observed PNe scales with the encompassed stellar mass. This means either disposing of wide-field observations or focusing on the bright central regions of galaxies. In this work we take this second approach and conduct a census of the PNe population in the central regions of galaxies in the Fornax cluster, using VLT/MUSE data for the early-type galaxies observed over the course of the Fornax3D survey. Using such integral-field spectroscopic observations to carefully separate the nebular emission from the stellar continuum, we isolated [OIII] 5007 {\AA} sources of interest, filtered out unresolved impostor sources or kinematic outliers, and present a catalogue of 1350 unique PNe sources across 21 early-type galaxies, which includes their positions, [OIII] 5007 {\AA} line magnitudes, and line-of-sight velocities. Using the PNe catalogued within each galaxy, we present independently derived distance estimates based on the fit to the entire observed PNLF observed while carefully accounting for the PNe detection incompleteness. With these individual measurements, we arrive at an average distance to the Fornax cluster itself of 19.86 $\pm$ 0.32 Mpc ($\mu_{PNLF}$ = 31.49 $\pm$ 0.04 mag). Our PNLF distance measurements agree well with previous distances based on surface brightness fluctuations, finding no significant systematic offsets between the two methods as otherwise reported in previous studies.


title: Multi-wavelength Observations of Sgr A*. I. 2019 July 18

abstract: We present and analyze ALMA submillimeter observations from a multi-wavelength campaign of Sgr A* during 18 July 2019. In addition to the submillimeter, we utilize concurrent mid-IR (Spitzer) and X-ray (Chandra) observations. The submillimeter emission lags less than $\delta t\approx30$ minutes behind the mid-IR data. However, the entire submillimeter flare was not observed, raising the possibility that the time delay is a consequence of incomplete sampling of the light curve. The decay of the submillimeter emission is not consistent with synchrotron cooling. Therefore, we analyze these data adopting an adiabatically expanding synchrotron source that is initially optically thick or thin in the submillimeter, yielding time-delayed or synchronous flaring with the IR, respectively. The time-delayed model is consistent with a plasma blob of radius $0.8~R_{\text{S}}$ (Schwarzschild radius), electron power-law index $p=3.5$ ($N(E)\propto E^{-p}$), equipartition magnetic field of $B_{\text{eq}}\approx90$ Gauss, and expansion velocity $v_{\text{exp}}\approx0.004c$. The simultaneous emission is fit by a plasma blob of radius $2~R_{\text{S}}$, $p=2.5$, $B_{\text{eq}}\approx27$ Gauss, and $v_{\text{exp}}\approx0.014c$. Since the submillimeter time delay is not completely unambiguous, we cannot definitively conclude which model better represents the data. This observation presents the best evidence for a unified flaring mechanism between submillimeter and X-ray wavelengths and places significant constraints on the source size and magnetic field strength. We show that concurrent observations at lower frequencies would be able to determine if the flaring emission is initially optically thick or thin in the submillimeter.


title: Laplace surface dynamics, revisited: satellites, exo-planets and debris with distant, eccentric companions

abstract: To date, studies of $\textit{Laplace Surface}$ dynamics have concerned themselves with test particle orbits of fixed shape and orientation in the combined field of an oblate central body (to which the particle is bound) and a distant, inclined, companion which is captured to quadrupolar order. While amply sufficient for satellites around planets on near-circular orbits, the quadrupolar approximation fails to capture essential dynamical features induced by a wide binary companion (be it a star, a planet or a black hole) on a fairly eccentric orbit. With similar such astronomical settings in mind, we extend the classical Laplace framework to higher multipoles, and map out the backbone of stationary orbits, now complexified by the broken axial symmetry. Eccentric and inclined Laplace equilibria, which had been presaged in systems of large enough mutual inclination, are here delineated over a broad range of mutually inclined perturbations. We recover them for test particles in the field of a hot Jupiter and a wide eccentric stellar binary, highlighting their relevance for the architecture of multi-planet systems in binaries. We then extend and deploy our machinery closer to home, as we consider the secular dynamics of Trans-Neptunian Objects (TNOs) in the presence of a putative ninth planet. We show how generalized Laplace equilibria seed islands for Trans-Neptunian objects to be sheltered around, islands within chaotic seas which we capture via Poincar$\'{e}$ sections, while highlighting a beautiful interplay between Laplace and Kozai-Lidov secular dynamical structures. An eminently classical tale revived for the exo-planetary 21st century!


title: Emergent parametric resonances and time-crystal phases in driven BCS systems

abstract: We study the out-of-equilibrium dynamics of a Bardeen-Cooper-Schrieffer condensate subject to a periodic drive. We demonstrate that the combined effect of drive and interactions results in emerging parametric resonances, analogous to a vertically driving pendulum. In particular, Arnold tongues appear when the driving frequency matches $2\Delta_0/n$, with $n$ a natural number, and $\Delta_0$ the equilibrium gap parameter. Inside the Arnold tongues we find a commensurate time-crystal condensate which retains the $U(1)$ symmetry breaking of the parent superfluid/superconducting phase and shows an additional time-translational symmetry breaking. Outside these tongues, the synchronized collective Higgs mode found in quench protocols is stabilized without the need of a strong perturbation. Our results are directly relevant to cold-atom and condensed-matter systems and do not require very long energy relaxation times to be observed.


title: Classification of Small Triorthogonal Codes

abstract: Triorthogonal codes are a class of quantum error correcting codes used in magic state distillation protocols. We classify all triorthogonal codes with $n+k \le 38$, where $n$ is the number of physical qubits and $k$ is the number of logical qubits of the code. We find $38$ distinguished triorthogonal subspaces and show that every triorthogonal code with $n+k\le 38$ descends from one of these subspaces through elementary operations such as puncturing and deleting qubits. Specifically, we associate each triorthogonal code with a Reed-Muller polynomial of weight $n+k$, and classify the Reed-Muller polynomials of low weight using the results of Kasami, Tokura, and Azumi and an extensive computerized search. In an appendix independent of the main text, we improve a magic state distillation protocol by reducing the time variance due to stochastic Clifford corrections.


title: On the possibility that PbZrO$_3$ not be antiferroelectric

abstract: Lead zirconate (PbZrO$_3$) is considered the prototypical antiferroelectric material with an antipolar ground state. Yet, several experimental and theoretical works hint at a partially polar behaviour in this compound, indicating that the polarization may not be completely compensated. In this work we propose a simple ferrielectric structure for lead zirconate. First-principles calculations reveal this state to be more stable than the commonly accepted antiferroelectric phase at low temperatures, possibly up to room temperature, suggesting that PbZrO$_3$ may not be antiferroelectric at ambient conditions. We discuss the implications of our discovery, how it can be reconciled with experimental observations and how the ferrielectric phase could be obtained in practice.


title: Enhanced Photonic Maxwell's Demon with Correlated Baths

abstract: Maxwell's Demon is at the heart of the interrelation between quantum information processing and thermodynamics. In this thought experiment, a demon generates a temperature gradient between two thermal baths initially at equilibrium by gaining information at the single-particle level and applying classical feed-forward operations, allowing for the extraction of work. Here we implement a photonic version of Maxwell's Demon with active feed-forward in a fibre-based system using ultrafast optical switches. We experimentally show that, if correlations exist between the two thermal baths, the Demon can generate a temperature difference over an order of magnitude larger than without correlations, and so extract more work. Our work demonstrates the great potential of photonic experiments -- which provide a unique degree of control on the system -- to access new regimes in quantum thermodynamics.


title: An Arecibo Search for Fast Radio Transients from M87

abstract: The possible origin of millisecond bursts from the giant elliptical galaxy M87 has been scrutinized since the earliest searches for extragalactic fast radio transients undertaken in the late 1970s. Motivated by rapid technological advancements in recent years, we conducted $\rm \simeq 10~hours$ of L-band ($\rm 1.15-1.75~GHz$) observations of the core of M87 with the Arecibo radio telescope in 2019. Adopting a matched filtering approach, we searched our data for single pulses using trial dispersion measures up to $\rm 5500~pc~cm^{-3}$ and burst durations between $\rm 0.3-123~ms$. We find no evidence of astrophysical bursts in our data above a 7$\sigma$ detection threshold. Our observations thus constrain the burst rate from M87 to $\rm \lesssim 0.1~bursts~hr^{-1}$ above $\rm 1.4~Jy~ms$, the most stringent upper limit obtained to date. Our non-detection of radio bursts is consistent with expectations of giant pulse emission from a Crab-like young neutron star population in M87. However, the dense, strongly magnetized interstellar medium surrounding the central $\sim 10^9 \ M_{\odot}$ supermassive black hole of M87 may potentially harbor magnetars that can emit detectable radio bursts during their flaring states.


title: Closing the window on WIMP Dark Matter

abstract: We study scenarios where Dark Matter is a weakly interacting particle (WIMP) embedded in an ElectroWeak multiplet. In particular, we consider real SU(2) representations with zero hypercharge, that automatically avoid direct detection constraints from tree-level Z-exchange. We compute for the first time all the calculable thermal masses for scalar and fermionic WIMPs, including Sommerfeld enhancement and bound states formation at leading order in gauge boson exchange and emission. WIMP masses of few hundred TeV are shown to be compatible both with s-wave unitarity of the annihilation cross-section, and perturbativity. We also provide theory uncertainties on the masses for all multiplets, which are shown to be significant for large SU(2) multiplets. We then outline a strategy to probe these scenarios at future experiments. Electroweak 3-plets and 5-plets have masses up to about 16 TeV and can efficiently be probed at a high energy muon collider. We study various experimental signatures, such as single and double gauge boson emission with missing energy, and disappearing tracks, and determine the collider energy and luminosity required to probe the thermal Dark Matter masses. Larger multiplets are out of reach of any realistic future collider, but can be tested in future gamma ray telescopes and possibly in large-exposure liquid Xenon experiments.


title: Sneutrino Tribrid Inflation, Metastable Cosmic Strings and Gravitational Waves

abstract: We present a successful realization of sneutrino tribrid inflation model based on a gauged $U(1)_{B-L}$ extension of Minimal Supersymmetric Standard Model (MSSM). A single interaction term involving the $B-L$ Higgs field and the right-handed neutrinos serves multiple purposes. These include the generation of heavy Majorana masses for the right-handed neutrinos to provide an explanation for the tiny neutrino masses via the seesaw mechanism, a realistic scenario for reheating and non-thermal leptogenesis with a reheat temperature as low as $10^6$ GeV, and a successful realization of inflation with right-handed sneutrino as the inflaton. The matter parity which helps avoid rapid proton decay survives as a $Z_{2}$ subgroup of a $U(1)$ $R$-symmetry. Depending on the choice of model parameters yields the following predicted range of the tensor to scalar ratio, $3 \times 10^{-11}\lesssim r\lesssim 7\times 10^{-4}$ ($ 6 \times 10^{-7} \lesssim r \lesssim 0.01 $), and the running of the scalar spectral index, $-0.00022 \lesssim dn_s/d\ln k \lesssim -0.0026$ ($-0.00014 \lesssim dn_s/d\ln k \lesssim 0.005$), along with the $B-L$ breaking scale, $ 3 \times 10^{14}\lesssim M/ \text{GeV}\lesssim 5 \times 10^{15}$ ($ 6 \times 10^{15}\lesssim M/ \text{GeV}\lesssim 2 \times 10^{16}$), calculated at the central value of the scalar spectral index, $n_s =0.966$, reported by Planck 2018. The possibility of realizing metastable cosmic strings in a grand unified theory setup is briefly discussed. The metastable cosmic string network admits string tension values in the range $10^{-8} \lesssim G\mu_s \lesssim 10^{-6}$, and predicts a stochastic gravitational wave background lying within the 2-$\sigma$ bounds of the recent NANOGrav 12.5-yr data.


title: Hilbert space fragmentation in a 2D quantum spin system with subsystem symmetries

abstract: We consider a 2D quantum spin model with ring-exchange interaction that has subsystem symmetries associated to conserved magnetization along rows and columns of a square lattice, which implies the conservation of the global dipole moment. In a certain regime, the model is non-integrable, but violates the eigenstate thermalization hypothesis through an extensive Hilbert space fragmentation, including an exponential number of inert subsectors with trivial dynamics, arising from kinetic constraints. While subsystem symmetries are quite restrictive for the dynamics, we show that they alone cannot account for such a number of inert states, even with infinite-range interactions. We present a procedure for constructing shielding structures that can separate and disentangle dynamically active regions from each other. Notably, subsystem symmetries allow the thickness of the shields to be dependent only on the interaction range rather than on the size of the active regions, unlike in the case of generic dipole-conserving systems.


title: Hodge classes on the moduli space of W(E_6)-covers and the geometry of A_6

abstract: In previous work we showed that the Hurwitz space of W(E_6)-covers of the projective line branched over 24 points dominates via the Prym-Tyurin map the moduli space A_6 of principally polarized abelian 6-folds. Here we determine the 25 Hodge classes on the Hurwitz space of W(E_6)-covers corresponding to the 25 irreducible representations of the Weyl group W(E_6). This result has direct implications to the intersection theory of the toroidal compactification A_6. In the final part of the paper, we present an alternative, elementary proof of our uniformization result on A_6 via Prym-Tyurin varieties of type W(E_6).


title: Measuring binary black hole orbital-plane spin orientations

abstract: Binary black hole spins are among the key observables for gravitational wave astronomy. Among the spin parameters, their orientations within the orbital plane, $\phi_1$, $\phi_2$ and $\Delta \phi=\phi_1-\phi_2$, are critical for understanding the prevalence of the spin-orbit resonances and merger recoils in binary black holes. Unfortunately, these angles are particularly hard to measure using current detectors, LIGO and Virgo. Because the spin directions are not constant for precessing binaries, the traditional approach is to measure the spin components at some reference stage in the waveform evolution, typically the point at which the frequency of the detected signal reaches 20 Hz. However, we find that this is a poor choice for the orbital-plane spin angle measurements. Instead, we propose measuring the spins at a fixed dimensionless time or frequency near the merger. This leads to significantly improved measurements for $\phi_1$ and $\phi_2$ for several gravitational wave events. Furthermore, using numerical relativity injections, we demonstrate that $\Delta \phi$ will also be better measured near the merger for louder signals expected in the future. Finally, we show that numerical relativity surrogate models are key for reliably measuring the orbital-plane spin orientations, even at moderate signal-to-noise ratios like $\sim 30-45$.


title: Hints of spin-orbit resonances in the binary black hole population

abstract: Binary black hole spin measurements from gravitational wave observations can reveal the binary's evolutionary history. In particular, the spin orientations of the component black holes within the orbital plane, $\phi_1$ and $\phi_2$, can be used to identify binaries caught in the so-called spin-orbit resonances. In a companion paper, we demonstrate that $\phi_1$ and $\phi_2$ are best measured near the merger of the two black holes. In this work, we use these spin measurements to provide the first constraints on the full six-dimensional spin distribution of merging binary black holes. In particular, we find that there is a preference for $\Delta \phi = \phi_1 - \phi_2 \sim \pm \pi$ in the population, which can be a signature of spin-orbit resonances. We also find a preference for $\phi_1 \sim -\pi/4$ with respect to the line of separation near merger, which has not been predicted for any astrophysical formation channel. However, the strength of these preferences depends on our prior choices, and we are unable to constrain the widths of the $\phi_1$ and $\Delta \phi$ distributions. Therefore, more observations are necessary to confirm the features we find. Finally, we derive constraints on the distribution of recoil kicks in the population, and use this to estimate the fraction of merger remnants retained by globular and nuclear star clusters. We make our spin and kick population constraints publicly available.


title: Integral decomposition for the solutions of the generalized Cattaneo equation

abstract: We present the integral decomposition for the fundamental solution of the generalized Cattaneo equation with both time derivatives smeared through convoluting them with some memory kernels. For power-law kernels $t^{-\alpha}$, $\alpha\in(0,1]$ this equation becomes the time fractional one governed by the Caputo derivatives which highest order is 2. To invert the solutions from the Fourier-Laplace domain to the space-time domain we use analytic methods based on the Efross theorem and find out that solutions looked for are represented by integral decompositions which tangle the fundamental solution of the standard Cattaneo equation with non-negative and normalizable functions being uniquely dependent on the memory kernels. Furthermore, the use of methodology arising from the theory of complete Bernstein functions allows us to assign such constructed integral decompositions the interpretation of subordination. This fact is preserved in two limit cases built into the generalized Cattaneo equations, i.e., either the diffusion or the wave equations. We point out that applying the Efross theorem enables us to go beyond the standard approach which usually leads to the integral decompositions involving the Gaussian distribution describing the Brownian motion. Our approach clarifies puzzling situation which takes place for the power-law kernels $t^{-\alpha}$ for which the subordination based on the Brownian motion does not work if $\alpha\in(1/2,1]$.


title: Using materials for quasiparticle engineering

abstract: The fundamental excitations in superconductors - Bogoliubov quasiparticles - can be either a resource or a liability in superconducting devices: they are what enables photon detection in microwave kinetic inductance detectors, but they are a source of errors in qubits and electron pumps. To improve operation of the latter devices, ways to mitigate quasiparticle effects have been devised; in particular, combining different materials quasiparticles can be trapped where they do no harm and their generation can be impeded. We review recent developments in these mitigation efforts and discuss open questions.


title: Phylogenetic consensus networks: Computing a consensus of 1-nested phylogenetic networks

abstract: An important and well-studied problem in phylogenetics is to compute a \emph{consensus tree} so as to summarize the common features within a collection of rooted phylogenetic trees, all whose leaf-sets are bijectively labeled by the same set~(X) of species. More recently, however, it has become of interest to find a consensus for a collection of more general, rooted directed acyclic graphs all of whose sink-sets are bijectively labeled by~(X), so called rooted \emph{phylogenetic networks}. These networks are used to analyse the evolution of species that cross with one another, such as plants and viruses. In this paper, we introduce an algorithm for computing a consensus for a collection of so-called 1-\emph{nested} phylogenetic networks. Our approach builds on a previous result by Rosell\'o et al. that describes an encoding for any 1-nested phylogenetic network in terms of a collection of ordered pairs of subsets of (X).More specifically, we characterize those collections of ordered pairs that arise as the encoding of some 1-nested phylogenetic network, and then use this characterization to compute a \emph{consensus network} for a collection of~$t$ 1-nested networks in $O(t|X|^2+|X|^3)$ time. Applying our algorithm to a collection of phylogenetic trees yields the well-known majority rule consensus tree. Our approach leads to several new directions for futurework, and we expect that it should provide a useful new tool to help understand complex evolutionary scenarios.


title: Phases and Quantum Phase Transitions in Anisotropic Antiferromagnetic Kitaev-Heisenberg-$\Gamma$ magnet

abstract: We study the Kitaev-Heisenberg-$\Gamma$ model with antiferromagnetic Kitaev exchanges in the strong anisotropic (toric code) limit to understand the phases and the intervening phase transitions between the gapped $Z_2$ quantum spin liquid and the spin-ordered (in the Heisenberg limit) as well as paramagnetic phases (in the pseudo-dipolar, $\Gamma$, limit). We find that the paramagnetic phase obtained in the large $\Gamma$ limit has no topological entanglement entropy and is proximate to a gapless critical point of a system described by an equal superposition of differently oriented stacked one-dimensional $Z_2\times Z_2$ symmetry protected topological phases. Using a combination of exact diagonalization calculations and field-theoretic analysis we map out the phases and phase transitions to reveal the complete phase diagram as a function of the Heisenberg, the Kitaev, and the pseudo-dipolar interactions. Our work shows a rich plethora of unconventional phases and phase transitions and provides a comprehensive understanding of the physics of anisotropic Kitaev-Heisenberg-$\Gamma$ systems along with our recent paper [Phys. Rev. B 102, 235124 (2020)] where the ferromagnetic Kitaev exchange was studied.


title: Mono2Micro: A Practical and Effective Tool for Decomposing Monolithic Java Applications to Microservices

abstract: In migrating production workloads to cloud, enterprises often face the daunting task of evolving monolithic applications toward a microservice architecture. At IBM, we developed a tool called Mono2Micro to assist with this challenging task. Mono2Micro performs spatio-temporal decomposition, leveraging well-defined business use cases and runtime call relations to create functionally cohesive partitioning of application classes. Our preliminary evaluation of Mono2Micro showed promising results. How well does Mono2Micro perform against other decomposition techniques, and how do practitioners perceive the tool? This paper describes the technical foundations of Mono2Micro and presents results to answer these two questions. To answer the first question, we evaluated Mono2Micro against four existing techniques on a set of open-source and proprietary Java applications and using different metrics to assess the quality of decomposition and tool's efficiency. Our results show that Mono2Micro significantly outperforms state-of-the-art baselines in specific metrics well-defined for the problem domain. To answer the second question, we conducted a survey of twenty-one practitioners in various industry roles who have used Mono2Micro. This study highlights several benefits of the tool, interesting practitioner perceptions, and scope for further improvements. Overall, these results show that Mono2Micro can provide a valuable aid to practitioners in creating functionally cohesive and explainable microservice decompositions.


title: Random Permutations -- A geometric point of view

abstract: We look at geometric limits of large random non-uniform permutations. We mainly consider two theories for limits of permutations: permuton limits, introduced by Hoppen, Kohayakawa, Moreira, Rath, and Sampaio to define a notion of scaling limits for permutations; and Benjamini-Schramm limits, introduced by the author to define a notion of local limits for permutations. The models of random permutations that we consider are mainly constrained models, that is, uniform permutations belonging to a given subset of the set of all permutations. We often identify this subset using pattern-avoidance, focusing on: permutations avoiding a pattern of length three, substitution-closed classes, (almost) square permutations, permutation families encoded by generating trees, and Baxter permutations. We explore some universal phenomena for the models mentioned above. For Benjamini-Schramm limits we explore a concentration phenomenon for the limiting objects. For permuton limits we deepen the study of some known universal permutons, called biased Brownian separable permutons, and we introduce some new ones, called Baxter permuton and skew Brownian permutons. In addition, for (almost) square permutations, we investigate the occurrence of a phase transition for the limiting permutons. On the way, we establish various combinatorial results both for permutations and other related objects. Among others, we give a complete description of the feasible region for consecutive patterns as the cycle polytope of a specific graph; and we find new bijections relating Baxter permutations, bipolar orientations, walks in cones, and a new family of discrete objects called coalescent-walk processes.


title: 3D-StyleGAN: A Style-Based Generative Adversarial Network for Generative Modeling of Three-Dimensional Medical Images

abstract: Image synthesis via Generative Adversarial Networks (GANs) of three-dimensional (3D) medical images has great potential that can be extended to many medical applications, such as, image enhancement and disease progression modeling. However, current GAN technologies for 3D medical image synthesis need to be significantly improved to be readily adapted to real-world medical problems. In this paper, we extend the state-of-the-art StyleGAN2 model, which natively works with two-dimensional images, to enable 3D image synthesis. In addition to the image synthesis, we investigate the controllability and interpretability of the 3D-StyleGAN via style vectors inherited form the original StyleGAN2 that are highly suitable for medical applications: (i) the latent space projection and reconstruction of unseen real images, and (ii) style mixing. We demonstrate the 3D-StyleGAN's performance and feasibility with ~12,000 three-dimensional full brain MR T1 images, although it can be applied to any 3D volumetric images. Furthermore, we explore different configurations of hyperparameters to investigate potential improvement of the image synthesis with larger networks. The codes and pre-trained networks are available online: https://github.com/sh4174/3DStyleGAN.


title: Well-posedness of Bayesian inverse problems for hyperbolic conservation laws

abstract: We study the well-posedness of the Bayesian inverse problem for scalar hyperbolic conservation laws where the statistical information about inputs such as the initial datum and (possibly discontinuous) flux function are inferred from noisy measurements. In particular, the Lipschitz continuity of the measurement to posterior map as well as the stability of the posterior to approximations, are established with respect to the Wasserstein distance. Numerical experiments are presented to illustrate the derived estimates.


title: Magnetic topology in coupled binaries, spin-orbital resonances, and flares

abstract: We consider topological configurations of the magnetically coupled spinning stellar binaries (e.g., merging neutron stars or interacting star-planet systems). We discuss conditions when the stellar spins and the orbital motion nearly `compensate' each other, leading to very {\it slow} overall winding of the coupled magnetic fields; slowly winding configurations allow gradual accumulation of magnetic energy, that is eventually released in a flare when the instability threshold is reached. We find that this slow winding can be global and/or local. We describe the topology of the relevant space $\mathbb{F}=T^1S^2$ as the unit tangent bundle of the two-sphere and find conditions for slowly winding configurations in terms of magnetic moments, spins and orbital momentum. These conditions become ambiguous near the topological bifurcation points; in certain cases they also depend on the relative phases of the spin and orbital motions. In the case of merging magnetized neutron stars, if one of the stars is a millisecond pulsar, spinning at $\sim$ 10 msec, the global resonance $\omega_1+\omega_2= 2 \Omega$ (spin-plus beat is two times the orbital period) occurs approximately a second before the merger; the total energy of the flare can be as large as $10\%$ of the total magnetic energy, producing bursts of luminosity $\sim 10^{44}$ erg s$^{-1}$. Higher order local resonances may have similar powers, since the amount of involved magnetic flux tubes may be comparable to the total connected flux.


title: Functional lower bounds for restricted arithmetic circuits of depth four

abstract: Recently, Forbes, Kumar and Saptharishi [CCC, 2016] proved that there exists an explicit $d^{O(1)}$-variate and degree $d$ polynomial $P_{d}\in VNP$ such that if any depth four circuit $C$ of bounded formal degree $d$ which computes a polynomial of bounded individual degree $O(1)$, that is functionally equivalent to $P_d$, then $C$ must have size $2^{\Omega(\sqrt{d}\log{d})}$. The motivation for their work comes from Boolean Circuit Complexity. Based on a characterization for $ACC^0$ circuits by Yao [FOCS, 1985] and Beigel and Tarui [CC, 1994], Forbes, Kumar and Saptharishi [CCC, 2016] observed that functions in $ACC^0$ can also be computed by algebraic $\Sigma\mathord{\wedge}\Sigma\Pi$ circuits (i.e., circuits of the form -- sums of powers of polynomials) of $2^{\log^{O(1)}n}$ size. Thus they argued that a $2^{\omega(\log^{O(1)}{n})}$ "functional" lower bound for an explicit polynomial $Q$ against $\Sigma\mathord{\wedge}\Sigma\Pi$ circuits would imply a lower bound for the "corresponding Boolean function" of $Q$ against non-uniform $ACC^0$. In their work, they ask if their lower bound be extended to $\Sigma\mathord{\wedge}\Sigma\Pi$ circuits. In this paper, for large integers $n$ and $d$ such that $\omega(\log^2n)\leq d\leq n^{0.01}$, we show that any $\Sigma\mathord{\wedge}\Sigma\Pi$ circuit of bounded individual degree at most $O\left(\frac{d}{k^2}\right)$ that functionally computes Iterated Matrix Multiplication polynomial $IMM_{n,d}$ ($\in VP$) over $\{0,1\}^{n^2d}$ must have size $n^{\Omega(k)}$. Since Iterated Matrix Multiplication $IMM_{n,d}$ over $\{0,1\}^{n^2d}$ is functionally in $GapL$, improvement of the afore mentioned lower bound to hold for quasipolynomially large values of individual degree would imply a fine-grained separation of $ACC^0$ from $GapL$.


title: Tight distance-dependent estimators for screening two-center and three-center short-range Coulomb integrals over Gaussian basis functions

abstract: We derive distance-dependent estimators for two-center and three-center electron repulsion integrals over a short-range Coulomb potential, $\textrm{erfc}(\omega r_{12})/r_{12}$. These estimators are much tighter than one based on the Schwarz inequality and can be viewed as a complement to the distance-dependent estimators for four-center short-range Coulomb integrals and for two-center and three-center full Coulomb integrals previously reported. Because the short-range Coulomb potential is commonly used in solid-state calculations, including those with the HSE functional and with our recently introduced range-separated periodic Gaussian density fitting, we test our estimators on a diverse set of periodic systems using a wide range of the range-separation parameter $\omega$. These tests demonstrate the robust tightness of our estimators, which are then used with integral screening to calculate periodic three-center short-range Coulomb integrals with linear scaling in system size.


title: MALS SALT-NOT survey of MIR-selected powerful radio-bright AGN at 0<z<3.5

abstract: We present results of an optical spectroscopic survey using SALT and NOT to build a WISE mid-infrared color-based, dust-unbiased sample of powerful radio-bright ($>$200 mJy at 1.4 GHz) AGN for the MeerKAT Absorption Line Survey (MALS). Our sample has 250 AGN (median $z=1.8$) showing emission lines, 26 with no emission lines, and 27 without optical counterparts. Overall, our sample is fainter ($\Delta i$=0.6 mag) and redder ($\Delta(g-i)$=0.2 mag) than radio-selected quasars, and representative of fainter quasar population detected in optical surveys. About 20% of the sources are narrow line AGN (NLAGN) $-$ 65% of these, at $z < 0.5$ are galaxies without strong nuclear emission, and 10% at $z>1.9$, have emission line ratios similar to radio galaxies. The farthest NLAGN in our sample is M1513$-$2524 ($z_{em}=3.132$), and the largest (size$\sim$330 kpc) is M0909$-$3133 ($z_{em}=0.884$). We discuss in detail 110 AGN at $1.9 < z < 3.5$. Despite representing the radio loudest quasars (median $R$=3685), their Eddington ratios are similar to the SDSS quasars having lower $R$. We detect 4 CIV BALQSOs, all among AGN with least $R$, and highest black hole masses and Eddington ratios. The BAL detection rate ($4^{+3}_{-2}$%) is consistent with that seen in extremely powerful ($L_{1.4GHz}>10^{25}$ WHz$^{-1}$) quasars. Using optical light-curves, radio polarization and $\gamma$-ray detections, we identify 7 high-probability BL Lacs. We also summarize the full MALS footprint to search for HI 21-cm and OH 18-cm lines at $z<2$.


title: Fiber decomposition of deterministic reaction networks with applications

abstract: Deterministic reaction networks (RNs) are tools to model diverse biological phenomena characterized by particle systems, when there are abundant number of particles. Examples include but are not limited to biochemistry, molecular biology, genetics, epidemiology, and social sciences. In this chapter we propose a new type of decomposition of RNs, called fiber decomposition. Using this decomposition, we establish lifting of mass-action RNs preserving stationary properties, including multistationarity and absolute concentration robustness. Such lifting scheme is simple and explicit which imposes little restriction on the reaction networks. We provide examples to illustrate how this lifting can be used to construct RNs preserving certain dynamical properties.


title: A Cooperative Optimal Mining Model for Bitcoin

abstract: We analyze Bitcoin mining from the perspective of a game and propose an optimal mining model that maximizes profits of pools and miners. The model is a two-stage Stackelberg game in which each stage forms a sub-game. In stage I, pools are the leaders who assign a computing power to be consumed by miners. In stage II, miners decide of their power consumption and distribution. They find themselves in a social dilemma in which they must choose between mining in solo, therefore prioritizing their individual preferences, and participating in a pool for the collective interest. The model relies on a pool protocol based on a simulated game in which the miners compete for the reward won by the pool. The solutions for the stage I sub-game and the simulated protocol game are unique and stable Nash equilibriums while the stage II sub-game leads to a stable cooperative equilibrium only when miners choose their strategies according to certain criteria. We conclude that the cooperative optimal mining model has the potential to favor Bitcoin decentralization and stability. Mainly, the social dilemma faced by miners together with the balance of incentives ensure a certain distribution of the network computing power between pools and solo miners, while equilibriums in the game solutions provide stability to the system.


title: Flavor-Violating Higgs Decays and Stellar Cooling Anomalies in Axion Models

abstract: We study a class of DFSZ-like models for the QCD axion that can address observed anomalies in stellar cooling. Stringent constraints from SN1987A and neutron stars are avoided by suppressed couplings to nucleons, while axion couplings to electrons and photons are sizable. All axion couplings depend on few parameters that also control the extended Higgs sector, in particular lepton flavor-violating couplings of the Standard Model-like Higgs boson $h$. This allows us to correlate axion and Higgs phenomenology, and we find that that ${\rm BR}(h \to \tau e)$ can be as large as the current experimental bound of 0.22%, while ${\rm BR} (h \to \mu \mu)$ can be larger than in the Standard Model by up to 70%. Large parts of the parameter space will be tested by the next generation of axion helioscopes such as the IAXO experiment.


title: Emerging mechanisms of magnetocaloric effect in phase-separated metals

abstract: We present a study of the magnetocaloric effect in metallic systems exhibiting first-order magnetic transitions and focus on consequences of magnetic phase separation. We account for ferrimagnetic, ferromagnetic, and Neel antiferromagnetic order. Based on the archetypal Hubbard model being treated within the mean-field approximation, we provide and explore its implications on the field-induced entropy change in metallic system with phase separation. Chosen framework allows us to properly analyze phase volumes' dependence on parameters of phase-separated (PS) system. Moreover, an account for phase separation boundaries as functions of magnetic field provides a natural splitting of the PS region, where each subregion corresponds to a different temperature dependence of entropy change: moving from one subregion to the other produces a kink, followed by a strong linear growth of entropy change. We encounter a second-order magnetic transition from paramagnetic to antiferromagnetic phase in PS region that occurs for particular parameter values. Despite the fact that both phases have zero total magnetization, the transition has a strong impact on entropy change.


title: TLA: Twitter Linguistic Analysis

abstract: Linguistics has been instrumental in developing a deeper understanding of human nature. Words are indispensable to bequeath the thoughts, emotions, and purpose of any human interaction, and critically analyzing these words can elucidate the social and psychological behavior and characteristics of these social animals. Social media has become a platform for human interaction on a large scale and thus gives us scope for collecting and using that data for our study. However, this entire process of collecting, labeling, and analyzing this data iteratively makes the entire procedure cumbersome. To make this entire process easier and structured, we would like to introduce TLA(Twitter Linguistic Analysis). In this paper, we describe TLA and provide a basic understanding of the framework and discuss the process of collecting, labeling, and analyzing data from Twitter for a corpus of languages while providing detailed labeled datasets for all the languages and the models are trained on these datasets. The analysis provided by TLA will also go a long way in understanding the sentiments of different linguistic communities and come up with new and innovative solutions for their problems based on the analysis.


title: Synchronization transition in the two-dimensional Kuramoto model with dichotomous noise

abstract: We numerically study the celebrated Kuramoto model of identical oscillators arranged on the sites of a two-dimensional periodic square lattice and subject to nearest neighbor interactions and dichotomous noise. In the nonequilibrium stationary state attained at long time, the model exhibits a Berezinskii-Kosterlitz-Thouless ($BKT$)-like transition between a phase at low noise amplitude characterized by quasi long-range order (critically ordered phase) and algebraic decay of correlations and a phase at high noise amplitude that is characterized by complete disorder and exponential decay of correlations. The interplay between the noise amplitude and the noise correlation time is investigated, and the complete, nonequilibrium stationary-state phase diagram of the model is obtained. We further study the dynamics of a single topological defect for various amplitude and correlation time of the noise. Our analysis reveals that a finite correlation time promotes vortex excitations, thereby lowering the critical noise amplitude of the transition with an increase in correlation time. In the suitable limit, the resulting phase diagram allows to estimate the critical temperature of the equilibrium $BKT$ transition, which is consistent with that obtained from the study of the dynamics in the Gaussian white noise limit.


title: Frequency-Domain Data-Driven Controller Synthesis for Unstable LPV Systems

abstract: Synthesizing controllers directly from frequency-domain measurement data is a powerful tool in the linear time-invariant framework. Ever-increasing performance requirements necessitate extending these approaches to account for plant variations. The aim of this paper is to develop frequency-domain analysis and synthesis conditions for local internal stability and $\mathcal{H}_\infty$-performance of single-input single-output linear parameter-varying systems. The developed synthesis procedure only requires frequency-domain measurement data of the system and does not need a parametric model of the plant. The capabilities of the synthesis procedure are demonstrated on an unstable nonlinear system.


title: $f({\sf R})$ Gravity Wormholes sourced by a Phantom Scalar Field

abstract: We derive an exact wormhole spacetime supported by a phantom scalar field in the context of $f({\sf R})$ gravity. Without specifying the form of the $f({\sf R})$ function, the scalar field self-interacts with a mass term potential that is derived from the scalar equation and the resulting $f({\sf R})$ model is purely supported by the scalar field and it is free of ghosts and avoids the tachyonic instability.


title: Physical limitations of the Hohenberg-Mermin-Wagner theorem

abstract: The Hohenberg-Mermin-Wagner (HMW) theorem states that infrared (IR) fluctuations prevent long-range order which breaks continuous symmetries in two dimensions (2D), at finite temperatures. We note that the theorem becomes physically effective for superconductivity (SC) only for astronomical sample sizes, so it does not prevent 2D SC in practice. We systematically explore the sensitivity of the magnetic and SC versions of the theorem to finite-size and disorder effects. For magnetism, finite-size effects, disorder, and perpendicular coupling can all restore the order parameter at a non-negligible value of $T_c$ equally well, making the physical reason for finite $T_c$ sample-dependent. For SC, an alternative version of the HMW theorem is presented, in which the temperature cutoff is set by Cooper pairing, in place of the Fermi energy in the standard version. It still allows 2D SC at $2$--$3$ times the room temperature when the interaction scale is large and Cooper pairs are small, the case with high-$T_c$ SC in the cuprates. Thus IR fluctuations do not prevent 2D SC at room temperatures in samples of any reasonable size, by any known version of the HMW argument. A possible approach to derive mechanism-dependent upper bounds for SC $T_c$ is pointed out.


title: Evolution of decaying particles and decay products in various scenarios for the future expansion of the universe

abstract: We examine nonrelativistic particles that decay into relativistic products in big rip, little rip, and pseudo-rip models for the future evolution of the universe. In contrast to decays that occur in standard $\Lambda$CDM, the evolution of the ratio $r$ of the energy density of the relativistic decay products to the energy density of the initially decaying particles can decrease with time in all of these models. In big rip and little rip models, $r$ always goes to zero asymptotically, while this ratio evolves to infinity or a constant in pseudo-rip models.


abstract: Traditionally, quantization is designed to minimize the reconstruction error of a data source. When considering downstream classification tasks, other measures of distortion can be of interest; such as the 0-1 classification loss. Furthermore, it is desirable that the performance of these quantizers not deteriorate once they are deployed into production, as relearning the scheme online is not always possible. In this work, we present a class of algorithms that learn distributed quantization schemes for binary classification tasks. Our method performs well on unseen data, and is faster than previous methods proportional to a quadratic term of the dataset size. It works by regularizing the 0-1 loss with the reconstruction error. We present experiments on synthetic mixture and bivariate Gaussian data and compare training, testing, and generalization errors with a family of benchmark quantization schemes from the literature. Our method is called Regularized Classification-Aware Quantization.


title: Next-to SV resummed Drell-Yan cross section beyond leading-logarithm

abstract: We present the resummed predictions for inclusive cross section for Drell-Yan (DY) production up to next-to-next-to leading logarithmic ($\rm \overline{ NNLL}$) accuracy taking into account both soft virtual (SV) and next-to SV (NSV) threshold logarithms. We restrict ourselves to resummed contributions only from quark anti-quark ($q \bar q$) initiated channels. The resummation is performed in Mellin-$N$ space. We derive the $N$-dependent coefficients and the $N$-independent constants to desired accuracy for our study. The resummed results are matched through the minimal prescription procedure with the fixed order results. We find that the resummation, taking into account the NSV terms, appreciably increases the cross section while decreasing the sensitivity to renormalisation scale. \textcolor{black}{We observe that, at 13 TeV LHC energies, the SV+NSV resummation at $\rm \overline{ NLL} (\rm \overline{NNLL})$ gives about 8\% (2\%) corrections respectively to the NLO (NNLO) results for the considered $Q$ range: 150-3500 GeV}. In addition, the absence of quark gluon initiated contributions to NSV part in the resummed terms leaves large factorisation scale dependence indicating their importance at NSV level. We also study the numerical impact of $N$-independent constants and explore the ambiguity involved in exponetiating them. Finally we present our predictions for the neutral Drell-Yan process at various center of mass of energies.


title: An Efficient Multi-objective Evolutionary Approach for Solving the Operation of Multi-Reservoir System Scheduling in Hydro-Power Plants

abstract: This paper tackles the short-term hydro-power unit commitment problem in a multi-reservoir system - a cascade-based operation scenario. For this, we propose a new mathematical modelling in which the goal is to maximize the total energy production of the hydro-power plant in a sub-daily operation, and, simultaneously, to maximize the total water content (volume) of reservoirs. For solving the problem, we discuss the Multi-objective Evolutionary Swarm Hybridization (MESH) algorithm, a recently proposed multi-objective swarm intelligence-based optimization method which has obtained very competitive results when compared to existing evolutionary algorithms in specific applications. The MESH approach has been applied to find the optimal water discharge and the power produced at the maximum reservoir volume for all possible combinations of turbines in a hydro-power plant. The performance of MESH has been compared with that of well-known evolutionary approaches such as NSGA-II, NSGA-III, SPEA2, and MOEA/D in a realistic problem considering data from a hydro-power energy system with two cascaded hydro-power plants in Brazil. Results indicate that MESH showed a superior performance than alternative multi-objective approaches in terms of efficiency and accuracy, providing a profit of \$412,500 per month in a projection analysis carried out.


title: Asymmetric One-Dimensional Slow Electron Holes

abstract: Slow solitary positive-potential peaks sustained by trapped electron deficit in a plasma with asymmetric ion velocity distributions are in principle asymmetric, involving a potential change across the hole. It is shown theoretically how to construct such asymmetric electron holes, thus providing fully consistent solutions of the one-dimensional Vlasov-Poisson equation for a wide variety of prescribed background ion velocity distributions. Because of ion reflection forces experienced by the hole, there is generally only one discrete slow hole velocity that is in equilibrium. Moreover the equilibrium is unstable unless there is a local minimum in the ion velocity distribution, in which the hole velocity then resides. For stable equilibria with Maxwellian electrons, the potential drop across the hole is shown to be $\Delta\phi\simeq{2\over 9}f'''{T_e\over e}({e\psi\over m_i})^2$, where $\psi$ is the hole peak potential, $f'''$ is the third derivative of the background ion velocity distribution function at the hole velocity, and $T_e$ the electron temperature. Potential asymmetry is small for holes of the amplitudes usually observed, $\psi\lesssim 0.5T_e/e$.


title: Spectrum of localized states in fermionic chains with defect and adiabatic charge pumping

abstract: In this paper, we study the localized states of a generic quadratic fermionic chain with finite-range couplings and an inhomogeneity in the hopping (defect) that breaks translational invariance. When the hopping of the defect vanishes, which represents an open chain, we obtain a simple bulk-edge correspondence: the zero-energy modes localized at the ends of the chain are related to the roots of a polynomial determined by the couplings of the Hamiltonian of the bulk. From this result, we define an index that characterizes the different topological phases of the system and can be easily computed by counting the roots of the polynomial. As the defect is turned on and varied adiabatically, the zero-energy modes may cross the energy gap and connect the valence and conduction bands. We analyze the robustness of the connection between bands against perturbations of the Hamiltonian. The pumping of states from one band to the other allows the creation of particle-hole pairs in the bulk. An important ingredient for our analysis is the transformation of the Hamiltonian under the standard discrete symmetries, $C$, $P$, $T$, as well as a fourth one, peculiar to our system, that is related to the existence of a gap and localized states.


title: Online Projected Gradient Descent for Stochastic Optimization with Decision-Dependent Distributions

abstract: This paper investigates the problem of tracking solutions of stochastic optimization problems with time-varying costs that depend on random variables with decision-dependent distributions. In this context, we propose the use of an online stochastic gradient descent method to solve the optimization, and we provide explicit bounds in expectation and in high probability for the distance between the optimizers and the points generated by the algorithm. In particular, we show that when the gradient error due to sampling is modeled as a sub-Weibull random variable, then the tracking error is ultimately bounded in expectation and in high probability. The theoretical findings are validated via numerical simulations in the context of charging optimization of a fleet of electric vehicles.


title: Wood Stove Combustion Modeling and Simulation: Technical Review and Recommendations

abstract: With the rise in wood stove use in recent years, a number of experimental research efforts have been undertaken with the objective of developing cleaner, and more efficient stove designs. However, numerical modeling of stoves is still in its nascent years, thus making it imperative to start adapting more mature computational techniques, established in other combustion-based applications (e.g., combustion engines), to wood stoves. This study focuses on a critical review of the state-of-the-art in wood combustion modeling techniques and on providing recommendations for needed enhancements. In order to develop an optimum conceptual numerical model for wood stoves, the basic processes inside the stove are broken down and isolated into generic computational domains: surface reactions, gas reactions, and the fluid domain. Then state-of-the-art numerical efforts for wood stove and similar applications are broken down in a similar manner to highlight points of strength and possible improvements. Finally, based on the comparison between the different models, an improved modeling approach is proposed as a road map for future research efforts to achieve a higher fidelity wood stove numerical model.


title: Superconductivity in disordered locally noncentrosymmetric materials: an application to CeRh$_2$As$_2$

abstract: Layered three-dimensional centrosymmetric crystals can exhibit characteristics of noncentrosymmetric materials. This happens when each individual layer alone lacks inversion, but, when combined, inversion symmetry is restored; hence the designation: locally noncentrosymmertic superconductors (LNCSs). In LNCSs, the effects of impurities and subdominant magnetic field induced pairing channels remain unexplored. Using a minimal model, we examine all pairing channels and show that there is always a subdominant superconducting instability that is favored at high magnetic fields, which can substantially alter the magnetic field -- temperature phase diagram. Also, we find that the phase diagram responds to disorder in a non-monotonic way, which can be subjected to experimental verification. We apply these ideas to the recently unveiled two-phase superconducting phase diagram of CeRh$_2$As$_2$. We identify the two phases as singlet-triplet mixed even- and odd-parity states at low and a high fields, respectively. Furthermore, we predict the presence of two superconducting phases also for in-plane magnetic fields in cleaner samples, since a high-field phase could have been so far hindered by impurity effects.


title: Quantification of Captopril using Ultra High Performance Liquid Chromatography

abstract: We develop an efficient, sensitive, accurate, cost effective, and greener liquid chromatography method for the determination of a quantification method concerning Captopril in pharmaceutical dosage forms by using ultra high performance liquid chromatography (UHPLC) equipment. Satisfactory resolution was able to be achieved using a Water Aquity UHPLC BEH C18 1.7 $\mu$m $\times$ 2.1 mm $\times$ 50 mm column and an isocratic mobile phase composition of Methanol, Milli-Q water and Trifluoroacetic acid in the ratio of 55:45:0.05 v/v/v respectively at a flow rate of 0.1mL/min and a wavelength detection of 220nm. The retention time for Captopril was found to be 1.744 minutes while for Captopril disulfide it was found to be 2.657 minutes. The method was validated based on the ICH and USP guidelines.


title: Registration of 3D Point Sets Using Correntropy Similarity Matrix

abstract: This work focuses on Registration or Alignment of 3D point sets. Although the Registration problem is a well established problem and it's solved using multiple variants of Iterative Closest Point (ICP) Algorithm, most of the approaches in the current state of the art still suffers from misalignment when the \textit{Source} and the \textit{Target} point sets are separated by large rotations and translation. In this work, we propose a variant of the Standard ICP algorithm, where we introduce a Correntropy Relationship Matrix in the computation of rotation and translation component which attempts to solve the large rotation and translation problem between \textit{Source} and \textit{Target} point sets. This matrix is created through correntropy criterion which is updated in every iteration. The correntropy criterion defined in this approach maintains the relationship between the points in the \textit{Source} dataset and the \textit{Target} dataset. Through our experiments and validation we verify that our approach has performed well under various rotation and translation in comparison to the other well-known state of the art methods available in the Point Cloud Library (PCL) as well as other methods available as open source. We have uploaded our code in the github repository for the readers to validate and verify our approach https://github.com/aralab-unr/CoSM-ICP.


title: The Foata-Fuchs proof of Cayley's formula, and its probabilistic uses

abstract: We present a very simple bijective proof of Cayley's formula due to Foata and Fuchs (1970). This bijection turns out to be very useful when seen through a probabilistic lens; we explain some of the ways in which it can be used to derive probabilistic identities, bounds, and growth procedures for random trees with given degrees, including random d-ary trees. We also introduce a partial order on the degree sequences of rooted trees, and conjecture that it induces a stochastic partial order on heights of random rooted trees with given degrees.


title: New searches for continuous gravitational waves from seven fast pulsars

abstract: We conduct searches for continuous gravitational waves from seven pulsars, that have not been targeted in continuous wave searches of Advanced LIGO data before. We target emission at exactly twice the rotation frequency of the pulsars and in a small band around such frequency. The former search assumes that the gravitational wave quadrupole is changing phase-locked with the rotation of the pulsar. The search over a range of frequencies allows for differential rotation between the component emitting the radio signal and the component emitting the gravitational waves, for example the crust or magnetosphere versus the core. Timing solutions derived from the Arecibo 327-MHz Drift-Scan Pulsar Survey (AO327) observations are used. No evidence of a signal is found and upper limits are set on the gravitational wave amplitude. For one of the pulsars we probe gravitational wave intrinsic amplitudes just a factor of 3.8 higher than the spin-down limit, assuming a canonical moment of inertia of $10^{38}$ kg m$^2$. Our tightest ellipticity constraint is $1.5 \times 10^{-8}$, which is a value well within the range of what a neutron star crust could support.


title: Machine Learning Approaches to Automated Flow Cytometry Diagnosis of Chronic Lymphocytic Leukemia

abstract: Flow cytometry is a technique that measures multiple fluorescence and light scatter-associated parameters from individual cells as they flow a single file through an excitation light source. These cells are labeled with antibodies to detect various antigens and the fluorescence signals reflect antigen expression. Interpretation of the multiparameter flow cytometry data is laborious, time-consuming, and expensive. It involves manual interpretation of cell distribution and pattern recognition on two-dimensional plots by highly trained medical technologists and pathologists. Using various machine learning algorithms, we attempted to develop an automated analysis for clinical flow cytometry cases that would automatically classify normal and chronic lymphocytic leukemia cases. We achieved the best success with the Gradient Boosting. The XGBoost classifier achieved a specificity of 1.00 and a sensitivity of 0.67, a negative predictive value of 0.75, a positive predictive value of 1.00, and an overall accuracy of 0.83 in prospectively classifying cases with malignancies.


title: What Do You Get When You Cross Beam Search with Nucleus Sampling?

abstract: We combine beam search with the probabilistic pruning technique of nucleus sampling to create two deterministic nucleus search algorithms for natural language generation. The first algorithm, p-exact search, locally prunes the next-token distribution and performs an exact search over the remaining space. The second algorithm, dynamic beam search, shrinks and expands the beam size according to the entropy of the candidate's probability distribution. Despite the probabilistic intuition behind nucleus search, experiments on machine translation and summarization benchmarks show that both algorithms reach the same performance levels as standard beam search.


title: A flexible approach for variable selection in large-scale healthcare database studies with missing covariate and outcome data

abstract: Prior work has shown that combining bootstrap imputation with tree-based machine learning variable selection methods can provide good performances achievable on fully observed data when covariate and outcome data are missing at random (MAR). This approach however is computationally expensive, especially on large-scale datasets. We propose an inference-based method, called RR-BART, which leverages the likelihood-based Bayesian machine learning technique, Bayesian additive regression trees, and uses Rubin's rule to combine the estimates and variances of the variable importance measures on multiply imputed datasets for variable selection in the presence of MAR data. We conduct a representative simulation study to investigate the practical operating characteristics of RR-BART, and compare it with the bootstrap imputation based methods. We further demonstrate the methods via a case study of risk factors for 3-year incidence of metabolic syndrome among middle-aged women using data from the Study of Women's Health Across the Nation (SWAN). The simulation study suggests that even in complex conditions of nonlinearity and nonadditivity with a large percentage of missingness, RR-BART can reasonably recover both prediction and variable selection performances, achievable on the fully observed data. RR-BART provides the best performance that the bootstrap imputation based methods can achieve with the optimal selection threshold value. In addition, RR-BART demonstrates a substantially stronger ability of detecting discrete predictors. Furthermore, RR-BART offers substantial computational savings. When implemented on the SWAN data, RR-BART adds to the literature by selecting a set of predictors that had been less commonly identified as risk factors but had substantial biological justifications.


title: A Gemini-NIFS view of the merger remnant NGC 34

abstract: The merger remnant NGC 34 is a local luminous infrared galaxy (LIRG) hosting a nuclear starburst and a hard X-ray source associated with a putative, obscured Seyfert~2 nucleus. In this work, we use adaptive optics assisted near infrared (NIR) integral field unit observations of this galaxy to map the distribution and kinematics of the ionized and molecular gas in its inner $\mathrm{1.2\,kpc \times 1.2\,kpc}$, with a spatial resolution of 70~pc. The molecular and ionized gas kinematics is consistent with a disc with projected major axis along a mean PA~=~$\mathrm{-9^{\circ}.2 \pm 0^{\circ}.9}$. Our main findings are that NGC~34 hosts an AGN and that the nuclear starburst is distributed in a circumnuclear star-formation ring with inner and outer radii of $\approx$~60 and 180~pc, respectively, as revealed by maps of the $\mathrm{[Fe II] / Pa\beta}$ and $\mathrm{H_{2} / Br\gamma}$ emission-line ratios, and corroborated by PCA Tomography analysis. The spatially resolved NIR diagnostic diagram of NGC~34 also identifies a circumnuclear structure dominated by processes related to the stellar radiation field and a nuclear region where $[Fe II]$ and H$_2$ emissions are enhanced relative to the hydrogen recombination lines. We estimate that the nuclear X-ray source can account for the central H$_2$ enhancement and conclude that $[Fe II]$ and H$_2$ emissions are due to a combination of photo-ionization by young stars, excitation by X-rays produced by the AGN and shocks. These emission lines show nuclear, broad, blue-shifted components that can be interpreted as nuclear outflows driven by the AGN.


title: An Analytic Solution to the Kozai-Lidov Evolution Equations

abstract: A test particle in a noncoplanar orbit about a member of a binary system can undergo Kozai-Lidov oscillations in which tilt and eccentricity are exchanged. An initially circular highly inclined particle orbit can reach high eccentricity. We consider the nonlinear secular evolution equations previously obtained in the quadrupole approximation. For the important case that the initial eccentricity of the particle orbit is zero, we derive an analytic solution for the particle orbital elements as a function of time that is exact within the quadrupole approximation. The solution involves only simple trigonometric and hyperbolic functions. It simplifies in the case that the initial particle orbit is close to being perpendicular to the binary orbital plane. The solution also provides an accurate description of particle orbits with nonzero but sufficiently small initial eccentricity. It is accurate over a range of initial eccentricity that broadens at higher initial inclinations. In the case of an initial inclination of pi/3, an error of 1% at maximum eccentricity occurs for initial eccentricities of about 0.1.


title: Stable and efficient FEM-BEM coupling with OSRC regularisation for acoustic wave transmission

abstract: The finite element method (FEM) and the boundary element method (BEM) can numerically solve the Helmholtz system for acoustic wave propagation. When an object with heterogeneous wave speed or density is embedded in an unbounded exterior medium, the coupled FEM-BEM algorithm promises to combine the strengths of each technique. The FEM handles the heterogeneous regions while the BEM models the homogeneous exterior. Even though standard FEM-BEM algorithms are effective, they do require stabilisation at resonance frequencies. One such approach is to add a regularisation term to the system of equations. This algorithm is stable at all frequencies but also brings higher computational costs. This study proposes a regulariser based on the on-surface radiation conditions (OSRC). The OSRC operators are also used to precondition the boundary integral operators and combined with incomplete LU factorisations for the volumetric weak formulation. The proposed preconditioning strategy improves the convergence of iterative linear solvers significantly, especially at higher frequencies.


title: Uncertainty Estimation and Out-of-Distribution Detection for Counterfactual Explanations: Pitfalls and Solutions

abstract: Whilst an abundance of techniques have recently been proposed to generate counterfactual explanations for the predictions of opaque black-box systems, markedly less attention has been paid to exploring the uncertainty of these generated explanations. This becomes a critical issue in high-stakes scenarios, where uncertain and misleading explanations could have dire consequences (e.g., medical diagnosis and treatment planning). Moreover, it is often difficult to determine if the generated explanations are well grounded in the training data and sensitive to distributional shifts. This paper proposes several practical solutions that can be leveraged to solve these problems by establishing novel connections with other research works in explainability (e.g., trust scores) and uncertainty estimation (e.g., Monte Carlo Dropout). Two experiments demonstrate the utility of our proposed solutions.


title: kNet: A Deep kNN Network To Handle Label Noise

abstract: Deep Neural Networks require large amounts of labeled data for their training. Collecting this data at scale inevitably causes label noise.Hence,the need to develop learning algorithms that are robust to label noise. In recent years, k Nearest Neighbors (kNN) emerged as a viable solution to this problem. Despite its success, kNN is not without its problems. Mainly, it requires a huge memory footprint to store all the training samples and it needs an advanced data structure to allow for fast retrieval of the relevant examples, given a query sample. We propose a neural network, termed kNet, that learns to perform kNN. Once trained, we no longer need to store the training data, and processing a query sample is a simple matter of inference. To use kNet, we first train a preliminary network on the data set, and then train kNet on the penultimate layer of the preliminary network.We find that kNet gives a smooth approximation of kNN,and cannot handle the sharp label changes between samples that kNN can exhibit. This indicates that currently kNet is best suited to approximate kNN with a fairly large k. Experiments on two data sets show that this is the regime in which kNN works best,and can therefore be replaced by kNet.In practice, kNet consistently improve the results of all preliminary networks, in all label noise regimes, by up to 3%.


title: Recent Developments in Inference: Practicalities for Applied Economics

abstract: We provide a review of recent developments in the calculation of standard errors and test statistics for statistical inference. While much of the focus of the last two decades in economics has been on generating unbiased coefficients, recent years has seen a variety of advancements in correcting for non-standard standard errors. We synthesize these recent advances in addressing challenges to conventional inference, like heteroskedasticity, clustering, serial correlation, and testing multiple hypotheses. We also discuss recent advancements in numerical methods, such as the bootstrap, wild bootstrap, and randomization inference. We make three specific recommendations. First, applied economists need to clearly articulate the challenges to statistical inference that are present in data as well as the source of those challenges. Second, modern computing power and statistical software means that applied economists have no excuse for not correctly calculating their standard errors and test statistics. Third, because complicated sampling strategies and research designs make it difficult to work out the correct formula for standard errors and test statistics, we believe that in the applied economics profession it should become standard practice to rely on asymptotic refinements to the distribution of an estimator or test statistic via bootstrapping. Throughout, we reference built-in and user-written Stata commands that allow one to quickly calculate accurate standard errors and relevant test statistics.


title: Geometry-based circulation of local thermal current in quantum harmonic and Bose Hubbard systems

abstract: A geometry-based mechanism for generating steady-state internal circulation of local thermal currents is demonstrated by harmonically coupled quantum oscillators formulated by the Redfield quantum master equation (RQME) and the Bose Hubbard model (BHM) of phonons formulated by the Lindblad quantum master equation (LQME) using the simple multi-path geometry of a triangle. Driven by two reservoirs at different temperatures, both systems can exhibit an atypical local thermal current flowing against the total current. However, the total thermal current behaves normally. While the RQME of harmonically coupled quantum oscillators allows an analytical solution, the LQME of the interacting BHM can be solved numerically. The emergence of the geometry-based circulation in both systems demonstrates the ubiquity and robustness of the mechanism. In the high-temperature limit, the results agree with the classical results, confirming the generality of the geometric-based circulation across the quantum and classical boundary. Possible experimental implications and applications are briefly discussed.


title: On anomalous optical beam shifts at near-normal incidence

abstract: We develop the theory of optical beam shifts (both Goos-Hanchen and Imbert-Fedorov) for the case of near-normal incidence, when the incident angle becomes comparable with the angular beam divergence. Such a situation naturally leads to strong enhancement of the shifts reported recently [ACS Photonics 6, 2530 (2019)]. Experimental results find complete and rigorous explanation in our generalized theory. In addition, the developed theory uncovers the unified origin of the anomalous beam shifts enhancement via the Berry phase singularity. We also propose a simple experimental scheme involving quarter-wave plate that allows to observe the giant transverse and longitudinal, spatial and angular beam shifts simultaneously. Our results can find applications in spin-orbit photonics, polarization optics, sensing applications, and quantum weak measurements.


title: Quasar clustering at redshift 6

abstract: Large-scale surveys over the last years have revealed about 300 QSOs at redshift above 6. Follow-up observations identified surprising properties, such as the very high black hole (BH) masses, spatial correlations with surrounding cold gas of the host galaxy, or high CIV-MgII velocity shifts. In particular, the discovery of luminous high-redshift quasars suggests that at least some black holes likely have large masses at birth and grow efficiently. We aim at quantifying quasar pairs at high redshift for a large sample of objects. This provides a new key constraint on a combination of parameters related to the origin and assembly for the most massive black holes: BH formation efficiency and clustering, growth efficiency and relative contribution of BH mergers. We observed 116 spectroscopically confirmed QSOs around redshift 6 with the simultaneous 7-channel imager GROND in order to search for companions. Applying identical colour-colour cuts as for those which led to the spectroscopically confirmed QSO, we perform LePHARE fits to the 26 best QSO pair candidates, and obtained spectroscopic observations for 11 of those. e do not find any QSO pair with a companion brighter than M1450(AB) < -26 mag within our 0.1-3.3 h^-1 cMpc search radius, in contrast to the serendipitous findings in the redshift range 4--5. However, a low fraction of such pairs at this luminosity and redshift is consistent with indications from present-day cosmological-scale galaxy evolution models. In turn, the incidence of L- and T-type brown dwarfs which occupy a similar colour space as z ~ 6 QSOs, is higher than expected, by a factor of 5 and 20, respectively.


title: Stray magnetic field and stability of time-dependent viscous electron flow

abstract: A hydrodynamic flow of electrons driven by an oscillating electric field in Dirac and Weyl semimetals is investigated. It is found that a double-peak profile of the electric current appears and is manifested in a stray magnetic field with peaks in one of the field components. The nontrivial current profile originates from the interplay of viscous and inertial properties of the electron fluid as well as the boundary conditions. Analytical results are supported by numerical calculations in samples of different geometries such as straight channels, nozzles, and cavities. The double-peak profile of the current is found to be qualitatively insensitive to a specific form of the time dependence of the oscillating electric field and is stable with respect to the sample geometry. A phase diagram and criteria for observing the double-peak structure are determined. In addition, it is shown that nozzles and cavities provide an efficient means to locally enhance or reduce the fluid velocity.


title: Which is the ontology of Dark Matter and Dark Energy?

abstract: We adopt in this work the idea that the building blocks of the visible Universe belong to a class of the irreducible representations of the Poincare group of transformations (the "things") endowed with classificatory quantum numbers ("the properties"). After a discussion of this fundamentality, the question of the nature of both "dark" components of the Universe which are deemed necessary, but have not been observed, is analyzed within this context. We broadly discuss the ontology of dark matter/dark energy in relation to the irreducible representations of the Poincar\'e group + quantum numbers, pointing out some cases in which the candidates can be associated to them, and others for which a reclassification of both the dark and visible (ordinary) components would be needed.


title: Applicability of the two-particle quantization condition to partially-quenched theories

abstract: Partial quenching allows one to consider correlation functions and amplitudes that do not arise in the corresponding unquenched theory. For example, physical $s$-wave pion scattering can be decomposed into $I=0$ and $2$ amplitudes, while, in a partially-quenched extension, the larger symmetry group implies that there are more than two independent scattering amplitudes. It has been proposed that the finite-volume quantization condition of L\"uscher holds for the correlation functions associated with each of the two-particle amplitudes that arise in partially-quenched theories. Using partially-quenched chiral perturbation theory, we show that this proposal fails for those correlation functions for which the corresponding one-loop amplitudes do not satisfy $s$-wave unitarity. For partially-quenched amplitudes that, while being unphysical, do satisfy one-loop $s$-wave unitarity, we argue that the proposal is plausible. Implications for previous work are discussed.


title: Complexity of Source-Sink Monotone 2-Parameter Min Cut

abstract: There are many applications of max flow with capacities that depend on one or more parameters. Many of these applications fall into the "Source-Sink Monotone" framework, a special case of Topkis's monotonic optimization framework, which implies that the parametric min cuts are nested. When there is a single parameter, this property implies that the number of distinct min cuts is linear in the number of nodes, which is quite useful for constructing algorithms to identify all possible min cuts. When there are multiple Source-Sink Monotone parameters and the vector of parameters are ordered in the usual vector sense, the resulting min cuts are still nested. However, the number of distinct min cuts was an open question. We show that even with only two parameters, the number of distinct min cuts can be exponential in the number of nodes.


title: Turbulent field fluctuations in gyrokinetic and fluid plasmas

abstract: A key uncertainty in the design and development of magnetic confinement fusion energy reactors is predicting edge plasma turbulence. An essential step in overcoming this uncertainty is the validation in accuracy of reduced turbulent transport models. Drift-reduced Braginskii two-fluid theory is one such set of reduced equations that has for decades simulated boundary plasmas in experiment, but significant questions exist regarding its predictive ability. To this end, using a novel physics-informed deep learning framework, we demonstrate the first ever direct quantitative comparisons of turbulent field fluctuations between electrostatic two-fluid theory and electromagnetic gyrokinetic modelling with good overall agreement found in magnetized helical plasmas at low normalized pressure. This framework is readily adaptable to experimental and astrophysical environments, and presents a new technique for the numerical validation and discovery of reduced global plasma turbulence models.

categories: physics.plasm-ph cs.LG physics.comp-ph stat.ML


title: Robust unrelated parallel machine scheduling problem with interval release dates

abstract: This paper presents a profound analysis of the robust job scheduling problem with uncertain release dates on unrelated machines. Our model involves minimizing the worst-case makespan and interval uncertainty where each release date belongs to a well-defined interval. Robust optimization requires scenario-based decision-making. A finite subset of feasible scenarios to determine the worst-case regret (a deviation from the optimal makespan) for a particular schedule is indicated. We formulate a mixed-integer nonlinear programming model to solve the underlying problem via three (constructive) greedy algorithms. Polynomial-time solvable cases are also discussed in detail. The algorithms solve the robust combinatorial problem using the makespan criterion and its non-deterministic counterpart. Computational testing compares both robust solutions and different decomposition strategies. Finally, the results confirm that a decomposition strategy applied to the makespan criterion is enough to create a competitive robust schedule.


title: An Exact Algorithm for Large-scale Non-convex Quadratic Facility Location

abstract: We study a general class of quadratic capacitated $p$-location problems facility location problems with single assignment where a non-separable, non-convex, quadratic term is introduced in the objective function to account for the interaction cost between facilities and customer assignments. This problem has many applications in the field of transportation and logistics where its most well-known special case is the single-allocation hub location problem and its many variants. The non-convex, binary quadratic program is linearized by applying a reformulation-linearization technique and the resulting continuous auxiliary variables are projected out using Benders decomposition. The obtained Benders reformulation is then solved using an exact branch-and-cut algorithm that exploits the underlying network flow structure of the decomposed separation subproblems to efficiently generate strong Pareto-optimal Benders cuts. Additional enhancements such as a matheuristic, a partial enumeration procedure, and variable elimination tests are also embedded in the proposed algorithmic framework. Extensive computational experiments on benchmark instances (with up to 500 nodes) and on a new set of instances (with up to 1,000 nodes) of four variants of single-allocation hub location problems confirm the algorithm's ability to scale to large-scale instances.


title: Around Hilbert's theorem: the center of a circle is not constructible by straightedge alone

abstract: In order to state the theorem in the title formally and to review its rigorous proof, we extend and make more precise the Uspenskiy-Shen-Akopyan-Fedorov model of Euclidean constructions with arbitrary points; we also introduce formalizations for infinite configurations and for the projective plane. We exemplify the proof method by simpler and not so well known results that it is impossible to construct the unit length, or a given point, by compass and straightedge from nothing by means of classical arbitrary points. On the other hand we construct any given point by compass and straightedge from nothing by means of arbitrary points determined by horizontal segments. We quote a "proof" of Hilbert's theorem from the literature and explain why it is problematic. We rigorously prove Hilbert's theorem and present three variants of it, the last one for the projective plane.


title: Surfaces close to the Severi lines

abstract: In this Thesis we study surfaces of general type with maximal Albanese dimension for which the quantity $K_X^2-4\chi(\mathcal{O}_X)-4(q-2)$ vanishes or is "small", that is surfaces close to the Severi lines. Over the complex numbers, it is known that a surface $X$, provided that $K_X^2<\frac{9}{2}\chi(\mathcal{O}_X)$, has to satisfy the inequality $K_X^2-4\chi(\mathcal{O}_X)-4(q-2)\geq 0$. We give a constructive and complete classification of surfaces for which equality holds: these are surfaces whose canonical model is a double cover of an Abelian surface ($q=2$) or of a product elliptic surface ($q\geq 3$) branched over an ample divisor with at most negligible singularities which intersects the elliptic fibre twice in the latter case. We also prove, in the same hypothesis, that a surface $X$ with $K_X^2\neq 4\chi(\mathcal{O}_X)+4(q-2)$ satisfies $K_X^2\geq 4\chi(\mathcal{O}_X)+8(q-2)$ and we give a characterization of surfaces for which the equality holds. These are surfaces whose canonical model is a double cover of an isotrivial smooth elliptic surface fibration branched over an ample divisor with at most negligible singularities whose intersection with the elliptic fibre is $4$. Because these results are intimately related to theory of double covers, we see that their proof extend almost step by step to the case of any algebraically closed field of characteristic different from $2$. We also give some partial results over algebraically closed fields of characteristic $2$ after a study of double covers in that case.


title: Evaluating Effectiveness of Public Health Intervention Strategies for Mitigating COVID-19 Pandemic

abstract: Coronavirus disease 2019 (COVID-19) pandemic is an unprecedented global public health challenge. In the United States (US), state governments have implemented various non-pharmaceutical interventions (NPIs), such as physical distance closure (lockdown), stay-at-home order, mandatory facial mask in public in response to the rapid spread of COVID-19. To evaluate the effectiveness of these NPIs, we propose a nested case-control design with propensity score weighting under the quasi-experiment framework to estimate the average intervention effect on disease transmission across states. We further develop a method to test for factors that moderate intervention effect to assist precision public health intervention. Our method takes account of the underlying dynamics of disease transmission and balance state-level pre-intervention characteristics. We prove that our estimator provides causal intervention effect under assumptions. We apply this method to analyze US COVID-19 incidence cases to estimate the effects of six interventions. We show that lockdown has the largest effect on reducing transmission and reopening bars significantly increase transmission. States with a higher percentage of non-white population are at greater risk of increased $R_t$ associated with reopening bars.


title: Modeling accretion shocks at the disk-envelope interface -- Sulfur chemistry

abstract: As material from an infalling protostellar envelope hits the forming disk, an accretion shock may develop which could (partially) alter the envelope material entering the disk. Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) indicate that emission originating from warm SO and SO$_2$ might be good tracers of such accretion shocks. The goal of this work is to test under what shock conditions the abundances of gas-phase SO and SO$_2$ increase in an accretion shock at the disk-envelope interface. Detailed shock models including gas dynamics are computed using the Paris-Durham shock code for non-magnetized J-type accretion shocks in typical inner envelope conditions. The effect of pre-shock density, shock velocity, and strength of the ultraviolet (UV) radiation field on the abundance of warm SO and SO$_2$ is explored. Warm gas-phase chemistry is efficient in forming SO under most J-type shock conditions considered. In lower-velocity (~3 km/s) shocks, the abundance of SO is increased through subsequent reactions starting from thermally desorbed CH$_4$ toward H$_2$CO and finally SO. In higher velocity (>4 km/s) shocks, both SO and SO$_2$ are formed through reactions of OH and atomic S. The strength of the UV radiation field is crucial for SO and in particular SO$_2$ formation through the photodissociation of H$_2$O. Thermal desorption of SO and SO$_2$ ice is only relevant in high-velocity (>5 km/s) shocks at high densities ($>10^7$ cm$^{-3}$). Warm emission from SO and SO$_2$ is a possible tracer of accretion shocks at the disk-envelope interface as long as a local UV field is present. Additional observations with ALMA at high-angular resolution could provide further constraints. Moreover, the James Webb Space Telescope will give access to other possible slow, dense shock tracers such as H$_2$, H$_2$O, and [S I] 25$\mu$m.


title: Local wellposedness for the free boundary incompressible Euler equations with interfaces that exhibit cusps and corners of nonconstant angle

abstract: We prove that free boundary incompressible Euler equations are locally well posed in a class of solutions in which the interfaces can exhibit corners and cusps. Contrary to what happens in all the previously known non-$C^1$ water waves, the angle of these crests can change in time.


title: Minute-timescale variability in the X-ray emission of the highest redshift blazar

abstract: We report on two Chandra observations of the quasar PSO J0309+27, the most distant blazar observed so far (z=6.1), performed eight months apart, in March and November 2020. Previous Swift-XRT observation showed that this object is one of the brightest X-ray sources beyond redshift 6.0 ever observed so far. This new data-set confirmed the high flux level and unveiled a spectral change occurred on a very short timescale (250s rest-frame), caused by a significant softening of the emission spectrum. This kind of spectral variability, on a such short interval, has never been reported in the X-ray emission of a flat spectrum radio quasar. A possible explanation is given by the emission produced by the inverse Compton scatter of the quasar UV photons by the cold electrons present in a fast shell moving along the jet. Although this bulk comptonization emission should be an unavoidable consequence of the standard leptonic jet model, this would be the first time that it is observed.


title: Empowering End-users with Object-aware Processes

abstract: Business process management systems from various vendors are used by companies around the globe. Most of these systems allow for the full or partial automation of business processes by ensuring that tasks and data are presented to the right person at the right time during process execution. However, almost all established BPMS employ the activity-centric process support paradigm, in which the various forms, i.e., the main way for users to input data into the process, have to be created by hand. Furthermore, traditional activity-centric process management systems are limited in their flexibility as all possible execution variants have to be taken into account by the process modeler. Therefore, large amounts of research have gone into developing alternative process support paradigms, with a large focus on enabling more flexibly executable processes. This article takes one of these paradigms, object-aware process management, and presents the concepts we developed while researching the possibility of bringing the power and flexibility of non-activity-centric process support paradigms to the people that matter: the end-users working with the processes. The contribution of this article are the concepts, ideas, and lessons learned during the development and evaluation of the PHILharmonicFlows runtime user interface, which allows for the generation of an entire user interface, complete with navigation and forms, based on an object-aware process model. This novel approach allows for the generation of entire information systems, complete with data storage, process logic, and now fully functional user interfaces in a fully generic fashion from data-centric object-aware process models.


title: The spectrum of tensor perturbations in warm inflation

abstract: We compute the spectrum of tensor perturbations in warm inflation. We find that the spectrum, besides the standard component $\propto {H^2}/{M_P^2}$ associated to the amplification of the tensor vacuum fluctuations, acquires a component $\propto {\ell_{\rm mfp}\,T^5}/{M_P^4}$, where $\ell_{\rm mfp}$ and $T$ are respectively the mean free path and the temperature of the thermal degrees of freedom. The new contribution is due to the direct production of gravitational waves by the thermal bath, and can exceed the standard one in a viable region of parameter space. This contribution is dominated by thermal fluctuations at scales longer than $\ell_{\rm mfp}$.


title: Assessing the Cost of Network Simplifications in Long-Term Hydrothermal Dispatch Planning Models

abstract: The sustainable utilization of hydro energy relies on accurate estimates of the opportunity cost of the water. This value is calculated through long-term hydrothermal dispatch problems (LTHDP), and the recent literature has raised awareness about the consequences of modeling simplifications in these problems. The inaccurate representation of Kirchhoff's voltage law under the premise of a DC power flow is an example. Under a non-linear AC model, however, the LTHDP becomes intractable, and the literature lacks an accurate evaluation method of different modeling alternatives. In this paper, we extend the state-of-the-art cost-assessment framework of network approximations for LTHDP and bring relevant and practical new insights. First, we increase the quality of the assessment by using an AC power flow to simulate and compare the performance of five policies based on different network approximations. Second, we find that the tightest network relaxation (based on semidefinite programming) is not the one exhibiting the best performance. Results show that the DC {power flow} with quadratic losses approximation exhibits the lowest expected cost and inconsistency gaps. Finally, its computational burden is lower than that exhibited by the semidefinite relaxation, whereas market distortions are significantly reduced in comparison to previously published benchmarks based on DC power flow.


title: Decomposition of cubic graphs with cyclic connectivity 5

abstract: Let $G$ be a cyclically $5$-connected cubic graph with a $5$-edge-cut separating $G$ into two cyclic components $G_1$ and $G_2$. We prove that each component $G_i$ can be completed to a cyclically $5$-connected cubic graph by adding three vertices, unless $G_i$ is a cycle of length five. Our work extends similar results by Andersen et al. for cyclic connectivity $4$ from 1988.


title: Log-symmetric models with cure fraction with application to leprosy reactions data

abstract: In this paper, we propose a log-symmetric survival model with cure fraction, considering that the distributions of lifetimes for susceptible individuals belong to the log-symmetric class of distributions. This class has continuous, strictly positive, and asymmetric distributions, including the log-normal, log-$t$-Student, Birnbaum-Saunders, log-logistic I, log-logistic II, log-normal-contaminated, log-exponential-power, and log-slash distributions. The log-symmetric class is quite flexible and allows for including bimodal distributions and outliers. This includes explanatory variables through the parameter associated with the cure fraction. We evaluate the performance of the proposed model through extensive simulation studies and consider a real data application to evaluate the effect of factors on the immunity to leprosy reactions in patients with Hansen's disease.


title: Efficient (j,k)-Domination in Regular Graphs

abstract: Rubalcaba and Slater (Robert R. Rubalcaba and Peter J. Slater. Efficient (j,k)-domination. Discuss. Math. Graph Theory, 27(3):409-423, 2007.) define a $(j,k)$-dominating function on graph $X$ as a function $f:V(X)\rightarrow \{0,\ldots,j\}$ so that for each $v\in V(X)$, $f(N[v])\geq k$, where $N[v]$ is the closed neighbourhood of $v$. Such a function is efficient if all of the vertex inequalities are met with equality. They give a simple necessary condition for efficient domination, namely: if $X$ is an $r$-regular graph on $n$ vertices that has an efficient $(1,k)$-dominating function, then the size of the corresponding dominating set divides $n\cdot k$. The Hamming graph $H(q,d)$ is the graph on the vectors $\mathbb{Z}_q^d$ where two vectors are adjacent if and only if they are at Hamming distance $1$. We show that if $q$ is prime, then the previous necessary condition is sufficient for $H(q,d)$ to have an efficient $(1,k)$-dominating function. This result extends a result of Lee (Jaeun Lee. Independent perfect domination sets in Cayley graphs. J. Graph Theory, 37(4):213-219, 2001.) on independent perfect domination in Cayley graphs. We mention difficulties that arise for $H(q,d)$ when $q$ is a prime power but not prime.