title: Freeform microfluidic networks encapsulated in laser printed three-dimensional macro-scale glass objects

abstract: Large-scale microfluidic microsystems with complex three-dimensional (3D) configurations are highly in demand by both fundamental research and industrial application, holding the potentials for fostering a wide range of innovative applications such as lab-on-a-chip and organ-on-a-chip as well as continuous-flow manufacturing of fine chemicals. However, freeform fabrication of such systems remains challenging for most of the current fabrication techniques in terms of fabrication resolution, flexibility, and achievable footprint size. Here, we report ultrashort pulse laser microfabrication of freeform microfluidic circuits with high aspect ratios and tunable diameters embedded in 3D printed glass objects. We achieve uniform microfluidic channel diameter by carefully distributing a string of extra access ports along the microfluidic channels for avoiding the over-etching in the thin microfluidic channels. After the chemical etching is completed, the extra access ports are sealed using carbon dioxide laser induced localized glass melting. We demonstrate a model hand of fused silica with a size of ~3 cm * 2.7 cm * 1.1 cm in which the whole blood vessel system is encapsulated.


title: A formula to calculate the invariant $J$ of a quasi-homogeneous map germ

abstract: In this work, we consider a quasi-homogeneous, corank $1$, finitely determined map germ $f$ from $(\mathbb{C}^2,0)$ to $(\mathbb{C}^3,0)$. We consider the invariants $m(f(D(f))$ and $J$, where $m(f(D(f))$ denotes the multiplicity of the image of the double point curve $D(f)$ of $f$ and $J$ denotes the number of tacnodes that appears in a stabilization of the transversal slice curve of $f(\mathbb{C}^2)$. We present formulas to calculate $m(f(D(f))$ and $J$ in terms of the weights and degrees of $f$.


title: Chance-constrained optimal inflow control in hyperbolic supply systems with uncertain demand

abstract: In this paper, we address the task of setting up an optimal production plan taking into account an uncertain demand. The energy system is represented by a system of hyperbolic partial differential equations (PDEs) and the uncertain demand stream is captured by an Ornstein-Uhlenbeck process. We determine the optimal inflow depending on the producer's risk preferences. The resulting output is intended to optimally match the stochastic demand for the given risk criteria. We use uncertainty quantification for an adaptation to different levels of risk aversion. More precisely, we use two types of chance constraints to formulate the requirement of demand satisfaction at a prescribed probability level. In a numerical analysis, we analyze the chance-constrained optimization problem for the Telegrapher's equation and a real-world coupled gas-to-power network.


title: Orbital Hall Insulating Phase in Transition Metal Dichalcogenide Monolayers

abstract: We show that H-phase transition metal dichalcogenides (TMDs) monolayers such as MoS$_2$ and WSe$_2$, are orbital Hall insulators. They present very large orbital Hall conductivity plateaus in their semiconducting gap, where the spin Hall conductivity vanishes. Our results open the possibility of using TMDs for orbital current injection and orbital torque transfers that surpass their spin-counterparts in spin-orbitronics devices. The orbital Hall effect (OHE) in TMD monolayers occurs even in the absence of spin-orbit coupling. It can be linked to exotic momentum-space Dresselhaus-like orbital textures, analogous to the spin-momentum locking in 2D Dirac fermions that arise from a combination of orbital attributes and lattice symmetry.


title: Authentication Against a Myopic Adversary

abstract: We consider keyless authentication for point-to-point communication in the presence of a myopic adversary. In particular, the adversary has access to a non-causal noisy version of the transmission and may use this knowledge to choose the channel state of an arbitrarily-varying channel between legitimate users; the receiver is successful if it either decodes to the correct message or correctly detects adversarial interference. We show that a channel condition called U-overwritability, which allows the adversary to make its false message appear legitimate and untampered with, is a sufficient condition for zero authentication capacity. We present a useful way to compare adversarial channels, and show that once an AVC becomes U-overwritable, it remains U-overwritable for all "less myopic" adversaries. Finally, we show that stochastic encoders are necessary for positive authentication capacity in some cases, and examine in detail a binary adversarial channel that illustrates this necessity. Namely, for this binary channel, we show that when the adversarial channel is degraded with respect to the main channel between users, the no-adversary capacity of the underlying channel is achievable with a deterministic encoder. Otherwise, provided the channel to the adversary is not perfect, a stochastic encoder is necessary for positive authentication capacity; if such an encoder is allowed, the no-adversary capacity is again achievable.


title: A broadband look of the Accreting Millisecond X-ray Pulsar SAX J1748.9-2021 using AstroSat and XMM-Newton

abstract: SAX J1748.9-2021 is a transient accretion powered millisecond X-ray pulsar located in the Globular cluster NGC 6440. We report on the spectral and timing analysis of SAX J1748.9-2021 performed on AstroSat data taken during its faint and short outburst of 2017. We derived the best-fitting orbital solution for the 2017 outburst and obtained an average local spin frequency of 442.361098(3) Hz. The pulse profile obtained from 3-7 keV and 7-20 keV energy bands suggest constant fractional amplitude ~0.5% for fundamental component, contrary to previously observed energy pulse profile dependence. Our AstroSat observations revealed the source to be in a hard spectral state. The 1-50 keV spectrum from SXT and LAXPC on-board AstroSat can be well described with a single temperature blackbody and thermal Comptonization. Moreover, we found that the combined spectra from XMM-Newton (EPIC-PN) and AstroSat (SXT+LAXPC) indicated the presence of reflection features in the form of iron (Fe K${\alpha}$) line that we modeled with the reflection model xillvercp. One of the two X-ray burst observed during the AstroSat/LAXPC observation showed hard X-ray emission (>30 keV) due to Compton up-scattering of thermal photons by the hot corona. Time resolved analysis performed on the bursts revealed complex evolution in emission radius of blackbody for second burst suggestive of mild photospheric radius expansion.


title: Nine New Repeating Fast Radio Burst Sources from CHIME/FRB

abstract: We report on the discovery and analysis of bursts from nine new repeating fast radio burst (FRB) sources found using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources span a dispersion measure (DM) range of 195 to 1380 pc cm$^{-3}$. We detect two bursts from three of the new sources, three bursts from four of the new sources, four bursts from one new source, and five bursts from one new source. We determine sky coordinates of all sources with uncertainties of $\sim$10$^\prime$. We detect Faraday rotation measures for two sources, with values $-20(1)$ and $-499.8(7)$ rad m$^{-2}$, that are substantially lower than the RM derived from bursts emitted by FRB 121102. We find that the DM distribution of our events, combined with the nine other repeaters discovered by CHIME/FRB, is indistinguishable from that of thus far non-repeating CHIME/FRB events. However, as previously reported, the burst widths appear statistically significantly larger than the thus far non-repeating CHIME/FRB events, further supporting the notion of inherently different emission mechanisms and/or local environments. These results are consistent with previous work, though are now derived from 18 repeating sources discovered by CHIME/FRB during its first year of operation. We identify candidate galaxies that may contain FRB 190303.J1353+48 (DM = 222.4 pc cm$^{-3}$).


title: Optimizing spontaneous parametric down-conversion sources for boson sampling

abstract: An important step for photonic quantum technologies is the demonstration of a quantum advantage through boson sampling. In order to prevent classical simulability of boson sampling, the photons need to be almost perfectly identical and almost without losses. These two requirements are connected through spectral filtering, improving one leads to a decrease of the other. A proven method of generating single photons is spontaneous parametric downconversion (SPDC). We show that an optimal trade-off between indistinguishability and losses can always be found for SPDC. We conclude that a 50-photon scattershot boson-sampling experiment using SPDC sources is possible from a computational complexity point of view. To this end, we numerically optimize SPDC sources under the regime of weak pumping and with a single spatial mode.


title: N\'eron blowups and low-degree cohomological applications

abstract: We define dilatations of general schemes and study their basic properties. Dilatations of group schemes are -- in favorable cases -- again group schemes, called N\'eron blowups. We give two applications to their cohomology in degree zero (integral points) and degree one (torsors): we prove a canonical Moy-Prasad isomorphism that identifies the graded pieces in the congruent filtration of $G$ with the graded pieces in its Lie algebra $\mathfrak g$, and we show that many level structures on moduli stacks of $G$-bundles are encoded in torsors under N\'eron blowups of $G$.


title: Guesswork with Quantum Side Information

abstract: What is the minimum number of guesses needed on average to correctly guess a realization of a random variable? The answer to this question led to the introduction of the notion of a quantity called guesswork by Massey in 1994, which can be viewed as an alternate security criterion to entropy. In this paper, we consider the guesswork in the presence of quantum side information, and show that a general sequential guessing strategy is equivalent to performing a single measurement and choosing a guessing strategy from the outcome. We use this result to deduce entropic one-shot and asymptotic bounds on the guesswork in the presence of quantum side information, and to formulate a semi-definite program (SDP) to calculate the quantity. We evaluate the guesswork for a simple example involving the BB84 states, both numerically and analytically, and prove a continuity result that certifies the security of slightly imperfect key states when the guesswork is used as the security criterion.


title: GeV-TeV Counterparts of SS 433/W50 from Fermi-LAT and HAWC Observations

abstract: The extended jets of the microquasar SS 433 have been observed in optical, radio, X-ray, and recently very-high-energy (VHE) $\gamma$-rays by HAWC. The detection of HAWC $\gamma$-rays with energies as great as 25 TeV motivates searches for high-energy $\gamma$-ray counterparts in the Fermi-LAT data in the 100 MeV--300 GeV band. In this paper, we report on the first-ever joint analysis of Fermi-LAT and HAWC observations to study the spectrum and location of $\gamma$-ray emission from SS~433. Our analysis finds common emission sites of GeV-to-TeV $\gamma$-rays inside the eastern and western lobes of SS 433. The total flux above 1 GeV is $\sim 1\times10^{-10}\,\rm cm^{-2}\,s^{-1}$ in both lobes. The $\gamma$-ray spectrum in the eastern lobe is consistent with inverse-Compton emission by an electron population that is accelerated by jets. To explain both the GeV and TeV flux, the electrons need to have a soft intrinsic energy spectrum, or undergo a quick cooling process due to synchrotron radiation in a magnetized environment.


title: Gamma-ray emission in alpha-particle reactions with C, Mg, Si, Fe

abstract: Cross sections for the strongest gamma-ray emission lines produced in alpha-particle reactions with C, Mg, Si, Fe have been measured in the range E_alpha = 50 - 90 MeV at the center for proton therapy at the Helmholtz-Zentrum Berlin. Data for more than 60 different gamma-ray lines were determined, with particular efforts for lines that are in cross section compilations/evaluations with astrophysical purpose, and where data exist at lower projectile energies. The data are compared with predictions of a modern nuclear reaction code and cross-section curves of the latest evaluation for gamma-ray line emission in accelerated-particle interactions in solar flares.


title: Quasilinear systems of Jordan block type and the mKP hierarchy

abstract: We demonstrate that commuting quasilinear systems of Jordan block type are parametrised by solutions of the modified KP hierarchy. Systems of this form naturally occur as hydrodynamic reductions of multi-dimensional linearly degenerate dispersionless integrable PDEs.


title: Investigating a Deep Learning Method to Analyze Images from Multiple Gamma-ray Telescopes

abstract: Imaging atmospheric Cherenkov telescope (IACT) arrays record images from air showers initiated by gamma rays entering the atmosphere, allowing astrophysical sources to be observed at very high energies. To maximize IACT sensitivity, gamma-ray showers must be efficiently distinguished from the dominant background of cosmic-ray showers using images from multiple telescopes. A combination of convolutional neural networks (CNNs) with a recurrent neural network (RNN) has been proposed to perform this task. Using CTLearn, an open source Python package using deep learning to analyze data from IACTs, with simulated data from the upcoming Cherenkov Telescope Array (CTA), we implement a CNN-RNN network and find no evidence that sorting telescope images by total amplitude improves background rejection performance.


title: Missing Mass Concentration for Markov Chains

abstract: The problem of missing mass in statistical inference (posed by McAllester and Ortiz, NIPS'02; most recently revisited by Changa and Thangaraj, ISIT'2019) seeks to estimate the weight of symbols that have not been sampled yet from a source. So far all the approaches have been focused on the IID model which, although overly simplistic, is already not straightforward to tackle. The non-trivial part is in handling correlated events and sums of variables with very different scales where classical concentration inequalities do not yield good bounds. In this paper we develop the research on missing mass further, solving the problem for Markov chains. We reduce the problem to studying the tails of hitting times and finding \emph{log-additive approximations} to them. More precisely, we combine the technique of majorization and certain estimates on set hitting times to show how the problem can be eventually reduced back to the IID case. Our contribution are a) new technique to obtain missing mass bounds - we replace traditionally used negative association by majorization which works for a wider class of processes b) first (exponential) concentration bounds for missing mass in Markov chain models c) simplifications of recent results on set hitting times and d) simplified derivation of missing mass estimates for memory-less sources.


