title: Approximate thermofield dynamics of interacting fermions

abstract: We analyze the many-particle Schrodinger equation for fermions in a thermal ensemble by introducing an exponential operator expansion, defined in the context of thermofield dynamics. The expansion is optimized variationally at each time step through changes in the basis of excitations, which leads to a method of generating approximate differential equations to solve the time dependent problem, and can also be used to cool the system in imaginary time. The method is applied for a specific set of basis transformations and truncation scheme, leading to an explicit set of differential equations that reduce to the Hartree Fock solution in the low temperature limit. This procedure can also be generalized to include quantum correlation, which will be pursued in a future publication.


title: On the confidentiality of controller states under sensor attacks

abstract: With the emergence of cyber-attacks on control systems it has become clear that improving the security of control systems is an important task in today's society. We investigate how an attacker that has access to the measurements transmitted from the plant to the controller can perfectly estimate the internal state of the controller. This attack on sensitive information of the control loop is, on the one hand, a violation of the privacy, and, on the other hand, a violation of the security of the closed-loop system if the obtained estimate is used in a larger attack scheme. Current literature on sensor attacks often assumes that the attacker has already access to the controller's state. However, this is not always possible. We derive conditions for when the attacker is able to perfectly estimate the controller's state. These conditions show that if the controller has unstable poles a perfect estimate of the controller state is not possible. Moreover, we propose a defence mechanism to render the attack infeasible. This defence is based on adding uncertainty to the controller dynamics. We also discuss why an unstable controller is only a good defence for certain plants. Finally, simulations with a three-tank system verify our results.


title: QED corrections to the $^2P_{1/2}-{}^2P_{3/2}$ fine structure in fluorinelike ions: model Lamb shift operator approach

abstract: In [Li ${\it et \, al.}$, Phys. Rev. A ${\bf 98}$, 020502(R) (2018)] it was claimed that the model-potential computations of the Lamb shift on the $^2P_{1/2}-{}^2P_{3/2}$ fine structure in fluorinelike uranium lead to a discrepancy between theory and experiment. Later, it was reported by [Volotka ${\it et \, al.}$, Phys. Rev. A ${\bf 100}$, 010502(R) (2019)] that {\it ab initio} QED calculation, including the first-order one-electron QED contributions and the related effects of two-electron screening, yields the result which restores the agreement between theory and experiment and strongly disagrees with the model-potential Lamb shift values. In the present paper, the model Lamb shift operator [Shabaev ${\it et \, al.}$, Phys. Rev. A ${\bf 88}$, 012513 (2013)] is used to evaluate the QED effects on the $^2P_{1/2}-{}^2P_{3/2}$ fine structure in F-like ions. The calculations are performed by incorporating this operator into the Dirac-Coulomb-Breit equation employing different methods. It is demonstrated that the methods, based on including the Lamb shift operator either into the Dirac-Fock equations or into the calculations by perturbation theory, lead to the theoretical results which are in good agreement with each other and with experiment. The restriction of these results to the first order in the QED effects leads to a value which agrees with the aforementioned ${\it ab\, initio}$ QED result.


title: Density-based clustering of crystal orientations and misorientations and the orix python library

abstract: Crystal orientation mapping experiments typically measure orientations that are similar within grains and misorientations that are similar along grain boundaries. Such (mis)orientation data will cluster in (mis)orientation space and clusters are more pronounced if preferred orientations or special orientation relationships are present. Here, cluster analysis of (mis)orientation data is described and demonstrated using distance metrics incorporating crystal symmetry and the density-based clustering algorithm DBSCAN. Frequently measured (mis)orientations are identified as corresponding to grains, grain boundaries or orientation relationships, which are visualised both spatially and in three-dimensional (mis)orientation spaces. A new open-source python library, orix, is also reported.


title: Frequency chirped continuous-wave sodium laser guide stars

abstract: We numerically study a method to increase the photon return flux of continuous-wave laser guide stars using one-dimensional atomic cooling principles. The method relies on chirping the laser towards higher frequencies following the change in velocity of sodium atoms due to recoil, which raises atomic populations available for laser excitation within the Doppler distribution. The efficiency of this effect grows with the average number of atomic excitations between two atomic collisions in the mesosphere. We find the parameters for maximizing the return flux and evaluate the performance of chirping for operation at La Palma. According to our simulations, the optimal chirp rate lies between 0.8-1.0 MHz/$\mu$s and an increase in the fluorescence of the sodium guide star up to 60% can be achieved with current 20 W-class guide star lasers.


title: Spectral optimization for Robin Laplacian on domains admitting parallel coordinates

abstract: In this paper we deal with spectral optimization for the Robin Laplacian on a family of planar domains admitting parallel coordinates, namely a fixed-width strip built over a smooth closed curve and the exterior of a convex set with a smooth boundary. We show that if the curve length is kept fixed, the first eigenvalue referring to the fixed-width strip is for any value of the Robin parameter maximized by a circular annulus. Furthermore, we prove that the second eigenvalue in the exterior of a convex domain $\Omega$ corresponding to a negative Robin parameter does not exceed the analogous quantity for a disk whose boundary has a curvature larger than or equal to the maximum of that for $\partial\Omega$.


title: A note on Horvitz-Thompson estimators for rare subgroup analysis in the presence of interference

abstract: When there is interference, a subject's outcome depends on the treatment of others and treatment effects may take on several different forms. This situation arises often, particularly in vaccine evaluation. In settings where interference is likely, two-stage cluster randomized trials have been suggested as a means of estimating some of the causal contrast of interest. Working in the finite population setting to investigate rare and unplanned subgroup analyses using some of the estimators that have been suggested in the literature, include Horvitz-Thompson, Hajek, and what might be called the natural extension of the marginal estimators suggested in Hudgens and Halloran 2008. I define the estimands of interest conditional on individual, group and both individual and group baseline variables, giving unbiased Horvitz-Thompson style estimates for each. I also provide variance estimators for several estimators. I show that the Horvitz-Thompson (HT) type estimators are always unbiased provided at least one subject within the group or population, whatever the level of interest for the estimator, is in the subgroup of interest. This is not true of the "natural" or the Hajek style estimators, which will often be undefined for rare subgroups.


title: Signatures of first galaxies at FIR/sub-mm wavelengths

abstract: We explore the possibility of detecting first galaxies with a generic far-infrared/sub-millimeter telescope by applying an analytical model of primordial dust emission. As shown in previous works, galaxies at redshifts $z>7$ experience a strong negative K-correction in such a way that systems of similar masses are brighter at higher $z$. In addition, at a given mass and $z$, our model predicts that luminosity fluxes increase proportionally to the dust-to-metal ratios ($D/M$) of primeval sources. We evaluate the observability of model sources at different $z>7$ as a function of the observed survey area ($\Delta \Omega$) and sensitivity ($S$) of a generic instrument. Assuming $\Delta \Omega \sim 10~{\rm deg}^2$ and a plausible $S \sim 1~{\mu}{\rm Jy}$ for a near future survey, we could assure the detection of at least one typical source with $D/M \sim 5\times10^{-3}$ at $z>7$. For $S\gtrsim 1~{\mu}{\rm Jy}$ and $\Delta \Omega \lesssim 10~{\rm deg}^2$, higher than typical $D/M$ are required to detect at least one individual source at $z>7$. The observability of model galaxies is also affected by the size distribution of dust grains, specially towards higher $z$.


title: Vibrational Coupling to Epsilon-Near-Zero Waveguide Modes

abstract: Epsilon near zero modes offer extreme field enhancement that can be utilized for developing enhanced sensing schemes. However, demonstrations of enhanced spectroscopies have largely exploited surface polaritons, mostly due to the challenges of coupling a vibrational transition to volume-confined epsilon near zero modes. Here we fabricate high aspect ratio gratings (up to 24.8 um height with greater than 5 {\mu}m pitch) of 4H-SiC, with resonant modes that couple to transverse magnetic and transverse electric incident fields. These correspond to metal-insulator-metal waveguide modes propagating downwards into the substrate. The cavity formed by the finite length of the waveguide allows for strong absorption of incident infrared light (>80%) with Q factors in excess of 90, including an epsilon near zero waveguide mode with {\epsilon}eff=0.0574+0.008i. The localization of the electromagnetic fields within the gap between the grating teeth suggests an opportunity to realize a new platform for studying vibrational coupling in liquid environments, with potential opportunities for enhanced spectroscopies. We show that these modes are supported in anhydrous and aqueous environments, and that high aspect ratio gratings coherently couple to the vibrational transition in the surrounding liquid.


title: Searching for optical companions to four binary millisecond pulsars with the Gran Telescopio Canarias

abstract: We report on multi-band photometric observations of four binary millisecond pulsars with the Gran Telescopio Canarias. The observations led to detection of binary companions to PSRs J1630+3734, J1741+1351 and J2042+0246 in the Sloan g', r' and i' bands. Their magnitudes in the r' band are $\approx$24.4, 24.4 and 24.0, respectively. We also set a 3$\sigma$ upper limit on the brightness of the PSR J0557+1550 companion in the r' band of $\approx$25.6 mag. Combining the optical data with the radio timing measurements and white dwarf cooling models, we show that the detected companions are cool low-mass white dwarfs with temperatures and ages in the respective ranges of (4-7)$\times 10^3$ K and 2-5 Gyr. All the detected white dwarfs are found to likely have either pure hydrogen or mixed helium-hydrogen atmospheres.


title: Optical polarization properties of AGN with significant VLBI-Gaia offsets

