title: Current-induced magnetization hysteresis defines atom trapping in a superconducting atomchip

abstract: The physics of superconducting films, and especially the role of remnant magnetization has a defining influence on the magnetic fields used to hold and manipulate atoms on superconducting atomchips. We magnetically trap ultracold ^{87}Rb atoms on a 200{\mu}m wide and 500nm thick cryogenically cooled niobium Z wire structure. By measuring the distance of the atomcloud to the trapping wire for different transport currents and bias fields, we probe the trapping characteristics of the niobium superconducting structure. At distances closer than the trapping wire width, we observe a different behaviour than that of normal conducting wire traps. Furthermore, we measure a stable magnetic trap at zero transport current. These observations point to the presence of a remnant magnetization in our niobium film which is induced by a transport current. This current-induced magnetization defines the trap close to the chip surface. Our measurements agree very well with an analytic prediction based on the critical state model (CSM). Our results provide a new tool to control atom trapping on superconducting atomchips by designing the current distribution through its current history.


title: The Linear BESS Model at the LHC

abstract: In this work we consider the Linear BESS model at the LHC. This model can be seen as an adequate benchmark for exploring the phenomenological consequences of a composite Higgs sector since its particle content is the one we would expect in a realistic low energy description of modern (Technicolor inspired) dynamical electroweak symmetry breaking scenarios. Additionally, the model exhibits the property of decoupling, producing a good ultraviolet behavior. We focus on the limits on the masses of the new heavy vector particles imposed by direct resonance searches, recent measurements of the decay of the Higgs boson into two photons and the electroweak precision tests. We found that the model is capable to accommodate the existing experimental constrains provided that the spin-1 resonances are heavier than 3.4 TeV.


title: Influence of outer-layer finite-size effects on the rupture kinetics of a thin polymer film embedded in an immiscible matrix

abstract: In capillary-driven fluid dynamics, simple departures from equilibrium offer the chance to quantitatively model the resulting relaxations. These dynamics in turn provide insight on both practical and fundamental aspects of thin-film hydrodynamics. In this work, we describe a model trilayer dewetting experiment elucidating the effect of solid, no-slip confining boundaries on the bursting of a liquid film in a viscous environment. This experiment was inspired by an industrial polymer processing technique, multilayer coextrusion, in which thousands of alternating layers are stacked atop one another. When pushed to the nanoscale limit, the individual layers are found to break up on time scales shorter than the processing time. To gain insight on this dynamic problem, we here directly observe the growth rate of holes in the middle layer of the trilayer films described above, wherein the distance between the inner film and solid boundary can be orders of magnitude larger than its thickness. In otherwise identical experimental conditions, thinner films break up faster than thicker ones. This observation is found to agree with a scaling model that balances capillary driving power and viscous dissipation with, crucially, a no-slip boundary condition at the solid substrate/viscous environment boundary. In particular, even for the thinnest middle-layers, no finite-size effect is needed to explain the data. The dynamics of hole growth is captured by a single master curve over four orders of magnitude in the dimensionless hole radius and time, and is found to agree well with predictions including analytic expressions for the dissipation.


title: Seglearn: A Python Package for Learning Sequences and Time Series

abstract: Seglearn is an open-source python package for machine learning time series or sequences using a sliding window segmentation approach. The implementation provides a flexible pipeline for tackling classification, regression, and forecasting problems with multivariate sequence and contextual data. This package is compatible with scikit-learn and is listed under scikit-learn Related Projects. The package depends on numpy, scipy, and scikit-learn. Seglearn is distributed under the BSD 3-Clause License. Documentation includes a detailed API description, user guide, and examples. Unit tests provide a high degree of code coverage.


title: Implications of Higgs Discovery for the Strong CP Problem and Unification

abstract: A $Z_2$ symmetry that extends the weak interaction, $SU(2)_L \rightarrow SU(2)_L \times SU(2)'$, and the Higgs sector, $H(2) \rightarrow H(2,1) + H'(1,2)$, yields a Standard Model quartic coupling that vanishes at scale $v' = <H'>~\gg~<H>$. Near $v'$, theories either have a "prime" sector, or possess "Left-Right" (LR) symmetry with $SU(2)' = SU(2)_R$. If the $Z_2$ symmetry incorporates spacetime parity, these theories can solve the strong CP problem. The LR theories have all quark and lepton masses arising from operators of dimension 5 or more, requiring Froggatt-Nielsen structures. Two-loop contributions to $\bar{\theta}$ are estimated and typically lead to a neutron electric dipole moment of order $10^{-27}$e cm that can be observed in future experiments. Minimal models, with gauge group $SU(3) \times SU(2)_L \times SU(2)_L \times U(1)_{B-L}$, have precise gauge coupling unification for $v' = 10^{10\pm1}$ GeV, successfully correlating gauge unification with the observed Higgs mass of $125$ GeV. With $SU(3) \times U(1)_{B-L}$ embedded in $SU(4)$, the central value of the unification scale is reduced from $10^{16-17}$ GeV to below $10^{16}$ GeV, improving the likelihood of proton decay discovery. Unified theories based on $SO(10) \times CP$ are constructed that have $H+H'$ in a ${\bf 16}$ or ${\bf 144}$ and generate higher-dimensional flavor operators, while maintaining perturbative gauge couplings.


title: Stochastic PDE Limit of the Six Vertex Model

abstract: We study the stochastic six vertex model and prove that under weak asymmetry scaling (i.e., when the parameter $\Delta\to 1^+$ so as to zoom into the ferroelectric/disordered phase critical point) its height function fluctuations converge to the solution to the KPZ equation. We also prove that the one-dimensional family of stochastic Gibbs states for the symmetric six vertex model converge under the same scaling to the stationary solution to the stochastic Burgers equation. Our proofs rely upon the Markov (self) duality of our model. The starting point is an exact microscopic Hopf-Cole transform for the stochastic six vertex model which follows from the model's known one-particle Markov self-duality. Given this transform, the crucial step is to establish self-averaging for specific quadratic function of the transformed height function. We use the model's two-particle self-duality to produce explicit expressions (as Bethe ansatz contour integrals) for conditional expectations from which we extract time-decorrelation and hence self-averaging in time. The crux of our Markov duality method is that the entire convergence result reduces to precise estimates on the one-particle and two-particle transition probabilities. Previous to our work, Markov dualities had only been used to prove convergence of particle systems to linear Gaussian SPDEs (e.g. the stochastic heat equation with additive noise).

categories: math.PR cond-mat.stat-mech math-ph math.MP


title: A Markov Chain Monte Carlo Approach to Cost Matrix Generation for Scheduling Performance Evaluation

abstract: In high performance computing, scheduling of tasks and allocation to machines is very critical especially when we are dealing with heterogeneous execution costs. Simulations can be performed with a large variety of environments and application models. However, this technique is sensitive to bias when it relies on random instances with an uncontrolled distribution. We use methods from the literature to provide formal guarantee on the distribution of the instance. In particular, it is desirable to ensure a uniform distribution among the instances with a given task and machine heterogeneity. In this article, we propose a method that generates instances (cost matrices) with a known distribution where tasks are scheduled on machines with heterogeneous execution costs.


title: Statistical diagonalization of a random biased Hamiltonian: the case of the eigenvectors

abstract: We present a non perturbative calculation technique providing the mixed moments of the overlaps between the eigenvectors of two large quantum Hamiltonians: $\hat{H}_0$ and $\hat{H}_0+\hat{W}$, where $\hat{H}_0$ is deterministic and $\hat{W}$ is random. We apply this method to recover the second order moments or Local Density Of States in the case of an arbitrary fixed $\hat{H}_0$ and a Gaussian $\hat{W}$. Then we calculate the fourth order moments of the overlaps in the same setting. Such quantities are crucial for understanding the local dynamics of a large composite quantum system. In this case, $\hat{H}_0$ is the sum of the Hamiltonians of the system subparts and $\hat{W}$ is an interaction term. We test our predictions with numerical simulations.


title: An Early Modeling Approach to Digital Electronics