title: Identification and nonlinearity compensation of hysteresis using NARX models

abstract: This paper deals with two problems: the identification and compensation of hysteresis nonlinearity in dynamical systems using nonlinear polynomial autoregressive models with exogenous inputs (NARX). First, based on gray-box identification techniques, some constraints on the structure and parameters of NARX models are proposed to ensure that the identified models display a key-feature of hysteresis. In addition, a more general framework is developed to explain how hysteresis occurs in such models. Second, two strategies to design hysteresis compensators are presented. In one strategy the compensation law is obtained through simple algebraic manipulations performed on the identified models. It has been found that the compensators based on gray-box models outperform the cases with models identified using black-box techniques. In the second strategy, the compensation law is directly identified from the data. Both numerical and experimental results are presented to illustrate the efficiency of the proposed procedures.


title: Real-Time Long Range Trajectory Replanning for MAVs in the Presence of Dynamic Obstacles

abstract: Real-time long-range local planning is a challenging task, especially in the presence of dynamics obstacles. We propose a complete system which is capable of performing the local replanning in real-time. Desired trajectory is needed in the system initialization phase; system starts initializing sub-components of the system including point cloud processor, trajectory estimator and planner. Afterwards, the multi-rotary aerial vehicle starts moving on the given trajectory. When it detects obstacles, it replans the trajectory from the current pose to pre-defined distance incorporating the desired trajectory. Point cloud processor is employed to identify the closest obstacles around the vehicle. For replanning, Rapidly-exploring Random Trees (RRT*) is used with two modifications which allow planning the trajectory in milliseconds scales. Once we replanned the desired path, velocity components(x,y and z) and yaw rate are calculated. Those values are sent to the controller at a constant frequency to maneuver the vehicle autonomously. Finally, we have evaluated each of the components separately and tested the complete system in the simulated and real environments.


title: R-motivic stable stems

abstract: We compute some R-motivic stable homotopy groups. For $s - w \leq 11$, we describe the motivic stable homotopy groups $\pi_{s,w}$ of a completion of the R-motivic sphere spectrum. We apply the $\rho$-Bockstein spectral sequence to obtain R-motivic Ext groups from the C-motivic Ext groups, which are well-understood in a large range. These Ext groups are the input to the R-motivic Adams spectral sequence. We fully analyze the Adams differentials in a range, and we also analyze hidden extensions by $\rho$, 2, and $\eta$. As a consequence of our computations, we recover Mahowald invariants of many low-dimensional classical stable homotopy elements.


title: Calculation of the Ion-Ion Recombination Rate Coefficient via a Hybrid Continuum- Molecular Dynamics Approach

abstract: Accurate calculation of the ion-ion recombination rate coefficient has been of long-standing interest, as it controls the ion concentration in gas phase systems and in aerosols. We describe the development of a hybrid continuum-molecular dynamics approach to determine the ion-ion recombination rate coefficient. The approach is based on the limiting sphere method classically used for transition regime collision phenomena in aerosols. When ions are sufficiently far from one another, ion-ion relative motion is described by diffusion equations while within a critical distance, molecular dynamics (MD) simulations are used to model ion-ion motion. MD simulations are parameterized using the AMBER force-field as well as by considering partial charges on atoms. Ion-neutral gas collisions are modeled in two mutually exclusive cubic domains composed of 103 gas atoms each, which remain centered on the recombining ions throughout calculations. Example calculations are reported for NH4+ recombination with NO2- in He, across a pressure range from 10 kPa to 10,000 kPa. Excellent agreement is found in comparison of calculations to literature values for the 100 kPa recombination rate coefficient (1.0 x 10-12 m3 s-1) in He. We also recover the experimentally observed increase in recombination rate coefficient with pressure at sub-atmospheric pressures, and the observed decrease in recombination rate coefficient in the high pressure continuum limit. We additionally find that non-dimensionalized forms of rate coefficients are consistent with recently developed equations for the dimensionless charged particle-ion collision rate coefficient based on Langevin dynamics simulations.


title: Solving inverse-PDE problems with physics-aware neural networks

abstract: We propose a novel composite framework to find unknown fields in the context of inverse problems for partial differential equations (PDEs). We blend the high expressibility of deep neural networks as universal function estimators with the accuracy and reliability of existing numerical algorithms for partial differential equations as custom layers in semantic autoencoders. Our design brings together techniques of computational mathematics, machine learning and pattern recognition under one umbrella to incorporate domain-specific knowledge and physical constraints to discover the underlying hidden fields. The network is explicitly aware of the governing physics through a hard-coded PDE solver layer in contrast to most existing methods that incorporate the governing equations in the loss function or rely on trainable convolutional layers to discover proper discretizations from data. This subsequently focuses the computational load to only the discovery of the hidden fields and therefore is more data efficient. We call this architecture Blended inverse-PDE networks (hereby dubbed BiPDE networks) and demonstrate its applicability for recovering the variable diffusion coefficient in Poisson problems in one and two spatial dimensions, as well as the diffusion coefficient in the time-dependent and nonlinear Burgers' equation in one dimension. We also show that this approach is robust to noise.


title: On linear stability and syzygy stability for rank 2 linear series

abstract: In previous works, the authors investigated the relationships between linear stability of a generated linear series $|V|$ on a curve $C$, and slope stabillity of the vector bundle $M_{V,L} := \ker (V \otimes \mathcal{O}_C \to L)$. In particular, the second named author and L. Stoppino conjecture that, for a complete linear system $|L|$, linear (semi)stability is equivalent to slope (semi)stability of $M_V$, and the first and third named authors proved that this conjecture holds for hyperelliptic and for generic curves. In this work we provide a counterexample to this conjecture on any smooth plane curve of degree $7$.


title: FBI Transform in Gevrey Classes and Anosov Flows

abstract: An analytic FBI transform is built on compact manifolds without boundary, that satisfies all the expected properties. It enables the study of microlocal analytic and Gevrey regularity on such manifolds. This tool is then used to study the Ruelle spectrum of Anosov flows with Gevrey coefficients. In particular, finite order for the associated dynamical determinant is proved.


title: Auction-based Charging Scheduling with Deep Learning Framework for Multi-Drone Networks

abstract: State-of-the-art drone technologies have severe flight time limitations due to weight constraints, which inevitably lead to a relatively small amount of available energy. Therefore, frequent battery replacement or recharging is necessary in applications such as delivery, exploration, or support to the wireless infrastructure. Mobile charging stations (i.e., mobile stations with charging equipment) for outdoor ad-hoc battery charging is one of the feasible solutions to address this issue. However, the ability of these platforms to charge the drones is limited in terms of the number and charging time. This paper designs an auction-based mechanism to control the charging schedule in multi-drone setting. In this paper, charging time slots are auctioned, and their assignment is determined by a bidding process. The main challenge in developing this framework is the lack of prior knowledge on the distribution of the number of drones participating in the auction. Based on optimal second-price-auction, the proposed formulation, then, relies on deep learning algorithms to learn such distribution online. Numerical results from extensive simulations show that the proposed deep learning-based approach provides effective battery charging control in multi-drone scenarios.


title: A Deep Learning Approach Towards Prediction of Faults in Wind Turbines

