title: On statistical models on super trees

abstract: We consider a particular example of interplay between statistical models related to CFT on one hand, and to the spectral properties of ODE, known as ODE/IS correspondence, on the other hand. We focus at the representation of wave functions of Schr\"odinger operators in terms of spectral properties of associated transfer matrices on "super trees" (the trees whose vertex degree changes with the distance from the root point). Such trees with varying branchings encode the structure of the Fock space of the model. We discuss basic spectral properties of "averaged random matrix ensembles" in terms of Hermite polynomials for the transfer matrix of super trees. At small "branching velocities" we have related the problem of paths counting on super trees to the statistics of area-weighted one-dimensional Dyck paths. We also discuss the connection of the spectral statistics of random walks on super trees with the Kardar-Parisi-Zhang scaling.

categories: math-ph cond-mat.stat-mech hep-th math.MP


title: Do Two Symmetry Breaking Transitions in Photosynthetic Light Harvesting Complexes Form One, Two or More Kibble Zurek Model Topological Defects?

abstract: Kibble and Zurek proposed that rapid symmetry breaking transitions in the hot, early universe could result in causally disconnected topological defects such as cosmic strings. This type of first order transition has analogues in certain second order transitions present in condensed matter such as liquid crystals, super fluids, and charge density waves in terms of flux tubes or vortices. Recently, we discovered that Rhodopseudomonas acidophilus photosynthetic light harvesting complex might have different types of coherent ground and excited states, suggesting that there are two different symmetry breaking transitions. The B 850 ground states comprise eight identical rings each containing 18 bacteriochlorophyll components, and each ring has undergone a Bose Einstein phase transition to a charge density wave that lowers the energy. The excited state coherence results from polariton formation from the non-crossing of bosons, here an extension of exciton theory. The result is short-lived quasi-particles with very low mass that can form an unusual BEC. We suggest the oriented, circular B 850 and enclosed singlet B 875 compounds create a new cavity structure with some attributes of a nano pillar. Since both the ground and excited states should contain solitons, we envisage three fast light pulse experiments could be able to map both the Kibble Zurek Model phase transitions and energy transfers as a function of light intensity and time in this complex at room temperature.


title: Open-Access Full-Duplex Wireless in the ORBIT Testbed

abstract: In order to support experimentation with full-duplex (FD) wireless, we recently integrated an open-access FD transceiver in the ORBIT testbed. In this report, we present the design and implementation of the FD transceiver and interfaces, and provide examples and guidelines for experimentation. In particular, an ORBIT node with a National Instruments (NI)/Ettus Research Universal Software Radio Peripheral (USRP) N210 software-defined radio (SDR) was equipped with the Columbia FlexICoN Gen-1 customized RF self-interference (SI) canceller box. The RF canceller box includes an RF SI canceller that is implemented using discrete components on a printed circuit board (PCB) and achieves 40dB RF SI cancellation across 5MHz bandwidth. We provide an FD transceiver baseline program and present two example FD experiments where 90dB and 85dB overall SI cancellation is achieved for a simple waveform and PSK modulated signals across both the RF and digital domains. We also discuss potential FD wireless experiments that can be conducted based on the implemented open-access FD transceiver and baseline program.


title: On Frobenius semisimplicity in Hida families

abstract: Let $p\geq 5$ be a prime and $\ell\neq p$ be a prime not dividing the tame level of a $p$-ordinary Hida family. We prove that the actions of the Frobenius element at $\ell$ on the Galois representations attached to almost all arithmetic specializations are semisimple and non-scalar. If $k_f$ denotes the weight of a cusp form $f(z)= \sum_{n\geq 1} a_\ell(f) e^{2\pi i n z}$, then the inequality $$|a_\ell(f) | \leq 2 \ell^{(k_f-1)/2},$$ predicted by the Ramanujan conjecture, is a strict inequality for almost all members $f$ of the family.


title: Bell monogamy relations in arbitrary qubit networks

abstract: Characterizing trade-offs between simultaneous violations of multiple Bell inequalities in a large network of qubits is computationally demanding. We propose a graph-theoretic approach to efficiently produce Bell monogamy relations in arbitrary arrangements of qubits. All the relations obtained for bipartite Bell inequalities are tight and leverage only a single Bell monogamy relation. This feature is unique to bipartite Bell inequalities, as we show that there is no finite set of such elementary monogamy relations for multipartite inequalities. Nevertheless, many tight monogamy relations for multipartite inequalities can be obtained with our method as shown in explicit examples.


title: A one-phase interior point method for nonconvex optimization

abstract: The work of Wachter and Biegler suggests that infeasible-start interior point methods (IPMs) developed for linear programming cannot be adapted to nonlinear optimization without significant modification, i.e., using a two-phase or penalty method. We propose an IPM that, by careful initialization and updates of the slack variables, is guaranteed to find a first-order certificate of local infeasibility, local optimality or unboundedness of the (shifted) feasible region. Our proposed algorithm differs from other IPM methods for nonconvex programming because we reduce primal feasibility at the same rate as the barrier parameter. This gives an algorithm with more robust convergence properties and closely resembles successful algorithms from linear programming. We implement the algorithm and compare with IPOPT on a subset of CUTEst problems. Our algorithm requires a similar median number of iterations, but fails on only 9% of the problems compared with 16% for IPOPT. Experiments on infeasible variants of the CUTEst problems indicate superior performance for detecting infeasibility. The code for our implementation can be found at https://github.com/ohinder/OnePhase .


title: XYZ - SU3 Breakings from Laplace Sum Rule at Higher Orders: Summary

abstract: This talk reviews and summarizes some of our results in [1] on XYZ- SU3 Breakings obtained from QCD Laplace Sum Rules (LSR) at next-to-next-leading order (N2LO) of perturbative (PT) theory and including next-to-leading order (NLO) SU3 breaking corrections and leading order (LO) contributions of dimensions d< (6 - 8) non-perturbative condensates. We conclude that the observed X states are good candidates for being 1^++ and 0^++ molecules states. We find that the SU3 breakings are relatively small for the masses (< 10 (resp. 3) %) for the charm (resp. bottom) channels while they are large (< 20 %) for the couplings. Like in the chiral limit case, the couplings decrease faster: 1/m_b^3/2 than 1/m_b^1/2 of HQET. Our approach cannot clearly separate ( within the errors ) some molecule states from the four-quark ones with the same quantum numbers.


title: Game of Drones - Detecting Streamed POI from Encrypted FPV Channel

abstract: Drones have created a new threat to people's privacy. We are now in an era in which anyone with a drone equipped with a video camera can use it to invade a subject's privacy by streaming the subject in his/her private space over an encrypted first person view (FPV) channel. Although many methods have been suggested to detect nearby drones, they all suffer from the same shortcoming: they cannot identify exactly what is being captured, and therefore they fail to distinguish between the legitimate use of a drone (for example, to use a drone to film a selfie from the air) and illegitimate use that invades someone's privacy (when the same operator uses the drone to stream the view into the window of his neighbor's apartment), a distinction that in some cases depends on the orientation of the drone's video camera rather than on the drone's location. In this paper we shatter the commonly held belief that the use of encryption to secure an FPV channel prevents an interceptor from extracting the POI that is being streamed. We show methods that leverage physical stimuli to detect whether the drone's camera is directed towards a target in real time. We investigate the influence of changing pixels on the FPV channel (in a lab setup). Based on our observations we demonstrate how an interceptor can perform a side-channel attack to detect whether a target is being streamed by analyzing the encrypted FPV channel that is transmitted from a real drone (DJI Mavic) in two use cases: when the target is a private house and when the target is a subject.


title: Rayleigh fractionation in high-Rayleigh-number solutal convection in porous media

abstract: We study the fractionation of two components between a well-mixed gas and a saturated convecting porous layer. Motivated by geological carbon dioxide (CO$_2$) storage we assume that convection is driven only by the dissolved concentration of the first component, while the second acts as a tracer with increased diffusivity. Direct numerical simulations for convection at high Rayleigh numbers reveal that the partitioning of the components, in general, does not follow a Rayleigh fractionation trend, as commonly assumed. Initially, increases in tracer diffusivity also increase its flux, because the diffusive boundary layer penetrates deeper into the flow. However, for $D_2\geq 10\, D_1$, where $D_1$ and $D_2$ are, respectively, the diffusion coefficients of CO$_2$ and the tracer in water, the transverse leakage of tracer between up- and down-welling plumes reduces the tracer flux. Rayleigh fractionation between components is only realized in the limit of two gases with very large differences in solubility and initial concentration in the gas.


title: Unidirectional spin density wave state in metallic (Sr1-xLax)2IrO4

abstract: Materials that exhibit both strong spin orbit coupling and electron correlation effects are predicted to host numerous new electronic states. One prominent example is the Jeff =1/2 Mott state in Sr2IrO4, where introducing carriers is predicted to manifest high temperature superconductivity analogous to the S=1/2 Mott state of La2CuO4. While bulk superconductivity currently remains elusive, anomalous quasi-particle behaviors paralleling those in the cuprates such as pseudogap formation and the formation of a d-wave gap are observed upon electron-doping Sr2IrO4. Here we establish a magnetic parallel between electron-doped Sr2IrO4 and hole-doped La2CuO4 by unveiling a spin density wave state in electron-doped Sr2IrO4. Our magnetic resonant x-ray scattering data reveal the presence of an incommensurate magnetic state reminiscent of the diagonal spin density wave state observed in the monolayer cuprate (La1-xSrx)2CuO4. This link supports the conjecture that the quenched Mott phases in electron-doped Sr2IrO4 and hole-doped La2CuO4 support common competing electronic phases.


title: The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample: Anisotropic Baryon Acoustic Oscillations measurements in Fourier-space with optimal redshift weights

abstract: We present a measurement of the anisotropic and isotropic Baryon Acoustic Oscillations (BAO) from the extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample with optimal redshift weights. Applying the redshift weights improves the constraint on the BAO dilation parameter $\alpha(z_{\rm eff})$ by 17\%. We reconstruct the evolution history of the BAO distance indicators in the redshift range of $0.8<z<2.2$. This paper is part of a set that analyses the eBOSS DR14 quasar sample.


title: A converse to Schreier's index-rank formula