abstract: Significant positional offsets of the value from 1 mas to more than 10 mas were found previously between radio (VLBI) and optical (Gaia) positions of active galactic nuclei (AGN). They happen preferentially parallel to the parsec-scale jet direction. AGN with VLBI-to-Gaia offsets pointed downstream the jet are found to have favourably higher optical polarization, as expected if extended optical jets dominate in the emission and shift the Gaia centroid away from the physical nucleus of the source. Upstream offsets with the suggested domination of accretion disks manifest themselves through the observed low optical polarization. Direction of linear optical polarization is confirmed to preferentially align with parsec-scale jets in AGN with dominant jets consistent with a toroidal magnetic field structure. Our findings support the disk-jet interpretation of the observed positional offsets. These results call on an intensification of AGN optical polarization monitoring programs in order to collect precious observational data. Taken together with the continued VLBI and Gaia observations, they will allow researchers to reconstruct detailed models of the disk-jet system in AGN on parsec scales.


title: Statistics of relative velocity for particles settling under gravity in a turbulent flow

abstract: We study the joint probability distributions of separation, $R$, and radial component of the relative velocity, $V_{\rm R}$, of particles settling under gravity in a turbulent flow. We also obtain the moments of these distributions and analyze their anisotropy using spherical harmonics. We find that the qualitative nature of the joint distributions remains the same as no gravity case. Distributions of $V_{\rm R}$ for fixed values of $R$ show a power-law dependence on $V_{\rm R}$ for a range of $V_{\rm R}$, exponent of the power-law depends on the gravity. Effects of gravity are also manifested in the following ways: (a) moments of the distributions are anisotropic; the degree of anisotropy depends on particle's Stokes number, but does not depend on $R$ for small values of $R$. (b) mean velocity of collision between two particles is decreased for particles having equal Stokes numbers but increased for particles having different Stokes numbers. For the later, collision velocity is set by the difference in their settling velocities.


title: A new and unusual LBV-like outburst from a Wolf-Rayet star in the outskirts of M33

abstract: MCA-1B (also called UIT003) is a luminous hot star in the western outskirts of M33, classified over 20yr ago with a spectral type of Ofpe/WN9 and identified then as a candidate luminous blue variable (LBV). Palomar Transient Factory data reveal that this star brightened in 2010, with a light curve resembling that of the classic LBV star AFAnd in M31. Other Ofpe/WN9 stars have erupted as LBVs, but MCA-1B was unusual because it remained hot. It showed a WN-type spectrum throughout its eruption, whereas LBVs usually get much cooler. MCA-1B showed an almost four-fold increase in bolometric luminosity and a doubling of its radius, but its temperature stayed around 29kK. As it faded, it shifted to even hotter temperatures, exhibiting a WN7/WN8-type spectrum, and doubling its wind speed. MCA-1B is reminiscent of some supernova impostors, and its location resembles the isolated environment of SN 2009ip. It is most similar to HD5980 (in the SMC) and GR 290 (also in M33). Whereas these two LBVs exhibited B-type spectra in eruption, MCA-1B is the first clear case where a Wolf-Rayet (WR) spectrum persisted at all times. Together, MCA-1B, HD 5980, and GR 290 constitute a class of WN-type LBVs, distinct from S Doradus LBVs. They are most interesting in the context of LBVs at low metallicity, a possible post-LBV/WR transition in binaries, and as likely Type~Ibn supernova progenitors.


title: Multispectral Photon-Counting for Medical Imaging and Beam Characterization

abstract: We present the current status of our project of developing a photon counting detector for medical imaging. An example motivation lays in producing a monitoring and dosimetry device for boron neutron capture therapy, currently not commercially available. Our approach combines in-house developed detectors based on cadmium telluride or thick silicon with readout chip technology developed for particle physics experiments at CERN. Here we describe the manufacturing process of our sensors as well as the processing steps for the assembly of first prototypes. The prototypes use currently the PSI46digV2.1-r readout chip. The accompanying readout electronics chain that was used for first measurements will also be discussed. Finally we present an advanced algorithm developed by us for image reconstruction using such photon counting detectors with focus on boron neutron capture therapy. This work is conducted within a consortium of Finnish research groups from Helsinki Institute of Physics, Aalto University, Lappeenranta-Lahti University of Technology LUT and Radiation and Nuclear Safety Authority (STUK) under the RADDESS program of Academy of Finland. Keywords: Solid state detectors, X-ray detectors, Gamma detectors, Neutron detectors, Instrumentation for hadron therapy, Medical-image reconstruction methods and algorithms.


title: Easy Capacitated Facility Location Problems, with Connections to Lot-Sizing

abstract: In this note, we consider the capacitated facility location problem when the transportation costs of the instance satisfy the Monge property. We show that a straightforward dynamic program finds the optimal solution when the demands are polynomially bounded. When demands are not polynomially bounded, we give a fully polynomial-time approximation scheme by adapting an algorithm and analysis of Van Hoesel and Wagelmans.


title: Learning Generative Models using Denoising Density Estimators

abstract: Learning probabilistic models that can estimate the density of a given set of samples, and generate samples from that density, is one of the fundamental challenges in unsupervised machine learning. We introduce a new generative model based on denoising density estimators (DDEs), which are scalar functions parameterized by neural networks, that are efficiently trained to represent kernel density estimators of the data. Leveraging DDEs, our main contribution is a novel technique to obtain generative models by minimizing the KL-divergence directly. We prove that our algorithm for obtaining generative models is guaranteed to converge to the correct solution. Our approach does not require specific network architecture as in normalizing flows, nor use ordinary differential equation solvers as in continuous normalizing flows. Experimental results demonstrate substantial improvement in density estimation and competitive performance in generative model training.


title: A Geometric Interpretation of the Quadratic Formula

abstract: This article provides a simple geometric interpretation of the quadratic formula. The geometry helps to demystify the formula's complex appearance and casts it into a much simpler existence, thus potentially benefits early algebra students.


title: Enhanced attraction between particles in a bidisperse mixture with random pair-wise interactions

abstract: Motivated by growing interests in multicomponent metallic alloys and complex fluids, we study a complex mixture with bidispersity in size and polydispersity in energy. The energy polydispersity in the system is introduced by considering random pair-wise interactions between the particles. Extensive molecular dynamics simulations are performed to compute potential energy and neighborhood identity ordering (NIO) parameter as a function of temperature for a wide range of parameters including size-ratio and concentration of the two species by quenching it from a high temperature fluid state to a crystalline state. Our findings demonstrate an enhancement of the neighborhood identity ordering on addition of particles of different sizes. Moreover, a comparatively higher increase in NIO parameter is achieved by tuning the size-ratio of the particles. We also propose NIO parameter to be a good marker to differentiate systems (below the liquid-to-solid transition temperature) having different values of size-ratio and concentrations. Effect of cooling rates on NIO parameter is also discussed.


title: SirenLess: reveal the intention behind news

abstract: News articles tend to be increasingly misleading nowadays, preventing readers from making subjective judgments towards certain events. While some machine learning approaches have been proposed to detect misleading news, most of them are black boxes that provide limited help for humans in decision making. In this paper, we present SirenLess, a visual analytical system for misleading news detection by linguistic features. The system features article explorer, a novel interactive tool that integrates news metadata and linguistic features to reveal semantic structures of news articles and facilitate textual analysis. We use SirenLess to analyze 18 news articles from different sources and summarize some helpful patterns for misleading news detection. A user study with journalism professionals and university students is conducted to confirm the usefulness and effectiveness of our system.


title: Exploring the acceleration mechanisms for particle injection and power-law formation during trans-relativistic magnetic reconnection

abstract: Magnetic reconnection in the relativistic and trans-relativistic regimes is able to accelerate particles to hard power law energy spectra $f \propto \gamma^{-p}$ (approaching $p=1$). The underlying acceleration mechanism that determines the spectral shape is currently a topic of intense investigation. By means of fully kinetic plasma simulations, we carry out a study of particle acceleration during magnetic reconnection in the trans-relativistic regime of a proton-electron plasma. While earlier work in this parameter regime has focused on the effects of electric field parallel to the local magnetic field on the particle injection (from thermal energy to the lower energy bound of the power-law spectrum), here we examine the roles of both parallel and perpendicular electric fields to gain a more complete understanding on the injection process and further development of a power-law spectrum. We show that the parallel electric field does contribute significantly to particle injection, and is more important in the initial phase of magnetic reconnection. However, as the simulation proceeds, the acceleration by the perpendicular electric field becomes more important for particle injection and completely dominates the acceleration responsible for the high-energy power-law spectrum. This holds robustly, in particular for longer reconnection times and larger systems, i.e. in simulations that are more indicative of the processes in astrophysical sources.


title: Practical method to reclassify Web of Science articles into unique subject categories and broad disciplines

abstract: Classification of bibliographic items into subjects and disciplines in large databases is essential for many quantitative science studies. The Web of Science classification of journals into ~250 subject categories, which has served as a basis for many studies, is known to have some fundamental problems and several practical limitations that may affect the results from such studies. Here we present an easily reproducible method to perform reclassification of the Web of Science into existing subject categories and into 14 broad areas. Our reclassification is at a level of articles, so it preserves disciplinary differences that may exist among individual articles published in the same journal. Reclassification also eliminates ambiguous (multiple) categories that are found for 50% of items, and assigns a discipline/field category to all articles that come from broad-coverage journals such as Nature and Science. The correctness of the assigned subject categories is evaluated manually and is found to be ~95%.


title: Gasoline Pricing Policies for Transportation Safety