abstract: With the rising costs of conventional sources of energy, the world is moving towards sustainable energy sources including wind energy. Wind turbines consist of several electrical and mechanical components and experience an enormous amount of irregular loads, making their operational behaviour at times inconsistent. Operations and Maintenance (O&M) is a key factor in monitoring such inconsistent behaviour of the turbines in order to predict and prevent any incipient faults which may occur in the near future. Machine learning has been applied to the domain of wind energy over the last decade for analysing, diagnosing and predicting wind turbine faults. In particular, we follow the idea of modelling a turbine's performance as a power curve where any power outputs that fall off the curve can be seen as performance errors. Existing work using this idea has used data from a turbine's Supervisory Control & Acquisition (SCADA) system to filter and analyse fault & alarm data using regression techniques. In contrast to previous work, we explore how deep learning can be applied to fault prediction from open access meteorological data only.


title: Raychaudhuri-based reconstruction of anisotropic Einstein-Maxwell equation in 1+3 covariant formalism of $f(R)$-gravity

abstract: Recently, a new strategy to the reconstruction of $f(R)$-gravity models based on the Raychaudhuri equation has been suggested by Choudhury et al. In this paper, utilizing this method, the reconstruction of anisotropic Einstein-Maxwell equation in the $1+3$ covariant formalism of $f(R)$-gravity is investigated in four modes: $i.$ Reconstruction from a negative constant deceleration parameter refereeing to an ever-accelerating universe; $ii.$ Reconstruction from a constant jerk parameter $j=1$ which recovers celebrated $\Lambda \text{CDM}$ mode of evolution; $iii.$ Reconstruction from a variable jerk parameter $j=Q(t)$; and $iv.$ Reconstruction from a slowly varying jerk parameter. Furthermore, two suggestions for enhancing the method are proposed.


title: Paradoxical Results of Long-Term Potentiation explained by Voltage-based Plasticity Rule

abstract: Experiments have shown that the same stimulation pattern that causes Long-Term Potentiation in proximal synapses, will induce Long-Term Depression in distal ones. In order to understand these, and other, surprising observations we use a phenomenological model of Hebbian plasticity at the location of the synapse. Our computational model describes the Hebbian condition of joint activity of pre- and post-synaptic neuron in a compact form as the interaction of the glutamate trace left by a presynaptic spike with the time course of the postsynaptic voltage. We test the model using experimentally recorded dendritic voltage traces in hippocampus and neocortex. We find that the time course of the voltage in the neighborhood of a stimulated synapse is a reliable predictor of whether a stimulated synapse undergoes potentiation, depression, or no change. Our model can explain the existence of different -- at first glance seemingly paradoxical -- outcomes of synaptic potentiation and depression experiments depending on the dendritic location of the synapse and the frequency or timing of the stimulation.


title: In Defense of Grid Features for Visual Question Answering

abstract: Popularized as 'bottom-up' attention, bounding box (or region) based visual features have recently surpassed vanilla grid-based convolutional features as the de facto standard for vision and language tasks like visual question answering (VQA). However, it is not clear whether the advantages of regions (e.g. better localization) are the key reasons for the success of bottom-up attention. In this paper, we revisit grid features for VQA, and find they can work surprisingly well - running more than an order of magnitude faster with the same accuracy (e.g. if pre-trained in a similar fashion). Through extensive experiments, we verify that this observation holds true across different VQA models (reporting a state-of-the-art accuracy on VQA 2.0 test-std, 72.71), datasets, and generalizes well to other tasks like image captioning. As grid features make the model design and training process much simpler, this enables us to train them end-to-end and also use a more flexible network design. We learn VQA models end-to-end, from pixels directly to answers, and show that strong performance is achievable without using any region annotations in pre-training. We hope our findings help further improve the scientific understanding and the practical application of VQA. Code and features will be made available.


title: L2-Stability Analysis of the SM-NLMS Algorithm

abstract: In this letter, we analyze, two properties, the local and the global robustness of the set-membership normalized least mean square (SM-NLMS) algorithm. We will show that the SM-NLMS algorithm has l2-stability. Indeed, the SM-NLMS algorithm never diverges; no matter how the parameters of the SM-NLMS algorithm has been selected. Ultimately, the numerical simulations corroborate the validity of our analysis.


title: Plasmon-polaritons and diffraction on the layer of asymmetric hyperbolic metamaterial

abstract: We consider the plasmon polaritons along a layer of hyperbolic metamaterial propagating in the plane of the anisotropy axis with an arbitrary its orientation. As a layer material, we use periodic plane-layered artificial medium or hyperbolic metamaterial of thin metal and dielectric layers and produce its homogenization. The conditions for the existence of fast, slow, leakage, gliding flowing, forward and backward plasmon-polaritons are found. The Fresnel formulas for the diffraction of a plane wave of arbitrary polarization on such a structure are obtained. The dispersion of plasmon-polaritons and plane wave diffraction are calculated.


title: Encode, Shuffle, Analyze Privacy Revisited: Formalizations and Empirical Evaluation

abstract: Recently, a number of approaches and techniques have been introduced for reporting software statistics with strong privacy guarantees. These range from abstract algorithms to comprehensive systems with varying assumptions and built upon local differential privacy mechanisms and anonymity. Based on the Encode-Shuffle-Analyze (ESA) framework, notable results formally clarified large improvements in privacy guarantees without loss of utility by making reports anonymous. However, these results either comprise of systems with seemingly disparate mechanisms and attack models, or formal statements with little guidance to practitioners. Addressing this, we provide a formal treatment and offer prescriptive guidelines for privacy-preserving reporting with anonymity. We revisit the ESA framework with a simple, abstract model of attackers as well as assumptions covering it and other proposed systems of anonymity. In light of new formal privacy bounds, we examine the limitations of sketch-based encodings and ESA mechanisms such as data-dependent crowds. We also demonstrate how the ESA notion of fragmentation (reporting data aspects in separate, unlinkable messages) improves privacy/utility tradeoffs both in terms of local and central differential-privacy guarantees. Finally, to help practitioners understand the applicability and limitations of privacy-preserving reporting, we report on a large number of empirical experiments. We use real-world datasets with heavy-tailed or near-flat distributions, which pose the greatest difficulty for our techniques; in particular, we focus on data drawn from images that can be easily visualized in a way that highlights reconstruction errors. Showing the promise of the approach, and of independent interest, we also report on experiments using anonymous, privacy-preserving reporting to train high-accuracy deep neural networks on standard tasks---MNIST and CIFAR-10.


title: Candidate Population III stellar complex at z=6.629 in the MUSE Deep Lensed Field

abstract: We discovered a strongly lensed (\mu >40) Lya emission at z=6.629 (S/N~18) in the MUSE Deep Lensed Field (MDLF) targeting the Hubble Frontier Field galaxy cluster MACS~J0416. Dedicated lensing simulations imply that the Lya emitting region necessarily crosses the caustic. The arc-like shape of the Lya extends 3 arcsec on the observed plane and is the result of two merged multiple images, each one with a de-lensed Lya luminosity L<~2.8 x 10^(40) erg/s arising from a confined region (< 150 pc effective radius). A spatially unresolved HST counterpart is barely detected at S/N~2 after stacking the near-infrared bands, corresponding to an observed(intrinsic) magnitude m_(1500)>~30.8(>~35.0). The inferred rest-frame Lya equivalent width is EWo > 1120 A if the IGM transmission is T(IGM)<0.5. The low luminosities and the extremely large Lya EWo match the case of a Population~III star complex made of several dozens stars (~ 10^4 Msun) which irradiate a HII region crossing the caustic. While the Lya and stellar continuum are among the faintest ever observed at this redshift, the continuum and the Lya emissions could be affected by differential magnification, possibly biasing the EWo estimate. The aforementioned tentative HST detection tend to favor a large EWo, making such a faint Pop~III candidate a key target for the James Webb Space Telescope and Extremely Large Telescopes.


title: Gravitational-wave detection and parameter estimation for accreting black-hole binaries and their electromagnetic counterpart

abstract: We study the impact of gas accretion on the orbital evolution of black-hole binaries initially at large separation in the band of the planned Laser Interferometer Space Antenna (LISA). We focus on two sources: (i)~stellar-origin black-hole binaries~(SOBHBs) that can migrate from the LISA band to the band of ground-based gravitational-wave observatories within weeks/months; and (ii) intermediate-mass black-hole binaries~(IMBHBs) in the LISA band only. Because of the large number of observable gravitational-wave cycles, the phase evolution of these systems needs to be modeled to great accuracy to avoid biasing the estimation of the source parameters. Accretion affects the gravitational-wave phase at negative ($-4$) post-Newtonian order, and is therefore dominant for binaries at large separations. If accretion takes place at the Eddington or at super-Eddington rate, it will leave a detectable imprint on the dynamics of SOBHBs. In optimistic astrophysical scenarios, a multiwavelength strategy with LISA and a ground-based interferometer can detect about $10$ (a few) SOBHB events for which the accretion rate can be measured at $50\%$ ($10\%$) level. In all cases the sky position can be identified within much less than $0.4\,{\rm deg}^2$ uncertainty. Likewise, accretion at $\gtrsim 10\%$ ($\gtrsim 100\%$) of the Eddington rate can be measured in IMBHBs up to redshift $z\approx 0.1$ ($z\approx 0.5$), and the position of these sources can be identified within less than $0.01\,{\rm deg}^2$ uncertainty. Altogether, a detection of SOBHBs or IMBHBs would allow for targeted searches of electromagnetic counterparts to black-hole mergers in gas-rich environments with future X-ray detectors (such as Athena) and radio observatories (such as SKA).


title: Evaluation of probabilistic photometric redshift estimation approaches for The Rubin Observatory Legacy Survey of Space and Time (LSST)