abstract: In "Subgroups of free profinite groups and large subfields of Q" (Israel J. Math. 39 (1981), no. 1-2, pages 25-45; MR 617288) A. Lubotzky and L. van den Dries raise the question whether a finitely generated, residually finite group is necessarily free if the rank function on its subgroups of finite index satisfies Schreier's well-known index rank relation (see Question 2 on p. 34). I answered this question in 1980 but, so far, I have not published my answer. This note fills the omission; it gives an amended and abridged version of my original proof.


title: Asymmetry Hurts: Private Information Retrieval Under Asymmetric Traffic Constraints

abstract: We consider the classical setting of private information retrieval (PIR) of a single message (file) out of $M$ messages from $N$ distributed databases under the new constraint of \emph{asymmetric traffic} from databases. In this problem, the \emph{ratios between the traffic} from the databases are constrained, i.e., the ratio of the length of the answer string that the user (retriever) receives from the $n$th database to the total length of all answer strings from all databases is constrained to be $\tau_n$. This may happen if the user's access to the databases is restricted due database availability, channel quality to the databases, and other factors. For this problem, for fixed $M$, $N$, we develop a general upper bound $\bar{C}(\boldsymbol{\tau})$, which generalizes the converse proof of Sun-Jafar, where database symmetry was inherently used. Our converse bound is a piece-wise affine function in the traffic ratio vector $\boldsymbol{\tau}=(\tau_1, \cdots, \tau_N)$. For the lower bound, we explicitly show the achievability of $\binom{M+N-1}{M}$ corner points. For the remaining traffic ratio vectors, we perform time-sharing between these corner points. The recursive structure of our achievability scheme is captured via a system of difference equations. The upper and lower bounds exactly match for $M=2$ and $M=3$ for any $N$ and any $\boldsymbol{\tau}$. The results show strict loss of PIR capacity due to the asymmetric traffic constraints compared with the symmetric case of Sun-Jafar which implicitly uses $\tau_n=\frac{1}{N}$ for all $n$.


title: Quadrature Compound: An approximating family of distributions

abstract: Compound distributions allow construction of a rich set of distributions. Typically they involve an intractable integral. Here we use a quadrature approximation to that integral to define the quadrature compound family. Special care is taken that this approximation is suitable for computation of gradients with respect to distribution parameters. This technique is applied to discrete (Poisson LogNormal) and continuous distributions. In the continuous case, quadrature compound family naturally makes use of parameterized transformations of unparameterized distributions (a.k.a "reparameterization"), allowing for gradients of expectations to be estimated as the gradient of a sample mean. This is demonstrated in a novel distribution, the diffeomixture, which is is a reparameterizable approximation to a mixture distribution.


title: A proof of Bertrand's theorem using the theory of isochronous potentials

abstract: We give an alternative proof for the celebrated Bertrand's theorem as a corollary of the isochronicity of a certain family of centers.


title: Rational points and prime values of polynomials in moderately many variables

abstract: We derive the Hasse principle and weak approximation for pencils of certain varieties in the spirit of work by Colliot-Th\'el\`ene,Sansuc and Harpaz-Skorobogatov-Wittenberg. Our varieties are defined through polynomials in many variables and part of our work is devoted to establishing Schinzel's hypothesis for polynomials of this kind. This last part is achieved by using arguments behind Birch's well-known result regarding the Hasse principle for complete intersections with the notable difference that we prove our result in 50% fewer variables than in the classical Birch setting. We also study the problem of square-free values of an integer polynomial with 66.6% fewer variables than in the Birch setting.


title: Smoluchowski's discrete coagulation equation with forcing

abstract: In this article we study an extension of Smoluchowski's discrete coagulation equation, where particle in- and output takes place. This model is frequently used to describe aggregation processes in combination with sedimentation of clusters. More precisely, we show that the evolution equation is well-posed for a large class of coagulation kernels and output rates. Additionally, in the long-time limit we prove that solutions converge to a unique equilibrium with exponential rate under a suitable smallness condition on the coefficients.


title: A Langevin approach to multi-scale modeling

abstract: In plasmas, distribution functions often demonstrate long anisotropic tails or otherwise significant deviations from local Maxwellians. The tails, especially if they are pulled out from the bulk, pose a serious challenge for numerical simulations as resolving both the bulk and the tail on the same mesh is often challenging. A multi-scale approach, providing evolution equations for the bulk and the tail individually, could offer a resolution in the sense that both populations could be treated on separate meshes, or different reduction techniques applied to the bulk and the tail population. In this letter, we propose a multi-scale method which allows us to split a distribution function into a bulk and a tail so that both populations remain genuine, non-negative distribution functions and may carry density, momentum, and energy. The proposed method is based on the observation that the motion of an individual test particle in a plasma obeys a stochastic differential equation, also referred to as a Langevin equation. This allows us to define transition probabilities between the bulk and the tail and to provide evolution equations for both populations separately.


title: The Mean Ultraviolet Spectrum of a Representative Sample of Faint z~3 Lyman Alpha Emitters

abstract: We discuss the rest-frame ultraviolet emission line spectra of a large (~100) sample of low luminosity redshift z~3.1 Lyman alpha emitters (LAEs) drawn from a Subaru imaging survey in the SSA22 survey field. Our earlier work based on smaller samples indicated that such sources have high [OIII]/[OII] line ratios possibly arising from a hard ionising spectrum that may be typical of similar sources in the reionisation era. With optical spectra secured from VLT/VIMOS, we re-examine the nature of the ionising radiation in a larger sample using the strength of the high ionisation diagnostic emission lines of CIII]1909, CIV1549, HeII1640, and OIII]1661,1666 in various stacked subsets. Our analysis confirms earlier suggestions of a correlation between the strength of Ly-alpha and CIII] emission and we find similar trends with broad band UV luminosity and rest-frame UV colour. Using various diagnostic line ratios and our stellar photoionisation models, we determine both the gas phase metallicity and hardness of the ionisation spectrum characterised by xi_ion - the number of Lyman continuum photons per UV luminosity. We confirm our earlier suggestion that xi_ion is significantly larger for LAEs than for continuum-selected Lyman break galaxies, particularly for those LAEs with the faintest UV luminosities. We briefly discuss the implications for cosmic reionisation if the metal-poor intensely star-forming systems studied here are representative examples of those at much higher redshift.


title: Precision determination of absolute neutron flux

abstract: A technique for establishing the total neutron rate of a highly-collimated monochromatic cold neutron beam was demonstrated using a method of an alpha-gamma counter. The method involves only the counting of measured rates and is independent of neutron cross sections, decay chain branching ratios, and neutron beam energy. For the measurement, a target of 10B-enriched boron carbide totally absorbed the neutrons in a monochromatic beam, and the rate of absorbed neutrons was determined by counting 478keV gamma rays from neutron capture on 10B with calibrated high-purity germanium detectors. A second measurement based on Bragg diffraction from a perfect silicon crystal was performed to determine the mean de Broglie wavelength of the beam to a precision of 0.024 %. With these measurements, the detection efficiency of a neutron monitor based on neutron absorption on 6Li was determined to an overall uncertainty of 0.058 %. We discuss the principle of the alpha-gamma method and present details of how the measurement was performed including the systematic effects. We also describe how this method may be used for applications in neutron dosimetry and metrology, fundamental neutron physics, and neutron cross section measurements.


title: Design Studies for a TPC Readout Plane Using Zigzag Patterns with Multistage GEM Detectors

abstract: A new Time Projection Chamber (TPC) is currently under development for the sPHENIX experiment at RHIC. The TPC will be read out using multistage GEM detectors on each end and will be divided into approximately 40 pad layers in radius. Each pad layer is required to provide a spatial resolution of ~250 microns, which must be achieved with a minimal channel count in order to minimize the overall cost of the detector. The current proposal is to make the pads into a zigzag shape in order to enhance charge sharing among neighboring pads. This will allow for the possibility to interpolate the hit position to high precision, resulting in a position resolution many times better than the 2mm pitch of the readout pads. This paper discusses various simulation studies that were carried out to optimize the size and shape of the zigzag pads for the readout board for the TPC, along with the technical challenges in fabricating it. It also describes the performance of the first prototype readout board obtained from measurements carried out in the laboratory using a highly collimated X-ray source.


title: Construction and Performance of the Barrel Electromagnetic Calorimeter for the GlueX Experiment

abstract: The barrel calorimeter is part of the new spectrometer installed in Hall D at Jefferson Lab for the GlueX experiment. The calorimeter was installed in 2013, commissioned in 2014 and has been operating routinely since early 2015. The detector configuration, associated Monte Carlo simulations, calibration and operational performance are described herein. The calorimeter records the time and energy deposited by charged and neutral particles created by a multi-GeV photon beam. It is constructed as a lead and scintillating-fiber calorimeter and read out with 3840 large-area silicon photomultiplier arrays. Particles impinge on the detector over a wide range of angles, from normal incidence at 90 degrees down to 11.5 degrees, which defines a geometry that is fairly unique among calorimeters. The response of the calorimeter has been measured during a running experiment and performs as expected for electromagnetic showers below 2.5 GeV. We characterize the performance of the BCAL using the energy resolution integrated over typical angular distributions for $\pi^0$ and $\eta$ production of $\sigma_E/E$=5.2\%/$\sqrt{E(\rm{GeV})} \oplus$ 3.6\% and a timing resolution of $\sigma$\,=\,150\,ps at 1\,GeV.


title: Prospects for indirect detection of frozen-in dark matter

abstract: We study observational consequences arising from dark matter (DM) of non-thermal origin, produced by dark freeze-out from a hidden sector heat bath. We assume this heat bath was populated by feebly-coupled mediator particles, produced via a Higgs portal interaction with the Standard Model (SM). The dark sector then attained internal equilibrium with a characteristic temperature different from the SM photon temperature. We find that even if the coupling between the DM and the SM sectors is very weak, the scenario allows for indirect observational signals. We show how the expected strength of these signals depends on the temperature of the hidden sector at DM freeze-out.


title: Measurement-based universal blind quantum computation with minor resources