abstract: An Early modeling approach of transistors characterized by simplicity and accuracy in representing intrinsic non-linearities is applied to the characterization of propagation delay and level transition switching properties of NPN and PNP small signal transistors. Eight types of devices were considered, each represented by 5 samples taken from the same lot, totaling 20 NPN and 20 PNP transistors. Four switching time measurements were experimentally obtained, and the transistors also had their Early parameters $V_a$ (the Early voltage) and $s$ (a proportionality constant such that $R_o = 1/tan(s I_B)$ accurately estimated by using an experimental-numeric procedure that involves Hough transform accumulation in order to identify the crossing of the base current ($I_B$) indexed characteristic isolines, yielding the respective $V_a$. The timing measurements exhibited strong positive Pearson correlations when taken pairwise. When these measurements were compared individually to the respective Early parameters, no significant Pearson correlation was obtained. However, a strong relationship was observed between the product of the two Early parameters and the ratio between the fall and rise time. A Pearson correlation coefficient of 0.78 was observed between these variables in the case of NPN devices. This suggests that transistors with larger average current gain tend to have more similar rise and fall times. The different relationship observed for PNP devices (Pearson 0.41) suggests some intrinsic difference in the way the Early parameters influence the rise and fall times of small signal transistors.


title: Photonic waveguide mode to free-space Gaussian beam extreme mode converter

abstract: Integration of photonic chips with atomic, micromechanical, chemical and biological systems can advance science and open many possibilities in chip-scale devices and technology. Compact photonic structures for direct coupling of light between high-index single-mode waveguides and arbitrary free-space modes spanning hundreds of waves in cross-section would eliminate bulky optical components and enable integration of photonics into many new applications requiring wide beams, structured light and centimeter-scale propagation distances with low diffraction-limited losses. Conventional fiber-coupling approaches do not scale well for accurate, low-loss coupling across the extremely large mode scale mismatch ($\approx10^6$ times in modal area). Here we present an extreme mode converter that can transform the photonic waveguide mode to the diffraction-limited, free-space Gaussian beam, with a beam waist of about $160~\mu$m. Using two identical converters, we demonstrate a grating-to-grating coupling that couples the radiating beam back to the chip through a mirror reflection in free-space. Operating at 780~nm for integration with chip-scale atomic vapor cell cavities, our design can be adapted for visible, telecommunication or other wavelengths. Furthermore, other types of beams can be implemented by using the 2-stage expansion approach presented in this paper.


title: Dimensionally regularized Tsallis' Statistical Mechanics and two-body Newton's gravitation

abstract: Typical Tsallis' statistical mechanics' quantifiers like the partition function and the mean energy exhibit poles. We are speaking of the partition function ${\cal Z}$ and the mean energy $<{\cal U}>$. The poles appear for distinctive values of Tsallis' characteristic real parameter $q$, at a numerable set of rational numbers of the $q-$line. These poles are dealt with dimensional regularization resources. The physical effects of these poles on the specific heats are studied here for the two-body classical gravitation potential.


title: Non-extensive statistical mechanics of a self-gravitating gas

abstract: The statistical mechanics of a cloud of particles interacting via their gravitational potentials is an old problem which encounters some issues when the traditional Boltzmann-Gibbs statistics is applied. In this article, we consider the generalized statistics of Tsallis and analyze the statistical and thermodynamical implications for a self-gravitating gas, obtaining analytical and convergent expressions for the equation of state and specific heat in the canonical as well as microcanonical ensembles. Although our results are comparable in both ensembles, it turns out that only in the canonical case the thermodynamic quantities depend explicitly on the non-extensivity parameter, indicating that the question of ensemble equivalence for Tsallis statistics must be further reviewed.


title: Spectral Statistics of Non-Hermitian Random Matrix Ensembles

abstract: Recently Burkhardt et. al. introduced the $k$-checkerboard random matrix ensembles, which have a split limiting behavior of the eigenvalues (in the limit all but $k$ of the eigenvalues are on the order of $\sqrt{N}$ and converge to semi-circular behavior, with the remaining $k$ of size $N$ and converging to hollow Gaussian ensembles). We generalize their work to consider non-Hermitian ensembles with complex eigenvalues; instead of a blip new behavior is seen, ranging from multiple satellites to annular rings. These results are based on moment method techniques adapted to the complex plane as well as analysis of singular values, and we further isolate the singular value joint density formula for the Complex Symmetric Gaussian Ensemble.


title: Modeling cure fraction with frailty term in latent risk: a Bayesian approach

abstract: In this paper, we propose a flexible cure rate model with frailty term in latent risk, which is obtained by incorporating a frailty term in risk function of latent competing causes. The number of competing causes of the event of interest follows negative binomial distribution and the frailty variable follows power variance function distribution, in which includes other frailty models such as gamma, positive stable and inverse Gaussian frailty models as special cases. The proposed model takes into account the presence of covariates and right-censored survival data suitable for populations with a cure rate. Besides, it allows quantifying the degree of unobserved heterogeneity induced by unobservable risk factors, in which is important to explain the survival time. Once the posterior distribution has not close form, Markov chain Monte Carlo simulations are considered for estimation procedure. We performed several simulation studies and the practical relevance of the proposed model is demonstrated in a real data set.


title: Analysis of Student Satisfaction Toward Quality of Service Facility

abstract: The development of higher education is very rapid rise to the tight competition both public universities and private colleges. XYZ University realized to win the competition, required continuous quality improvement, including the quality of existing service facilities. Amenities quality services is believed to support the success of the learning activities and improve user satisfaction. This study aims to determine the extent to which the quality of the services effect on user satisfaction. The research method used is survey-based questionnaire that measure perception and expectation. The results showed a gap between perception and expectations of the respondents have a negative value for each item. This means XYZ service facility at the university is not currently meet the expectations of society members. Three service facility that has the lowest index is based on the perception of respondents is a laboratory (2.56), computer and multimedia (2.63) as well as wifi network (2.99). The magnitude of the correlation between satisfaction with the quality of service facilities is 0.725 which means a strong and positive relationship. The influence of the quality of service facilities to the satisfaction of the students is 0.525 meaning that the variable quality of the services facility can explain 52.5% of the variable satisfaction. The study provided recommendations for improvements to enhance the quality of services facility at the XYZ university facilities.


title: Non-Markovian quantum-classical ratchet for ultrafast long-range electron-hole separation in condensed phases

abstract: In organic photovoltaic systems, a photogenerated molecular exciton in the donor domain dissociates into a hole and an electron at the donor/acceptor heterojunction, and subsequently separate into free charge carriers that can be extracted as photocurrents. The recombination of the once-separated electron and hole is a major loss mechanism in photovoltaic systems, which controls their performance. Hence, efficient photovoltaic systems need built-in ratchet mechanisms, namely, ultrafast charge separation and retarded charge recombination. In order to obtain insight into the internal working of the experimentally observed ultrafast long-range charge separation and protection against charge recombination, we theoretically investigate a potential ratchet mechanism arising from the combination of quantum delocalization and its destruction by performing numerically accurate quantum-dynamics calculations on a model system. It is demonstrated that the non-Markovian effect originating from the slow polaron formation strongly suppresses the electron transfer reaction back to the interfacial charge-transfer state stabilized at the donor/accepter interface and that it plays a critical role in maintaining the long-range electron--hole separation.


title: Markovian robots: minimal navigation strategies for active particles

abstract: We explore minimal navigation strategies for active particles in complex, dynamical, external fields, introducing a class of autonomous, self-propelled particles which we call Markovian robots (MR). These machines are equipped with a navigation control system (NCS) that triggers random changes in the direction of self-propulsion of the robots. The internal state of the NCS is described by a Boolean variable that adopts two values. The temporal dynamics of this Boolean variable is dictated by a closed Markov chain -- ensuring the absence of fixed points in the dynamics -- with transition rates that may depend exclusively on the instantaneous, local value of the external field. Importantly, the NCS does not store past measurements of this value in continuous, internal variables. We show that, despite the strong constraints, it is possible to conceive closed Markov chain motifs that lead to nontrivial motility behaviors of the MR in one, two and three dimensions. By analytically reducing the complexity of the NCS dynamics, we obtain an effective description of the long-time motility behavior of the MR that allows us to identify the minimum requirements in the design of NCS motifs and transition rates to perform complex navigation tasks such as adaptive gradient following, detection of minima or maxima, or selection of a desired value in a dynamical, external field. We put these ideas in practice by assembling a robot that operates by the proposed minimalistic NCS to evaluate the robustness of MR, providing a proof-of-concept that is possible to navigate through complex information landscapes with such a simple NCS whose internal state can be stored in one bit. These ideas may prove useful for the engineering of miniaturized robots.


title: Acoustic wave polarization and energy flow in periodic beam lattice materials

abstract: The free propagation of acoustic plane waves through cellular periodic materials is generally accompanied by a flow of mechanical energy across the adjacent cells. The paper focuses on the energy transport related to dispersive waves propagating through nondissipative microstructured materials. The generic microstructure of the periodic cell is described by a beam lattice model, suitably reduced to the minimal space of dynamic degrees-of-freedom. The linear eigenproblem governing the wave propagation is stated and the complete eigensolution is considered to study both the real-valued dispersion functions and the complex-valued waveforms of the propagating elastic waves. First, a complete family of nondimensional quantities (polarization factors) is proposed to quantify the linear polarization or quasi-polarization, according to a proper energetic criterion. Second, a vector variable related to the periodic cell is introduced to assess the directional flux of mechanical energy, in analogy to the Umov-Poynting vector related to the material point in solid mechanics. The physical-mathematical relation between the energy flux and the velocity of the energy transport is recognized. The formal equivalence between the energy and the group velocity is pointed out, according to the mechanical assumptions. Finally, all the theoretical developments are successfully applied to the prototypical beam lattice material characterized by a periodic tetrachiral microstructure.


title: Comments on frequency dependent ac conductivity in polymeric materials at low frequency regime

abstract: The AC conductivity response in a broad frequency range of disordered materials is of great interest not only for technological applications, but also from a theoretical point of view. The Jonscher power exponent value, and its temperature dependence, is a very important parameter in dielectric data analysis as well as the physical interpretation of conduction mechanisms in disordered materials. In some cases the power exponent of AC conductivity has been reported to be greater than 1 at the low frequency regime. This fact seems to contradict the universal dynamic response. The present work focuses on the analysis of dielectric spectroscopy measurements in polymeric materials, below 100 MHz. The apparent power exponent n gets values in the range (0,1) and is directly related to the characteristics of mobile charges at shorter time scales, in the case of the occurrence of DC conduction and the slowest polarization mechanism that is due to the charge motions within sort length scales, in log(epsilon'')-log(frequenvy) plot. The emergence of apparent n values in the range [1,2], for a relatively narrow frequency range, may be attributed to an additional molecular dipolar relaxation contribution at higher frequencies, in log(epsilon'')-log(frequency) plot. The appearance of apparent n values in the range (1,2], can be assigned to the existence of a well defined minimum between DC conductivity contribution and a molecular dipolar dispersion or between two well separated dielectric loss mechanisms, in log(epsilon'')-log(frequency) plots, above the crossover frequency. In these latter cases, the apparent power exponent n is merely related to the Havriliak-Negami equation shape parameters of the higher frequencies molecular dipolar relaxations.


title: Task dependent Deep LDA pruning of neural networks

abstract: With deep learning's success, a limited number of popular deep nets have been widely adopted for various vision tasks. However, this usually results in unnecessarily high complexities and possibly many features of low task utility. In this paper, we address this problem by introducing a task-dependent deep pruning framework based on Fisher's Linear Discriminant Analysis (LDA). The approach can be applied to convolutional, fully-connected, and module-based deep network structures, in all cases leveraging the high decorrelation of neuron motifs found in the pre-decision space and cross-layer deconv dependency. Moreover, we examine our approach's potential in network architecture search for specific tasks and analyze the influence of our pruning on model robustness to noises and adversarial attacks. Experimental results on datasets of generic objects (ImageNet, CIFAR100) as well as domain specific tasks (Adience, and LFWA) illustrate our framework's superior performance over state-of-the-art pruning approaches and fixed compact nets (e.g. SqueezeNet, MobileNet). The proposed method successfully maintains comparable accuracies even after discarding most parameters (98%-99% for VGG16, up to 82% for the already compact InceptionNet) and with significant FLOP reductions (83% for VGG16, up to 64% for InceptionNet). Through pruning, we can also derive smaller, but more accurate and more robust models suitable for the task.


title: Information Geometry of the Gaussian Space

abstract: We discuss the Pistone-Sempi exponential manifold on the finite-dimensional Gaussian space. We consider the role of the entropy, the continuity of translations, Poincar\'e-type inequalities, the generalized differentiability of probability densities of the Gaussian space.


title: On The Correct Thermo-dynamic Potential for Electro-static Dielectric Energy

abstract: Various types of equilibrium processes involve electric fields. In some cases, the electrical energy appears to be negative (e.g. if the voltage is fixed by an external source). This paper explains how to derive the correct thermo-dynamic potential for electro-static phenomena, whether the voltage is fixed, or the charge is fixed, or some combination is fixed. In particular, we explain, in complete detail, why fixing the voltage introduces "a minus sign" in the electrical energy. Two explanations are given. The first explanation is based on a lumped-parameter argument (i.e. a lumped-capacitor model). The second explanation uses a distributed parameter model (i.e. a partial differential equation (PDE) model) of a dielectric medium, in this case, we allow for non-linearity and external polarization effects. Connections with Legendre (duality) transforms are also discussed.


title: Robust Blind Deconvolution via Mirror Descent

abstract: We revisit the Blind Deconvolution problem with a focus on understanding its robustness and convergence properties. Provable robustness to noise and other perturbations is receiving recent interest in vision, from obtaining immunity to adversarial attacks to assessing and describing failure modes of algorithms in mission critical applications. Further, many blind deconvolution methods based on deep architectures internally make use of or optimize the basic formulation, so a clearer understanding of how this sub-module behaves, when it can be solved, and what noise injection it can tolerate is a first order requirement. We derive new insights into the theoretical underpinnings of blind deconvolution. The algorithm that emerges has nice convergence guarantees and is provably robust in a sense we formalize in the paper. Interestingly, these technical results play out very well in practice, where on standard datasets our algorithm yields results competitive with or superior to the state of the art. Keywords: blind deconvolution, robust continuous optimization


title: Extended depth-of-field in holographic image reconstruction using deep learning based auto-focusing and phase-recovery

abstract: Holography encodes the three dimensional (3D) information of a sample in the form of an intensity-only recording. However, to decode the original sample image from its hologram(s), auto-focusing and phase-recovery are needed, which are in general cumbersome and time-consuming to digitally perform. Here we demonstrate a convolutional neural network (CNN) based approach that simultaneously performs auto-focusing and phase-recovery to significantly extend the depth-of-field (DOF) in holographic image reconstruction. For this, a CNN is trained by using pairs of randomly de-focused back-propagated holograms and their corresponding in-focus phase-recovered images. After this training phase, the CNN takes a single back-propagated hologram of a 3D sample as input to rapidly achieve phase-recovery and reconstruct an in focus image of the sample over a significantly extended DOF. This deep learning based DOF extension method is non-iterative, and significantly improves the algorithm time-complexity of holographic image reconstruction from O(nm) to O(1), where n refers to the number of individual object points or particles within the sample volume, and m represents the focusing search space within which each object point or particle needs to be individually focused. These results highlight some of the unique opportunities created by data-enabled statistical image reconstruction methods powered by machine learning, and we believe that the presented approach can be broadly applicable to computationally extend the DOF of other imaging modalities.


title: Calculation of the Mean Strain of Smooth Non-uniform Strain Fields Using Conventional FBG Sensors

abstract: In the past few decades, fibre Bragg grating (FBG) sensors have gained a lot of attention in the field of distributed point strain measurement. One of the most interesting properties of these sensors is the presumed linear relationship between the strain and the peak wavelength shift of the FBG reflected spectra. However, subjecting sensors to a non-uniform stress field will in general result in a strain estimation error when using this linear relationship. In this paper we propose a new strain estimation algorithm that accurately estimates the mean strain value in the case of smooth non-uniform strain distributions. To do so, we first introduce an approximation of the classical transfer matrix model, which we will refer to as the approximated transfer matrix model (ATMM). This model facilitates the analysis of FBG reflected spectra under arbitrary strain distributions, particularly by providing a closed-form approximation of the side-lobes of the reflected spectra. Based on this new formulation, we derive a maximum likelihood estimator of the mean strain value. The algorithm is validated using both computer simulations and experimental FBG measurements. Compared to state-of-the-art methods, which typically introduce errors of tens of microstrains, the proposed method is able to compensate for this error. In the typical examples that were analysed in this study, mean strain errors of around 60 microstrains were compensated.


title: On locally repeated values of arithmetic functions over $\mathbb F_q[T]$

abstract: The frequency of occurrence of "locally repeated" values of arithmetic functions is a common theme in analytic number theory, for instance in the Erd\H{o}s-Mirsky problem on coincidences of the divisor function at consecutive integers, the analogous problem for the Euler totient function, and the quantitative conjectures of Erd\H{o}s, Pomerance and Sark\H{o}zy and of Graham, Holt and Pomerance on the frequency of occurrences. In this paper we introduce the corresponding problems in the setting of polynomials over a finite field, and completely solve them in the large finite field limit.


title: Efficient Search of QC-LDPC Codes with Girths 6 and 8 and Free of Elementary Trapping Sets with Small Size

abstract: One of the phenomena that influences significantly the performance of low-density parity-check codes is known as trapping sets. An $(a,b)$ elementary trapping set, or simply an ETS where $a$ is the size and $b$ is the number of degree-one check nodes and $\frac{b}{a}<1$, causes high decoding failure rate and exert a strong influence on the error floor. In this paper, we provide sufficient conditions for exponent matrices to have fully connected $(3,n)$-regular QC-LDPC codes with girths 6 and 8 whose Tanner graphs are free of small ETSs. Applying sufficient conditions on the exponent matrix to remove some 8-cycles results in removing all 4-cycles, 6-cycles as well as some small elementary trapping sets. For each girth we obtain a lower bound on the lifting degree and present exponent matrices with column weight three whose corresponding Tanner graph is free of certain ETSs.


title: Giant negative electrostriction and dielectric tunability in a van der Waals layered ferroelectric

abstract: The interest in ferroelectric van der Waals crystals arises from the potential to realize ultrathin ferroic systems owing to the reduced surface energy of these materials and the layered structure that allows for exfoliation. Here, we quantitatively unravel giant negative electrostriction of van der Waals layered copper indium thiophosphate (CIPS), which exhibits an electrostrictive coefficient Q33 as high as -3.2 m4/C2 and a resulting bulk piezoelectric coefficient d33 up to -85 pm/V. As a result, the electromechanical response of CIPS is comparable in magnitude to established perovskite ferroelectrics despite possessing a much smaller spontaneous polarization of only a few uC/cm2. In the paraelectric state, readily accessible owing to low transition temperatures, CIPS exhibits large dielectric tunability, similar to widely-used barium strontium titanate, and as a result both giant and continuously tunable electromechanical response. The persistence of electrostrictive and tunable responses in the paraelectric state indicates that even few layer films or nanoparticles will sustain significant electromechanical functionality, offsetting the inevitable suppression of ferroelectric properties in the nanoscale limit. These findings can likely be extended to other ferroelectric transition metal thiophosphates and (quasi-) two-dimensional materials and might facilitate the quest towards novel ultrathin functional devices incorporating electromechanical response.


title: Theory of pore-driven and end-pulled polymer translocation dynamics through a nanopore: An overview

abstract: We review recent progress on the theory of dynamics of polymer translocation through a nanopore based on the iso-flux tension propagation (IFTP) theory. We investigate both pore-driven translocation of flexible and a semi-flexible polymers, and the end-pulled case of flexible chains by means of the IFTP theory and extensive molecular dynamics (MD) simulations. The validity of the IFTP theory can be quantified by the waiting time distributions of the monomers which reveal the details of the dynamics of the translocation process. The IFTP theory allows a parameter-free description of the translocation process and can be used to derive exact analytic scaling forms in the appropriate limits, including the influence due to the pore friction that appears as a finite-size correction to asymptotic scaling. We show that in the case of pore-driven semi-flexible and end-pulled polymer chains the IFTP theory must be augmented with an explicit {\it trans} side friction term for a quantitative description of the translocation process.


title: Constraints on mass, spin and magnetic field of microquasar H~1743-322 from observations of QPOs

abstract: The study of quasi-periodic oscillations (QPOs) of X-ray flux observed in many microquasars can provide a powerful tool for testing of the phenomena occurring in strong gravity regime. QPOs phenomena can be well related to the oscillations of charged particles in accretion disks orbiting Kerr black holes immersed in external large-scale magnetic fields. In the present paper we study the model of magnetic relativistic precession and provide estimations of the mass and spin of the central object of the microquasar H~1743-322 which is a candidate for a black hole. Moreover, we discuss the possible values of external magnetic field and study its influence on the motion of charged particles around rotating black hole.


title: Circular repetition thresholds on some small alphabets: Last cases of Gorbunova's conjecture

abstract: A word is called $\beta$-free if it has no factors of exponent greater than or equal to $\beta$. The repetition threshold $\mathrm{RT}(k)$ is the infimum of the set of all $\beta$ such that there are arbitrarily long $k$-ary $\beta$-free words (or equivalently, there are $k$-ary $\beta$-free words of every sufficiently large length, or even every length). These three equivalent definitions of the repetition threshold give rise to three natural definitions of a repetition threshold for circular words. The infimum of the set of all $\beta$ such that - there are arbitrarily long $k$-ary $\beta$-free circular words is called the weak circular repetition threshold, denoted $\mathrm{CRT}_{\mathrm{W}}(k)$; - there are $k$-ary $\beta$-free circular words of every sufficiently large length is called the intermediate circular repetition threshold, denoted $\mathrm{CRT}_{\mathrm{I}}(k)$; - there are $k$-ary $\beta$-free circular words of every length is called the strong circular repetition threshold, denoted $\mathrm{CRT}_{\mathrm{S}}(k)$. We prove that $\mathrm{CRT}_{\mathrm{S}}(4)=\tfrac{3}{2}$ and $\mathrm{CRT}_{\mathrm{S}}(5)=\tfrac{4}{3}$, confirming a conjecture of Gorbunova and providing the last unknown values of the strong circular repetition threshold. We also prove that $\mathrm{CRT}_{\mathrm{I}}(3)=\mathrm{CRT}_{\mathrm{W}}(3)=\mathrm{RT}(3)=\tfrac{7}{4}$.


title: Mapping The Neutrino Floor For Dark Matter-Electron Direct Detection Experiments

abstract: We study the discovery reach of future Dark Matter (DM) Direct Detection experiments using DM-electron scattering in the presence of the solar neutrino background. At these low energies traditional methods for nuclear and electronic recoil discrimination fail, implying that the neutrino-{\it nucleus} scattering background can be sizable. We calculate discovery limits based on ionization values of signal and background, and quantify the dependence on the ionization model. Moreover, we explore how the dependence of the DM cross section discovery limits vary with exposure, electronic/nuclear recoil discrimination, DM form factors, and DM astrophysical uncertainties.


title: The Pontrjagin Dual of 4-Dimensional Spin Bordism

abstract: The goal of this paper is to study the Pontrjagin dual of (reduced) 4-dimensional Spin bordism. That is to say, we consider the functor from the category of topological spaces to the category of compact abelian groups that associates to each space X the compact group of homomorphisms from the reduced 4-dimensional Spin bordism of X to the circle. In a previous paper, we studied the analogous problem for 3-dimensional Spin bordism. Our work was motivated by some questions from physics. The physicists are primarily interested in the case when X is the classifying space of a finite group, but our arguments are valid for general X. We describe the dual group, G(X), as equivalence classes of triples of cochains (w,p,a) on X, triples satisfying certain relations with a product. We also describe the pairing between such triples and a closed 4-dimensional Spin manifold mapping to X, the pairing that produces the identification of G(X) with the Pontrjagin dual of the reduced 4-dimensional Spin bordism of X.


title: Large anomalous Nernst coefficient in an oxide Skyrmion crystal Chern insulator

abstract: A sizable transverse thermoelectric coefficient N , large to the extent that it potentially serves applications, is predicted to arise, by means of first-principles calculations, in a Skyrmion crystal assumed on EuO monolayer where carrier electrons are introduced upon a quantum anomalous Hall insulating phase of Chern number C = 2. This encourages future experiments to pursue such an effect.


title: Spinning boson stars and hairy black holes with non-minimal coupling

abstract: We obtain spinning boson star solutions and hairy black holes with synchronised hair in the Einstein-Klein-Gordon model, wherein the scalar field is massive, complex and with a non-minimal coupling to the Ricci scalar. The existence of these hairy black holes in this model provides yet another manifestation of the universality of the synchronisation mechanism to endow spinning black holes with hair. We study the variation of the physical properties of the boson stars and hairy black holes with the coupling parameter between the scalar field and the curvature, showing that they are, qualitatively, identical to those in the minimally coupled case. By discussing the conformal transformation to the Einstein frame, we argue that the solutions herein provide new rotating boson star and hairy black hole solutions in the minimally coupled theory, with a particular potential, and that no spherically symmetric hairy black hole solutions exist in the non-minimally coupled theory, under a condition of conformal regularity.


title: Generic Zero-Cost Reuse for Dependent Types

abstract: Dependently typed languages are well known for having a problem with code reuse. Traditional non-indexed algebraic datatypes (e.g. lists) appear alongside a plethora of indexed variations (e.g. vectors). Functions are often rewritten for both non-indexed and indexed versions of essentially the same datatype, which is a source of code duplication. We work in a Curry-style dependent type theory, where the same untyped term may be classified as both the non-indexed and indexed versions of a datatype. Many solutions have been proposed for the problem of dependently typed reuse, but we exploit Curry-style type theory in our solution to not only reuse data and programs, but do so at zero-cost (without a runtime penalty). Our work is an exercise in dependently typed generic programming, and internalizes the process of zero-cost reuse as the identity function in a Curry-style theory.


title: A Theorem for Secrecy in Tagged Protocols Using the Theory of Witness-Functions

abstract: In this paper, we enunciate the theorem of secrecy in tagged protocols using the theory of witness-functions and we run a formal analysis on a new tagged version of the Needham-Schroeder public-key protocol using this theorem. We discuss the significance of tagging in securing cryptographic protocols as well.


title: Connectivity-Preserving Coordination Control of Multi-Agent Systems with Time-Varying Delays

abstract: This paper presents a distributed position synchronization strategy that also preserves the initial communication links for single-integrator multi-agent systems with time-varying delays. The strategy employs a coordinating proportional control derived from a specific type of potential energy, augmented with damping injected through a dynamic filter. The injected damping maintains all agents within the communication distances of their neighbours, and asymptotically stabilizes the multi-agent system, in the presence of time delays. Regarding the closed-loop single-integrator multi-agent system as a double-integrator system suggests an extension of the proposed strategy to connectivity-preserving coordination of Euler-Lagrange networks with time-varying delays. Lyapunov stability analysis and simulation results validate the two designs.


title: Learning the Localization Function: Machine Learning Approach to Fingerprinting Localization

abstract: Considered as a data-driven approach, Fingerprinting Localization Solutions (FPSs) enjoy huge popularity due to their good performance and minimal environment information requirement. This papers addresses applications of artificial intelligence to solve two problems in Received Signal Strength Indicator (RSSI) based FPS, first the cumbersome training database construction and second the extrapolation of fingerprinting algorithm for similar buildings with slight environmental changes. After a concise overview of deep learning design techniques, two main techniques widely used in deep learning are exploited for the above mentioned issues namely data augmentation and transfer learning. We train a multi-layer neural network that learns the mapping from the observations to the locations. A data augmentation method is proposed to increase the training database size based on the structure of RSSI measurements and hence reducing effectively the amount of training data. Then it is shown experimentally how a model trained for a particular building can be transferred to a similar one by fine tuning with significantly smaller training numbers. The paper implicitly discusses the new guidelines to consider about deep learning designs when they are employed in a new application context.


title: Network and Panel Quantile Effects Via Distribution Regression

abstract: This paper provides a method to construct simultaneous confidence bands for quantile functions and quantile effects in nonlinear network and panel models with unobserved two-way effects, strictly exogenous covariates, and possibly discrete outcome variables. The method is based upon projection of simultaneous confidence bands for distribution functions constructed from fixed effects distribution regression estimators. These fixed effects estimators are debiased to deal with the incidental parameter problem. Under asymptotic sequences where both dimensions of the data set grow at the same rate, the confidence bands for the quantile functions and effects have correct joint coverage in large samples. An empirical application to gravity models of trade illustrates the applicability of the methods to network data.


title: Fast Bayesian inference in large Gaussian graphical models

abstract: Despite major methodological developments, Bayesian inference for Gaussian graphical models remains challenging in high dimension due to the tremendous size of the model space. This article proposes a method to infer the marginal and conditional independence structures between variables by multiple testing of hypotheses. Specifically, we introduce closed-form Bayes factors under the Gaussian conjugate model to evaluate the null hypotheses of marginal and conditional independence between variables. Their computation for all pairs of variables is shown to be extremely efficient, thereby allowing us to address large problems with thousands of nodes. Moreover, we derive exact tail probabilities from the null distributions of the Bayes factors. These allow the use of any multiplicity correction procedure to control error rates for incorrect edge inclusion. We demonstrate the proposed approach to graphical model selection on various simulated examples as well as on a large gene expression data set from The Cancer Genome Atlas.


title: Quantum glass of interacting bosons with off-diagonal disorder

abstract: We study disordered interacting bosons described by the Bose-Hubbard model with Gaussian-distributed random tunneling amplitudes. It is shown that the off-diagonal disorder induces a spin-glass-like ground state, characterized by randomly frozen quantum-mechanical U(1) phases of bosons. To access criticality, we employ the "$n$-replica trick", as in the spin-glass theory, and the Trotter-Suzuki method for decomposition of the statistical density operator, along with numerical calculations. The interplay between disorder, quantum and thermal fluctuations leads to phase diagrams exhibiting a glassy state of bosons, which are studied as a function of model parameters. The considered system may be relevant for quantum simulators of optical-lattice bosons, where the randomness can be introduced in a controlled way. The latter is supported by a proposition of experimental realization of the system in question.


title: On the Parameterized Computation of Minimum Volume Outer Ellipsoid of Minkowski Sum of Ellipsoids

abstract: We consider the problem of computing certain parameterized minimum volume outer ellipsoidal (MVOE) approximation of the Minkowski sum of a finite number of ellipsoids. We clarify connections among several parameterizations available in the literature, obtain novel analysis results regarding the conditions of optimality, and based on the same, propose two new algorithms for computing the parameterized MVOE. Numerical results reveal faster runtime for the proposed algorithms than the state-of-the-art semidefinite programming approach of computing the same.


title: Comparative study of structural and electronic properties of GaSe and InSe polytypes

abstract: Equilibrium crystal structures, electron band dispersions and band gap values of layered GaSe and InSe semiconductors, each being represented by four polytypes, are studied via first-principles calculations within the density functional theory (DFT). A number of practical algorithms to take into account dispersion interactions are tested, from empirical Grimme corrections to many-body dispersion schemes. Due to the utmost technical accuracy achieved in the calculations, nearly degenerate energy-volume curves of different polytypes are resolved, and the conclusions concerning the relative stability of competing polytypes drawn. The predictions are done as for how the equilibrium between different polytypes will be shifted under the effect of hydrostatic pressure. The band structures are inspected under the angle of identifying features specific for different polytypes, and with respect to modifications of the band dispersions brought about by the use of modified Becke-Johnson (mBJ) scheme for the exchange-correlation (XC) potential. As another way to improve the predictions of band gaps values, hybrid functional calculations according to the HSE06 scheme are performed for the band structures, and the relation with the mBJ results discussed. Both methods nicely agree with experimental results and with state-of-the-art GW calculations. Some discrepancies are identified in cases of close competition between the direct and indirect gap (e.g., in GaSe); moreover, the accurate placement of bands revealing relatively localized states is slightly different according to mBJ and HSE06 schemes.


title: Globally Stable Output Feedback Synchronization of Teleoperation with Time-Varying Delays

abstract: This paper presents a globally stable teleoperation control strategy for systems with time-varying delays that eliminates the need for velocity measurements through novel augmented Immersion and Invariance velocity observers. The new observers simplify a recent constructive Immersion and Invariance velocity observer to achieve globally convergent velocity estimation with only $n+2$ states, where $n$ is the number of degrees of freedom of the master and slave robots. They introduce dynamic scaling factors to accelerate the speed of convergence of the velocity estimates and, thus, to limit the energy generated by the velocity estimation errors and to guarantee sufficient estimate-based damping injection to dissipate the energy generated by the time-varying delays. The paper shows that Proportional plus damping control with the simplified and augmented Immersion and Invariance-based velocity observers can synchronize the free master and slave motions in the presence of time-varying delays without using velocity measurements. Numerical results illustrate the estimation performance of the new observers and the stability of a simulated two degrees-of-freedom nonlinear teleoperation system with time-varying delays under the proposed output feedback Proportional plus damping control.


title: An Economic Bubble Model and Its First Passage Time

abstract: We introduce a new diffusion process Xt to describe asset prices within an economic bubble cycle. The main feature of the process, which differs from existing models, is the drift term where a mean-reversion is taken based on an exponential decay of the scaled price. Our study shows the scaling factor on Xt is crucial for modelling economic bubbles as it mitigates the dependence structure between the price and parameters in the model. We prove both the process and its first passage time are well-defined. An efficient calibration scheme, together with the probability density function for the process are given. Moreover, by employing the perturbation technique, we deduce the closed-form density for the downward first passage time, which therefore can be used in estimating the burst time of an economic bubble. The object of this study is to understand the asset price dynamics when a financial bubble is believed to form, and correspondingly provide estimates to the bubble crash time. Calibration examples on the US dot-com bubble and the 2007 Chinese stock market crash verify the effectiveness of the model itself. The example on BitCoin prediction confirms that we can provide meaningful estimate on the downward probability for asset prices.


title: Entropy-based closure for probabilistic learning on manifolds

abstract: In a recent paper, the authors proposed a general methodology for probabilistic learning on manifolds. The method was used to generate numerical samples that are statistically consistent with an existing dataset construed as a realization from a non-Gaussian random vector. The manifold structure is learned using diffusion manifolds and the statistical sample generation is accomplished using a projected Ito stochastic differential equation. This probabilistic learning approach has been extended to polynomial chaos representation of databases on manifolds and to probabilistic nonconvex constrained optimization with a fixed budget of function evaluations. The methodology introduces an isotropic-diffusion kernel with hyperparameter {\epsilon}. Currently, {\epsilon} is more or less arbitrarily chosen. In this paper, we propose a selection criterion for identifying an optimal value of {\epsilon}, based on a maximum entropy argument. The result is a comprehensive, closed, probabilistic model for characterizing data sets with hidden constraints. This entropy argument ensures that out of all possible models, this is the one that is the most uncertain beyond any specified constraints, which is selected. Applications are presented for several databases.


title: Large Thermal Hall Effect in $\alpha$-RuCl$_3$: Evidence for Heat Transport by Kitaev-Heisenberg Paramagnons

abstract: The honeycomb Kitaev model in a magnetic field is a source of a topological quantum spin liquid with Majorana fermions and gauge flux excitations as fractional quasiparticles. We present experimental results for the thermal Hall effect of the material $\alpha$-RuCl$_{3}$ which recently emerged as a prime candidate for realizing such physics. At temperatures above long-range magnetic ordering $T\gtrsim T_N\approx8$ K, we observe with an applied magnetic field $B$ perpendicular to the honeycomb layers a sizeable positive transversal heat conductivity $\kappa_{xy}$ which increases linearly with $B$. Upon raising the temperature, $\kappa_{xy}(T)$ increases strongly, exhibits a broad maximum at around 30 K, and eventually becomes negligible at $T\gtrsim 125$ K. Remarkably, the longitudinal heat conductivity $\kappa_{xx}(T)$ exhibits a sizeable positive thermal magnetoresistance effect. Thus, our findings provide clear-cut evidence for longitudinal and transverse magnetic heat transport and underpin the unconventional nature of the quasiparticles in the paramagnetic phase of $\alpha$-RuCl$_{3}$.


title: Detection of the closest Jovian exoplanet in the Epsilon Indi triple system

abstract: We confirm the trend in the radial velocity data for Epsilon Indi A suggesting a long-period planetary companion and find significant curvature is present, sufficient to quantify Epsilon Indi Ab as a cold Jupiter with a minimum mass of $2.71_{-0.44}^{+2.19}~M_{\rm Jup}$ on a nearly circular orbit with a semi-major axis of $12.82_{-0.71}^{+4.18}$ au and an orbital period of $52.62_{-4.12}^{+27.70}$ yr. We also identify other significant signals in the radial velocity data. We investigate a variety of spectral diagnostics and interpret these signals as arising from activity-induced radial velocity variations. In particular, the 2500 and 278 d signals are caused by magnetic cycles. While a planetary signal might be present in the 17.8 d signal, the origin of 17.8 and 11 d signals are most easily interpreted as arising in the rotation of the star with a period of about 35 d. We find that traditional activity indicators have a variety of sensitivities. In particular, the sodium lines and CaHK index are sensitive to all activity-induced signals. The line bisector measurement is sensitive to stellar rotation signal while H$\alpha$ is sensitive to the secondary magnetic cycle. In general, because of their different sensitivities these activity indicators introduce extra noise if included in the noise model whereas differential RVs provide a robust proxy to remove wavelength-dependent noise efficiently. Based on these analyses, we propose an activity diagnostics procedure for the detection of low amplitude signals in high precision radial velocity data. Thus the Epsilon Indi system comprises of at least Epsilon Indi A, Ab as well as a long period brown dwarf binary Ba and Bb; so it provides a benchmark case for our understanding of the formation of gas giants and brown dwarfs.


title: Boltzmann approach to spin-orbit-induced transport in effective quantum theories

abstract: In model studies of the spin/anomalous Hall effect, effective Hamiltonians often serve as the starting point. However, a complete effective quantum theory contains not only the effective Hamiltonian but also the relation linking the physical observables to the canonical ones. We construct the semiclassical Boltzmann (SB) transport framework in the weak disorder-potential regime directly in the level of the effective quantum theory, and confirm this construction by formulating a generalized Kohn-Luttinger density matrix transport theory also in this level. The link and difference between the present SB theory and previous phenomenological Boltzmann, quantum kinetic and usual Kubo-Streda theories are clarified. We also present the slightly generalized Kubo-Streda formula in the level of the effective quantum theory. In this level, it is the generalized Kubo-Streda formula rather than the usual one that leads to the same physical interpretations as the present SB theory. In the application to a Rashba 2D effective model, a nonzero spin Hall effect important in the case of strong Rashba coupling but neglected in previous theories is found.


title: Comparing Fixed and Adaptive Computation Time for Recurrent Neural Networks

abstract: Adaptive Computation Time for Recurrent Neural Networks (ACT) is one of the most promising architectures for variable computation. ACT adapts to the input sequence by being able to look at each sample more than once, and learn how many times it should do it. In this paper, we compare ACT to Repeat-RNN, a novel architecture based on repeating each sample a fixed number of times. We found surprising results, where Repeat-RNN performs as good as ACT in the selected tasks. Source code in TensorFlow and PyTorch is publicly available at https://imatge-upc.github.io/danifojo-2018-repeatrnn/


title: Optimal price management in retail energy markets: an impulse control problem with asymptotic estimates

abstract: We consider a retailer who buys energy in the wholesale market and resells it to final consumers. The retailer has to decide when to intervene to change the price he asks to his customers, in order to maximize his income. We model the problem as an infinite-horizon stochastic impulse control problem. We characterize an optimal price strategy and provide analytical existence results for the equations involved. We then investigate the dependence on the intervention cost. In particular, we prove that the measure of the continuation region is asymptotic to the fourth root of the cost. Finally, we provide some numerical results and consider a suitable extension of the model.


title: A Staggered Explicit-Implicit Finite Element Formulation for Electroactive Polymers

abstract: Electroactive polymers such as dielectric elastomers (DEs) have attracted significant attention in recent years. Computational techniques to solve the coupled electromechanical system of equations for this class of materials have universally centered around fully coupled monolithic formulations, which while generating good accuracy requires significant computational expense. However, this has significantly hindered the ability to solve large scale, fully three-dimensional problems involving complex deformations and electromechanical instabilities of DEs. In this work, we provide theoretical basis for the effectiveness and accuracy of staggered explicit-implicit finite element formulations for this class of electromechanically coupled materials, and elicit the simplicity of the resulting staggered formulation. We demonstrate the stability and accuracy of the staggered approach by solving complex electromechanically coupled problems involving electroactive polymers, where we focus on problems involving electromechanical instabilities such as creasing, wrinkling, and bursting drops. In all examples, essentially identical results to the fully monolithic solution are obtained, showing the accuracy of the staggered approach at a significantly reduced computational cost.


title: Neutron Energy Dependence of Delayed Neutron Yields and its Assessments

abstract: Incident neutron energy dependence of delayed neutron yields of uranium and plutonium isotopes is investigated. A summation calculation of decay and fission yield data is employed, and the energy dependence of the latter part is considered in a phenomenological way. Our calculation systematically reproduces the energy dependence of delayed neutron yields by introducing an energy dependence of the most probable charge and the odd-even effect. The calculated fission yields are assessed by comparison with JENDL/FPY-2011, delayed neutron activities, and decay heats. Although the fission yields in this work are optimized to delayed neutron yields, the calculated decay heats are in good agreement with the experimental data. Comparison of the fission yields calculated in this work and JENDL/FPY-2011 gave an important insight for the evaluation of the next JENDL nuclear data.


title: Financial Contagion in a Generalized Stochastic Block Model

abstract: One of the most defining features of the global financial network is its inherent complex and intertwined structure. From the perspective of systemic risk it is important to understand the influence of this network structure on default contagion. Using sparse random graphs to model the financial network, asymptotic methods turned out powerful to analytically describe the contagion process and to make statements about resilience. So far, however, they have been limited to so-called {\em rank one} models in which informally the only network parameter is the degree sequence (see (Amini et. al. 2016) and (Detering et. al. 2019) for example) and the contagion process can be described by a one dimensional fix-point equation. These networks fail to account for a pronounced block structure such as core/periphery or a network composed of different connected blocks for different countries. We present a much more general model here, where we distinguish vertices (institutions) of different types and let edge probabilities and exposures depend on the types of both, the receiving and the sending vertex plus additional parameters. Our main result allows to compute explicitly the systemic damage caused by some initial local shock event, and we derive a complete characterisation of resilient respectively non-resilient financial systems. This is the first instance that default contagion is rigorously studied in a model outside the class of rank one models and several technical challenges arise. Moreover, in contrast to previous work, in which networks could be classified as resilient or non resilient, independent of the distribution of the shock, information about the shock becomes important in our model and a more refined resilience condition arises. Among other applications of our theory we derive resilience conditions for the global network based on subnetwork conditions only.


title: Mislearning from Censored Data: The Gambler's Fallacy and Other Correlational Mistakes in Optimal-Stopping Problems

abstract: I study endogenous learning dynamics for people who misperceive intertemporal correlations in random sequences. Biased agents face an optimal-stopping problem. They are uncertain about the underlying distribution and learn its parameters from predecessors. Agents stop when early draws are "good enough," so predecessors' experiences contain negative streaks but not positive streaks. When agents wrongly expect systematic reversals (the "gambler's fallacy"), they understate the likelihood of consecutive below-average draws, converge to over-pessimistic beliefs about the distribution's mean, and stop too early. Agents uncertain about the distribution's variance overestimate it to an extent that depends on predecessors' stopping thresholds. I also analyze how other misperceptions of intertemporal correlation interact with endogenous data censoring.


title: Emotional Attachment Framework for People-Oriented Software

abstract: In organizational and commercial settings, people often have clear roles and workflows against which functional and non-functional requirements can be extracted. However, in more social settings, such as platforms for enhancing social interaction, successful applications are driven more by using emotional engagement than functionality, and the drivers of user engagement are difficult to identify. A key challenge is to understand people's emotional goals so that they can be incorporated into the design. This paper proposes a novel framework called the Emotional Attachment Framework, which is based on existing models and theories of emotional attachment. Its aim is to facilitate the process of getting a deeper insight into emotional goals in software engineering. To demonstrate the framework in use, emotional goals are elicited for a software application that aims to provide help for homeless people. To measure the effectiveness and efficiency of the proposed technique in this study, a series of evaluations are undertaken: a semi-controlled experiment, a comparison analysis, and domain expert and end-user evaluation. The results indicate that the Emotional Attachment Framework has the potential to give better insight during analysis of emotional goals.


title: Coexistence of Triple Nodal Points, Nodal Links, and Unusual Flat Bands in intermetallic ${\cal A}$Pd$_3$ (${\cal A}$=Pb, Sn)

abstract: We investigate the electronic structure and several properties, and topological character, of the cubic time-reversal invariant intermetallic compounds PbPd$_3$ and SnPd$_3$ using density functional theory based methods. These compounds have a dispersionless band along the $\Gamma-X$ line, forming the top of the Pd $4d$ bands and lying within a few meV of the Fermi level $E_F$. Effects of the flat band on transport and optical properties have been inspected by varying the doping concentration treated with the virtual crystal approximation for substitution on the Pb site. In the absence of spin-orbit coupling (SOC), we find triple nodal points and three-dimensional nodal loops, which are known to lead to surface bands and drumhead states, respectively, which we discuss for PbPd$_3$. SOC removes degeneracy in most of the zone, providing a topological index $Z_2$=1 on the $k_z=0$ plane that indicates a topological character on that plane. The isovalent and isostructural compound SnPd$_3$ shows only minor differences in its electronic structures, so it is expected to display similar electronic, transport, and topological properties.


title: Clear and Cloudy Exoplanet Forecasts for JWST: Maps, Retrieved Composition and Constraints on Formation with MIRI and NIRCam

abstract: The James Webb Space Telescope (JWST) will measure exoplanet transmission and eclipse spectroscopy at un-precedented precisions to better understand planet structure, dynamics, chemistry and formation. These are essential tools on the march towards biosignature searches on potentially habitable planets. We explore a range of exoplanet atmospheric conditions and forecast the expected results with JWST. We take realistic CHIMERA models that match existing Spitzer and HST results and simulate the spectra achievable with the JWST MIRI + NIRCam Guaranteed Time Observations (GTO) survey. We then retrieve atmospheric parameters from these spectra to estimate the precision to which the planets' atmospheric compositions can be measured. We find that emission spectra have well-constrained unimodal solutions but transmission spectra near 10X solar abundance and solar C/O ratios can suffer from bimodal solutions. Broad wavelength coverage as well as higher precision data can resolve bimodal solutions and provide dramatically better atmospheric parameter constraints. We find that metallicities can be measured to within 20% to 170%, which approaches the precisions on Solar System planets, and C/O ratios can be constrained to ~10% to 60%, assuming that observers can leverage short wavelength data to select the correct solution from the bimodal posteriors. These compositional precisions are sufficient to validate or refute predictions from disk formation models on final atmospheric abundances as long as their history is not erased by planet evolution processes. We also show the extent to which eclipse mapping with JWST is possible on our brightest system HD 189733 b.


title: Development of a Mach-Zehnder Modulator Photonic Local Oscillator Source

abstract: This paper describes the development of a photonic local oscillator (LO) source based on a 3-stage Mach-Zehnder modulator (MZM) device. The MZM laser synthesizer demonstrates the feasibility of providing the photonic reference LO for the Atacama Large Millimeter Array telescope located in Chile. This MZM approach to generating an LO by radio RF modulation of a monochromatic optical source provides the merits of wide frequency coverage of 4-130 GHz, tuning speed of about 0.2 seconds, and residual integrated phase noise performance of 0.3 degrees RMS at 100 GHz.


title: Limits on turbulent propagation of energy in cool-core clusters of galaxies

abstract: We place constraints on the propagation velocity of bulk turbulence within the intracluster medium of three clusters and an elliptical galaxy. Using Reflection Grating Spectrometer measurements of turbulent line broadening, we show that for these clusters, the 90% upper limit on turbulent velocities when accounting for instrumental broadening is too low to propagate energy radially to the cooling radius of the clusters within the required cooling time. In this way, we extend previous Hitomi-based analysis on the Perseus cluster to more clusters, with the intention of applying these results to a future, more extensive catalog. These results constrain models of turbulent heating in AGN feedback by requiring a mechanism which can not only provide sufficient energy to offset radiative cooling, but resupply that energy rapidly enough to balance cooling at each cluster radius.


title: Moduli identification methods in Type II compactifications

abstract: Recent work on four dimensional effective descriptions of the heterotic string has identified the moduli of such systems as being given by kernels of maps between ordinary Dolbeault cohomology groups. The maps involved are defined by the supergravity data of the background solutions. Such structure is seen both in the case of Calabi-Yau compactifications with non-trivial constraints on moduli arising from the gauge bundle and in the case of some non-Kahler compactifications of the theory. This description of the moduli has allowed the explicit computation of the moduli stabilization effects of a wide range of non-trivial gauge bundles on Calabi-Yau three-folds. In this paper we examine to what extent the ideas and techniques used in this work can be extended to the case of Type IIB string theory. Certain simplifications arise in the Type IIB case in comparison to the heterotic situation. However, complications also arise due to the richer supergravity data of the theory inducing a more involved map structure. We illustrate our discussion with several concrete examples of compactification of Type IIB string theory on conformal CICY three-folds with flux.


title: Viewport-Driven Rate-Distortion Optimized 360{\deg} Video Streaming

abstract: The growing popularity of virtual and augmented reality communications and 360{\deg} video streaming is moving video communication systems into much more dynamic and resource-limited operating settings. The enormous data volume of 360{\deg} videos requires an efficient use of network bandwidth to maintain the desired quality of experience for the end user. To this end, we propose a framework for viewport-driven rate-distortion optimized 360{\deg} video streaming that integrates the user view navigation pattern and the spatiotemporal rate-distortion characteristics of the 360{\deg} video content to maximize the delivered user quality of experience for the given network/system resources. The framework comprises a methodology for constructing dynamic heat maps that capture the likelihood of navigating different spatial segments of a 360{\deg} video over time by the user, an analysis and characterization of its spatiotemporal rate-distortion characteristics that leverage preprocessed spatial tilling of the 360{\deg} view sphere, and an optimization problem formulation that characterizes the delivered user quality of experience given the user navigation patterns, 360{\deg} video encoding decisions, and the available system/network resources. Our experimental results demonstrate the advantages of our framework over the conventional approach of streaming a monolithic uniformly encoded 360{\deg} video and a state-of-the-art reference method. Considerable video quality gains of 4 - 5 dB are demonstrated in the case of two popular 4K 360{\deg} videos.


title: Boosted Density Estimation Remastered

abstract: There has recently been a steady increase in the number iterative approaches to density estimation. However, an accompanying burst of formal convergence guarantees has not followed; all results pay the price of heavy assumptions which are often unrealistic or hard to check. The Generative Adversarial Network (GAN) literature --- seemingly orthogonal to the aforementioned pursuit --- has had the side effect of a renewed interest in variational divergence minimisation (notably $f$-GAN). We show that by introducing a weak learning assumption (in the sense of the classical boosting framework) we are able to import some recent results from the GAN literature to develop an iterative boosted density estimation algorithm, including formal convergence results with rates, that does not suffer the shortcomings other approaches. We show that the density fit is an exponential family, and as part of our analysis obtain an improved variational characterisation of $f$-GAN.


title: Protocol architectures for IoT domains

abstract: In this work we discuss proxy architectures which interconnect IoT domains running CoAP with the rest of Internet including micro datacenters and other domains building scalable hierarchical architectures. We assume that CoAP domain is terminated by an IoT proxy with cache, and we investigate several design issues with respect to successful data transmission, round trip delay and energy consumption. We present performance data for the case when proxy autonomously maintains data freshness which clearly point to efficient design choices.


title: A polynomial Ansatz for Norm-conserving Pseudopotentials

abstract: We show that efficient norm-conserving pseudopotentials for electronic structure calculations can be obtained from a polynomial Ansatz for the potential. Our pseudopotential is a polynomial of degree ten in the radial variable and fulfills the same smoothness conditions imposed by the Troullier-Martins method [Phys. Rev. B 43, 1993 (1991)] where pseudopotentials are represented by a polynomial of degree twenty-two. We compare our method to the Troullier-Martins approach in electronic structure calculations for diamond and iron in the bcc structure and find that the two methods perform equally well in calculations of the total energy. However, first and second derivatives of the total energy with respect to atomic coordinates converge significantly faster with the plane wave cutoff if the standard Troullier-Martins potentials are replaced by the pseudopotentials introduced here.


title: CalibNet: Geometrically Supervised Extrinsic Calibration using 3D Spatial Transformer Networks

abstract: 3D LiDARs and 2D cameras are increasingly being used alongside each other in sensor rigs for perception tasks. Before these sensors can be used to gather meaningful data, however, their extrinsics (and intrinsics) need to be accurately calibrated, as the performance of the sensor rig is extremely sensitive to these calibration parameters. A vast majority of existing calibration techniques require significant amounts of data and/or calibration targets and human effort, severely impacting their applicability in large-scale production systems. We address this gap with CalibNet: a self-supervised deep network capable of automatically estimating the 6-DoF rigid body transformation between a 3D LiDAR and a 2D camera in real-time. CalibNet alleviates the need for calibration targets, thereby resulting in significant savings in calibration efforts. During training, the network only takes as input a LiDAR point cloud, the corresponding monocular image, and the camera calibration matrix K. At train time, we do not impose direct supervision (i.e., we do not directly regress to the calibration parameters, for example). Instead, we train the network to predict calibration parameters that maximize the geometric and photometric consistency of the input images and point clouds. CalibNet learns to iteratively solve the underlying geometric problem and accurately predicts extrinsic calibration parameters for a wide range of mis-calibrations, without requiring retraining or domain adaptation. The project page is hosted at https://epiception.github.io/CalibNet


title: Enforcing constraints for interpolation and extrapolation in Generative Adversarial Networks

abstract: We suggest ways to enforce given constraints in the output of a Generative Adversarial Network (GAN) generator both for interpolation and extrapolation (prediction). For the case of dynamical systems, given a time series, we wish to train GAN generators that can be used to predict trajectories starting from a given initial condition. In this setting, the constraints can be in algebraic and/or differential form. Even though we are predominantly interested in the case of extrapolation, we will see that the tasks of interpolation and extrapolation are related. However, they need to be treated differently. For the case of interpolation, the incorporation of constraints is built into the training of the GAN. The incorporation of the constraints respects the primary game-theoretic setup of a GAN so it can be combined with existing algorithms. However, it can exacerbate the problem of instability during training that is well-known for GANs. We suggest adding small noise to the constraints as a simple remedy that has performed well in our numerical experiments. The case of extrapolation (prediction) is more involved. During training, the GAN generator learns to interpolate a noisy version of the data and we enforce the constraints. This approach has connections with model reduction that we can utilize to improve the efficiency and accuracy of the training. Depending on the form of the constraints, we may enforce them also during prediction through a projection step. We provide examples of linear and nonlinear systems of differential equations to illustrate the various constructions.


title: Activation cross-section measurement of proton induced reactions on cerium

abstract: In the framework of a systematic study of proton induced nuclear reactions on lanthanides we have measured the excitation functions on natural cerium for the production of 142,139,138m,137Pr, 141,139,137m,137g,135Ce and 133La up to 65 MeV proton energy using the activation method with stacked-foil irradiation technique and high-resolution gamma-ray spectrometry. The cross-sections of the investigated reactions were compared with the data retrieved from the TENDL-2014 and TENDL-2015 libraries, based on the latest version of the TALYS code system. No earlier experimental data were found in the literature. The measured cross-section data are important for further improvement of nuclear reaction models and for practical applications in nuclear medicine, other labeling and activation studies.


title: Generalized Optimization of High Capacity Compressive Imaging Systems

abstract: One of the greatest challenges in applying compressive sensing (CS) signal processing techniques to electromagnetic imaging applications is designing a sensing matrix that has good reconstruction capabilities. Compressive reflector antennas (CRA) are a class of antennas that have been shown to provide enhanced image reconstruction performance over traditional reflector antennas (TRA) when CS techniques are employed. In this paper, we present a unified CRA design method, which considers both the sensing capacity and efficiency of the antenna, and can be used for both compressive imaging and multiple-input multiple-output (MIMO) communication applications. The unified design method is assessed for a CRA configuration in which dielectric scatterers are added to the surface of a TRA. The design results demonstrate the ability of the unified design method to enhance the CS reconstruction capabilities of the CRA.


title: The Zeldovich approximation and wide-angle redshift-space distortions

abstract: The contribution of line-of-sight peculiar velocities to the observed redshift of objects breaks the translational symmetry of the underlying theory, modifying the predicted 2-point functions. These `wide angle effects' have mostly been studied using linear perturbation theory in the context of the multipoles of the correlation function and power spectrum. In this work we present the first calculation of wide angle terms in the Zeldovich approximation, which is known to be more accurate than linear theory on scales probed by the next generation of galaxy surveys. We present the exact result for dark matter and perturbatively biased tracers as well as the small angle expansion of the configuration- and Fourier-space two-point functions and the connection to the multi-frequency angular power spectrum. We compare different definitions of the line-of-sight direction and discuss how to translate between them. We show that wide angle terms can reach tens of percent of the total signal in a measurement at low redshift in some approximations, and that a generic feature of wide angle effects is to slightly shift the Baryon Acoustic Oscillation scale.


title: Sensing Matrix Design via Capacity Maximization for Block Compressive Sensing Applications

abstract: It is well established in the compressive sensing (CS) literature that sensing matrices whose elements are drawn from independent random distributions exhibit enhanced reconstruction capabilities. In many CS applications, such as electromagnetic imaging, practical limitations on the measurement system prevent one from generating sensing matrices in this fashion. Although one can usually randomized the measurements to some degree, these sensing matrices do not achieve the same reconstruction performance as the truly randomized sensing matrices. In this paper, we present a novel method, based upon capacity maximization, for designing sensing matrices with enhanced block-sparse signal reconstruction capabilities. Through several numerical examples, we demonstrate how our method significantly enhances reconstruction performance.


title: Nonmonotonic bias dependence of local spin accumulation signals in ferromagnet/semiconductor lateral spin-valve devices

abstract: We find extraordinary behavior of the local two-terminal spin accumulation signals in ferromagnet (FM)/semiconductor (SC) lateral spin-valve devices. With respect to the bias voltage applied between two FM/SC Schottky tunnel contacts, the local spin-accumulation signal can show nonmonotonic variations, including a sign inversion. A part of the nonmonotonic features can be understood qualitatively by considering the rapid reduction in the spin polarization of the FM/SC interfaces with increasing bias voltage. In addition to the sign inversion of the FM/SC interface spin polarization, the influence of the spin-drift effect in the SC layer and the nonlinear electrical spin conversion at a biased FM/SC contact are discussed.


title: Using mm-Waves for Secret Key Establishment

abstract: The fact that Millimeter Wave (mmWave) communication needs to be directional is usually perceived as a challenge; in this paper we argue that it enables efficient secret key sharing that are unconditionally secure from passive eavesdroppers, by building on packet erasures. We showcase the potential of our approach in two setups: mmWave-based WiFi networks and vehicle platooning. We show that in the first case, we can establish a few hundred secret bits with minimal changes to standard communication protocol; while in both cases, with the right choice of parameters, we can potentially establish keys in the order of tenths of Mbps. These first results are based on some simplifying assumptions, yet we believe they give incentives to further explore such techniques.


title: Can Decentralized Status Update Achieve Universally Near-Optimal Age-of-Information in Wireless Multiaccess Channels?

abstract: In an Internet-of-Things system where status data are collected from sensors and actuators for time-critical applications, the freshness of data is vital and can be quantified by the recently proposed age-of-information (AoI) metric. In this paper, we first consider a general scenario where multiple terminals share a common channel to transmit or receive randomly generated status packets. The optimal scheduling problem to minimize AoI is formulated as a restless multi-armed bandit problem. To solve the problem efficiently, we derive the Whittle's index in closed-form and establish the indexability thereof. Compared with existing work, we extend the index policy for AoI optimization to incorporate stochastic packet arrivals and optimal packet management (buffering the latest packet). Inspired by the index policy which has near-optimal performance but is centralized by nature, a decentralized status update scheme, i.e., the index-prioritized random access policy (IPRA), is further proposed, achieving universally near-optimal AoI performance and outperforming state-of-the-arts in the literature.


title: Deep Pose Consensus Networks

abstract: In this paper, we address the problem of estimating a 3D human pose from a single image, which is important but difficult to solve due to many reasons, such as self-occlusions, wild appearance changes, and inherent ambiguities of 3D estimation from a 2D cue. These difficulties make the problem ill-posed, which have become requiring increasingly complex estimators to enhance the performance. On the other hand, most existing methods try to handle this problem based on a single complex estimator, which might not be good solutions. In this paper, to resolve this issue, we propose a multiple-partial-hypothesis-based framework for the problem of estimating 3D human pose from a single image, which can be fine-tuned in an end-to-end fashion. We first select several joint groups from a human joint model using the proposed sampling scheme, and estimate the 3D poses of each joint group separately based on deep neural networks. After that, they are aggregated to obtain the final 3D poses using the proposed robust optimization formula. The overall procedure can be fine-tuned in an end-to-end fashion, resulting in better performance. In the experiments, the proposed framework shows the state-of-the-art performances on popular benchmark data sets, namely Human3.6M and HumanEva, which demonstrate the effectiveness of the proposed framework.


title: Hyperbolic Blockade: Suppression of the Photonic Density of States and the Spontaneous Emission Rate at the Interface with Conducting Medium

abstract: Surface scattering of free electrons strongly modifies the electromagnetic response near the interface. Due to the inherent anisotropy of the surface scattering that necessarily reverses the normal the interface component of the electron velocity while its tangential component may remain the same, a thin layer near a high-quality interface shows strong dielectric anisotropy. The formation of the resulting hyperbolic dispersion layers near the metal-dielectric interface strongly modifies the local density of states, and leads to orders of magnitude changes in all associated phenomena.


title: Relative singularity categories II: DG models

abstract: We study the relationship between singularity categories and relative singularity categories and discuss constructions of differential graded algebras of relative singularity categories. As consequences, we obtain structural results, which are known or generalise known results, on singularity categories of algebras with radical square zero, of non-commutative deformations of Kleinian singularities, of $SL_3(\mathbb{C})$-quotient singularities and of Gorenstein toric threefolds.


title: The Epistemology of Nondeterminism

abstract: This paper proposes new semantics for nondeterministic program execution, replacing the standard relational semantics for propositional dynamic logic (PDL). Under these new semantics, program execution is represented as fundamentally deterministic (i.e., functional), while nondeterminism emerges as an epistemic relationship between the agent and the system: intuitively, the nondeterministic outcomes of a given process are precisely those that cannot be ruled out in advance. We formalize these notions using topology and the framework of dynamic topological logic (DTL). We show that DTL can be used to interpret the language of PDL in a manner that captures the intuition above, and moreover that continuous functions in this setting correspond exactly to deterministic processes. We also prove that certain axiomatizations of PDL remain sound and complete with respect to the corresponding classes of dynamic topological models. Finally, we extend the framework to incorporate knowledge using the machinery of subset space logic, and show that the topological interpretation of public announcements coincides exactly with a natural interpretation of test programs.


title: Activation cross sections of Proton Induced Nuclear Reactions on neodymium up to 65 MeV

abstract: In the frame of a systematic study of the activation cross sections of charged particle induced nuclear reactions on rare earths for production of therapeutic radionuclides, proton induced reactions on neodymium were measured up to 65 MeV energy, above 45 MeV for the first time. The excitation functions of the natNd(p,x)150,149,148m,148g,146,144,143,141Pm, 149,147,141,140,139m,138Nd, 142,138mPr and 139gCe nuclear reactions were assessed by using stacked foil activation technique and high resolution gamma-spectrometry. The excitation functions were compared to the theoretical predictions, available in the TENDL-2015 library based on latest version of the TALYS code. The application of the data for medical isotope production is shortly discussed.


title: Extrapolating quantum observables with machine learning: Inferring multiple phase transitions from properties of a single phase

abstract: We present a machine-learning method for predicting sharp transitions in a Hamiltonian phase diagram by extrapolating the properties of quantum systems. The method is based on Gaussian Process regression with a combination of kernels chosen through an iterative procedure maximizing the predicting power of the kernels. The method is capable of extrapolating across the transition lines. The calculations within a given phase can be used to predict not only the closest sharp transition, but also a transition removed from the available data by a separate phase. This makes the present method particularly valuable for searching phase transitions in the parts of the parameter space that cannot be probed experimentally or theoretically.


title: Distortion of Magnetic Fields in a Starless Core III: Polarization--Extinction Relationship in FeSt 1-457

abstract: The relationship between dust polarization and extinction was determined for the cold dense starless molecular cloud core FeSt 1-457 based on the background star polarimetry of dichroic extinction at near-infrared wavelengths. Owing to the known (three-dimensional) magnetic field structure, the observed polarizations from the core were corrected by considering (a) the subtraction of the ambient polarization component, (b) the depolarization effect of inclined distorted magnetic fields, and (c) the magnetic inclination angle of the core. After these corrections, a linear relationship between polarization and extinction was obtained for the core in the range up to $A_V \approx 20$ mag. The initial polarization vs. extinction diagram changed dramatically after the corrections of (a) to (c), with the correlation coefficient being refined from 0.71 to 0.79. These corrections should affect the theoretical interpretation of the observational data. The slope of the finally obtained polarization--extinction relationship is $P_H / E_{H-K_s} = 11.00 \pm 0.72$ $\%$ ${\rm mag}^{-1}$, which is close to the statistically estimated upper limit of the interstellar polarization efficiency (Jones 1989). This consistency suggests that the upper limit of interstellar polarization efficiency might be determined by the observational viewing angle toward polarized astronomical objects.


title: Study of activation cross sections of deuteron induced reactions on erbium in the 32-50 MeV energy range

abstract: Activation cross sections of the natEr(d,x)163,165,166,167,168,170Tm and natEr(d,x)171,161Er nuclear reactions have been measured in the 32-50 MeV energy range, above 40 MeV for the first time. The activation method with stacked foil irradiation technique and gamma-ray spectroscopy were used. The experimental cross sections were compared with the theoretical predictions in the TENDL-2015 library.


title: SUCAG: Stochastic Unbiased Curvature-aided Gradient Method for Distributed Optimization

abstract: We propose and analyze a new stochastic gradient method, which we call Stochastic Unbiased Curvature-aided Gradient (SUCAG), for finite sum optimization problems. SUCAG constitutes an unbiased total gradient tracking technique that uses Hessian information to accelerate con- vergence. We analyze our method under the general asynchronous model of computation, in which each function is selected infinitely often with possibly unbounded (but sublinear) delay. For strongly convex problems, we establish linear convergence for the SUCAG method. When the initialization point is sufficiently close to the optimal solution, the established convergence rate is only dependent on the condition number of the problem, making it strictly faster than the known rate for the SAGA method. Furthermore, we describe a Markov-driven approach of implementing the SUCAG method in a distributed asynchronous multi-agent setting, via gossiping along a random walk on an undirected communication graph. We show that our analysis applies as long as the graph is connected and, notably, establishes an asymptotic linear convergence rate that is robust to the graph topology. Numerical results demonstrate the merits of our algorithm over existing methods.


title: Low-Frequency Carbon Recombination Lines in the Orion Molecular Cloud Complex

abstract: We detail tentative detections of low-frequency carbon radio recombination lines from within the Orion molecular cloud complex observed at 99-129 MHz. These tentative detections include one alpha transition and one beta transition over three locations and are located within the diffuse regions of dust observed in the infrared at 100umm, the Halpha emission detected in the optical, and the synchrotron radiation observed in the radio. With these observations, we are able to study the radiation mechanism transition from collisionally pumped to radiatively pumped within the HII regions of the Orion molecular cloud complex.


title: Randomness and Permutations in Coordinate Descent Methods

abstract: We consider coordinate descent (CD) methods with exact line search on convex quadratic problems. Our main focus is to study the performance of the CD method that use random permutations in each epoch and compare it to the performance of the CD methods that use deterministic orders and random sampling with replacement. We focus on a class of convex quadratic problems with a diagonally dominant Hessian matrix, for which we show that using random permutations instead of random with-replacement sampling improves the performance of the CD method in the worst-case. Furthermore, we prove that as the Hessian matrix becomes more diagonally dominant, the performance improvement attained by using random permutations increases. We also show that for this problem class, using any fixed deterministic order yields a superior performance than using random permutations. We present detailed theoretical analyses with respect to three different convergence criteria that are used in the literature and support our theoretical results with numerical experiments.


title: $L(\mathbb{R})$ with Determinacy Satisfies the Suslin Hypothesis

abstract: The Suslin hypothesis states that there are no nonseparable complete dense linear orderings without endpoints which have the countable chain condition. $\mathsf{ZF + AD^+ + V = L(\mathscr{P}(\mathbb{R}))}$ proves the Suslin hypothesis. In particular, if $L(\mathbb{R}) \models \mathsf{AD}$, then $L(\mathbb{R})$ satisfies the Suslin hypothesis, which answers a question of Foreman.