abstract: Economic factors can have substantial effects on transportation crash trends. This study makes a comprehensive examination of the relationship between the retail gasoline price (including state and federal fuel taxes) and transportation fatal crashes from 2007 to 2016 in the US. Data on motor vehicle, bicycle and pedestrian fatal crashes come from Fatality Analysis Reporting System (FARS) provided by the National Highway Safety Administration (NHTSA) and the gasoline price data is from U.S. Energy Information Administration (EIA). Random effect negative binomial regression models are used to estimate the impact of inflation-adjusted gasoline prices on trends of transportation fatal crashes. Initial results combined with results of previous studies showed that gender and transportation mean type (motorcycle, non-motorcycle, bicycle and pedestrian) play prominent roles in interpreting the final model, so by using random effect negative binomial regression, seven models are developed to evaluate the effects of gasoline price changes on total population, male, female, motorcyclists, non-motorcyclists, bicyclists and pedestrians separately. Our findings suggest that increasing the gasoline prices will not significantly alter the number of total fatal crashes. However, by looking at different vehicle types, it is estimated that one dollar increase in adjusted gasoline price is associated with 24.2% increase in the number of motorcycle fatal crashes, 1.9% decrease in the number of non-motorcycle fatal crashes, and 0.7% decrease in the number of pedestrian fatal crashes. Also, there is no noticeable difference between male and female in response to the gasoline price changes.


title: Complex Solutions to Bessel SDEs and SLEs

abstract: We consider a variant of Bessel SDE by allowing the solution to be complex valued. Such SDEs appear naturally while studying the trace of Schramm-Loewner-Evolutions (SLE). We establish the existence and uniqueness of the strong solution to such SDEs when the dimension is negative. We also consider the stochastic flow associated to such SDEs and prove that it is almost surely continuous. Our proofs are based on an improvement of the derivative estimate of Rohde-Schramm \cite{RS05}. We finally show the connection between such stochastic flows and SLE$_{\kappa}$ for $\kappa <4$.


title: LHC dijet angular distributions as a probe for the dimension-six triple gluon vertex

abstract: In the absence of a concrete discovery of new physics at the LHC, global analyses of the standard model effective field theory (SMEFT) are important to find and describe the impact of new physics beyond the energy reach of the LHC. Among the SMEFT operators that can be constrained via various measurements, the dimension six triple gluon operator involves neither the Higgs boson nor the top quark, yet its variation can have measurable effects on top and Higgs production. Without independent constraints on its impact, the sensitivity of measurements in the top and Higgs sectors to new physics is reduced. We show that the dijet angular distribution is a powerful observable for probing the triple gluon operator. We set the most stringent limit on the triple gluon effective coupling by reinterpreting the results of a search for new phenomena in dijet events using 35.9 fb$^{-1}$ of pp collision data collected at $\sqrt{s}$ = 13 TeV performed by the CMS collaboration. The obtained limit on the strength of the triple gluon operator is far below the sensitivity of that which can be derived from top quark and Higgs measurements and thus this operator can be neglected in global SMEFT analyses.


title: Shadowing and Stability in p-adic dynamics

abstract: In this paper, we study dynamical properties as shadowing and structural stability for a class of dynamics on $\mathbb{Z}_p$ and $\mathbb{Q}_p$, where $p \geq 2$ is a prime number. In particular, we prove that if $f: \mathbb{Z}_p \to \mathbb{Z}_p$ is a $(p^{-k},p^{m})$ ( $0 < m \leq k$ integers ) locally scaling map then $f$ is shadowing and structurally stable. We also study the number of conjugacy classes of these maps and we consider the above properties for $1$-Lipschitz maps of $\mathbb{Z}_p$ and for extensions of the shift map, contractions and dilatations on $\mathbb{Q}_p$.


title: Liquid-liquid microphase separation leads to formation of membraneless organelles

abstract: Proteins and nucleic acids can spontaneously self-assemble into membraneless droplet-like compartments, both in vitro and in vivo. A key component of these droplets are multi-valent proteins that possess several adhesive domains with specific interaction partners (whose number determines total valency of the protein) separated by disordered regions. Here, using multi-scale simulations we show that such proteins self-organize into micro-phase separated droplets of various sizes as opposed to the Flory-like macro-phase separated equilibrium state of homopolymers or equilibrium physical gels. We show that the micro-phase separated state is a dynamic outcome of the interplay between two competing processes: a diffusion-limited encounter between proteins, and the dynamics within small clusters that results in exhaustion of available valencies whereby all specifically interacting domains find their interacting partners within smaller clusters, leading to arrested phase separation. We first model these multi-valent chains as bead-spring polymers with multiple adhesive domains separated by semi-flexible linkers and use Langevin Dynamics (LD) to assess how key timescales depend on the molecular properties of associating polymers. Using the time-scales from LD simulations, we develop a coarse-grained kinetic model to study this phenomenon at longer times. Consistent with LD simulations, the macro-phase separated state was only observed at high concentrations and large interaction valencies. Further, in the regime where cluster sizes approach macro-phase separation, the condensed phase becomes dynamically solid-like, suggesting that it might no longer be biologically functional. Therefore, the micro-phase separated state could be a hallmark of functional droplets formed by proteins with the sticker-spacer architecture.


title: High-precision four-loop mass and wave function renormalization in QED

abstract: The 4-loop QED mass and wave function renormalization constants $Z_2$ and $Z_m$ have been evaluated in the on-shell subtraction scheme with 1100 digits of precision. We also worked out the coefficients of the five color structures of the QCD renormalization constants $Z_2^{OS}$ and $Z_m^{OS}$ which can obtained from QED-like diagrams. The results agree with lower precision results available in the literature. Analytical fits were also obtained for all these quantities.


title: Wind profiling in the lower atmosphere from wind-induced perturbations to multirotor UAS

abstract: We present a model-based approach to wind velocity profiling using motion perturbations of a multirotor unmanned aircraft system (UAS) in both hovering and steady ascending flight. A state estimation framework was adapted to a set of closed-loop rigid body models identified for an off-the-shelf quadrotor. The quadrotor models used for wind estimation were characterized for hovering and steady ascending flight conditions ranging between 0 and 2 m/s. The closed-loop models were obtained using system identification algorithms to determine model structures and estimate model parameters. The wind measurement method was validated experimentally above the Virginia Tech Kentland Experimental Aircraft Systems Laboratory by comparing quadrotor and independent sensor measurements from a sonic anemometer and two SoDARs. Comparison results demonstrated quadrotor wind estimation in close agreement with the independent wind velocity measurements. Wind velocity profiles were difficult to validate using time-synchronized SoDAR measurements, however. Analysis of the noise intensity and signal-to-noise ratio of the SoDARs proved that close-proximity quadrotor operations can corrupt wind measurement from SoDARs.


title: Toward Generalized Clustering through an One-Dimensional Approach

abstract: After generalizing the concept of clusters to incorporate clusters that are linked to other clusters through some relatively narrow bridges, an approach for detecting patches of separation between these clusters is developed based on an agglomerative clustering, more specifically the single-linkage, applied to one-dimensional slices obtained from respective feature spaces. The potential of this method is illustrated with respect to the analyses of clusterless uniform and normal distributions of points, as well as a one-dimensional clustering model characterized by two intervals with high density of points separated by a less dense interstice. This partial clustering method is then considered as a means of feature selection and cluster identification, and two simple but potentially effective respective methods are described and illustrated with respect to some hypothetical situations.


title: Self-Consistent Quantum-Field Theory for the Characterization of Complex Random Media by Short Laser Pulses

abstract: We present a quantum field theoretical method for the characterization of disordered complex media with short laser pulses in an optical coherence tomography setup (OCT). We solve this scheme of coherent transport in space and time with weighted essentially nonoscillatory methods (WENO). WENO is preferentially used for the determination of highly nonlinear and discontinuous processes including interference effects and phase transitions like Anderson localization of light. The theory determines spatiotemporal characteristics of the scattering mean free path and the transmission cross section that are directly measurable in time-of-flight (ToF) and pump-probe experiments. The results are a measure of the coherence of multiple scattering photons in passive as well as in optically soft random media. Our theoretical results of ToF are instructive in spectral regions where material characteristics such as the scattering mean free path and the diffusion coefficient are methodologically almost insensitive to gain or absorption and to higher-order nonlinear effects. Our method is applicable to OCT and other advanced spectroscopy setups including samples of strongly scattering mono- and polydisperse complex nano- and microresonators.

categories: cond-mat.dis-nn cond-mat.mes-hall cond-mat.soft physics.comp-ph physics.optics


title: Rank-one Detector for Kronecker-Structured Constant Modulus Constellations

abstract: To achieve a reliable communication with short data blocks, we propose a novel decoding strategy for Kronecker-structured constant modulus signals that provides low bit error ratios (BERs) especially in the low energy per bit to noise power spectral density ratio $(E_b/N_0)$. The encoder exploits the fact that any M-PSK constellation can be factorized as Kronecker products of lower or equal order PSK constellation sets. A construction of two types of schemes is first derived. For such Kronecker-structured schemes, a conceptually simple decoding algorithm is proposed, referred to as Kronecker-RoD (rank-one detector). The decoder is based on a rank-one approximation of the "tensorized" received data block, has a built-in noise rejection capability and a smaller implementation complexity than state-of-the-art detectors. Compared with convolutional codes with hard and soft Viterbi decoding, Kronecker-RoD outperforms the latter in BER performance at same spectral efficiency.


title: Contribution of high risk groups' unmet needs may be underestimated in epidemic models without risk turnover: a mechanistic modelling analysis

abstract: BACKGROUND. Epidemic models of STIs are often used to characterize the contribution of risk groups to overall transmission by projecting the transmission population attributable fraction (tPAF) of unmet prevention and treatment needs within risk groups. However, evidence suggests that STI risk is dynamic over an individual's sexual life course, which manifests as turnover between risk groups. We sought to examine the mechanisms by which turnover influences modelled projections of the tPAF of high risk groups. METHODS. We developed a unifying, data-guided framework to simulate risk group turnover in deterministic, compartmental transmission models. We applied the framework to an illustrative model of an STI and examined the mechanisms by which risk group turnover influenced equilibrium prevalence across risk groups. We then fit a model with and without turnover to the same risk-stratified STI prevalence targets and compared the inferred level of risk heterogeneity and tPAF of the highest risk group projected by the two models. RESULTS. The influence of turnover on group-specific prevalence was mediated by three main phenomena: movement of infectious individuals between risk groups; changes to herd immunity; and changes in discordant partnerships. Faster turnover led to a smaller ratio of STI prevalence between the highest and lowest risk groups. Compared to the fitted model without turnover, the fitted model with turnover inferred greater risk heterogeneity and consistently projected a larger tPAF of the highest risk group over time. IMPLICATIONS. If turnover is not captured in epidemic models, the projected contribution of high risk groups, and thus, the potential impact of prioritizing interventions to address their needs, could be underestimated. To aid the next generation of tPAF models, data collection efforts to parameterize risk group turnover should be prioritized.