abstract: Blind quantum computation (BQC) enables a client with less quantum computational ability to delegate her quantum computation to a server with strong quantum computational power while preserving the client's privacy. Generally, many-qubit entangled states are often used to complete BQC tasks. But for a large-scale entangled state, it is difficult to be described since its Hilbert space dimension is increasing exponentially. Furthermore, the number of entangled qubits is limited in experiment of existing works. To tackle this problem, in this paper we propose a universal BQC protocol based on measurement with minor resources, where the trap technology is adopted to verify correctness of the server's measurement outcomes during computation and testing process. In our model there are two participants, a client who prepares initial single-qubit states and a server that performs universal quantum computation. The client is almost classical since she does not require any quantum computational power, quantum memory. To realize the client's universal BQC, we construct an $m\times n$ latticed state composed of six-qubit cluster states and eight-qubit cluster states, which needs less qubits than the brickwork state. Finally, we analyze and prove the blindness, correctness, universality and verifiability of our proposed BQC protocol.


title: A time-dependent Schr\"odinger equation for molecular core-hole dynamics

abstract: X-ray spectroscopy is an important tool for the investigation of matter. X rays primarily interact with inner-shell electrons creating core (inner-shell) holes that will decay on the time scale of attoseconds to few femtoseconds through electron relaxations involving the emission of a photon or an electron. The advent of femtosecond x-ray pulses expands x-ray spectroscopy to the time domain and will eventually allow the control of core-hole population on timescales comparable to core-vacancy lifetimes. For both cases, a theoretical approach that accounts for the x-ray interaction while the electron relaxations occur is required. Here we describe a time-dependent framework, based on solving the time-dependent Schr\"odinger equation, that is suitable for describing the induced electron and nuclear dynamics.


title: Ultrafast x-ray-induced nuclear dynamics in diatomic molecules using femtosecond x-ray/x-ray pump-probe spectroscopy

abstract: The capability of generating two intense, femtosecond x-ray pulses with controlled time delay opens the possibility of performing time-resolved experiments for x-ray induced phenomena. We have applied this capability to study the photoinduced dynamics in diatomic molecules. In molecules composed of low-Z elements, \textit{K}-shell ionization creates a core-hole state in which the main decay mode is an Auger process involving two electrons in the valence shell. After Auger decay, the nuclear wavepackets of the transient two-valence-hole states continue evolving on the femtosecond timescale, leading either to separated atomic ions or long-lived quasi-bound states. By using an x-ray pump and an x-ray probe pulse tuned above the \textit{K}-shell ionization threshold of the nitrogen molecule, we are able to observe ion dissociation in progress by measuring the time-dependent kinetic energy releases of different breakup channels. We simulated the measurements on N$_2$ with a molecular dynamics model that accounts for \textit{K}-shell ionization, Auger decay, and the time evolution of the nuclear wavepackets. In addition to explaining the time-dependent feature in the measured kinetic energy release distributions from the dissociative states, the simulation also reveals the contributions of quasi-bound states.


title: End-of-Use Core Triage in Extreme Scenarios Based on a Threshold Approach

abstract: Remanufacturing is a significant factor in securing sustainability through a circular economy. Sorting plays a significant role in remanufacturing pre-processing inspections. Its significance can increase when remanufacturing facilities encounter extreme situations, such as abnormally huge core arrivals. Our main objective in this work is switching from less efficient to a more efficient model and to characterize extreme behavior of core arrival in remanufacturing and applying the developed model to triage cores. Central tendency core flow models are not sufficient to handle extreme situations, however, complementary Extreme Value (EV) approaches have shown to improve model efficiency. Extreme core flows to remanufacturing facilities are rare but still likely and can adversely affect remanufacturing business operations. In this investigation, extreme end-of-use core flow is modelled by a threshold approach using the Generalized Pareto Distribution (GPD). It is shown that GPD has better performance than its maxima-block GEV counterpart from practical and data efficiency perspectives. The model is validated by a synthesized big dataset, tested by sophisticated statistical Anderson Darling (AD) test, and is applied to a case of extreme flow to a valve shop in order to predict probability of over-capacity arrivals that is critical in remanufacturing business management. Finally, the GPD model combined with triage strategies is used to initiate investigations into the efficacy of different triage methods in remanufacturing operations.


title: Long-term radio and X-ray evolution of the tidal disruption event ASASSN-14li

abstract: We report on late time radio and X-ray observations of the tidal disruption event candidate ASASSN-14li, covering the first 1000 days of the decay phase. For the first $\sim200$ days the radio and X-ray emission fade in concert. This phase is better fit by an exponential decay at X-ray wavelengths, while the radio emission is well described by either an exponential or the canonical $t^{-5/3}$ decay assumed for tidal disruption events. The correlation between radio and X-ray emission during this period can be fit as $L_{R}\propto L_{X}^{1.9\pm0.2}$. After 400 days the radio emission at $15.5\,\textrm{GHz}$ has reached a plateau level of $244\pm8\,\mu\textrm{Jy}$ which it maintains for at least the next 600 days, while the X-ray emission continues to fade exponentially. This steady level of radio emission is likely due to relic radio lobes from the weak AGN-like activity implied by historical radio observations. We note that while most existing models are based upon the evolution of ejecta which are decoupled from the central black hole, the radio : X-ray correlation during the declining phase is also consistent with core jet emission coupled to a radiatively efficient accretion flow.


title: Evidence for feedback and stellar-dynamically regulated bursty star cluster formation: the case of the Orion Nebula Cluster

abstract: (abridged) A scenario for the formation of multiple co-eval populations separated in age by about 1 Myr in very young clusters (VYCs, ages less than 10 Myr) and with masses in the range 600-20000 Msun is outlined. It rests upon a converging inflow of molecular gas building up a first population of pre-main sequence stars. The associated just-formed O stars ionise the inflow and suppress star formation in the embedded cluster. However, they typically eject each other out of the embedded cluster within 10^6 yr, that is before the molecular cloud filament can be ionised entirely. The inflow of molecular gas can then resume forming a second population. This sequence of events can be repeated multiply. This model is applied to the Orion Nebula Cluster (ONC), in which three well-separated pre-main sequences in the color-magnitude diagram of the cluster have recently been discovered. The mass-inflow history is constrained using this model and the number of OB stars ejected from each population are estimated for verification using Gaia data. As a consequence of the proposed model, the three runaway O star systems, AE Aur, mu Col and iota Ori, are considered as significant observational evidence for stellar-dynamical ejections of massive stars from the oldest population in the ONC.


title: Accretion of satellites onto central galaxies in clusters: merger mass ratios and orbital parameters

abstract: We study the statistical properties of mergers between central and satellite galaxies in galaxy clusters in the redshift range $0<z<1$, using a sample of dark-matter only cosmological N-body simulations from Le SBARBINE dataset. Using a spherical overdensity algorithm to identify dark-matter haloes, we construct halo merger trees for different values of the over-density $\Delta_c$. While the virial overdensity definition allows us to probe the accretion of satellites at the cluster virial radius $r_{vir}$, higher overdensities probe satellite mergers in the central region of the cluster, down to $\approx 0.06 r_{vir}$, which can be considered a proxy for the accretion of satellite galaxies onto central galaxies. We find that the characteristic merger mass ratio increases for increasing values of $\Delta_c$: more than $60\%$ of the mass accreted by central galaxies since $z\approx 1$ comes from major mergers. The orbits of satellites accreting onto central galaxies tend to be more tangential and more bound than orbits of haloes accreting at the virial radius. The obtained distributions of merger mass ratios and orbital parameters are useful to model the evolution of the high-mass end of the galaxy scaling relations without resorting to hydrodynamic cosmological simulations.


title: Stellar Streams Discovered in the Dark Energy Survey

abstract: We perform a search for stellar streams around the Milky Way using the first three years of multi-band optical imaging data from the Dark Energy Survey (DES). We use DES data covering $\sim 5000$ sq. deg. to a depth of $g > 23.5$ with a relative photometric calibration uncertainty of $< 1 \%$. This data set yields unprecedented sensitivity to the stellar density field in the southern celestial hemisphere, enabling the detection of faint stellar streams to a heliocentric distance of $\sim 50$ kpc. We search for stellar streams using a matched-filter in color-magnitude space derived from a synthetic isochrone of an old, metal-poor stellar population. Our detection technique recovers four previously known thin stellar streams: Phoenix, ATLAS, Tucana III, and a possible extension of Molonglo. In addition, we report the discovery of eleven new stellar streams. In general, the new streams detected by DES are fainter, more distant, and lower surface brightness than streams detected by similar techniques in previous photometric surveys. As a by-product of our stellar stream search, we find evidence for extra-tidal stellar structure associated with four globular clusters: NGC 288, NGC 1261, NGC 1851, and NGC 1904. The ever-growing sample of stellar streams will provide insight into the formation of the Galactic stellar halo, the Milky Way gravitational potential, as well as the large- and small-scale distribution of dark matter around the Milky Way.


title: $\mathcal{N}=2$ supersymmetry and anisotropic scale invariance

abstract: We find a class of scale-anomaly-free $\mathcal{N}=2$ supersymmetric quantum systems with non-vanishing potential terms where space and time scale with distinct exponents. Our results generalise the known case of the supersymmetric inverse square potential to a larger class of scaling symmetries.


title: Predicting the binary black hole population of the Milky Way with cosmological simulations

abstract: Binary black holes are the primary endpoint of massive stellar evolution. Their properties provide a unique opportunity to constrain binary evolution, which is still poorly understood. In this paper, we predict the inventory of binary black holes and their merger products in/around the Milky Way, and detail their main properties. We present the first combination of a high-resolution cosmological simulation of a Milky Way-mass galaxy with a binary population synthesis model. The hydrodynamic simulation, taken from the FIRE project, provides a cosmologically realistic star formation history for the galaxy and its stellar halo and satellites. We apply a metallicity-dependent evolutionary model to the star particles to produce individual binary black holes. We find that a million binary black holes have merged in the model Milky Way, and 3 million binaries are still present, with an average mass of 28 Msun per binary. Because the black hole progenitors are biased towards low metallicity stars, half reside in the stellar halo and satellites and 40 per cent of the binaries were formed outside the main galaxy. This trend increases with the masses of the black holes. The numbers and mass distribution of the merged systems is compatible with the LIGO/Virgo detections. Observations of these black holes will be challenging, both with electromagnetic methods and LISA. We find that a cosmologically realistic star formation history, with self-consistent metal enrichment and Galactic accretion history, are key ingredients for determining binary black hole rates that can be compared with observations to constrain massive binary evolution.