abstract: Many scientific investigations of photometric galaxy surveys require redshift estimates, whose uncertainty properties are best encapsulated by photometric redshift (photo-z) posterior probability density functions (PDFs). A plethora of photo-z PDF estimation methodologies abound, producing discrepant results with no consensus on a preferred approach. We present the results of a comprehensive experiment comparing twelve photo-z algorithms applied to mock data produced for The Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). By supplying perfect prior information, in the form of the complete template library and a representative training set as inputs to each code, we demonstrate the impact of the assumptions underlying each technique on the output photo-z PDFs. In the absence of a notion of true, unbiased photo-z PDFs, we evaluate and interpret multiple metrics of the ensemble properties of the derived photo-z PDFs as well as traditional reductions to photo-z point estimates. We report systematic biases and overall over/under-breadth of the photo-z PDFs of many popular codes, which may indicate avenues for improvement in the algorithms or implementations. Furthermore, we raise attention to the limitations of established metrics for assessing photo-z PDF accuracy; though we identify the conditional density estimate (CDE) loss as a promising metric of photo-z PDF performance in the case where true redshifts are available but true photo-z PDFs are not, we emphasize the need for science-specific performancemetrics.


title: Quantum embeddings for machine learning

abstract: Quantum classifiers are trainable quantum circuits used as machine learning models. The first part of the circuit implements a quantum feature map that encodes classical inputs into quantum states, embedding the data in a high-dimensional Hilbert space; the second part of the circuit executes a quantum measurement interpreted as the output of the model. Usually, the measurement is trained to distinguish quantum-embedded data. We propose to instead train the first part of the circuit -- the embedding -- with the objective of maximally separating data classes in Hilbert space, a strategy we call quantum metric learning. As a result, the measurement minimizing a linear classification loss is already known and depends on the metric used: for embeddings separating data using the l1 or trace distance, this is the Helstrom measurement, while for the l2 or Hilbert-Schmidt distance, it is a simple overlap measurement. This approach provides a powerful analytic framework for quantum machine learning and eliminates a major component in current models, freeing up more precious resources to best leverage the capabilities of near-term quantum information processors.


title: Quantifying many-body learning far from equilibrium with representation learning

abstract: Far-from-equilibrium many-body systems, from soap bubbles to suspensions to polymers, learn the drives that push them. This learning has been observed via thermodynamic properties, such as work absorption and strain. We move beyond these macroscopic properties that were first defined for equilibrium contexts: We quantify statistical mechanical learning with machine learning. Our toolkit relies on a structural parallel that we identify between far-from-equilibrium statistical mechanics and representation learning, which is undergone by neural networks that contain bottlenecks, including variational autoencoders. We train a variational autoencoder, via unsupervised learning, on configurations assumed by a many-body system during strong driving. We analyze the neural network's bottleneck to measure the many-body system's classification ability, memory capacity, discrimination ability, and novelty detection. Numerical simulations of a spin glass illustrate our technique. This toolkit exposes self-organization that eludes detection by thermodynamic measures, more reliably and more precisely identifying and quantifying learning by matter.


title: The Expansion of the Universe is Faster than Expected

abstract: The present rate of the expansion of our Universe, the Hubble constant, can be predicted from the cosmological model using measurements of the early Universe, or more directly measured from the late Universe. But as these measurements improved, a surprising disagreement between the two appeared. In 2019, a number of independent measurements of the late Universe using different methods and data provided consistent results making the discrepancy with the early Universe predictions increasingly hard to ignore. We review key advances realized by 2019: -- The local or late Universe measurement of the Hubble constant improved from 10% uncertainty twenty years ago to under 2% by the end of 2019. -- In 2019, multiple independent teams presented measurements with different methods and different calibrations to produce consistent results. -- These late Universe estimations disagree at 4$\sigma$ to 6$\sigma$ with predictions made from the Cosmic Microwave Background in conjunction with the standard cosmological model, a disagreement that is hard to explain or ignore.


title: Stratified disc wind models for the AGN broad-line region: ultraviolet, optical and X-ray properties

abstract: The origin, geometry and kinematics of the broad line region (BLR) gas in quasars and active galactic nuclei (AGN) are uncertain. We demonstrate that clumpy biconical disc winds illuminated by an AGN continuum can produce BLR-like spectra. We first use a simple toy model to illustrate that disc winds make quite good BLR candidates, because they are self-shielded flows and can cover a large portion of the ionizing flux-density ($\phi_H$-$n_H$) plane. We then conduct Monte Carlo radiative transfer and photoionization calculations, which fully account for self-shielding and multiple scattering in a non-spherical geometry. The emergent model spectra show broad emission lines with equivalent widths and line ratios comparable to those observed in AGN, provided that the wind has a volume filling factor of $f_V\lesssim0.1$. Similar emission line spectra are produced for a variety of wind geometries (polar or equatorial) and for launch radii that differ by an order of magnitude. The line emission arises almost exclusively from plasma travelling below the escape velocity, implying that `failed winds' are important BLR candidates. The behaviour of a line-emitting wind (and possibly any `smooth flow' BLR model) is similar to that of the locally optimally-emitting cloud (LOC) model originally proposed by Baldwin et al (1995), except that the gradients in ionization state and temperature are large-scale and continuous, rather than within or between distinct clouds. Our models also produce UV absorption lines and X-ray absorption features, and the stratified ionization structure can partially explain the different classes of broad absorption line quasars.


title: Nuclear Star Clusters

abstract: We review the current knowledge about nuclear star clusters (NSCs), the spectacularly dense and massive assemblies of stars found at the centers of most galaxies. Recent observational and theoretical work suggest that many NSC properties, including their masses, densities, and stellar populations vary with the properties of their host galaxies. Understanding the formation, growth, and ultimate fate of NSCs therefore is crucial for a complete picture of galaxy evolution. Throughout the review, we attempt to combine and distill the available evidence into a coherent picture of NSC evolution. Combined, this evidence points to a clear transition mass in galaxies of ~10^9 solar masses where the characteristics of nuclear star clusters change. We argue that at lower masses, NSCs are formed primarily from globular clusters that inspiral into the center of the galaxy, while at higher masses, star formation within the nucleus forms the bulk of the NSC. We also discuss the coexistence of NSCs and central black holes, and how their growth may be linked. The extreme densities of NSCs and their interaction with massive black holes lead to a wide range of unique phenomena including tidal disruption and gravitational wave events. Lastly, we review the evidence that many NSCs end up in the halos of massive galaxies stripped of the stars that surrounded them, thus providing valuable tracers of the galaxies' accretion histories.


title: High reddening patches in Gaia DR2. Possible artifacts or indication of star formation at the edge of the Galactic disk

abstract: Context: Deep GALEX UV data show that the extreme outskirts of some spiral galaxies are teeming with star formation. Such young stellar populations evolving so far away from the bulk of their host galaxies challenge our overall understanding of how star formation proceeds at galactic scales. It is at present unclear whether our own Milky Way may also exhibit ongoing and recent star formation beyond the conventional edge of the disk ($\sim 15$ kpc). Aims: Using \textit{Gaia} DR2 data, we aim to determine if such a population is present in the Galactic halo, beyond the nominal radius of the Milky Way disk. Methods: We studied the kinematics of \textit{Gaia} DR2 sources with parallax values between 1/60 and 1/30 milliarcseconds towards two regions that show abnormally high values of extinction and reddening; the results are compared with predictions from GALAXIA Galactic model. We also plotted the color-magnitude (CM) diagrams with heliocentric distances computed inverting the parallaxes, and studied the effects of the large parallax errors by Monte Carlo sampling. Results: The kinematics point towards a Galactic origin for one of the regions, while the provenance of the stars in the other is not clear. A spectroscopic analysis of some of the sources in the first region confirms that they are located in the halo. The CM diagram of the sources suggests that some of them are young.


title: Evidence for $\gtrsim{4}$ Gyr timescales of neutron star mergers from Galactic archaeology

abstract: The nucleosynthetic site of the rapid ($r$) neutron-capture process is currently being debated. The direct detection of the neutron star merger GW170817, through gravitational waves and electromagnetic radiation, has confirmed such events as important sources of the $r$-process elements. However, chemical evolution models are not able to reproduce the observed chemical abundances in the Milky Way when neutron star mergers are assumed to be the only $r$-process site and realistic time distributions of such events are taken into account. Now for the first time, we combine all the available observational evidence of the Milky Way and its dwarf galaxy satellites to show that the data can only be explained if there are (at least) two distinct $r$-process sites: a quick source with timescales comparable to core-collapse supernovae, $t_{quick}\lesssim10^8$ yr, and a delayed source with characteristic timescales $t_{delayed}\gtrsim4$ Gyr. The delayed $r$-process source most probably originates in neutron star mergers, as the timescale fits well with that estimated for GW170817. Given the short timescales of the quick source, it is likely associated with massive stars, though a specific fast-track channel for compact object mergers cannot be excluded at this point. Our approach demonstrates that only by looking at all the available data will we be able to solve the puzzle that is the $r$-process.


title: A fractional corner anomaly reveals higher-order topology

abstract: Spectral measurements of boundary localized in-gap modes are commonly used to identify topological insulators via the bulk-boundary correspondence. This can be extended to high-order topological insulators for which the most striking feature is in-gap modes at boundaries of higher co-dimension, e.g. the corners of a 2D material. Unfortunately, this spectroscopic approach is not always viable since the energies of the topological modes are not protected and they can often overlap the bulk bands, leading to potential misidentification. Since the topology of a material is a collective product of all its eigenmodes, any conclusive indicator of topology must instead be a feature of its bulk band structure, and should not rely on specific eigen-energies. For many topological crystalline insulators the key topological feature is fractional charge density arising from the filled bulk bands, but measurements of charge distributions have not been accessible to date. In this work, we experimentally measure boundary-localized fractional charge density of two distinct 2D rotationally-symmetric metamaterials, finding 1/4 and 1/3 fractionalization. We then introduce a new topological indicator based on collective phenomenology that allows unambiguous identification of higher-order topology, even in the absence of in-gap states. Finally, we demonstrate the higher-order bulk-boundary correspondence associated with this fractional feature by using boundary deformations to spectrally isolate localized corner modes where they were previously unobservable.


title: Reducible Fibers of Polynomial Maps