title: Sums of Squared Distances between Points on a Unit $n$-Sphere

abstract: In this paper, we prove two theorems concerning the sums of squared distances between points on a unit $n$-sphere that generalize two facts previously known about the case where the points are the vertices of a regular polygon. The first theorem is that, given a multiset of $V$ points on a unit $n$-sphere, the sum of the squared distances between these points is $V^2 ( 1 - d^2 )$ where $d$ is the distance between the centroid of the points and the center of the unit $n$-sphere (for any $n \geq 2$). The second is that, given a finite set of points on the unit $n$-sphere centered at the origin such that the point set is symmetric about the origin and the symmetry group of the point set acts transitively on the set, the sum of the squared distinct distances between these points is $2k + 2$ where $k$ is the number of distinct distances between the points (for any $n \geq 2$). Using the first theorem, we find a new way to calculate the potential energy function of a finite normalized frame.


title: Spinneret: Aiding Creative Ideation through Non-Obvious Concept Associations

abstract: Mind mapping is a popular way to explore a design space in creative thinking exercises, allowing users to form associations between concepts. Yet, most existing digital tools for mind mapping focus on authoring and organization, with little support for addressing the challenges of mind mapping such as stagnation and design fixation. We present Spinneret, a functional approach to aid mind mapping by providing suggestions based on a knowledge graph. Spinneret uses biased random walks to explore the knowledge graph in the neighborhood of an existing concept node in the mind map, and provides "suggestions" for the user to add to the mind map. A comparative study with a baseline mind-mapping tool reveals that participants created more diverse and distinct concepts with Spinneret, and reported that the suggestions inspired them to think of ideas they would otherwise not have explored.


title: Spectroscopy of mesons with bottom quarks

abstract: Preliminary results for the spectra of excited and exotic bottom mesons are presented. The calculation on a dynamical anisotropic ensemble exploits distillation, enabling the use of a large basis of interpolating operators including those proportional to the gluonic field strength which are relevant for hybrid states. A comparison of results with similar calculations in the light and charm sectors is discussed.


title: Rate-Constrained Shaping Codes for Structured Sources

abstract: Shaping codes are used to encode information for use on channels with cost constraints. Applications include data transmission with a power constraint and, more recently, data storage on flash memories with a constraint on memory cell wear. In the latter application, system requirements often impose a rate constraint. In this paper, we study rate-constrained fixed-to-variable length shaping codes for noiseless, memoryless costly channels and general i.i.d. sources. The analysis relies on the theory of word-valued sources. We establish a relationship between the code expansion factor and minimum average symbol cost. We then determine the expansion factor that minimizes the average cost per source symbol (total cost), corresponding to a conventional optimal source code with cost. An equivalence is established between codes minimizing average symbol cost and codes minimizing total cost, and a separation theorem is proved, showing that optimal shaping can be achieved by a concatenation of optimal compression and optimal shaping for a uniform i.i.d. source. Shaping codes often incorporate, either explicitly or implicitly, some form of non-equiprobable signaling. We use our results to further explore the connections between shaping codes and codes that map a sequence of i.i.d. source symbols into an output sequence of symbols that are approximately independent and distributed according to a specified target distribution, such as distribution matching (DM) codes. Optimal DM codes are characterized in terms of a new performance measure - generalized expansion factor (GEF) - motivated by the costly channel perspective. The GEF is used to study DM codes that minimize informational divergence and normalized informational divergence.


title: Giant shift upon strain on the fluorescence spectrum of V$_{\rm N}$N$_{\rm B}$ color centers in $h$-BN

abstract: We study the effect of strain on the physical properties of the nitrogen antisite-vacancy pair in hexagonal boron nitride ($h$-BN), a color center that may be employed as a quantum bit in a two-dimensional material. With group theory and ab-initio analysis we show that strong electron-phonon coupling plays a key role in the optical activation of this color center. We find a giant shift on the zero-phonon-line (ZPL) emission of the nitrogen antisite-vacancy pair defect upon applying strain that is typical of $h$-BN samples. Our results provide a plausible explanation for the experimental observation of quantum emitters with similar optical properties but widely scattered ZPL wavelengths and the experimentally observed dependence of the ZPL on the strain.


title: Effects of fluctuations on correlation functions in inhomogeneous mixtures

abstract: Approximate expressions for correlation functions in binary inhomogeneous mixtures are derived in a framework of the mesoscopic theory [Ciach A., Mol. Phys., 2011, {\textbf{109}}, 1101]. Fluctuation contribution is taken into account in a Brazovskii-type approximation. Explicit results are obtained for two model systems. In the two models, the diameters of the hard cores of particles are equal, and the interactions favour a periodic arrangement of alternating species A and B. However, the optimal distance between the species A and B is much different in the two models. Theoretical results for different temperature and volume fractions of the two components are compared with the results of Monte Carlo simulations, and the structure is illustrated by simulation snapshots. Despite different interaction potentials and different length scale of the local ordering, properties of the correlation functions in the two models are very similar.


title: Complete regularity of Ellis semigroups of $\mathbb Z$-actions

abstract: It is shown that the Ellis semigroup of a $\mathbb Z$-action on a compact totally disconnected space is completely regular if and only if forward proximality coincides with forward asymptoticity and backward proximality coincides with backward asymptoticity. Furthermore, the Ellis semigroup of a $\mathbb Z$- or $\mathbb R$-action for which forward proximality and backward proximality are transitive relations is shown to have at most two left minimal ideals. Finally, the notion of near simplicity of the Ellis semigroup is introduced and related to the above.


title: Global Prompt Proton Sensor Network: Monitoring Solar Energetic Protons based on GPS Satellite Constellation

abstract: Energetic particle instruments on board GPS satellites form a powerful global prompt proton sensor network (GPPSn) that provides an unprecedented opportunity to monitor and characterize solar energetic protons targeting the Earth. The medium-Earth-orbits of the GPS constellation have the unique advantage of allowing solar energetic protons to be simultaneously measured from multiple points in both open- and closed-field line regions. Examining two example intervals of solar proton events, we showcase in this study how GPS proton data are prepared, calibrated and utilized to reveal important features of solar protons, including their source, acceleration/scattering by interplanetary shocks, the relative position of Earth when impinged by these shocks, the shape of solar particle fronts, the access of solar protons inside the dynamic geomagnetic field, as well temporally-varying proton distributions in both energy and space. By comparing to Van Allen Probes data, GPS proton observations are further demonstrated not only to be useful for qualitatively monitoring the dynamics of solar protons, but also for quantitative scientific research including determining cutoff L-shells. Our results establish that this GPPSn can join forces with other existing solar proton monitors and contribute to observing, warning, understanding and ultimately forecasting the incoming solar energetic proton events.


title: Locating conical degeneracies in the spectra of parametric self-adjoint matrices

abstract: A simple iterative scheme is proposed for locating the parameter values for which a 2-parameter family of real symmetric matrices has a double eigenvalue. The convergence is proved to be quadratic. An extension of the scheme to complex Hermitian matrices (with 3 parameters) and to location of triple eigenvalues (5 parameters for real symmetric matrices) is also described. Algorithm convergence is illustrated in several examples: a real symmetric family, a complex Hermitian family, a family of matrices with an "avoided crossing" (no covergence) and a 5-parameter family of real symmetric matrices with a triple eigenvalue.


title: Anisotropic elliptic equations with gradient-dependent lower order terms and $L^1$ data

abstract: We prove the existence of a weak solution for a general class of Dirichlet anisotropic elliptic problems such as $\mathcal Au+\Phi(x,u,\nabla u)=\mathfrak{B}u+f$ in $\Omega$, where $\Omega$ is a bounded open subset of $\mathbb R^N$ and $f\in L^1(\Omega)$ is arbitrary. The principal part is a divergence-form nonlinear anisotropic operator $\mathcal A$, the prototype of which is $\mathcal A u=-\sum_{j=1}^N \partial_j(|\partial_j u|^{p_j-2}\partial_j u)$ with $p_j>1$ for all $1\leq j\leq N$ and $\sum_{j=1}^N (1/p_j)>1$. As a novelty in this paper, our lower order terms involve a new class of operators $\mathfrak B$ such that $\mathcal{A}-\mathfrak{B}$ is bounded, coercive and pseudo-monotone from $W_0^{1,\overrightarrow{p}}(\Omega)$ into its dual, as well as a gradient-dependent nonlinearity $\Phi$ with an "anisotropic natural growth" in the gradient and a good sign condition.


title: Coherent States on the Circle: Semiclassical Matrix Elements in the Context of Kummer Functions and the Zak transformation

abstract: We extend former results for coherent states on the circle in the literature in two ways. On the one hand, we show that expectation values of fractional powers of momentum operators can be computed exactly analytically by means of Kummer's confluent hypergeometric functions. Earlier, these expectation values have only been obtained by using suitable estimates. On the other hand, we consider the Zak transformation not only to map harmonic oscillator coherent states to coherent states on the circle as it has been discussed before, but we also use the properties of the Zak transformation to derive a relation between matrix elements with respect to coherent states in L2(R) and L2(S1). This provides an alternative way for computing semiclassical matrix elements for coherent states on the circle. In certain aspects, this method simplifies the semiclassical computations in particular if one is only interested in the classical limit, that is the zeroth order term in the semiclassical expansion.