title: Multi-timescale X-ray reverberation mapping of accreting black holes

abstract: Accreting black holes show characteristic reflection features in their X-ray spectrum, including an iron K$\alpha$ line, resulting from hard X-ray continuum photons illuminating the accretion disk. The reverberation lag resulting from the path length difference between direct and reflected emission provides a powerful tool to probe the innermost regions around both stellar-mass and supermassive black holes. Here, we present for the first time a reverberation mapping formalism that enables modeling of energy dependent time lags and variability amplitude for a wide range of variability timescales, taking the complete information of the cross-spectrum into account. We use a pivoting power-law model to account for the spectral variability of the continuum that dominates over the reverberation lags for longer time scale variability. We use an analytic approximation to self-consistently account for the non-linear effects caused by this continuum spectral variability, which have been ignored by all previous reverberation studies. We find that ignoring these non-linear effects can bias measurements of the reverberation lags, particularly at low frequencies. Since our model is analytic, we are able to fit simultaneously for a wide range of Fourier frequencies without prohibitive computational expense. We also introduce a formalism of fitting to real and imaginary parts of our cross-spectrum statistic, which naturally avoids some mistakes/inaccuracies previously common in the literature. We perform proof-of-principle fits to Rossi X-ray Timing Explorer data of Cygnus X-1.


title: Many-Body Spectral Reflection Symmetry and Protected Infinite-Temperature Degeneracy

abstract: Protected zero modes in quantum physics traditionally arise in the context of ground states of many-body Hamiltonians. Here we study the case where zero modes exist in the center of a reflection-symmetric many-body spectrum, giving rise to the notion of a protected "infinite-temperature" degeneracy. For a certain class of nonintegrable spin chains, we show that the number of zero modes is determined by a chiral index that grows exponentially with system size. We propose a dynamical protocol, feasible in ongoing experiments in Rydberg atom quantum simulators, to detect these many-body zero modes and their protecting spectral reflection symmetry. Finally, we consider whether the zero energy states obey the eigenstate thermalization hypothesis, as is expected of states in the middle of the many-body spectrum. We find intriguing differences in their eigenstate properties relative to those of nearby nonzero-energy eigenstates at finite system sizes.

categories: cond-mat.str-el cond-mat.quant-gas cond-mat.stat-mech quant-ph


title: Advanced Environmentally Resistant Lithium Fluoride Mirror Coatings for the Next-Generation of Broadband Space Observatories

abstract: Recent advances in the physical vapor deposition (PVD) of protective fluoride films have raised the far-ultraviolet (FUV: 912-1600 {\AA}) reflectivity of aluminum-based mirrors closer to the theoretical limit. The greatest gains, at more than 20%, have come for lithium fluoride-protected aluminum, which has the shortest wavelength cutoff of any conventional overcoat. Despite the success of the NASA FUSE mission, the use of lithium fluoride (LiF)-based optics is rare, as LiF is hygroscopic and requires handling procedures that can drive risk. With NASA now studying two large mission concepts for astronomy, Large UV-Optical-IR Surveyor (LUVOIR) and the Habitable Exoplanet Imaging Mission (HabEx), which mandate throughput down to 1000 {\AA}, the development of LiF-based coatings becomes crucial. This paper discusses steps that are being taken to qualify these new enhanced LiF-protected aluminum (eLiF) mirror coatings for flight. In addition to quantifying the hygroscopic degradation, we have developed a new method of protecting eLiF with an ultrathin (10-20 {\AA}) capping layer of a non-hygroscopic material to increase durability. We report on the performance of eLiF-based optics and assess the steps that need to be taken to qualify such coatings for LUVOIR, HabEx, and other FUV-sensitive space missions.


title: A Candidate $z\sim10$ Galaxy Strongly Lensed into a Spatially Resolved Arc

abstract: The most distant galaxies known are at z~10-11, observed 400-500 Myr after the Big Bang. The few z~10-11 candidates discovered to date have been exceptionally small- barely resolved, if at all, by the Hubble Space Telescope. Here we present the discovery of SPT0615-JD, a fortuitous z~10 (z_phot=9.9+/-0.6) galaxy candidate stretched into an arc over ~2.5" by the effects of strong gravitational lensing. Discovered in the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury program and companion S-RELICS Spitzer program, this candidate has a lensed H-band magnitude of 25.7+/-0.1 AB mag. With a magnification of \mu~4-7 estimated from our lens models, the de-lensed intrinsic magnitude is 27.6+/-0.3 AB mag, and the half-light radius is r_e<0.8 kpc, both consistent with other z>9 candidates. The inferred stellar mass (log [M* /M_Sun]=9.7^{+0.7}_{-0.5}) and star formation rate (\log [SFR/M_Sun yr^{-1}]=1.3^{+0.2}_{-0.3}) indicate that this candidate is a typical star-forming galaxy on the z>6 SFR-M* relation. We note that three independent lens models predict two counterimages, at least one of which should be of a similar magnitude to the arc, but these counterimages are not yet detected. Counterimages would not be expected if the arc were at lower redshift. However, the only spectral energy distributions capable of fitting the Hubble and Spitzer photometry well at lower redshifts require unphysical combinations of z~2 galaxy properties. The unprecedented lensed size of this z~10 candidate offers the potential for the James Webb Space Telescope to study the geometric and kinematic properties of a galaxy observed 500 Myr after the Big Bang.


title: Physical properties and scaling relations of molecular clouds: the effect of stellar feedback

abstract: Using hydrodynamical simulations of entire galactic discs similar to the Milky Way, reaching 4.6pc resolution, we study the origins of observed physical properties of giant molecular clouds (GMCs). We find that efficient stellar feedback is a necessary ingredient in order to develop a realistic interstellar medium (ISM), leading to molecular cloud masses, sizes, velocity dispersions and virial parameters in excellent agreement with Milky Way observations. GMC scaling relations observed in the Milky Way, such as the mass-size ($M$--$R$), velocity dispersion-size ($\sigma$--$R$), and the $\sigma$--$R\Sigma$ relations, are reproduced in a feedback driven ISM when observed in projection, with $M\propto R^{2.3}$ and $\sigma\propto R^{0.56}$. When analysed in 3D, GMC scaling relations steepen significantly, indicating potential limitations of our understanding of molecular cloud 3D structure from observations. Furthermore, we demonstrate how a GMC population's underlying distribution of virial parameters can strongly influence the scatter in derived scaling relations. Finally, we show that GMCs with nearly identical global properties exist in different evolutionary stages, where a majority of clouds being either gravitationally bound or expanding, but with a significant fraction being compressed by external ISM pressure, at all times.


title: Full-disc $^{13}$CO(1-0) mapping across nearby galaxies of the EMPIRE survey and the CO-to-H$_2$ conversion factor

abstract: Carbon monoxide (CO) provides crucial information about the molecular gas properties of galaxies. While $^{12}$CO has been targeted extensively, isotopologues such as $^{13}$CO have the advantage of being less optically thick and observations have recently become accessible across full galaxy discs. We present a comprehensive new dataset of $^{13}$CO(1-0) observations with the IRAM 30-m telescope of the full discs of 9 nearby spiral galaxies from the EMPIRE survey at a spatial resolution of $\sim$1.5kpc. $^{13}$CO(1-0) is mapped out to $0.7-1r_{25}$ and detected at high signal-to-noise throughout our maps. We analyse the $^{12}$CO(1-0)-to-$^{13}$CO(1-0) ratio ($\Re$) as a function of galactocentric radius and other parameters such as the $^{12}$CO(2-1)-to-$^{12}$CO(1-0) intensity ratio, the 70-to-160$\mu$m flux density ratio, the star-formation rate surface density, the star-formation efficiency, and the CO-to-H$_2$ conversion factor. We find that $\Re$ varies by a factor of 2 at most within and amongst galaxies, with a median value of 11 and larger variations in the galaxy centres than in the discs. We argue that optical depth effects, most likely due to changes in the mixture of diffuse/dense gas, are favored explanations for the observed $\Re$ variations, while abundance changes may also be at play. We calculate a spatially-resolved $^{13}$CO(1-0)-to-H$_2$ conversion factor and find an average value of $1.0\times10^{21}$ cm$^{-2}$ (K.km/s)$^{-1}$ over our sample with a standard deviation of a factor of 2. We find that $^{13}$CO(1-0) does not appear to be a good predictor of the bulk molecular gas mass in normal galaxy discs due to the presence of a large diffuse phase, but it may be a better tracer of the mass than $^{12}$CO(1-0) in the galaxy centres where the fraction of dense gas is larger.


title: Why informatics and general science need a conjoint basic definition of information

abstract: First the basic definition of information as a selection from a set of possibilities resp. domain is recalled. This also applies to digital information. The bits of digital information are parts of number sequences which represent a selection from a set of possibilities resp. domain. For faultless conversation sender and receiver of information must have the same definition of the domain (e.g. of language vocabulary). Up to now the definition of the domain and of its elements is derived from context and knowledge. The internet provides an additional important possibility: A link to a conjoint uniform definition of the domain at unique location on the internet. The associated basic information structure is called "Domain Vector" (DV) and has the structure "UL (of the domain definition) plus sequence of numbers". The "UL" is not only "Uniform Locator" of the domain definition. It also identifies a certain kind of information for later comparison and search. It can be a Uniform Resource Locator (URL) or an abbreviated equivalent, e.g. a hierarchic numeric pointer or a short local pointer to a table with global internet pointers. The DV structure can be used as general carrier of information which is language independent and more precise than language. A domain which contains DVs is called "Domain Space" (DS) and is defined as metric space. This allows similarity search according to user defined criteria, so that any kind of definable information can be made comparable and searchable according to user selected (relevant) and objectifiable (globally uniform) criteria. DS definitions can be reused in new DS definitions. Their elements, the DVs, are automatically globally uniformly identified and defined. Obviously such conjoint definition of comparable information has great potential. It also can avoid interoperability problems and redundant programming and so save high costs.