abstract: For a degree $n$ polynomial $f$ over the rationals, the elements in the fiber $f^{-1}(a)$ are of degree $n$ over $\mathbb Q$ for most rational values $a$ by Hilbert's irreducibility theorem. Determining the set of exceptional $a$'s without this property is a long standing open problem that is closely related to the Davenport--Lewis--Schinzel problem (1959) on reducibility of separated polynomials. As opposed to previous work which mostly concerns indecomposable $f$, we answer both problems for decomposable $f=f_1\circ\cdots\circ f_r$, as long as the indecomposable factors $f_i\in\mathbb Q[x]$ are of degree at least $5$ and are not $x^n$ or a Chebyshev polynomial composed with linear polynomials.


title: Deconfinement on axion domain walls

abstract: We study $SU(N_c)$ gauge theories with Dirac fermions in representations ${\cal{R}}$ of nonzero $N$-ality, coupled to axions. These theories have an exact discrete chiral symmetry, which has a mixed 't Hooft anomaly with general baryon-color-flavor backgrounds, called the "BCF anomaly" in arXiv:1909.09027. The infrared theory also has an emergent $\mathbb Z_{N_c}^{(1)}$ $1$-form center symmetry. We show that the BCF anomaly is matched in the infrared by axion domain walls. We argue that $\mathbb Z_{N_c}^{(1)}$ is spontaneously broken on axion domain walls, so that nonzero $N$-ality Wilson loops obey the perimeter law and probe quarks are deconfined on the walls. We give further support to our conclusion by using a calculable small-circle compactification to study the multi-scale structure of the axion domain walls and the microscopic physics of deconfinement on their worldvolume.


title: Does syntax need to grow on trees? Sources of hierarchical inductive bias in sequence-to-sequence networks

abstract: Learners that are exposed to the same training data might generalize differently due to differing inductive biases. In neural network models, inductive biases could in theory arise from any aspect of the model architecture. We investigate which architectural factors affect the generalization behavior of neural sequence-to-sequence models trained on two syntactic tasks, English question formation and English tense reinflection. For both tasks, the training set is consistent with a generalization based on hierarchical structure and a generalization based on linear order. All architectural factors that we investigated qualitatively affected how models generalized, including factors with no clear connection to hierarchical structure. For example, LSTMs and GRUs displayed qualitatively different inductive biases. However, the only factor that consistently contributed a hierarchical bias across tasks was the use of a tree-structured model rather than a model with sequential recurrence, suggesting that human-like syntactic generalization requires architectural syntactic structure.


title: The Cosmological Heavy Ion Collider: Fast Thermalization after Cosmic Inflation

abstract: Heavy-ion colliders have revealed the process of "fast thermalization". This experimental breakthrough has led to new theoretical tools to study the thermalization process at both weak and strong coupling. We apply this to the reheating epoch of inflationary cosmology, and the formation of a cosmological quark gluon plasma (QGP). We compute the thermalization time of the QGP at reheating, and find it is determined by the energy scale of inflation and the shear viscosity to entropy ratio $\eta/s$; or equivalently, the tensor-to-scalar ratio and the strong coupling constant at the epoch of thermalization. Thermalization is achieved near-instantaneously in low-scale inflation and in strongly coupled systems, and takes of order or less than a single e-fold of expansion for weakly-coupled systems or after high-scale inflation. We demonstrate that the predictions of inflation are robust to the physics of thermalization, and find a stochastic background of gravitational waves at frequencies accessible by interferometers, albeit with a small amplitude.


title: A Ray Tracing Technique for the Navigation on a Non-convex Pareto Front

abstract: A new interactive approach to navigate on approximations of in general non-convex but connected Pareto fronts is introduced. Given a finite number of precalculated representative Pareto-efficient solutions, an adapted Delaunay triangulation is generated. Based on interpolation and ray tracing techniques, real time navigation in the vicinity of Pareto-optimal solutions is made possible.


title: Further applications of the G function integral method

abstract: In our recent work we proposed a generalization of the beta integral method for derivation of the hypergeometric identities which can by analogy be termed "the G function integral method". In this paper we apply this technique to the cubic and the degenerate Miller-Paris transformations to get several new transformation and summation formulas for the generalized hypergeometric functions at a fixed argument. We further present an alternative approach for reducing the right hand sides resulting from our method to a single hypergeometric function which does not require the use of summation formulas.


title: Pinned QMA: The power of fixing a few qubits in proofs

abstract: What could happen if we pinned a single qubit of a system and fixed it in a particular state? First, we show that this can greatly increase the complexity of static questions -- ground state properties of local Hamiltonian problems with restricted types of terms. In particular, we show that the Pinned commuting and Pinned Stoquastic Local Hamiltonian problems are QMA complete. Second, we show that pinning a single qubit via often repeated measurements also results in universal quantum computation already with commuting and stoquastic Hamiltonians. Finally, we discuss variants of the Ground State Connectivity (GSCON) problem in light of pinning, and show that Stoquastic GSCON is QCMA complete. We hence identify a comprehensive picture of the computational power of pinning, reminiscent of the power of the one clean qubit model.


title: Can Giraffes Become Birds? An Evaluation of Image-to-image Translation for Data Generation

abstract: There is an increasing interest in image-to-image translation with applications ranging from generating maps from satellite images to creating entire clothes' images from only contours. In the present work, we investigate image-to-image translation using Generative Adversarial Networks (GANs) for generating new data, taking as a case study the morphing of giraffes images into bird images. Morphing a giraffe into a bird is a challenging task, as they have different scales, textures, and morphology. An unsupervised cross-domain translator entitled InstaGAN was trained on giraffes and birds, along with their respective masks, to learn translation between both domains. A dataset of synthetic bird images was generated using translation from originally giraffe images while preserving the original spatial arrangement and background. It is important to stress that the generated birds do not exist, being only the result of a latent representation learned by InstaGAN. Two subsets of common literature datasets were used for training the GAN and generating the translated images: COCO and Caltech-UCSD Birds 200-2011. To evaluate the realness and quality of the generated images and masks, qualitative and quantitative analyses were made. For the quantitative analysis, a pre-trained Mask R-CNN was used for the detection and segmentation of birds on Pascal VOC, Caltech-UCSD Birds 200-2011, and our new dataset entitled FakeSet. The generated dataset achieved detection and segmentation results close to the real datasets, suggesting that the generated images are realistic enough to be detected and segmented by a state-of-the-art deep neural network.


title: End-user experience with the SCoTSS Compton imager and directional survey spectrometer

abstract: The Silicon photomultiplier-based Compton Telescope for Safety and Security (SCoTSS) has been developed incorporating end-user requirements into the design process. The end-user group includes those responsible for mobile survey in the event of a radiological or nuclear accident, those responsible for radiation survey support to security operations at major events and at Canadian borders, as well as some of those responsible for Canadian defense applications. The SCoTSS development program has reached a technology readiness level of eight, and we are proceeding with field trials of the instrument in high-fidelity operational environments. Prospective end users have been involved in trial set up and execution, assuring applicability in their mission spaces. SCoTSS has been subject to trials involving hidden sources, heavily shielded sources, imager moving with respect to source, and complicated man-made surroundings. Our operators value high sensitivity for anomaly geolocation and mapping. End users also require an instrument which is capable of direction reconstruction in motion, as well as rapid imaging of a field of view. We have developed a "time to image" measure which allows for quantitative comparison of imagers of fundamentally different technology, where one design may have an advantage in terms of energy resolution and compactness and another design may have an advantage in terms of efficiency and cost effectiveness. We present here the performance of the SCoTSS imager in rapid direction finding. As well, we compare the time to image quantity for the SCoTSS imager and the H3D Polaris-H Quad imager where the data were taken under equivalent conditions. This quantitative measure of imaging performance can allow operators to make an informed choice of the design that meets their needs taking into consideration also weight and size as well as budgetary constraints.


title: Harry Kesten (1931-2019), a personal and scientific tribute

abstract: The mathematical achievements of Harry Kesten since the mid-1950s have revolutionized probability theory as a subject in its own right and in its associations with aspects of algebra, analysis, geometry, and statistical physics. Through his personality and scientific ability, he has framed the modern subject to a degree exceeded by no other. The impact of his work and personality is summarised in this memoir.


title: Unsupervised multi-modal Styled Content Generation

abstract: The emergence of deep generative models has recently enabled the automatic generation of massive amounts of graphical content, both in 2D and in 3D. Generative Adversarial Networks (GANs) and style control mechanisms, such as Adaptive Instance Normalization (AdaIN), have proved particularly effective in this context, culminating in the state-of-the-art StyleGAN architecture. While such models are able to learn diverse distributions, provided a sufficiently large training set, they are not well-suited for scenarios where the distribution of the training data exhibits a multi-modal behavior. In such cases, reshaping a uniform or normal distribution over the latent space into a complex multi-modal distribution in the data domain is challenging, and the generator might fail to sample the target distribution well. Furthermore, existing unsupervised generative models are not able to control the mode of the generated samples independently of the other visual attributes, despite the fact that they are typically disentangled in the training data. In this paper, we introduce UMMGAN, a novel architecture designed to better model multi-modal distributions, in an unsupervised fashion. Building upon the StyleGAN architecture, our network learns multiple modes, in a completely unsupervised manner, and combines them using a set of learned weights. We demonstrate that this approach is capable of effectively approximating a complex distribution as a superposition of multiple simple ones. We further show that UMMGAN effectively disentangles between modes and style, thereby providing an independent degree of control over the generated content.


title: A hyper luminous starburst at z=4.72 magnified by a lensing galaxy pair at z=1.48