title: George Sudarshan and Quantum Dynamics

abstract: These are my recollections of working with George Sudarshan from 2002 to 2008 when I was a PhD student in his group. During these years I learnt a lot of physics and also witness to some remarkable occurrences.


title: On the Performance of PDCCH in LTE and 5G New Radio

abstract: 5G New Radio (NR) Release 15 has been specified in June 2018. It introduces numerous changes and potential improvements for physical layer data transmissions, although only point-to-point (PTP) communications are considered. In order to use physical data channels such as the Physical Downlink Shared Channel (PDSCH), it is essential to guarantee a successful transmission of control information via the Physical Downlink Control Channel (PDCCH). Taking into account these two aspects, in this paper, we first analyze the PDCCH processing chain in NR PTP as well as in the state-of-the-art Long Term Evolution (LTE) point-to-multipoint (PTM) solution, i.e., evolved Multimedia Broadcast Multicast Service (eMBMS). Then, via link level simulations, we compare the performance of the two technologies, observing the Bit/Block Error Rate (BER/BLER) for various scenarios. The objective is to identify the performance gap brought by physical layer changes in NR PDCCH as well as provide insightful guidelines on the control channel configuration towards NR PTM scenarios.


title: Pioneering Studies on LTE eMBMS: Towards 5G Point-to-Multipoint Transmissions

abstract: The first 5G (5th generation wireless systems) New Radio Release-15 was recently completed. However, the specification only considers the use of unicast technologies and the extension to point-to-multipoint (PTM) scenarios is not yet considered. To this end, we first present in this work a technical overview of the state-of-the-art LTE (Long Term Evolution) PTM technology, i.e., eMBMS (evolved Multimedia Broadcast Multicast Services), and investigate the physical layer performance via link-level simulations. Then based on the simulation analysis, we discuss potential improvements for the two current eMBMS solutions, i.e., MBSFN (MBMS over Single Frequency Networks) and SC-PTM (Single-Cell PTM). This work explicitly focus on equipping the current eMBMS solutions with 5G candidate techniques, e.g., multiple antennas and millimeter wave, and its potentials to meet the requirements of next generation PTM transmissions.


title: Collaborative SLAM based on Wifi Fingerprint Similarity and Motion Information

abstract: Simultaneous localization and mapping (SLAM) has been extensively researched in past years particularly with regard to range-based or visual-based sensors. Instead of deploying dedicated devices that use visual features, it is more pragmatic to exploit the radio features to achieve this task, due to their ubiquitous nature and the widespread deployment of Wi-Fi wireless network. This paper presents a novel approach for collaborative simultaneous localization and radio fingerprint mapping (C-SLAM-RF) in large unknown indoor environments. The proposed system uses received signal strengths (RSS) from Wi-Fi access points (AP) in the existing infrastructure and pedestrian dead reckoning (PDR) from a smart phone, without a prior knowledge about map or distribution of AP in the environment. We claim a loop closure based on the similarity of the two radio fingerprints. To further improve the performance, we incorporate the turning motion and assign a small uncertainty value to a loop closure if a matched turning is identified. The experiment was done in an area of 130 meters by 70 meters and the results show that our proposed system is capable of estimating the tracks of four users with an accuracy of 0.6 meters with Tango-based PDR and 4.76 meters with a step counter-based PDR.


title: Heuristic Approach for Jointly Optimizing FeICIC and UAV Locations in Multi-Tier LTE-Advanced Public Safety HetNet

abstract: UAV enabled communications and networking can enhance wireless connectivity and support emerging services. However, this would require system-level understanding to modify and extend the existing terrestrial network infrastructure. In this paper, we integrate UAVs both as user equipment and base stations into existing LTE-Advanced heterogeneous network (HetNet) and provide system-level insights of this three-tier LTE-Advanced air-ground HetNet (AG-HetNet). This AG-HetNet leverages cell range expansion (CRE), ICIC, 3D beamforming, and enhanced support for UAVs. Using system-level understanding and through brute-force technique and heuristics algorithms, we evaluate the performance of AG-HetNet in terms of fifth percentile spectral efficiency (5pSE) and coverage probability. We compare 5pSE and coverage probability, when aerial base-stations (UABS) are deployed on a fixed hexagonal grid and when their locations are optimized using genetic algorithm (GA) and elitist harmony search algorithm based on genetic algorithm (eHSGA). Our simulation results show the heuristic algorithms outperform the brute-force technique and achieve better peak values of coverage probability and 5pSE. Simulation results also show that trade-off exists between peak values and computation time when using heuristic algorithms. Furthermore, the three-tier hierarchical structuring of FeICIC provides considerably better 5pSE and coverage probability than eICIC.


title: QoS Provision for Controlling Energy Consumption in Ad-hoc Wireless Sensor Networks

abstract: Adhoc wireless sensor network is an architecture of connected nodes, each node has its main elements such as sensors, computation and communications capabilities. Adhoc WSNs restrained energy sources result in a shorter lifetime of the sensor network and inefficient topology. In this paper, a new approach for saving and energy controlling is introduced using quality of service, the main reason is to reduce the nodes energy through discovering the best optimum route that meets quality of service requirements; quality of service technique is used to find the optimum methodology for nodes packets transmission and energy consumption.


title: HPGe detector field calculation methods demonstrated with an educational program, GeFiCa

abstract: A review of tools and methods to calculate electrostatic potentials and fields inside high-purity germanium detectors in various configurations is given. The methods are illustrated concretely with a new educational program named GeFiCa - Germanium detector Field Calculator. Demonstrated in GeFiCa are generic numerical calculations based on the successive over-relaxation method as well as analytic ones whenever simplification is possible due to highly symmetric detector geometries. GeFiCa is written in C++, and provided as an extension to the CERN ROOT libraries widely used in the particle physics community. Calculation codes for individual detectors, provided as ROOT macros and python scripts, are distributed along with the GeFiCa core library, serving as both examples showing the usage of GeFiCa and starting points for customized calculations. They can be run without compilation in a ROOT interactive session or directly from a Linux shell. The numerical results are saved in a ROOT tree, making full use of the I/O optimization and plotting functionalities in ROOT. The speed and precision of the calculation are comparable to other commonly used packages, which qualifies GeFiCa as a scientific research tool. However, the main focus of GeFiCa is to clearly explain and demonstrate the analytic and numeric methods to solve Poisson's equation, practical coding considerations and visualization methods, with intensive documentation and example macros. It serves as a one-stop resource for people who want to understand the operating mechanism of such a package under the hood.


title: An Impedance-Modulated Code-Division Microwave SQUID Multiplexer

abstract: Large arrays of cryogenic detectors, including transition-edge sensors (TESs) or magnetic micro-calorimeters (MMCs), are needed for future experiments across a wide range of applications. Complexities in integration and cryogenic wiring have driven efforts to develop cryogenic readout technologies with large multiplexing factors while maintaining minimal readout noise. One such example is the microwave SQUID multiplexer ($\mu$mux), which couples an incoming TES or magnetic calorimeter signal to a unique GHz-frequency resonance that is modulated in frequency. Here, we present a hybrid scheme combining the microwave SQUID multiplexer with code division multiplexing: the impedance-modulated code-division multiplexer (Z-CDM), which may enable an order of magnitude increase in multiplexing factor particularly for low-bandwidth signal applications.


title: Extended hydrodynamics of the degenerate partially spin polarized fermions with the short-range interaction up to the third order by the interaction radius approximation

abstract: A possibility of the hydrodynamic description of ultracold fermions via the microscopic derivation of the model is described. Differently truncated hydrodynamic models are derived and compared. All models are based on the microscopic many-particle Schr\"{o}dinger equation. Minimal coupling model based on the continuity and Euler equations are considered. The extended hydrodynamic model including the independent dynamics of the momentum flux (the pressure evolution) is derived. Influence of the spin polarization is described. The short-range interaction is considered in the isotropic limit. The interaction is considered up to the third order by the interaction radius. Therefore, the single fluid model of spin-1/2 fermions and the two fluid model of spin-1/2 fermions are under consideration in this paper. Spectra of bulk collective excitations are derived and compared in terms of different models.


title: Direct Calculation of Self-Gravitational Force for Infinitesimally Thin Gaseous Disks Using Adaptive Mesh Refinement

abstract: Yen et al. (2012) advanced a direct approach for the calculation of self-gravitational force to second order accuracy based on uniform grid discretization. This method improves the accuracy of N-body calculation by using exact integration of kernel functions and employing the Fast Fourier Transform (FFT) to reduce complexity of computation to nearly linear. This direct approach is free of artificial boundary conditions, however, the applicability is limited by the uniform discretization of grids. We report here an advancement in the direct method with the implementation of adaptive mesh refinement (AMR) and maintaining second-order accuracy, which breaks the barrier set by uniform grid discretization. The adoption of graphic process units (GPUs) can significantly speed up the computation and make application of this method possible to astrophysical systems of gaseous disk galaxies and protoplanetary disks.


title: Partonic transport model application to heavy flavor

abstract: Heavy-flavor particles are excellent probes of the properties of the hot and dense nuclear medium created in the relativistic heavy-ion collisions. Heavy-flavor transport coefficients in the quark-gluon plasma (QGP) stage of the collisions are particularly interesting, as they contain important information on the strong interaction at finite temperatures. Studying the heavy-flavor evolution in a dynamically evolving medium requires a comprehensive multi-stage modeling approach of both the medium and the probes, with an accurate implementation of the physical ingredients to be tested. For this purpose, I have developed a new partonic transport model (Linear-Boltzmann-plus-Diffusion-Transport-Model) LIDO and applied it to heavy quark propagation inside a QGP. The model has an improved implementation of parton in-medium bremsstrahlung and a flexible treatment of the probe-medium interactions, combining both large angle scatterings and diffusion processes. The model is then coupled to a high-energy event-generator, a hydrodynamic medium evolution and a hadronic transport model. Finally, applying a Bayesian analysis, I extract the heavy quark transport coefficients from a model-to-data comparison. The results, with uncertainty quantification, are found to be consistent with earlier extraction of the light-quark transport coefficients at high momentum and with first-principle calculations of the heavy flavor diffusion constant at low momentum.