title: An excess of massive stars in the local 30 Doradus starburst

abstract: The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analogue of large star-formation events in the distant Universe. We determine the recent formation history and the initial mass function (IMF) of massive stars in 30 Doradus based on spectroscopic observations of 247 stars more massive than 15 solar masses ($\mathrm{M}_\odot$). The main episode of massive star formation started about $8\,\mathrm{Myr}$ ago and the star-formation rate seems to have declined in the last $1\,\mathrm{Myr}$. The IMF is densely sampled up to $200\,\mathrm{M}_\odot$ and contains $32\pm12\%$ more stars above $30\,\mathrm{M}_\odot$ than predicted by a standard Salpeter IMF. In the mass range $15-200\,\mathrm{M}_\odot$, the IMF power-law exponent is $1.90^{+0.37}_{-0.26}$, shallower than the Salpeter value of 2.35.


title: Enforcing the Courant-Friedrichs-Lewy Condition in Explicitly Conservative Local Time Stepping Schemes

abstract: An optimally efficient explicit numerical scheme for solving fluid dynamics equations, or any other parabolic or hyperbolic system of partial differential equations, should allow local regions to advance in time with their own, locally constrained time steps. However, such a scheme can result in violation of the Courant-Friedrichs-Lewy (CFL) condition, which is manifestly non-local. Although the violations can be considered to be "weak" in a certain sense and the corresponding numerical solution may be stable, such calculation does not guarantee the correct propagation speed for arbitrary waves. We use an experimental fluid dynamics code that allows cubic "patches" of grid cells to step with independent, locally constrained time steps to demonstrate how the CFL condition can be enforced by imposing a condition on the time steps of neighboring patches. We perform several numerical tests that illustrate errors introduced in the numerical solutions by weak CFL condition violations and show how strict enforcement of the CFL condition eliminates these errors. In all our tests the strict enforcement of the CFL condition does not impose a significant performance penalty.


title: An operator-valued Lyapunov theorem

abstract: We generalize Lyapunov's convexity theorem for classical (scalar-valued) measures to quantum (operator-valued) measures. In particular, we show that the range of a nonatomic quantum probability measure is a weak*-closed convex set of quantum effects (positive operators bounded above by the identity operator) under a sufficient condition on the non-injectivity of integration. To prove the operator-valued version of Lyapunov's theorem, we must first define the notions of essentially bounded, essential support, and essential range for quantum random variables (Borel measurable functions from a set to the bounded linear operators acting on a Hilbert space).


title: QCD thermodynamics from lattice calculations with non-equilibrium methods: The SU(3) equation of state

abstract: A precise lattice determination of the equation of state in SU(3) Yang-Mills theory is carried out by means of a simulation algorithm, based on Jarzynski's theorem, that allows one to compute physical quantities in thermodynamic equilibrium, by driving the field configurations of the system out of equilibrium. The physical results and the computational efficiency of the algorithm are compared with other state-of-the-art lattice calculations, and the extension to full QCD with dynamical fermions and to other observables is discussed.


title: Lambert function methods for laser dynamics with time-delayed feedback

abstract: Time-delayed differential equations arise frequently in the study of nonlinear dynamics of lasers with optical feedback. Traditionally, one has resorted to numerical methods because the analytical solution of such equations are intractable. In this manuscript, we show that under some conditions, the rate equations model that is used to model semiconductor lasers with feedback can be analytically solved by using the Lambert W function. In particular, we discuss the conditions under which the coupled rate equations for the intra-cavity electric field and excess carrier inversion can be reduced to a single equation for the field, and how this single rate equation can be cast in a form that is amenable to the use of the Lambert W function. We conclude the manuscript with a similar discussion for two lasers coupled via time-delayed feedbacks.


title: Fundamental Solutions and Decay of Fully Non-local Problems

abstract: In this paper, we study a fully non-local reaction-diffusion equation which is non-local both in time and space. We apply subordination principles to construct the fundamental solutions of this problem, which we use to find a representation of the mild solutions. Moreover, using techniques of Harmonic Analysis and Fourier Multipliers, we obtain the temporal decay rates for the mild solutions.


title: Influence of image forces on the interlayer exchange interaction in magnetic tunnel junctions with ferroelectric barrier

abstract: We study interlayer exchange interaction in magnetic tunnel junctions with ferroelectric barrier. We focus on the influence of image forces on the voltage dependence of the interlayer magnetic interaction (magneto-electric effect). The influence of the image forces is twofold: 1) they significantly enforce magneto-electric effect occurring due to the surface charges at the interface between ferroelectric and ferromagnets; 2) in combination with voltage dependent dielectric constant of the ferroelectric barrier image forces cause an additional contribution to the magneto-electric effect in magnetic tunnel junctions. This contribution can exceed the one coming from surface charges. We compare the interlayer exchange coupling voltage variation with spin transfer torque effect and show that for half-metallic electrodes the interlayer exchange coupling variation is dominant and defines the magnetic state and dynamics of magnetization in the tunnel junction.


title: Magneto-electric effect in doped magnetic ferroelectrics

abstract: We propose a model of magneto-electric effect in doped magnetic ferroelectrics. This magneto-electric effect does not involve the spin-orbit coupling and is based purely on the Coulomb interaction. We calculate magnetic phase diagram of doped magnetic ferroelectrics. We show that magneto-electric coupling is pronounced only for ferroelectrics with low dielectric constant. We find that magneto-electric coupling leads to modification of magnetization temperature dependence in the vicinity of ferroelectric phase transition. A peak of magnetization appears. We find that magnetization of doped magnetic ferroelectrics strongly depends on applied electric field.


title: Topological Aspect of Graphene Patchwork with Regular Arrays of Nano Holes

abstract: Triangular and honeycomb lattices are dual to each other -- if we puncture holes into a featureless plane in a regular triangular alignment, the remaining body looks like a honeycomb lattice, and vice versa, if the holes are in a regular honeycomb alignment, the remaining body has a feature of triangular lattice. In this work, we reveal that the electronic states in graphene sheets with nano-sized holes in triangular and honeycomb alignments are also dual to each other in a topological sense. Namely, a regular hole array perforated in graphene can open a band gap in the energy-momentum dispersion of relativistic electrons in the pristine graphene, and the insulating states induced by triangular and honeycomb hole arrays are distinct in topology. In a graphene patchwork with regions of these two hole arrays put side by side counterpropagating topological currents emerge at the domain wall. This observation indicates that the cerebrated atomically thin sheet is where topological physics and nanotechnology meet.


title: Existence and continuity of solution trajectories of generalized equations with application in electronics

abstract: We consider a special form of parametric generalized equations arising from electronic circuits with AC sources and study the effect of perturbing the input signal on solution trajectories. Using methods of variational analysis and strong metric regularity property of an auxiliary map, we are able to prove the regularity properties of the solution trajectories inherited by the input signal. Furthermore, we establish the existence of continuous solution trajectories for the perturbed problem. This can be achieved via a result of uniform strong metric regularity for the auxiliary map. Key words and phrases: generalized equations, electronic circuits, strong metric regularity, uniform strong metric regularity, perturbations.


title: Magnetic Monopole Supercurrent through a Quantum Spin Ice Tunnel Junction

abstract: Magnetic monopoles are hypothetical particles that may exist as quantized sources and sinks of the magnetic field. In materials, they may appear in an emergent quantum electrodynamics described by a U(1) lattice gauge theory. Particularly, quantum spin ice hosts monopoles as bosonic spinons coupled to emergent gauge fields in a U(1) quantum spin liquid, namely, a deconfined Coulomb phase. When monopoles are condensed to form a long-range order, monopoles and gauge fields are screened and confined. Here we show, however, that monopole supercurrent flows across a junction of two ferromagnets that are weakly linked through and placed on top of the U(1) QSL, when a gauge-invariant phase difference of spinons across the junction is generated by quenching or an applied electric voltage parallel to the junction. This novel phenomenon paves the way to a new paradigm of spinonics for a dissipationless control of magnetism.


title: Nonperturbative renormalization-group approach preserving the momentum dependence of correlation functions