abstract: [Abridged] We discovered in the Herschel Reference Survey an extremely bright IR source with $S_{500}$~120mJy (Red Virgo 4 - RV4). Based on IRAM/EMIR and IRAM/NOEMA detections of the CO(5-4), CO(4-3), and [CI] lines, RV4 is located at z=4.724, yielding a total observed L$_{IR}$ of 1.1+/-0.6x0$^{14}$L$_{\odot}$. At the position of the Herschel emission, three blobs are detected with the VLA at 10cm. The CO(5-4) line detection of each blob confirms that they are at the same redshift with the same line width, indicating that they are multiple images of the same source. In Spitzer and deep optical observations, two sources, High-z Lens 1 (HL1) West and HL1 East, are detected at the center of the three VLA/NOEMA blobs. These two sources are placed at z=1.48 with XSHOOTER spectra, suggesting that they could be merging and gravitationally lensing the emission of RV4. HL1 is the second most distant lens known to date in strong lensing systems. The Einstein radius of the lensing system is 2.2"+/-0.2 (20kpc). The high redshift of HL1 and the large Einstein radius are highly unusual for a strong lensing system. We present the ISM properties of the background source RV4. Different estimates of the gas depletion time yield low values suggesting that RV4 is a SB galaxy. Among all high-z SMGs, this source exhibits one of the lowest L$_{[CI]}$ to L$_{IR}$ ratios, 3.2+/-0.9x10$^{-6}$, suggesting an extremely short gas tdepl of only 14+/-5Myr. It also shows a relatively high L$_{[CI]}$ to L$_{CO(4-3)}$ ratio (0.7+/-0.2) and low L$_{CO(5-4)}$ to L$_{IR}$ ratio (only ~50% of the value expected for normal galaxies) hinting a low density of gas. Finally, we discuss that the short tdepl of RV4 can be explained by either a very high SFE, which is difficult to reconcile with major mergers simulations of high-z galaxies, or a rapid decrease of SF, which would bias the estimate of tdepl toward low value.


title: Nested polyhedra and indices of orbits of Coxeter groups of non-crystallographic type

abstract: The invariants of finite-dimensional representations of simple Lie algebras, such as even-degree indices and anomaly numbers, are considered in the context of the non-crystallographic finite reflection groups $H_2$, $H_3$ and $H_4$. Using a representation-orbit replacement, the definitions and properties of the indices are formulated for individual orbits of the examined groups. The indices of orders two and four of the tensor product of $k$ orbits are determined. Using the branching rules for the non-crystallographic Coxeter groups, the embedding index is defined similarly to the Dynkin index of a representation. Moreover, since the definition of the indices can be applied to any orbit of non-crystallographic type, the algorithm allowing to search for the orbits of smaller radii contained within any considered one is presented for the Coxeter groups $H_2$ and $H_3$. The geometrical structures of nested polytopes are exemplified.


title: Understanding and Mitigating the Security Risks of Content Inclusion in Web Browsers

abstract: Thanks to the wide range of features offered by web browsers, modern websites include various types of content such as JavaScript and CSS in order to create interactive user interfaces. Browser vendors also provided extensions to enhance web browsers with additional useful capabilities that are not necessarily maintained or supported by default. However, included content can introduce security risks to users of these websites, unbeknownst to both website operators and users. In addition, the browser's interpretation of the resource URLs may be very different from how the web server resolves the URL to determine which resource should be returned to the browser. The URL may not correspond to an actual server-side file system structure at all, or the web server may internally rewrite parts of the URL. This semantic disconnect between web browsers and web servers in interpreting relative paths (path confusion) could be exploited by Relative Path Overwrite (RPO). On the other hand, even tough extensions provide useful additional functionality for web browsers, they are also an increasingly popular vector for attacks. Due to the high degree of privilege extensions can hold, extensions have been abused to inject advertisements into web pages that divert revenue from content publishers and potentially expose users to malware. In this thesis, I propose novel research into understanding and mitigating the security risks of content inclusion in web browsers to protect website publishers as well as their users.


title: A Simple Multiple Integral Solution to the Broken Stick Problem

abstract: Regard the closed interval $[0,1]$ as a stick. Partition $[0,1]$ into $n+1$ different intervals $I_1, \ \dots \ , I_{n+1},$ where $n \geq 2,$ which represent smaller sticks. The classical Broken Stick problem asks to find the probability that the lengths of these smaller sticks can be the side lengths of a polygon with $n+1$ sides. We will show that this probability is $1-\frac{n+1}{2^{n}}$ by using multiple integration.


title: Extreme Software Defined Radio -- GHz in Real Time

abstract: Software defined radio is a widely accepted paradigm for design of reconfigurable modems. The continuing march of Moore's law makes real-time signal processing on general purpose processors feasible for a large set of waveforms. Data rates in the low Mbps can be processed on low-power ARM processors, and much higher data rates can be supported on large x86 processors. The advantages of all-software development (vs. FPGA/DSP/GPU) are compelling: much wider pool of talent, lower development time and cost, and easier maintenance and porting. However, very high-rate systems (above 100 Mbps) are still firmly in the domain of custom and semi-custom hardware (mostly FPGAs). In this paper we describe an architecture and testbed for an SDR that can be easily scaled to support over 3 GHz of bandwidth and data rate up to 10 Gbps. The paper covers a novel technique to parallelize typically serial algorithms for phase and symbol tracking, followed by a discussion of data distribution for a massively parallel architecture. We provide a brief description of a mixed-signal front end and conclude with measurement results. To the best of the author's knowledge, the system described in this paper is an order of magnitude faster than any prior published result.


title: Commuting Service Platform: Concept and Analysis

abstract: We propose and investigate the concept of commuting service platforms (CSP) that leverage emerging mobility services to provide commuting services and connect directly commuters (employees) and their worksites (employers). By applying the two-sided market analysis framework, we show under what conditions a CSP may present the two-sidedness. Both the monopoly and duopoly CSPs are then analyzed. We showhowthe price allocation, i.e., the prices charged to commuters and worksites, can impact the participation and profit of the CSPs. We also add demand constraints to the duopoly model so that the participation rates ofworksites and employees are (almost) the same. With demand constraints, the competition between the two CSPs becomes less intense in general. Discussions are presented on how the results and findings in this paper may help build CSP in practice and how to develop new, CSP-based travel demand management strategies.


title: Protective measurement of open quantum systems

abstract: We study protective quantum measurements in the presence of an environment and decoherence. We consider the model of a protectively measured qubit that also interacts with a spin environment during the measurement. We investigate how the coupling to the environment affects the two characteristic properties of a protective measurement, namely, (i) the ability to leave the state of the system approximately unchanged and (ii) the transfer of information about expectation values to the apparatus pointer. We find that even when the interaction with the environment is weak enough not to lead to appreciable decoherence of the initial qubit state, it causes a significant broadening of the probability distribution for the position of the apparatus pointer at the conclusion of the measurement. This washing out of the pointer position crucially diminishes the accuracy with which the desired expectation values can be measured from a readout of the pointer. We additionally show that even when the coupling to the environment is chosen such that the state of the system is immune to decoherence, the environment may still detrimentally affect the pointer readout.


title: Two-color flat-top solitonic pulses in $\chi^{(2)}$ optical microresonators via second harmonic generation

abstract: We studied numerically the generation of the coherent frequency combs at second harmonic generation in $\chi^{(2)}$ microresonators via conventional frequency scan method. It was demonstrated for the first time that under particular conditions it is possible to generate two-color flat-top solitonic pulses, platicons, using pump amplitude modulation or controllable mode interaction approach, if the signs of the group velocity coefficients at pump frequency and its second harmonic are opposite but absolute values of these coefficients are rather close. It was revealed that platicons may be observed on both sides of the linear microresonator resonance (at positive, as well as negative pump frequency detunings). For the efficient platicon excitation, one needs simultaneous accurate matching of both microresonator free spectral ranges and resonant eigenfrequencies. Platicon generation processes were simulated numerically, excitation conditions and platicon generation domains were found for different generation methods, and the properties of generated platicons were studied for the different combinations of the medium parameters.


title: Log-Linear Dynamical Systems

abstract: We present log-linear dynamical systems, a dynamical system model for positive quantities. We explain the connection to linear dynamical systems and show how convex optimization can be used to identify and control log-linear dynamical systems. We illustrate system identification and control with an example from predator-prey dynamics. We conclude by discussing potential applications of the proposed model.


title: Ein optimiertes Gl\"attungsverfahren motiviert durch eine technische Fragestellung (An optimized smoothing method motivated by a technical problem)

abstract: Ausgehend von einer konkreten technischen Fragestellung diskutieren wir in dieser Notiz die Anwendung verschiedener Gl\"attungsverfahren auf Datens\"atze mit vorgegebener Struktur. Wir stellen die Verfahren im Detail vor und besprechen die Vorund Nachteile. Insbesondere stellen wir hier die symmetrisierte exponentielle Gl\"attung vor, die ein sehr gutes Gl\"attungsverhalten mit einem hohen Ma\ss~an Symmetrieerhaltung kombiniert. (Based on a specific technical question, we discuss the application of various smoothing methods to data sets with a particular structure. We present the methods in detail and discuss their advantages and disadvantages. In particular, we present the symmetrized exponential smoothing that combines very good smoothing behavior with a high degree of symmetry-conservation)


title: Modifications of Prony's Method for the Recovery and Sparse Approximation of Generalized Exponential Sums

abstract: In this survey we describe some modifications of Prony's method. In particular, we consider the recovery of general expansions into eigenfunctions of linear differential operators of first order and show, how these expansions can be recovered from function samples using generalized shift operators. We derive an ESPRIT-like algorithm for the generalized recovery method and show, that this approach can be directly used to reconstruct classical exponential sums from non-equispaced data. Furthermore, we derive a modification of Prony's method for sparse approximation with exponential sums which leads to a non-linear least-squares problem.


title: Time delay of slow electrons by a diatomic molecule described by non-overlapping atomic potentials model

abstract: We study the elastic scattering of slow electrons by two-atomic molecule in the frame of non-overlapping atomic potentials model. The molecular continuum wave function is represented as a combination of a plane wave and two spherical s-waves, generated by the centers of atomic spheres. The asymptotic of this function determines in closed form the amplitude of elastic electron scattering. We show that this amplitude cannot be represented as a series of spherical functions. Therefore, it is impossible to use straightly the usual S-matrix methods to determine the scattering phases for non-spherical targets. We show that far from molecule the continuum wave function can be presented as an expansion in other than spherical orthonormal functions. The coefficients of this expansion determine the molecular scattering phases for non-spherical molecular systems. In such an approach, we calculate the Wigner times delay for slow electron scattered by two-atomic target.


title: Towards GAN Benchmarks Which Require Generalization