title: Explainable Deep Convolutional Candlestick Learner

abstract: Candlesticks are graphical representations of price movements for a given period. The traders can discovery the trend of the asset by looking at the candlestick patterns. Although deep convolutional neural networks have achieved great success for recognizing the candlestick patterns, their reasoning hides inside a black box. The traders cannot make sure what the model has learned. In this contribution, we provide a framework which is to explain the reasoning of the learned model determining the specific candlestick patterns of time series. Based on the local search adversarial attacks, we show that the learned model perceives the pattern of the candlesticks in a way similar to the human trader.


title: Supersymmetric Yang-Mills Theories: not quite the usual perspective

abstract: In this paper, we take up an old thread of development concerning the characterization of supersymmetric theories without any use of anticommuting variables that goes back to one of the authors' very early work [1]. Our special focus here will be on the formulation of supersymmetric Yang-Mills theories, extending previous results beyond $D=4$ dimensions. This perspective is likely to provide new insights into these theories, and in particular the maximally extended $N=4$ theory. As a new result we present a novel derivation of the admissible dimensions for interacting (pure) super-Yang-Mills theories to exist. This article is dedicated to the memory of Peter Freund, amongst many other things an early contributor to supersymmetry, and an author of one of the very first papers on superconformal gauge theories [2]. The final section contains some personal reminiscences of H.N.'s encounters with Peter Freund.


title: Incompressible viscous fluids in $\mathbb{R}^2$ and SPDEs on graphs, in presence of fast advection and non smooth noise

abstract: The asymptotic behavior of a class of stochastic reaction-diffusion-advection equations in the plane is studied. We show that as the divergence-free advection term becomes larger and larger, the solutions of such equations converge to the solution of a suitable stochastic PDE defined on the graph associated with the Hamiltonian. Firstly, we deal with the case that the stochastic perturbation is given by a singular spatially homogeneous Wiener process taking values in the space of Schwartz distributions. As in previous works, we assume here that the derivative of the period of the motion on the level sets of the Hamiltonian does not vanish. Then, in the second part, without assuming this condition on the derivative of the period, we study a weaker type of convergence for the solutions of a suitable class of linear SPDEs.


title: Fubini theorems for analytic Yeh-Feynman integrals associated with Gaussian processes with applications

abstract: In this paper we study an analytic Yeh--Feynman integral and an analytic Yeh--Fourier--Feynman transform associated with Gaussian processes. Fubini theorems involving the generalized analytic Yeh--Feynman integrals are established. The Fubini theorems investigated in this paper are to express the iterated generalized Yeh--Feynman integrals associated with Gaussian processes as a single generalized Yeh--Feynman integral. Using our Fubini theorems, we next examined fundamental relationships (with extended versions) between generalized Yeh--Fourier--Feynman transforms and convolution products (with respect to Gaussian processes) of functionals on Yeh--Wiener space.


title: Global Sensitivity Analysis in Load Modeling via Low-rank Tensor

abstract: Growing model complexities in load modeling have created high dimensionality in parameter estimations, and thereby substantially increasing associated computational costs. In this paper, a tensor-based method is proposed for identifying composite load modeling (CLM) parameters and for conducting a global sensitivity analysis. Tensor format and Fokker-Planck equations are used to estimate the power output response of CLM in the context of simultaneously varying parameters under their full parameter distribution ranges. The proposed tensor structured is shown as effective for tackling high-dimensional parameter estimation and for improving computational performances in load modeling through global sensitivity analysis.


title: DeepRecSys: A System for Optimizing End-To-End At-scale Neural Recommendation Inference

abstract: Neural personalized recommendation is the corner-stone of a wide collection of cloud services and products, constituting significant compute demand of the cloud infrastructure. Thus, improving the execution efficiency of neural recommendation directly translates into infrastructure capacity saving. In this paper, we devise a novel end-to-end modeling infrastructure, DeepRecInfra, that adopts an algorithm and system co-design methodology to custom-design systems for recommendation use cases. Leveraging the insights from the recommendation characterization, a new dynamic scheduler, DeepRecSched, is proposed to maximize latency-bounded throughput by taking into account characteristics of inference query size and arrival patterns, recommendation model architectures, and underlying hardware systems. By doing so, system throughput is doubled across the eight industry-representative recommendation models. Finally, design, deployment, and evaluation in at-scale production datacenter shows over 30% latency reduction across a wide variety of recommendation models running on hundreds of machines.


title: Lifted Hybrid Variational Inference

abstract: A variety of lifted inference algorithms, which exploit model symmetry to reduce computational cost, have been proposed to render inference tractable in probabilistic relational models. Most existing lifted inference algorithms operate only over discrete domains or continuous domains with restricted potential functions, e.g., Gaussian. We investigate two approximate lifted variational approaches that are applicable to hybrid domains and expressive enough to capture multi-modality. We demonstrate that the proposed variational methods are both scalable and can take advantage of approximate model symmetries, even in the presence of a large amount of continuous evidence. We demonstrate that our approach compares favorably against existing message-passing based approaches in a variety of settings. Finally, we present a sufficient condition for the Bethe approximation to yield a non-trivial estimate over the marginal polytope.


title: Stability of Sampling for CUR Decompositions

abstract: This article studies how to form CUR decompositions of low-rank matrices via primarily random sampling, though deterministic methods due to previous works are illustrated as well. The primary problem is to determine when a column submatrix of a rank $k$ matrix also has rank $k$. For random column sampling schemes, there is typically a tradeoff between the number of columns needed to be chosen and the complexity of determining the sampling probabilities. We discuss several sampling methods and their complexities as well as stability of the method under perturbations of both the probabilities and the underlying matrix. As an application, we give a high probability guarantee of the exact solution of the Subspace Clustering Problem via CUR decompositions when columns are sampled according to their Euclidean lengths.


title: stratamatch: Prognostic ScoreStratification using a Pilot Design

abstract: Optimal propensity score matching has emerged as one of the most ubiquitous approaches for causal inference studies on observational data; However, outstanding critiques of the statistical properties of propensity score matching have cast doubt on the statistical efficiency of this technique, and the poor scalability of optimal matching to large data sets makes this approach inconvenient if not infeasible for sample sizes that are increasingly commonplace in modern observational data. The stratamatch package provides implementation support and diagnostics for `stratified matching designs,' an approach which addresses both of these issues with optimal propensity score matching for large-sample observational studies. First, stratifying the data enables more computationally efficient matching of large data sets. Second, stratamatch implements a `pilot design' approach in order to stratify by a prognostic score, which may increase the precision of the effect estimate and increase power in sensitivity analyses of unmeasured confounding.


title: The topology, geometry and dynamics of free groups. Part I: Outer space, fold paths, and the Nielsen/Whitehead problems

abstract: The object of this expository work is to try to unveil the topological/geometric intuition behind the theory of free groups and their automorphism and outer automorphism groups. The method we follow is to focus on a series of problems in the study of free groups, and use the solutions of those problems to motivate topological/geometric tools. We do not aim to write down proofs which minimize the number of alphanumeric characters. We instead strive to write down proofs which maximize the development of broadly applicable geometric tools. In Part I we study problems solved by Nielsen and Whitehead in the 1920's and 1930's, but we approach these problems from a modern topological/geometric viewpoint, and we formulate their solutions so as to motivate modern tools, including marked graphs, the outer space of a free group, and fold paths in outer space.


title: Sagnac Interferometer and the Quantum Nature of Gravity

abstract: We use a quantum variant of the Sagnac interferometer to argue for the quantum nature of gravity as well as to formulate a quantum version of the equivalence principle. We first present an original derivation of the phase acquired in the conventional Sagnac matter-wave interferometer, within the Hamiltonian formalism. Then we modify the interferometer in two crucial respects. The interfering matter wave is interfered along two different distances from the centre and the interferometer is prepared in a superposition of two different angular velocities. We argue that if the radial and angular degrees of freedom of the matter wave become entangled through this experiment, then, via the equivalence principle, the gravitational field must be non-classical.


title: Distinguished curves and first integrals on Poincar\'e-Einstein and other conformally singular geometries

abstract: We treat the problem of defining, and characterising in a practical way, an appropriate class of distinguished curves for Poincar\'e-Einstein manifolds, and other conformally singular geometries. These "generalised geodesics" agree with geodesics away from the conformal singularity set and are shown to satisfy natural "boundary conditions" at points where they meet or cross the metric singularity set. We also characterise when they coincide with conformal circles. In the case of (Poincar\'e-)Einstein manifolds, we are able to provide a very general theory of first integrals for these distinguished curves. As well as the general procedure outlined, a specific example is given.


title: Robustly decorrelating errors with mixed quantum gates

abstract: Coherent errors in quantum operations are ubiquitous. Whether arising from spurious environmental couplings or errors in control fields, such errors can accumulate rapidly and degrade the performance of a quantum circuit significantly more than an average gate fidelity may indicate. As Hastings [1] and Campbell [2] have recently shown, by replacing the deterministic implementation of a quantum gate with a randomized ensemble of implementations, on can dramatically suppress coherent errors. Our work begins by reformulating the results of Hastings and Campbell as a quantum optimal control problem. We then discuss a family of convex programs designed to improve the performance, implementability, and robustness of the resulting mixed quantum gates. Finally, we implement these mixed quantum gates on a superconducting qubit and discuss randomized benchmarking results consistent with a marked reduction in the coherent error. [1] M. B. Hastings, Quantum Information & Computation 17, 488 (2017). [2] E. Campbell, Physical Review A 95, 042306 (2017).