abstract: We present an approximation scheme of the nonperturbative renormalization group that preserves the momentum dependence of correlation functions. This approximation scheme can be seen as a simple improvement of the local potential approximation (LPA) where the derivative terms in the effective action are promoted to arbitrary momentum-dependent functions. As in the LPA the only field dependence comes from the effective potential, which allows us to solve the renormalization-group equations at a relatively modest numerical cost (as compared, e.g., to the Blaizot--Mend\'ez-Galain--Wschebor approximation scheme). As an application we consider the two-dimensional quantum O($N$) model at zero temperature. We discuss not only the two-point correlation function but also higher-order correlation functions such as the scalar susceptibility (which allows for an investigation of the "Higgs" amplitude mode) and the conductivity. In particular we show how, using Pad\'e approximants to perform the analytic continuation $i\omega_n\to\omega+i0^+$ of imaginary frequency correlation functions $\chi(i\omega_n)$ computed numerically from the renormalization-group equations, one can obtain spectral functions in the real-frequency domain.

categories: cond-mat.stat-mech cond-mat.quant-gas cond-mat.str-el hep-th


title: Giulio Fermi's contributions to biophysics and molecular biology

abstract: This paper presents a comprehensive list of the scientific articles of Giulio Fermi (1936-1997), son of the Italian-American physicist Enrico Fermi, published between 1962 and 1997. The initial research activity of Giulio was concerned with virology and biological cybernetics while, from 1975 onward, his work was completely devoted to protein crystallography. The crystallographic research was carried out in collaboration with Nobel laureate Max Perutz at the Medical Research Council (MRC) Laboratory of Molecular Biology in Cambridge (United Kingdom). A short biography of Giulio (Judd) Fermi appears inside John Finch's book A Nobel Fellow on Every Floor: A History of the Medical Research Council Laboratory of Molecular Biology published by the MRC in 2008.


title: Limiting curves for the dyadic odometer and the generalized Trollope-Delange formula

abstract: We study limiting curves resulting from deviations in partial sums in the ergodic theorem for the dyadic odometer and non-cylindric functions. In particular, we generalize the Trollope-Delange formula for the case of the weighted sum-of-binary-digits function and show that the Takagi-Landsberg curve arises.


title: Merging massive black holes: the right place and the right time

abstract: The LIGO/Virgo detections of gravitational waves from merging black holes of $\simeq$ 30 solar mass suggest progenitor stars of low metallicity (Z/Z$_{\odot} \lesssim 0.3$). In this talk I will provide constrains on where the progenitors of GW150914 and GW170104 may have formed, based on advanced models of galaxy formation and evolution combined with binary population synthesis models. First I will combine estimates of galaxy properties (star-forming gas metallicity, star formation rate and merger rate) across cosmic time to predict the low redshift BBH merger rate as a function of present day host galaxy mass, formation redshift of the progenitor system and different progenitor metallicities. I will show that the signal is dominated by binaries formed at the peak of star formation in massive galaxies with and binaries formed recently in dwarf galaxies. Then, I will present what very high resolution hydrodynamic simulations of different galaxy types can learn us about their black hole populations.


title: A detailed treatment of Doob's theorem

abstract: Doob's theorem provides guarantees of consistent estimation and posterior consistency under very general conditions. Despite the limitation that it only guarantees consistency on a set with prior probability 1, for many models arising in practice, Doob's theorem is an easy way of showing that consistency will hold almost everywhere. In this article, we give a detailed proof of Doob's theorem.


title: Level-dependent interpolatory Hermite subdivision schemes and wavelets

abstract: We study many properties of level-dependent Hermite subdivision, focusing on schemes preserving polynomial and exponential data. We specifically consider interpolatory schemes, which give rise to level-dependent multiresolution analyses through a prediction-correction approach. A result on the decay of the associated multiwavelet coefficients, corresponding to a uniformly continuous and differentiable function, is derived. It makes use of the approximation of any such function with a generalized Taylor formula expressed in terms of polynomials and exponentials.


title: Is 4C+29.48 a gamma-ray source?

abstract: The Fermi Large Area Telescope revealed that the extragalactic gamma-ray sky is dominated by blazars, active galactic nuclei (AGN) whose jet is seen at very small angle to the line of sight. To associate and then classify the gamma-ray sources, data have been collected from lower frequency surveys and observations. The gamma-ray source 3FGL J1323.0+2942 is associated with the radio source 4C+29.48 and classified as a blazar of unknown type, lacking optical spectrum and redshift. The higher-resolution radio data showed that 4C+29.48 comprises three bright radio-emitting features located within a ~1'-diameter area. We (re-)analyzed archival Very Large Array and unpublished very long baseline interferometry (VLBI) observations conducted by the Very Long Baseline Array and the European VLBI Network of 4C+29.48. We also collected data form optical, infrared and X-ray surveys. The northernmost complex of 4C+29.48 contains a blazar with a high brightness temperature compact core and a steep-spectrum jet feature. The blazar is positionally coincident with an optical source at a redshift of 1.142. Its mid-infrared colors also support its association with a gamma-ray emitting blazar. The two other radio complexes have steep radio spectra and do not have optical or infrared counterparts in currently available surveys. Based on the radio morphology, they are unlikely to be related to the blazar. We discuss the possibilities whether the two radio features are lobes of a radio galaxy, or gravitationally lensed images of a background source. We propose to associate 3FGL J1323.0+2942 in subsequent versions of the Fermi catalog with the blazar residing in northernmost complex. We suggest naming this source J1323+2941A to avoid misinterpretation arising from the fact that the coordinates of the currently listed radio counterpart 4C+29.48 is closer to a most probably unrelated radio source. (abridged)


title: Mathematical modeling of cholera epidemics in South Sudan

abstract: In this work, we analyze and model the cholera epidemics that affected South Sudan, the newest country in the world, during 2014 and 2015. South Sudan possibly represents one of the most difficult context in which adapt the deterministic mathematical cholera model, due to the unstable social and political situation that clearly affects the fluxes of people and the sanitary conditions, increasing the risk of large outbreaks. Despite the limitation of a static gravity model in describing the chaotic human mobility of South Sudan, the SIRB model, calibrated with a data assimilation technique (Ensemble Kalman Filter), retrieves the epidemic dynamics in the counties with the largest number of infected cases, showing the potentiality of the methodology in forecasting future outbreaks.


title: Schramm's formula for multiple loop-erased random walks

abstract: We revisit the computation of the discrete version of Schramm's formula for the loop-erased random walk derived by Kenyon. The explicit formula in terms of the Green function relies on the use of a complex connection on a graph, for which a line bundle Laplacian is defined. We give explicit results in the scaling limit for the upper half-plane, the cylinder and the Moebius strip. Schramm's formula is then extended to multiple loop-erased random walks.


title: Recognizing Material Properties from Images

abstract: Humans rely on properties of the materials that make up objects to guide our interactions with them. Grasping smooth materials, for example, requires care, and softness is an ideal property for fabric used in bedding. Even when these properties are not visual (e.g. softness is a physical property), we may still infer their presence visually. We refer to such material properties as visual material attributes. Recognizing these attributes in images can contribute valuable information for general scene understanding and material recognition. Unlike well-known object and scene attributes, visual material attributes are local properties with no fixed shape or spatial extent. We show that given a set of images annotated with known material attributes, we may accurately recognize the attributes from small local image patches. Obtaining such annotations in a consistent fashion at scale, however, is challenging. To address this, we introduce a method that allows us to probe the human visual perception of materials by asking simple yes/no questions comparing pairs of image patches. This provides sufficient weak supervision to build a set of attributes and associated classifiers that, while unnamed, serve the same function as the named attributes we use to describe materials. Doing so allows us to recognize visual material attributes without resorting to exhaustive manual annotation of a fixed set of named attributes. Furthermore, we show that this method may be integrated in the end-to-end learning of a material classification CNN to simultaneously recognize materials and discover their visual attributes. Our experimental results show that visual material attributes, whether named or automatically discovered, provide a useful intermediate representation for known material categories themselves as well as a basis for transfer learning when recognizing previously-unseen categories.


title: Dust arcs in the region of Jupiter's Trojan asteroids

abstract: Aims. The surfaces of the Trojan asteroids are steadily bombarded by interplanetary micrometeoroids, which releases ejecta of small dust particles. These particles form the faint dust arcs that are associated with asteroid clouds. Here we analyze the particle dynamics and structure of the arc in the region of the $L_4$ Trojan asteroids. Methods. We calculate the total cross section of the $L_4$ Trojan asteroids and the production rate of dust particles. The motion of the particles is perturbed by a variety of forces. We simulate the dynamical evolution of the dust particles, and explore the overall features of the Trojan dust arc. Results. The simulations show that the arc is mainly composed of grains in the size range 4-10 microns. Compared to the $L_4$ Trojan asteroids, the dust arc is distributed more widely in the azimuthal direction, extending to a range of [30,120] degrees relative to Jupiter. The peak number density does not develop at $L_4$. There exist two peaks that are azimuthally displaced from $L_4$.


title: Phases of superclimbing dislocation with long-range interaction between jogs

abstract: The main candidate for the superfluid pathways in solid Helium-4 are dislocations with Burgers vector along the hcp symmetry axis. Here we focus on quantum behavior of a generic edge dislocation which can perform superclimb -- climb supported by the superflow along its core. The role of the long range elastic interactions between jogs is addressed by Monte Carlo simulations. It is found that such interactions do not change qualitatively the phase diagram found without accounting for such forces. Their main effect consists of renormalizing the effective scale determining compressibility of the dislocation in the Tomonaga-Luttinger Liquid phase. It is also found that the quantum rough phase of the dislocation can be well described within the gaussian approximation which features off-diagonal long range order in 1D for the superfluid order parameter along the core.


title: Isovector and flavor-diagonal charges of the nucleon

abstract: We present an update on the status of the calculations of isovector and flavor-diagonal charges of the nucleon. The calculations of the isovector charges are being done using ten $2+1+1$-flavor HISQ ensembles generated by the MILC collaboration covering the range of lattice spacings $a \approx$ 0.12, 0.09, 0.06 fm and pion masses $M_\pi \approx$ 310, 220, 130 MeV. Excited-states contamination is controlled by using four-state fits to two-point correlators and three-states fits to the three-point correlators. The calculations of the disconnected diagrams needed to estimate flavor-diagonal charges are being done on a subset of six ensembles using the stocastic method. Final results are obtained using a simultaneous fit in $M_\pi^2$, the lattice spacing $a$ and the finite volume parameter $M_\pi L$ keeping only the leading order corrections.


title: Lifted graphene nanoribbons on gold: from smooth sliding to multiple stick-slip regimes

abstract: Graphene nanoribbons (GNRs) physisorbed on a Au(111) surface can be picked up, lifted at one end, and made slide by means of the tip of an atomic-force microscope. The dynamical transition from smooth sliding to multiple stick-slip regimes, the pushing/pulling force asymmetry, the presence of pinning, and its origin are real frictional processes in a nutshell, in need of a theoretical description. To this purpose, we conduct classical simulations of frictional manipulations for GNRs up to 30 nm in length, one end of which is pushed or pulled horizontally while held at different heights above the Au surface. These simulations allow us to clarify theoretically the emergence of stick-slip originating from the short 1D edges rather than the 2D "bulk", the role of adhesion, of lifting, and of graphene bending elasticity in determining the GNR sliding friction. The understanding obtained in this simple context is of additional value for more general cases.


title: Robust Propensity Score Computation Method based on Machine Learning with Label-corrupted Data

abstract: In biostatistics, propensity score is a common approach to analyze the imbalance of covariate and process confounding covariates to eliminate differences between groups. While there are an abundant amount of methods to compute propensity score, a common issue of them is the corrupted labels in the dataset. For example, the data collected from the patients could contain samples that are treated mistakenly, and the computing methods could incorporate them as a misleading information. In this paper, we propose a Machine Learning-based method to handle the problem. Specifically, we utilize the fact that the majority of sample should be labeled with the correct instance and design an approach to first cluster the data with spectral clustering and then sample a new dataset with a distribution processed from the clustering results. The propensity score is computed by Xgboost, and a mathematical justification of our method is provided in this paper. The experimental results illustrate that xgboost propensity scores computing with the data processed by our method could outperform the same method with original data, and the advantages of our method increases as we add some artificial corruptions to the dataset. Meanwhile, the implementation of xgboost to compute propensity score for multiple treatments is also a pioneering work in the area.


title: Dynamics and distribution of Jovian dust ejected from the Galilean satellites

abstract: In this paper, the dynamical analysis of the Jovian dust originating from the four Galilean moons is presented. High accuracy orbital integrations of dust particles are used to determine their dynamical evolution. A variety of forces are taken into account, including the Lorentz force, solar radiation pressure, Poynting-Robertson drag, solar gravity, the satellites' gravity, plasma drag, and gravitational effects due to non-sphericity of Jupiter. More than 20,000 dust particles from each source moon in the size range from 0.05 micron to 1 cm are simulated over 8,000 (Earth) years until each dust grain hits a sink (moons, Jupiter, or escape from the system). Configurations of dust number density in the Jovicentric equatorial inertial frame are calculated and shown. In a Jovicentric frame rotating with the Sun the dust distributions are found to be asymmetric. For certain small particle sizes, the dust population is displaced towards the Sun, while for certain larger sizes, the dust population is displaced away from the Sun. The average lifetime as a function of particle size for ejecta from each source moon is derived, and two sharp jumps in the average lifetime are analyzed. Transport of dust between the Galilean moons and to Jupiter is investigated. Most of the orbits for dust particles from Galilean moons are prograde, while, surprisingly, a small fraction of orbits are found to become retrograde mainly due to solar radiation pressure and Lorentz force. The distribution of orbital elements is also analyzed.


title: A $q$-analog of Jacobi's two squares formula and its applications

abstract: We consider a $q$-analog $r_2(n, q)$ of the number of representations of an integer as a sum of two squares $r_2(n)$. This $q$-analog is generated by the expansion of a product that was studied by Kronecker and Jordan. We generalize Jacobi's two squares formula from $r_2(n)$ to $r_2(n, q)$. We characterize the signs in the coefficients of $r_2(n, q)$ using the prime factors of $n$. We use $r_2(n, q)$ to characterize the integers which are the length of the hypotenuse of a primitive Pythagorean triangle.


title: Mass-deformed ABJM and Black Holes in AdS$_4$

abstract: We find a class of new supersymmetric dyonic black holes in four-dimensional maximal gauged supergravity which are asymptotic to the ${\rm SU(3)}\times {\rm U(1)}$ invariant AdS$_4$ Warner vacuum. These black holes can be embedded in eleven-dimensional supergravity where they describe the backreaction of M2-branes wrapped on a Riemann surface. The holographic dual description of these supergravity backgrounds is given by a partial topological twist on a Riemann surface of a three-dimensional $\mathcal{N}=2$ SCFT that is obtained by a mass-deformation of the ABJM theory. We compute explicitly the topologically twisted index of this SCFT and show that it accounts for the entropy of the black holes.


title: Numerical simulations of ICME-ICME interactions

abstract: We present hydrodynamical simulations of interacting Coronal Mass Ejections in the Interplanetary medium (ICMEs). In these events, two consecutive CMEs are launched from the Sun in similar directions within an interval of time of a few hours. In our numerical model, we assume that the ambient solar wind is characterized by its velocity and mass-loss rate. Then, the CMEs are generated when the flow velocity and mass-loss rate suddenly change, with respect to the ambient solar wind conditions during two intervals of time, which correspond to the durations of the CMEs. After their interaction, a merged region is formed and evolve as a single structure into the interplanetary medium. In this work, we are interested in the general morphology of this merged region, which depends on the initial parameters of the ambient solar wind and each of the CMEs involved. In order to understand this morphology, we have performed a parametric study in which we characterize the effects of the initial parameters variations on the density and velocity profiles at 1 AU, using as reference the well-documented event of July 25th, 2004. Based on this parametrization we were able to reproduce with a high accuracy the observed profiles. Then, we apply the parametrization results to the interaction events of May 23, 2010; August 1, 2010; and November 9, 2012. With this approach and using values for the input parameters within the CME observational errors, our simulated profiles reproduce the main features observed at 1 AU. Even though we do not take into account the magnetic field, our models give a physical insight into the propagation and interaction of ICMEs.


title: Convergence Analysis of Gradient Descent Algorithms with Proportional Updates

abstract: The rise of deep learning in recent years has brought with it increasingly clever optimization methods to deal with complex, non-linear loss functions. These methods are often designed with convex optimization in mind, but have been shown to work well in practice even for the highly non-convex optimization associated with neural networks. However, one significant drawback of these methods when they are applied to deep learning is that the magnitude of the update step is sometimes disproportionate to the magnitude of the weights (much smaller or larger), leading to training instabilities such as vanishing and exploding gradients. An idea to combat this issue is gradient descent with proportional updates. Gradient descent with proportional updates was introduced in 2017. It was independently developed by You et al (Layer-wise Adaptive Rate Scaling (LARS) algorithm) and by Abu-El-Haija (PercentDelta algorithm). The basic idea of both of these algorithms is to make each step of the gradient descent proportional to the current weight norm and independent of the gradient magnitude. It is common in the context of new optimization methods to prove convergence or derive regret bounds under the assumption of Lipschitz continuity and convexity. However, even though LARS and PercentDelta were shown to work well in practice, there is no theoretical analysis of the convergence properties of these algorithms. Thus it is not clear if the idea of gradient descent with proportional updates is used in the optimal way, or if it could be improved by using a different norm or specific learning rate schedule, for example. Moreover, it is not clear if these algorithms can be extended to other problems, besides neural networks. We attempt to answer these questions by establishing the theoretical analysis of gradient descent with proportional updates, and verifying this analysis with empirical examples.


title: Deep In-GPU Experience Replay

abstract: Experience replay allows a reinforcement learning agent to train on samples from a large amount of the most recent experiences. A simple in-RAM experience replay stores these most recent experiences in a list in RAM, and then copies sampled batches to the GPU for training. I moved this list to the GPU, thus creating an in-GPU experience replay, and a training step that no longer has inputs copied from the CPU. I trained an agent to play Super Smash Bros. Melee, using internal game memory values as inputs and outputting controller button presses. A single state in Melee contains 27 floats, so the full experience replay fits on a single GPU. For a batch size of 128, the in-GPU experience replay trained twice as fast as the in-RAM experience replay. As far as I know, this is the first in-GPU implementation of experience replay. Finally, I note a few ideas for fitting the experience replay inside the GPU when the environment state requires more memory.


title: Experimental Observation of Strong Coupling Between an Epsilon-Near-Zero Mode in a Deep Subwavelength Nanofilm and a Gap Plasmon Mode

abstract: Strong coupling is a phenomenon which occurs when the interaction between two resonance systems is so strong that the oscillatory energy exchange between them exceeds all dissipative loss channels. Each resonance can then no longer be described individually but only as a part of the coupled, hybrid system. Here, we show that strong coupling can occur in a deep subwavelength nanofilm supporting an epsilon near zero mode which is integrated into a metal-insulator-metal gap plasmon structure. To generate an epsilon near zero mode resonance in the short-wave infrared region, an indium tin oxide nanofilm of ~lambda/100 thickness is used. A polariton splitting value of 27%, corresponding to a normalized coupling rate of 0.135, is experimentally demonstrated. Simulations indicate that much larger coupling rates, well within the ultra-strong regime where the energy exchange rate is comparable with the frequency of light, are possible.


title: Detailed abundance analysis of globular clusters in the Local Group: NGC 147, NGC 6822, and Messier 33

abstract: We present new abundance measurements for eleven GCs in the Local Group galaxies NGC 147, NGC 6822, and Messier 33. These are combined with previously published observations of four GCs in the Fornax and WLM galaxies. The abundances were determined from analysis of integrated-light spectra, obtained with HIRES on the Keck I telescope and with UVES on the VLT. We find that the clusters with [Fe/H]<-1.5 are all alpha-enhanced at about the same level as Milky Way GCs. Their Na abundances are also generally enhanced relative to Milky Way halo stars, suggesting that these extragalactic GCs resemble their Milky Way counterparts in containing significant fractions of Na-rich stars. For [Fe/H]>-1.5, the GCs in M33 are also alpha-enhanced, while the GCs that belong to dwarfs (NGC 6822 SC7 and Fornax 4) have closer to Solar-scaled alpha-element abundances, thus mimicking the abundance trends observed in field stars in nearby dwarf galaxies. The abundance patterns in SC7 are remarkably similar to those in the Galactic GC Ruprecht 106, including significantly sub-solar [Na/Fe] and [Ni/Fe] ratios. In NGC 147, the GCs with [Fe/H]<-2.0 account for about 6% of the total luminosity of stars in the same metallicity range, a lower fraction than those previously found in the Fornax and WLM galaxies, but substantially higher than in the Milky Way halo.


title: Anatomy of triply-periodic network assemblies: Characterizing skeletal and inter-domain surface geometry of block copolymer gyroids

abstract: Triply-periodic networks (TPNs), like the well-known gyroid and diamond network phases, abound in soft matter assemblies, from block copolymers (BCPs), lyotropic liquid crystals and surfactants to functional architectures in biology. While TPNs are, in reality, volume-filling patterns of spatially-varying molecular composition, physical and structural models most often reduce their structure to lower-dimensional geometric objects: the {\it 2D interfaces} between chemical domains; and the {\it 1D skeletons} that thread through inter-connected, tubular domains. These lower-dimensional structures provide a useful basis of comparison to idealized geometries based on triply-periodic minimal, or constant-mean curvature surfaces, and shed important light on the spatially heterogeneous packing of molecular constituents that form the networks. Here, we propose a simple, efficient and flexible method to extract a 1D skeleton from 3D volume composition data of self-assembled networks. We apply this method to both self-consistent field theory predictions as well as experimental electron microtomography reconstructions of the double-gyroid phase of an ABA triblock copolymer. We further demonstrate how the analysis of 1D skeleton, 2D inter-domain surfaces, and combinations therefore, provide physical and structural insight into TPNs, across multiple length scales. Specifically, we propose and compare simple measures of {\it network chirality} as well as {\it domain thickness}, and analyze their spatial and statistical distributions in both ideal (theoretical) and non-ideal (experimental) double gyroid assemblies.


title: Topological freeness for $C^*$-correspondences

abstract: We study conditions that ensure uniqueness theorems of Cuntz-Krieger type for relative Cuntz-Pimsner algebras $\mathcal{O}(J,X)$ associated to a $C^*$-correspondence $X$ over a $C^*$-algebra $A$. We give general sufficient conditions phrased in terms of a multivalued map $\widehat{X}$ acting on the spectrum $\widehat{A}$ of $A$. When $X(J)$ is of Type I we construct a directed graph dual to $X$ and prove a uniqueness theorem using this graph. When $X(J)$ is liminal, we show that topological freeness of this graph is equivalent to the uniqueness property for $\mathcal{O}(J,X)$, as well as to an algebraic condition, which we call $J$-acyclicity of $X$. As an application we improve the Fowler-Raeburn uniqueness theorem for the Toeplitz algebra $\mathcal{T}_X$. We give new simplicity criteria for $\mathcal{O}_X$. We generalize and enhance uniqueness results for relative quiver $C^*$-algebras of Muhly and Tomforde. We also discuss applications to crossed products by endomorphisms.


title: Comparing heterogeneous entities using artificial neural networks of trainable weighted structural components and machine-learned activation functions

abstract: To compare entities of differing types and structural components, the artificial neural network paradigm was used to cross-compare structural components between heterogeneous documents. Trainable weighted structural components were input into machine-learned activation functions of the neurons. The model was used for matching news articles and videos, where the inputs and activation functions respectively consisted of term vectors and cosine similarity measures between the weighted structural components. The model was tested with different weights, achieving as high as 59.2% accuracy for matching videos to news articles. A mobile application user interface for recommending related videos for news articles was developed to demonstrate consumer value, including its potential usefulness for cross-selling products from unrelated categories.


title: Reconstruction of piecewise smooth wave speeds using multiple scattering

abstract: Let $c$ be a piecewise smooth wave speed on $\mathbb R^n$, unknown inside a domain $\Omega$. We are given the solution operator for the scalar wave equation $(\partial_t^2-c^2\Delta)u=0$, but only outside $\Omega$ and only for initial data supported outside $\Omega$. Using our recently developed scattering control method, we prove that piecewise smooth wave speeds are uniquely determined by this map, and provide a reconstruction formula. In other words, the wave imaging problem is solvable in the piecewise smooth setting under mild conditions. We also illustrate a separate method, likewise constructive, for recovering the locations of interfaces in broken geodesic normal coordinates using scattering control.


title: Visual and Semantic Knowledge Transfer for Large Scale Semi-supervised Object Detection

abstract: Deep CNN-based object detection systems have achieved remarkable success on several large-scale object detection benchmarks. However, training such detectors requires a large number of labeled bounding boxes, which are more difficult to obtain than image-level annotations. Previous work addresses this issue by transforming image-level classifiers into object detectors. This is done by modeling the differences between the two on categories with both image-level and bounding box annotations, and transferring this information to convert classifiers to detectors for categories without bounding box annotations. We improve this previous work by incorporating knowledge about object similarities from visual and semantic domains during the transfer process. The intuition behind our proposed method is that visually and semantically similar categories should exhibit more common transferable properties than dissimilar categories, e.g. a better detector would result by transforming the differences between a dog classifier and a dog detector onto the cat class, than would by transforming from the violin class. Experimental results on the challenging ILSVRC2013 detection dataset demonstrate that each of our proposed object similarity based knowledge transfer methods outperforms the baseline methods. We found strong evidence that visual similarity and semantic relatedness are complementary for the task, and when combined notably improve detection, achieving state-of-the-art detection performance in a semi-supervised setting.


title: Venus: The Making of an Uninhabitable World

abstract: The goals of the astrobiology community are focussed on developing a framework for the detection of biosignatures, or evidence thereof, on objects inside and outside of our solar system. A fundamental aspect of understanding the limits of habitable environments and detectable signatures is the study of where the boundaries of such environments can occur. Thus, the need to study the creation, evolution, and frequency of hostile environments for habitability is an integral part of the astrobiology story. These provide the opportunity to understand the bifurcation, between habitable and uninhabitable. The archetype of such a planet is the Earth's sister planet, Venus, and provides a unique opportunity to explore the processes that created a completely uninhabitable environment and thus define the conditions that can rule out bio-related signatures. We advocate a continued comprehensive study of our sister planet, including models of early atmospheres, compositional abundances, and Venus-analog frequency analysis from current and future exoplanet data. Moreover, new missions to Venus to provide in-situ data are necessary.


title: "Robust-squared" Imputation Models Using BART

abstract: Examples of "doubly robust" estimator for missing data include augmented inverse probability weighting (AIPWT) models (Robins et al., 1994) and penalized splines of propensity prediction (PSPP) models (Zhang and Little, 2009). Doubly-robust estimators have the property that, if either the response propensity or the mean is modeled correctly, a consistent estimator of the population mean is obtained. However, doubly-robust estimators can perform poorly when modest misspecification is present in both models (Kang and Schafer, 2007). Here we consider extensions of the AIPWT and PSPP models that use Bayesian Additive Regression Trees (BART; Chipman et al., 2010) to provide highly robust propensity and mean model estimation. We term these "robust-squared" in the sense that the propensity score, the means, or both can be estimated with minimal model misspecification, and applied to the doubly-robust estimator. We consider their behavior via simulations where propensities and/or mean models are misspecified. We apply our proposed method to impute missing instantaneous velocity (delta-v) values from the 2014 National Automotive Sampling System Crashworthiness Data System dataset and missing Blood Alcohol Concentration values from the 2015 Fatality Analysis Reporting System dataset. We found that BART applied to PSPP and AIPWT, provides a more robust and efficient estimate compared to PSPP and AIPWT, with the BART-estimated propensity score combined with PSPP providing the most efficient estimator with close to nominal coverage.


title: $\tau \to f_1 (1285)\pi^{-} \nu_{\tau}$ decay in the extended Nambu -- Jona-Lasinio model

abstract: Within the framework of the extended Nambu -- Jona-Lasinio model, we calculate the matrix element of the $\tau \to f_1(1285) \pi^{-} \nu_{\tau}$ decay, obtain the invariant mass distribution of the $f_1\pi$ -system and estimate the branching ratio Br$(\tau \to f_1 \pi^{-}\nu_{\tau})=4.0\times 10^{-4}$. The two types of contributions are considered: the contact interaction, and the axial-vector $I^G(J^{PC})=1^-(1^{++})$ resonance exchange. The latter includes the ground $a_1(1260)$ state, and its first radially excited state, $a_1(1640)$. The corrections caused by the $\pi -a_1$ transitions are taken into account. Our estimate is in a good agreement with the latest empirical result Br$(\tau \to f_1 \pi^{-} \nu_{\tau})=(3.9\pm 0.5)\times 10^{-4}$. The distribution function obtained for the decay $\tau \to f_1(1285) \pi^{-} \nu_{\tau}$ shows a clear signal of $a_1(1640)$ resonance which should be compared with future experimental data including our estimate of the decay width $\Gamma (a_1(1640) \to f_1 \pi )=14.1\,\mbox{MeV}$.


title: An overview of deep learning based methods for unsupervised and semi-supervised anomaly detection in videos

abstract: Videos represent the primary source of information for surveillance applications and are available in large amounts but in most cases contain little or no annotation for supervised learning. This article reviews the state-of-the-art deep learning based methods for video anomaly detection and categorizes them based on the type of model and criteria of detection. We also perform simple studies to understand the different approaches and provide the criteria of evaluation for spatio-temporal anomaly detection.


title: Moments in Time Dataset: one million videos for event understanding

abstract: We present the Moments in Time Dataset, a large-scale human-annotated collection of one million short videos corresponding to dynamic events unfolding within three seconds. Modeling the spatial-audio-temporal dynamics even for actions occurring in 3 second videos poses many challenges: meaningful events do not include only people, but also objects, animals, and natural phenomena; visual and auditory events can be symmetrical in time ("opening" is "closing" in reverse), and either transient or sustained. We describe the annotation process of our dataset (each video is tagged with one action or activity label among 339 different classes), analyze its scale and diversity in comparison to other large-scale video datasets for action recognition, and report results of several baseline models addressing separately, and jointly, three modalities: spatial, temporal and auditory. The Moments in Time dataset, designed to have a large coverage and diversity of events in both visual and auditory modalities, can serve as a new challenge to develop models that scale to the level of complexity and abstract reasoning that a human processes on a daily basis.


title: Domain-dependent stability analysis and parameter classification of a reaction-diffusion model on spherical geometries

abstract: In this work an activator-depleted reaction-diffusion system is investigated on polar coordinates with the aim of exploring the relationship and the corresponding influence of domain size on the types of possible diffusion-driven instabilities. Quantitative relationships are found in the form of necessary conditions on the area of a disk-shape domain with respect to the diffusion and reaction rates for certain types of diffusion-driven instabilities to occur. Robust analytical methods are applied to find explicit expressions for the eigenvalues and eigenfunctions of the diffusion operator on a disk-shape domain with homogenous Neumann boundary conditions in polar coordinates. Spectral methods are applied using chebyshev non-periodic grid for the radial variable and Fourier periodic grid on the angular variable to verify the nodal lines and eigen-surfaces subject to the proposed analytical findings. The full classification of the parameter space in light of the bifurcation analysis is obtained and numerically verified by finding the solutions of the partitioning curves inducing such a classification. Furthermore, analytical results are found relating the area of (disk-shape) domain with reaction-diffusion rates in the form of necessary conditions for the different types of bifurcations. These results are on one hand presented in the form of mathematical theorems with rigorous proofs, and, on the other hand using finite element method, each claim of the corresponding theorems are verified by obtaining the theoretically predicted behaviour of the dynamics in the numerical simulations.


title: Invariant measures for the periodic derivative nonlinear Schr\"odinger equation

abstract: We construct invariant measures associated to the integrals of motion of the periodic derivative nonlinear Schr\"odinger equation (DNLS) for small data in $L^2$ and we show these measures to be absolutely continuous with respect to the Gaussian measure. The key ingredient of the proof is the analysis of the gauge group of transformations associated to DNLS. As an intermediate step for our main result, we prove quasi-invariance with respect to the gauge maps of Gaussian measures on $L^2$.