abstract: For many evaluation metrics commonly used as benchmarks for unconditional image generation, trivially memorizing the training set attains a better score than models which are considered state-of-the-art; we consider this problematic. We clarify a necessary condition for an evaluation metric not to behave this way: estimating the function must require a large sample from the model. In search of such a metric, we turn to neural network divergences (NNDs), which are defined in terms of a neural network trained to distinguish between distributions. The resulting benchmarks cannot be "won" by training set memorization, while still being perceptually correlated and computable only from samples. We survey past work on using NNDs for evaluation and implement an example black-box metric based on these ideas. Through experimental validation we show that it can effectively measure diversity, sample quality, and generalization.


title: Relativistic three-body effects in hierarchical triples

abstract: The hierarchical three-body problem has many applications in relativistic astrophysics, and can play an important role in the formation of the binary black hole mergers detected by LIGO/Virgo. However, many studies have only included relativistic corrections responsible for the precession of pericenter of the inner and outer binaries, neglecting relativistic interactions between the three bodies. We revisit this problem and develop a fully consistent derivation of the secular three-body problem to first post-Newtonian order. We start with the Einstein-Infeld-Hoffman equations for a three-body system and expand the accelerations as a power series in the ratio of the semi-major axes of the inner ($a_1$) and outer ($a_2$) binary. We then perform a post-Keplerian, two-parameter expansion of the single-orbit-averaged Lagrange planetary equations in $\delta = v^2/c^2$ and $\epsilon = a_1/a_2$ using the method of multiple scales. Using this method, we derive previously-indentified secular effects at $\delta \epsilon^{5/2}$ order that arise directly from the equations of motion. We also calculate new secular effects through $\delta \epsilon^4$ order that can lead to eccentricity growth over many Lidov-Kozai cycles when the tertiary is much more massive than the inner binary. In such cases, inclusion of these effects can substantially alter the evolution of three-body systems as compared to an analysis in which they are neglected. Careful analysis of post-Newtonian three-body effects will be important to understand the formation and properties of coalescing binaries that form via three-body dynamical processes.


title: Lorentz invariance relations for twist-3 quark distributions

abstract: We calculate twist-3 parton ditribution functions (PDFs) using cut and uncut diagrams. Uncut diagrams lead to a Dirac delta function term. No such term appears when cut diagrams are used. We show that a $\delta(x)$ is necessary to satisfy the Lorentz invariance relations of twist-3 PDFs, except for the Burkhardt-Cottingham sum rule in QCD.


title: Temporal patterns of dispersal-induced synchronization in population dynamics

abstract: The mechanisms and properties of synchronization of oscillating ecological populations attract attention because it is a fairly common phenomenon and because spatial synchrony may elevate a risk of extinction and may lead to other environmental impacts. Conditions for stable synchronization in a system of linearly coupled predator-prey oscillators have been considered in the past. However, the spatial dispersion coupling may be relatively weak and may not necessarily lead to a stable, complete synchrony. If the coupling between oscillators is too weak to induce a stable synchrony, oscillators may be engaged into intermittent synchrony, when episodes of synchronized dynamics are interspersed with the episodes of desynchronized dynamics. In the present study we consider the temporal patterning of this kind of intermittent synchronized dynamics in a system of two dispersal-coupled Rosenzweig-MacArthur predator-prey oscillators. We consider the properties of the distributions of durations of desynchronized intervals and their dependence on the model parameters. We show that the temporal patterning of synchronous dynamics (an ecological network phenomenon) may depend on the properties of individual predator-prey patch (individual oscillator) and may vary independently of the strength of dispersal. We also show that if the dynamics of predator is slow relative to the dynamics of the prey (a situation that may promote brief but large outbreaks), dispersal-coupled predator-prey oscillating populations exhibit numerous short desynchronizations (as opposed to few long desynchronizations) and may require weaker dispersal in order to reach strong synchrony.


title: Dominance Move calculation using a MIP approach for comparison of multi and many-objective optimization solution sets

abstract: Dominance move (DoM) is a binary quality indicator that can be used in multiobjective optimization. It can compare solution sets while representing some important features such as convergence, spread, uniformity, and cardinality. DoM has an intuitive concept and considers the minimum move of one set needed to weakly Pareto dominate the other set. Despite the aforementioned properties, DoM is hard to calculate. The original formulation presents an efficient and exact method to calculate it in a biobjective case only. This work presents a new approach to calculate and extend DoM to deal with three or more objectives. The idea is to use a mixed integer programming (MIP) approach to calculate DoM. Some initial experiments, in the biobjective space, were done to verify the model correctness. Furthermore, other experiments, using three, five, and ten objective functions were done to show how the model behaves in higher dimensional cases. Algorithms such as IBEA, MOEAD, NSGAIII, NSGAII, and SPEA2 were used to generate the solution sets, however any other algorithms could be used with DoM indicator. The results have confirmed the effectiveness of the MIP DoM in problems with more than three objective functions. Final notes, considerations, and future research are discussed to exploit some solution sets particularities and improve the model and its use for other situations.


title: Forecasting multiple functional time series in a group structure: an application to mortality

abstract: When modeling sub-national mortality rates, we should consider three features: (1) how to incorporate any possible correlation among sub-populations to potentially improve forecast accuracy through multi-population joint modeling; (2) how to reconcile sub-national mortality forecasts so that they aggregate adequately across various levels of a group structure; (3) among the forecast reconciliation methods, how to combine their forecasts to achieve improved forecast accuracy. To address these issues, we introduce an extension of grouped univariate functional time series method. We first consider a multivariate functional time series method to jointly forecast multiple related series. We then evaluate the impact and benefit of using forecast combinations among the forecast reconciliation methods. Using the Japanese regional age-specific mortality rates, we investigate one-step-ahead to 15-step-ahead point and interval forecast accuracies of our proposed extension and make recommendations.


title: Network of Steel: Neural Font Style Transfer from Heavy Metal to Corporate Logos

abstract: We introduce a method for transferring style from the logos of heavy metal bands onto corporate logos using a VGG16 network. We establish the contribution of different layers and loss coefficients to the learning of style, minimization of artefacts and maintenance of readability of corporate logos. We find layers and loss coefficients that produce a good tradeoff between heavy metal style and corporate logo readability. This is the first step both towards sparse font style transfer and corporate logo decoration using generative networks. Heavy metal and corporate logos are very different artistically, in the way they emphasize emotions and readability, therefore training a model to fuse the two is an interesting problem.


title: On stochastic It\^o processes with drift in $L_{d}$

abstract: For It\^o stochastic processes in $\mathbb{R}^{d}$ with drift in $L_{d}$ Aleksandrov's type estimates are established in the elliptic and parabolic settings. They are applied to estimating the resolvent operators of the corresponding elliptic and parabolic operators in $L_{p}$ and $L_{p+1}$, respectively, where $p\geq d$.


title: Topological Insulator Laser Using Valley-Hall Photonic Crystals

abstract: Topological photonics has recently been proved a robust framework for manipulating light. Active topological photonic systems, in particular, enable richer fundamental physics by employing nonlinear light-matter interactions, thereby opening a new landscape for applications such as topological lasing. Here we report an all-dielectric topological insulator laser scheme based on semiconductor cavities formed by topologically distinct Kagome photonic crystals. The proposed planar semiconductor Kagome lattice allows broadband edge states below the light line due to photonic valley hall effect in telecommunication region, which provides a new route to retrieve nontrivial photonic topology and to develop integrated topological systems for robust light generation and transport.


title: ReluDiff: Differential Verification of Deep Neural Networks

abstract: As deep neural networks are increasingly being deployed in practice, their efficiency has become an important issue. While there are compression techniques for reducing the network's size, energy consumption and computational requirement, they only demonstrate empirically that there is no loss of accuracy, but lack formal guarantees of the compressed network, e.g., in the presence of adversarial examples. Existing verification techniques such as Reluplex, ReluVal, and DeepPoly provide formal guarantees, but they are designed for analyzing a single network instead of the relationship between two networks. To fill the gap, we develop a new method for differential verification of two closely related networks. Our method consists of a fast but approximate forward interval analysis pass followed by a backward pass that iteratively refines the approximation until the desired property is verified. We have two main innovations. During the forward pass, we exploit structural and behavioral similarities of the two networks to more accurately bound the difference between the output neurons of the two networks. Then in the backward pass, we leverage the gradient differences to more accurately compute the most beneficial refinement. Our experiments show that, compared to state-of-the-art verification tools, our method can achieve orders-of-magnitude speedup and prove many more properties than existing tools.


title: The existence of a transverse universal knot

abstract: We prove that there is a knot $K$ transverse to $\xi_{std}$, the tight contact structure of $S^3$, such that every contact 3-manifold $(M, \xi)$ can be obtained as a contact covering branched along $K$. By contact covering we mean a map $\varphi: M \to S^3$ branched along $K$ such that $\xi$ is contact isotopic to the lifting of $\xi_{std}$ under $\varphi$.


title: Obtaining a Canonical Polygonal Schema from a Greedy Homotopy Basis with Minimal Mesh Refinement

abstract: Any closed manifold of genus g can be cut open to form a topological disk and then mapped to a regular polygon with 4g sides. This construction is called the canonical polygonal schema of the manifold, and is a key ingredient for many applications in graphics and engineering, where a parameterization between two shapes with same topology is often needed. The sides of the 4g-gon define on the manifold a system of loops, which all intersect at a single point and are disjoint elsewhere. Computing a shortest system of loops of this kind is NP-hard. A computationally tractable alternative consists in computing a set of shortest loops that are not fully disjoint in polynomial time, using the greedy homotopy basis algorithm proposed by Erickson and Whittlesey, and then detach them in post processing via mesh refinement. Despite this operation is conceptually simple, known refinement strategies do not scale well for high genus shapes, triggering a mesh growth that may exceed the amount of memory available in modern computers, leading to failures. In this paper we study various local refinement operators to detach cycles in a system of loops, and show that there are important differences between them, both in terms of mesh complexity and preservation of the original surface. We ultimately propose two novel refinement approaches: the former minimizes the number of new elements in the mesh, possibly at the cost of a deviation from the input geometry. The latter allows to trade mesh complexity for geometric accuracy, bounding deviation from the input surface. Both strategies are trivial to implement, and experiments confirm that they allow to realize canonical polygonal schemas even for extremely high genus shapes where previous methods fail.