title: Large Scale Structure Reconstruction with Short-Wavelength Modes

abstract: Large scale density modes are difficult to measure because they are sensitive to systematic observational errors in galaxy surveys, but we can study them indirectly by observing their impact on small scale perturbations. Cosmological perturbation theory predicts that second-order density inhomogeneities are a convolution of a short- and a long-wavelength mode. This arises physically because small scale structures grow at different rates depending on the large scale environment in which they reside. This induces an off-diagonal term in the two-point statistics in Fourier space that we use as the basis for a quadratic estimator for the large scale field. We demonstrate that this quadratic estimator works well on an N-body simulation of size (2.5 h^{-1} Gpc)^3. In particular, the quadratic estimator successfully reconstructs the long-wavelength modes using only small-scale information. This opens up novel opportunities to study structure on the largest observable scales.


title: Homology operations revisited

abstract: The mod p homology of E-infinity spaces is a classical topic in algebraic topology traditionally approached in terms of Dyer--Lashof operations. In this paper, we offer a new perspective on the subject by providing a detailed investigation of an alternative family of homology operations equivalent to, but distinct from, the Dyer--Lashof operations. Among other things, we will relate these operations to the Dyer--Lashof operations, describe the algebra generated by them, and use them to describe the homology of free E-infinity spaces. We will also investigate the relationship between the operations arising from the additive and multiplicative E-infinity structures on an E-infinity ring space. The operations have especially good properties in this context, allowing for a simple and conceptual formulation of "mixed Adem relations" describing how the operations arising from the two different E-infinity structures interact.


title: 1+1+2 covariant formulation of light propagation in spacetime

abstract: We present a covariant approach to the problem of light beam propagation in a spacetime. We develop our considerations within the framework of classical geometric optics in general relativity. Using the concept of a screen surface orthogonal to the observer velocity and to the bundle of geodesics, we introduce covariant four-dimensional definitions for Sachs and Jacobi optical fields and for the area distance. Then we give relationships between them and derive their propagation equations together with initial conditions for these equations. Ultimately, for practical use, we transform the resulting formulas into the redshift-dependent form.


title: If the Prospect of Some Occupations Are Stagnating With Technological Advancement? A Task Attribute Approach to Detect Employment Vulnerability

abstract: Two distinct trends can prove the existence of technological unemployment in the US. First, there are more open jobs than the number of unemployed persons looking for a job, and second, the shift of the Beveridge curve. There have been many attempts to find the cause of technological unemployment. However, all of these approaches fail when it comes to evaluating the impact of modern technologies on employment future. This study hypothesizes that rather than looking into skill requirement or routine non-routine discrimination of tasks, a holistic approach is required to predict which occupations are going to be vulnerable with the advent of this 4th industrial revolution, i.e., widespread application of AI, ML algorithms, and Robotics. Three critical attributes are considered: bottleneck, hazardous, and routine. Forty-five relevant attributes are chosen from the O*NET database that can define these three types of tasks. Performing Principal Axis Factor Analysis, and K-medoid clustering, the study discovers a list of 367 vulnerable occupations. The study further analyzes the last nine years of national employment data and finds that over the previous four years, the growth of vulnerable occupations is only half than that of non-vulnerable ones despite the long rally of economic expansion.


title: Deciphering exciton-generation processes in quantum-dot electroluminescence

abstract: Electroluminescence (EL) of colloidal nanocrystals promises a new generation of high-performance and solution-processable light-emitting diodes (LEDs). The operation of nanocrystal-based LEDs relies on the recombination of electrically-generated excitons. However, a fundamental question, i.e, how excitons are electrically generated in individual nanocrystals, remains unanswered. Here, we reveal a molecular mechanism of sequential electron-hole injection for the exciton generation in nanocrystal-based EL devices. To decipher the corresponding elementary processes, we develop electrically-pumped single-nanocrystal spectroscopy. While hole injection into neutral quantum dots (QDs) is generally-considered to be inefficient, we find that the intermediate negatively-charged state of QD triggers confinement-enhanced Coulomb interactions, which simultaneously accelerate hole injection and hinder excessive electron injection. In-situ/operando spectroscopy on state-of-the-art QD-LEDs demonstrate that exciton generation at the ensemble level is consistent with the charge-confinement-enabled sequential electron-hole injection mechanism revealed at the single-nanocrystal level. Our findings provide a universal mechanism for enhancing charge balance in nanocrystal-based EL devices.


title: Interaction of forced Orr-Sommerfeld and Squire modes in a low-order representation of turbulent channel flow

abstract: A resolvent-based reduced-order representation is used to capture time-averaged second-order statistics in turbulent channel flow. The recently-proposed decomposition of the resolvent operator into two distinct families related to the Orr-Sommerfeld and Squire operators [K. Rosenberg and B. J. McKeon, Efficient representation of exact coherent states of the Navier-Stokes equations using resolvent analysis, Fluid Dynamics Research 51, 011401 (2019)] results in dramatic improvement in the ability to match all components of the energy spectra and the $uv$ cospectrum. The success of the new representation relies on the ability of the Squire modes to compete with the vorticity generated by Orr-Sommerfeld modes, which is demonstrated by decomposing the statistics into contributions from each family. It is then shown that this competition can be used to infer a phase relationship between the two sets of modes. Additionally, the relative Reynolds number scalings for the two families of resolvent weights are derived for the universal classes of resolvent modes presented by Moarref et al. [R. Moarref, A. S. Sharma, J. A. Tropp, and B. J. McKeon, Model-based scaling of the streamwise energy density in high-Reynolds-number turbulent channels, Journal of Fluid Mechanics 734, 275 (2013)]. These developments can be viewed as a starting point for further modeling efforts to quantify nonlinear interactions in wall-bounded turbulence.


title: Least squares binary quantization of neural networks

abstract: Quantizing weights and activations of deep neural networks results in significant improvement in inference efficiency at the cost of lower accuracy. A source of the accuracy gap between full precision and quantized models is the quantization error. In this work, we focus on the binary quantization, in which values are mapped to -1 and 1. We provide a unified framework to analyze different scaling strategies. Inspired by the pareto-optimality of 2-bits versus 1-bit quantization, we introduce a novel 2-bits quantization with provably least squares error. Our quantization algorithms can be implemented efficiently on the hardware using bitwise operations. We present proofs to show that our proposed methods are optimal, and also provide empirical error analysis. We conduct experiments on the ImageNet dataset and show a reduced accuracy gap when using the proposed least squares quantization algorithms.


title: The Hodge ring of varieties in positive characteristic

abstract: Let $k$ be a field of positive characteristic. We prove that the only linear relations between the Hodge numbers $h^{i,j}(X) = \dim H^j(X,\Omega_X^i)$ that hold for every smooth proper variety $X$ over $k$ are the ones given by Serre duality. We also show that the only linear combinations of Hodge numbers that are birational invariants of $X$ are given by the span of the $h^{i,0}(X)$ and the $h^{0,j}(X)$ (and their duals $h^{i,n}(X)$ and $h^{n,j}(X)$). The corresponding statements for compact K\"ahler manifolds were proven by Kotschick and Schreieder.


title: Numerical evolution of the interior geometry of charged black holes

abstract: Previously, we developed a late time approximation scheme to study the interior geometry of black holes. In the present paper we test this scheme with numerical relativity simulations. In particular, we present numerical relativity simulations of the interior geometry of charged spherically symmetric two-sided black holes with a spacelike singularity at $r = 0$. Our numerics are in excellent agreement with the late time approximation. We also demonstrate that the geometry near $r = 0$ is a scalarized Kasner geometry and compute the associated Kasner exponents.


title: Gallai Ramsey number for double stars

abstract: Given a graph $G$ and a positive integer $k$, the \emph{Gallai-Ramsey number} is defined to be the minimum number of vertices $n$ such that any $k$-edge coloring of $K_n$ contains either a rainbow (all different colored) copy of $G$ or a monochromatic copy of $G$. In this paper, we obtain general upper and lower bounds on the Gallai-Ramsey numbers for double stars $S(n,m)$, where $S(n,m)$ is the graph obtained from the union of two stars $K_{1,n}$ and $K_{1,m}$ by adding an edge between their centers. We also provide the sharp result in some cases.


title: The Dynamic-Mode Decomposition and Optimal Prediction

abstract: The Dynamic-Mode Decomposition (DMD) is a well established data-driven method of finding temporally evolving linear-mode decompositions of nonlinear time series. Traditionally, this method presumes that all relevant dimensions are sampled through measurement. To address dynamical systems in which the data may be incomplete or represent only partial observation of a more complex system, we extend the DMD algorithm by including a Mori-Zwanzig Decomposition to derive memory kernels that capture the averaged dynamics of the unresolved variables as projected onto the resolved dimensions. From this, we then derive what we call the Memory-Dependent Dynamic Mode Decomposition (MDDMD). Through numerical examples, the MDDMD method is shown to produce reasonable approximations of the ensemble-averaged dynamics of the full system given a single time series measurement of the resolved variables.


title: The 3.4 {\mu}m absorption of the Titan's stratosphere: contribution of ethane, propane, butane and complex hydrogenated organics

abstract: The complex organic chemistry harbored by the atmosphere of Titan has been investigated in depth by Cassini observations. Among them, a series of solar occultations performed by the VIMS instrument throughout the 13 years of Cassini revealed a strong absorption centered at 3.4 $\mu$m. Several molecules present in Titan's atmosphere create spectral features in that wavelength region, but their individual contributions are difficult to disentangle. In this work, we quantify the contribution of the various molecular species to the 3.4 $\mu$m band using a radiative transfer model. Ethane and propane are a significant component of the band but they are not enough to fit the shape perfectly, then we need something else. Polycyclic Aromatic Hydrocarbons (PAHs) and more complex polyaromatic hydrocarbons like Hydrogenated Amorphous Carbons (HACs) are the most plausible candidates because they are rich in C-H bonds. PAHs signature have already been detected above ~900 km, and they are recognized as aerosols particles precursors. High similarities between individual spectra impede abundances determinations.