title: Classification of Traffic Using Neural Networks by Rejecting: a Novel Approach in Classifying VPN Traffic

abstract: In this paper, we introduce a novel end-to-end traffic classification method to distinguish between traffic classes including VPN traffic in three layers of the Open Systems Interconnection (OSI) model. Classification of VPN traffic is not trivial using traditional classification approaches due to its encrypted nature. We utilize two well-known neural networks, namely multi-layer perceptron and recurrent neural network to create our cascade neural network focused on two metrics: class scores and distance from the center of the classes. Such approach combines extraction, selection, and classification functionality into a single end-to-end system to systematically learn the non-linear relationship between input and predicted performance. Therefore, we could distinguish VPN traffics from non-VPN traffics by rejecting the unrelated features of the VPN class. Moreover, we obtain the application type of non-VPN traffics at the same time. The approach is evaluated using the general traffic dataset ISCX VPN-nonVPN, and an acquired dataset. The results demonstrate the efficacy of the framework approach for encrypting traffic classification while also achieving extreme accuracy, $95$ percent, which is higher than the accuracy of the state-of-the-art models, and strong generalization capabilities.


title: Bound fermion states in pinned vortices in the surface states of a superconducting topological insulator: The Majorana bound state

abstract: By analytically solving the Bogoliubov-de Gennes equations we study the fermion bound states at the center of the core of a vortex in a two-dimensional superconductor. We consider three kinds of 2D superconducting models: (a) a standard type II superconductor in the mixed state with low density of vortex lines, (b) a superconductor with strong spin-orbit coupling locking the spin parallel to the momentum and (c) a superconductor with strong spin-orbit coupling locking the spin perpendicular to the momentum. The 2D superconducting states are induced via proximity effect between an $s$-wave superconductor and the surface states of a strong topological insulator. In case (a) the energy gap for the excitations is of order $\Delta_{\infty}^2/(2E_F)$, while for cases (b) and (c) a zero-energy Majorana state arises together with an equally spaced ($\Delta^2_{\infty}/E_F$) sequence of fermion excitations. The spin-momentum locking is key to the formation of the Majorana state. We present analytical expressions for the energy spectrum and the wave functions.


title: Classical-trajectory Monte Carlo calculations of differential electron emission in fast heavy-ion collisions with water molecules

abstract: A classical description of electron emission differential ionization cross sections for highly-charged high-velocity ions ($\sim$ 10 a.u.) impinging on water molecules is presented. We investigate the validity of the classical statistical mechanics description of ionization ($\hbar=0$ limit of quantum mechanics) in different ranges of electron emission energy and solid angle, where mechanisms such as soft and binary collisions are expected to contribute. The classical-trajectory Monte Carlo method is employed to calculate doubly and singly differential cross sections for C$^{6+}$, O$^{8+}$ and Si$^{13+}$ projectiles, and comparisons with Continuum Distorted Wave Eikonal Initial State theoretical results and with experimental data are presented. We implement a time-dependent screening effect in our model, in the spirit of mean-field theory to investigate its effect for highly charged projectiles. We also focus on the role of an accurate description of the molecular target by means of a three-center potential to show its effect on differential cross sections. Very good agreement with experiments is found at medium to high electron emission energies.


title: Understanding the atmospheric properties and chemical composition of the ultra-hot Jupiter HAT-P-7b II. Mapping the effects of gas kinetics

abstract: The atmospheres of ultra-hot Jupiters are commonly considered to be at thermochemical equilibrium. We aim to provide disequilibrium chemistry maps for a global understanding of the chemistry in HAT-P-7b's atmosphere and assess the importance of disequilibrium chemistry on UHJs. We apply a hierarchical modelling approach utilising 97 1D atmospheric profiles from 3D GCM of HAT-P-7b. For each 1D profile, we evaluate our kinetic cloud formation model consistently with the local gas-phase composition in chemical equilibrium. We then evaluate quenching results from a zeroth-order approximation in comparison to a kinetic gas-phase approach. We find that the zeroth-order approach of estimating quenching points agrees well with the full gas-kinetic modeling results. Chemical disequilibrium has the greatest effect on the nightside and morning abundance of species such as H, H$_2$O, CH$_4$, CO$_2$, HCN, and all C$_n$H$_m$ molecules; heavier C$_n$H$_m$ molecules are more affected by disequilibrium processes. CO abundance, however, is affected only marginally. While dayside abundances also notably change, those around the evening terminator of HAT-P-7b are the least affected by disequilibrium processes. The latter finding may partially explain the consistency of observed transmission spectra of UHJs with atmospheres in thermochemical equilibrium. Photochemistry only negligibly affects molecular abundances and quenching levels. In general, the quenching points of HAT-P-7b's atmosphere are at much lower pressures in comparison to the cooler hot-jupiters. We propose several avenues to look for the effect of disequilibrium processes on UHJs that are, in general, based on abundance and opacity measurements at different local times. It remains a challenge to completely disentangle this from the chemical effects of clouds and that of a primordial non-solar abundance.


title: Bayesian inference of the incompressibility, skewness and kurtosis of nuclear matter from empirical pressures in relativistic heavy-ion collisions

abstract: Within the Bayesian statistical framework we infer the incompressibility $K_0$, skewness $J_0$ and kurtosis $Z_0$ parameters of symmetric nuclear matter (SNM) at its saturation density $\rho_0$ using the constraining bands on the pressure in cold SNM in the density range of 1.3$\rho_0$ to 4.5$\rho_0$ from transport model analyses of kaon production and nuclear collective flow in relativistic heavy-ion collisions. As the default option assuming the $K_0$, $J_0$ and $Z_0$ have Gaussian prior probability distribution functions (PDFs) with the means and variances of $235\pm 30$, $-200\pm 200$ and $-146\pm 1728$ MeV, their posterior most probable values are narrowed down to 192$^{+12}_{-16}$ MeV, -180$^{+100}_{-110}$ MeV and 200$^{+250}_{-250}$ at 68\% confidence level, respectively. The results are largely independent of the prior PDFs of $J_0$ and $Z_0$ used. However, if one adopts the strong belief that the incompressibility $K_0$ has a uniform prior PDF within its absolute boundary of 220-260 MeV as one can find easily in the literature, the posterior most probable values of $K_0$, $J_0$ and $Z_0$ shift to $K_0=220^{+6}_{-0}$ MeV, $J_0=-390^{+60}_{-70}$ MeV and $Z_0=600^{+200}_{-200}$ MeV, respectively. While the posterior PDFs of the SNM EOS parameters depend somewhat on the prior PDF of $K_0$ used, the results from using different prior PDFs are qualitatively consistent. The uncertainties of all three parameters are significantly reduced especially for the $J_0$ and $Z_0$ parameters compared to their current values.


title: The bounds for the number of linear extensions via chain and antichain coverings

abstract: Let $(\mathcal{P},\leqslant)$ be a finite poset. Define the numbers $a_1,a_2,\ldots$ (respectively, $c_1,c_2,\ldots$) so that $a_1+\ldots+a_k$ (respectively, $c_1+\ldots+c_k$) is the maximal number of elements of $\mathcal{P}$ which may be covered by $k$ antichains (respectively, $k$ chains.) Then the number $e(\mathcal{P})$ of linear extensions of poset $\mathcal{P}$ is not less than $\prod a_i!$ and not more than $n!/\prod c_i!$. A corollary: if $\mathcal{P}$ is partitioned onto disjoint antichains of size $b_1,b_2, \ldots$, then $e(\mathcal{P})\geqslant \prod b_i!$.


title: Retouchdown: Adding Touchdown to StreetLearn as a Shareable Resource for Language Grounding Tasks in Street View

abstract: The Touchdown dataset (Chen et al., 2019) provides instructions by human annotators for navigation through New York City streets and for resolving spatial descriptions at a given location. To enable the wider research community to work effectively with the Touchdown tasks, we are publicly releasing the 29k raw Street View panoramas needed for Touchdown. We follow the process used for the StreetLearn data release (Mirowski et al., 2019) to check panoramas for personally identifiable information and blur them as necessary. These have been added to the StreetLearn dataset and can be obtained via the same process as used previously for StreetLearn. We also provide a reference implementation for both of the Touchdown tasks: vision and language navigation (VLN) and spatial description resolution (SDR). We compare our model results to those given in Chen et al. (2019) and show that the panoramas we have added to StreetLearn fully support both Touchdown tasks and can be used effectively for further research and comparison.


title: Mean Field Quarks on the Light Front

abstract: We present a new approach for the calculation of the valence quark distributions in the nucleon based on the scenario in which the spectrum of the valence quarks at x>0.05 is generated through three main mechanisms: interaction of valence quarks with the mean field generated by the residual nucleon system, two and three quark short range interactions through gluon exchanges. In the current report we present the first phase of the project in which we develop a non-perturbative model for valence quark interaction in the mean field of the nucleonic interior to describe their distribution in the moderate x region (0.05 < x < 0.4). The short range quark-quark interaction effects in our approach generate the high x tail of valence quark distributions. The presented non-perturbative model is based on the picture in which three relativistic valence quarks occupy the nucleon core at distances of $\le 0.5$ Fm while interacting in the mean field generated by the residual nucleon system. The calculations are based on the assumption of the a factorization of the internal interaction of short-range three valence quarks with the long-range interaction of these quarks with the residual system. The theoretical approach is based on effective light-front diagrammatic approach which allows us to introduce the valence quark and residual system wave functions in a consistent way The parameters of these wave functions are fixed by the position of the peak of the xf_q(x) distribution of valence quarks at Q_0 corresponding to the charm-quark mass. With few parameters we achieved a very reasonable description of the up and down valence quark distributions in the moderate x region (x < 0.4), where one expects the mean field dynamics to dominate. The model, however, systematically underestimates the high $x$ region where enhanced contributions from partonic short-range correlations are expected.