title: On Computation and Generalization of Generative Adversarial Imitation Learning

abstract: Generative Adversarial Imitation Learning (GAIL) is a powerful and practical approach for learning sequential decision-making policies. Different from Reinforcement Learning (RL), GAIL takes advantage of demonstration data by experts (e.g., human), and learns both the policy and reward function of the unknown environment. Despite the significant empirical progresses, the theory behind GAIL is still largely unknown. The major difficulty comes from the underlying temporal dependency of the demonstration data and the minimax computational formulation of GAIL without convex-concave structure. To bridge such a gap between theory and practice, this paper investigates the theoretical properties of GAIL. Specifically, we show: (1) For GAIL with general reward parameterization, the generalization can be guaranteed as long as the class of the reward functions is properly controlled; (2) For GAIL, where the reward is parameterized as a reproducing kernel function, GAIL can be efficiently solved by stochastic first order optimization algorithms, which attain sublinear convergence to a stationary solution. To the best of our knowledge, these are the first results on statistical and computational guarantees of imitation learning with reward/policy function approximation. Numerical experiments are provided to support our analysis.


title: Central Limit Theorems on Compact Metric Spaces

abstract: We produce a series of Central Limit Theorems (CLTs) associated to compact metric measure spaces $(K,d,\eta)$, with $\eta$ a reasonable probability measure. For the first CLT, we can ignore $\eta$ by isometrically embedding $K$ into ${\mathcal C}(K)$, the space of continuous functions on $K$ with the sup norm, and then applying known CLTs for sample means on Banach spaces (Theorem 3.1). However, the sample mean makes no sense back on $K$, so using $\eta$ we develop a CLT for the sample Fr\'echet mean (Corollary 4.1). This involves working on the closed convex hull of the embedded image of $K$. To work in the easier Hilbert space setting of $L^2(K,\eta)$, we have to modify the metric $d$ to a related metric $d_\eta$. We obtain an $L^2$-CLT for both the sample mean and the sample Fr\'echet mean (Theorem 5.1), and we relate the Fr\'echet sample and population means on the closed convex hull to the Fr\'echet means on the image of $K$. Since the $L^2$ and $L^\infty$ norms play important roles, in Section 6 we develop a metric-measure criterion relating $d$ and $\eta$ under which all $L^p$ norms are equivalent.


title: On the construction of non-Hermitian Hamiltonians with all-real spectra through supersymmetric algorithms

abstract: The energy spectra of two different quantum systems are paired through supersymmetric algorithms. One of the systems is Hermitian and the other is characterized by a complex-valued potential, both of them with only real eigenvalues in their spectrum. The superpotential that links these systems is complex-valued, parameterized by the solutions of the Ermakov equation, and may be expressed either in nonlinear form or as the logarithmic derivative of a properly chosen complex-valued function. The non-Hermitian systems can be constructed to be either parity-time-symmetric or non-parity-time-symmetric.


title: Enhancing Dynamic-Mode Decomposition via Multi-Scale Analysis

abstract: Through the use of wavelet based Besov norms, we compute nontrivial multiscale nonlinear features of a given data set so as to enhance the standard Dynamic-Mode Decomposition algorithm. Thus we are able to build sophisticated observables which enhance algorithm performance without placing undue computational burdens on the user.


title: Moir\'e magnets

abstract: We introduce a general framework to study moir\'e structures of two-dimensional Van der Waals magnets using continuum field theory. The formalism eliminates quasiperiodicity and allows a full understanding of magnetic structures and their excitations. In particular, we analyze in detail twisted bilayers of N\'eel antiferromagnets on the honeycomb lattice. A rich phase diagram with non-collinear twisted phases is obtained, and spin waves are further calculated. Direct extensions to zig-zag antiferromagnets and ferromagnets are also presented. We anticipate the results and formalism presented to lead to a broad range of applications to both fundamental research and experiments.


title: Approximation and Convergence of Large Atomic Congestion Games

abstract: We consider the question of whether, and in what sense, Wardrop equilibria provide a good approximation for Nash equilibria in atomic unsplittable congestion games with a large number of small players. We examine two different definitions of small players. In the first setting, we consider games where each player's weight is small. We prove that when the number of players goes to infinity and their weights to zero, the random flows in all (mixed) Nash equilibria for the finite games converge in distribution to the set of Wardrop equilibria of the corresponding nonatomic limit game. In the second setting, we consider an increasing number of players with a unit weight that participate in the game with a decreasingly small probability. In this case, the Nash equilibrium flows converge in total variation towards Poisson random variables whose expected values are Wardrop equilibria of a different nonatomic game with suitably-defined costs. The latter can be viewed as symmetric equilibria in a Poisson game in the sense of Myerson, establishing a plausible connection between the Wardrop model for routing games and the stochastic fluctuations observed in real traffic. In both settings we provide explicit approximation bounds, and we study the convergence of the price of anarchy. Beyond the case of congestion games, we prove a general result on the convergence of large games with random players towards Poisson games.


title: Self-guided Approximate Linear Programs

abstract: Approximate linear programs (ALPs) are well-known models based on value function approximations (VFAs) to obtain policies and lower bounds on the optimal policy cost of discounted-cost Markov decision processes (MDPs). Formulating an ALP requires (i) basis functions, the linear combination of which defines the VFA, and (ii) a state-relevance distribution, which determines the relative importance of different states in the ALP objective for the purpose of minimizing VFA error. Both these choices are typically heuristic: basis function selection relies on domain knowledge while the state-relevance distribution is specified using the frequency of states visited by a heuristic policy. We propose a self-guided sequence of ALPs that embeds random basis functions obtained via inexpensive sampling and uses the known VFA from the previous iteration to guide VFA computation in the current iteration. Self-guided ALPs mitigate the need for domain knowledge during basis function selection as well as the impact of the initial choice of the state-relevance distribution, thus significantly reducing the ALP implementation burden. We establish high probability error bounds on the VFAs from this sequence and show that a worst-case measure of policy performance is improved. We find that these favorable implementation and theoretical properties translate to encouraging numerical results on perishable inventory control and options pricing applications, where self-guided ALP policies improve upon policies from problem-specific methods. More broadly, our research takes a meaningful step toward application-agnostic policies and bounds for MDPs.


title: Neural Data Server: A Large-Scale Search Engine for Transfer Learning Data

abstract: Transfer learning has proven to be a successful technique to train deep learning models in the domains where little training data is available. The dominant approach is to pretrain a model on a large generic dataset such as ImageNet and finetune its weights on the target domain. However, in the new era of an ever-increasing number of massive datasets, selecting the relevant data for pretraining is a critical issue. We introduce Neural Data Server (NDS), a large-scale search engine for finding the most useful transfer learning data to the target domain. NDS consists of a dataserver which indexes several large popular image datasets, and aims to recommend data to a client, an end-user with a target application with its own small labeled dataset. The dataserver represents large datasets with a much more compact mixture-of-experts model, and employs it to perform data search in a series of dataserver-client transactions at a low computational cost. We show the effectiveness of NDS in various transfer learning scenarios, demonstrating state-of-the-art performance on several target datasets and tasks such as image classification, object detection and instance segmentation. Neural Data Server is available as a web-service at http://aidemos.cs.toronto.edu/nds/.


title: Complex Variability of Kepler AGN Revealed by Recurrence Analysis

abstract: The advent of new time domain surveys and the imminent increase in astronomical data expose the shortcomings in traditional time series analysis (such as power spectra analysis) in characterising the abundantly varied, complex and stochastic light curves of Active Galactic Nuclei (AGN). Recent applications of novel methods from non-linear dynamics have shown promise in characterising higher modes of variability and time-scales in AGN. Recurrence analysis in particular can provide complementary information about characteristic time-scales revealed by other methods, as well as probe the nature of the underlying physics in these objects. Recurrence analysis was developed to study the recurrences of dynamical trajectories in phase space, which can be constructed from one-dimensional time series such as light curves. We apply the methods of recurrence analysis to two optical light curves of Kepler-monitored AGN. We confirm the detection and period of an optical quasi-periodic oscillation in one AGN, and confirm multiple other time-scales recovered from other methods ranging from 5 days to 60 days in both objects. We detect regions in the light curves that deviate from regularity, provide evidence of determinism and non-linearity in the mechanisms underlying one light curve (KIC 9650712), and determine a linear stochastic process recovers the dominant variability in the other light curve (Zwicky 229--015). We discuss possible underlying processes driving the dynamics of the light curves and their diverse classes of variability.


title: Learning landmark guided embeddings for animal re-identification

abstract: Re-identification of individual animals in images can be ambiguous due to subtle variations in body markings between different individuals and no constraints on the poses of animals in the wild. Person re-identification is a similar task and it has been approached with a deep convolutional neural network (CNN) that learns discriminative embeddings for images of people. However, learning discriminative features for an individual animal is more challenging than for a person's appearance due to the relatively small size of ecological datasets compared to labelled datasets of person's identities. We propose to improve embedding learning by exploiting body landmarks information explicitly. Body landmarks are provided to the input of a CNN as confidence heatmaps that can be obtained from a separate body landmark predictor. The model is encouraged to use heatmaps by learning an auxiliary task of reconstructing input heatmaps. Body landmarks guide a feature extraction network to learn the representation of a distinctive pattern and its position on the body. We evaluate the proposed method on a large synthetic dataset and a small real dataset. Our method outperforms the same model without body landmarks input by 26% and 18% on the synthetic and the real datasets respectively. The method is robust to noise in input coordinates and can tolerate an error in coordinates up to 10% of the image size.