title: Constraining nuclear star cluster formation using MUSE-AO observations of the early-type galaxy FCC47

abstract: Nuclear star clusters (NSCs) are found in at least 70% of all galaxies, but their formation path is still unclear. In the most common scenarios, NSCs form in-situ from the galaxy's central gas reservoir, through merging of globular clusters (GCs), or through a combination of the two. As the scenarios pose different expectations for angular momentum and stellar population properties of the NSC in comparison to the host galaxy and the GC system, it is necessary to characterise the stellar light, NSC and GCs simultaneously. The large NSC (r$_\rm{eff} = 66$ pc) and rich GC system of the early-type Fornax cluster galaxy FCC47 (NGC1336) render this galaxy an ideal laboratory to constrain NSC formation. Using MUSE science verification data assisted by adaptive optics, we obtained maps for the stellar kinematics and for stellar-population properties of FCC47. We extracted the spectra of the central NSC and determined line-of-sight velocities of 24 GCs and metallicities of five. FCC47 shows two decoupled components (KDCs): a rotating disk and the NSC. Our orbit-based dynamical Schwarzschild model revealed that the NSC is a distinct kinematic feature and it constitutes the peak of metallicity and old ages in the galaxy. The main body consists of two counter-rotating populations and is dominated by a more metal-poor population. The GC system is bimodal with a dominant metal-poor population and the total GC system mass is $\sim 17\%$ of the NSC mass ($\sim$ 7 $\times$ $10^8 M_\odot$). The rotation, high metallicity and high mass of the NSC cannot be uniquely explained by GC-inspiral and most likely requires additional, but quickly quenched, in-situ formation. The presence of two KDCs most probably are evidence of a major merger that has altered the structure of FCC47 significantly, indicating the important role of galaxy mergers in forming the complex kinematics in the galaxy-NSC system.


title: PAGAN: Video Affect Annotation Made Easy

abstract: How could we gather affect annotations in a rapid, unobtrusive, and accessible fashion? How could we still make sure that these annotations are reliable enough for data-hungry affect modelling methods? This paper addresses these questions by introducing PAGAN, an accessible, general-purpose, online platform for crowdsourcing affect labels in videos. The design of PAGAN overcomes the accessibility limitations of existing annotation tools, which often require advanced technical skills or even the on-site involvement of the researcher. Such limitations often yield affective corpora that are restricted in size, scope and use, as the applicability of modern data-demanding machine learning methods is rather limited. The description of PAGAN is accompanied by an exploratory study which compares the reliability of three continuous annotation tools currently supported by the platform. Our key results reveal higher inter-rater agreement when annotation traces are processed in a relative manner and collected via unbounded labelling.


title: Finite free convolutions via Weingarten calculus

abstract: We consider the three finite free convolutions for polynomials studied in a recent paper by Marcus, Spielman, and Srivastava. Each can be described either by direct explicit formulae or in terms of operations on randomly rotated matrices. We present an alternate approach to the equivalence between these descriptions, based on combinatorial Weingarten methods for integration over the unitary and orthogonal groups. A key aspect of our approach is to identify a certain \emph{quadrature property}, which is satisfied by some important series of subgroups of the unitary groups (including the groups of unitary, orthogonal, and signed permutation matrices), and which yields the desired convolution formulae.


title: Folding pathways to crumpling in thermalized elastic frames

abstract: The mechanical properties of thermally excited two-dimensional crystalline membranes can depend dramatically on their geometry and topology. A particularly relevant example is the effect on the crumpling transition of holes in the membrane. Here we use molecular dynamics simulations to study the case of elastic frames (sheets with a single large hole in the center) and find that the system approaches the crumpled phase through a sequence of origami-like folds at decreasing length scales when temperature is increased. We use normal-normal correlation functions to quantify the temperature-dependent number of folds.


title: Learning Representations from Imperfect Time Series Data via Tensor Rank Regularization

abstract: There has been an increased interest in multimodal language processing including multimodal dialog, question answering, sentiment analysis, and speech recognition. However, naturally occurring multimodal data is often imperfect as a result of imperfect modalities, missing entries or noise corruption. To address these concerns, we present a regularization method based on tensor rank minimization. Our method is based on the observation that high-dimensional multimodal time series data often exhibit correlations across time and modalities which leads to low-rank tensor representations. However, the presence of noise or incomplete values breaks these correlations and results in tensor representations of higher rank. We design a model to learn such tensor representations and effectively regularize their rank. Experiments on multimodal language data show that our model achieves good results across various levels of imperfection.


title: Broad-band spectral study of X-ray transient MAXI J1820+070 using Swift/XRT and NuSTAR

abstract: We report on a \textit{NuSTAR} and \textit{Swift}/XRT observation of the newly discovered X-ray transient MAXI J1820+070. \textit{Swift}/XRT and \textit{NuSTAR} have concurrently observed the newly detected source on 14 March 2018. We have simultaneously fitted the broad-band spectra obtained from \textit{Swift}/XRT and \textit{NuSTAR}. The observed joint spectra in the energy range 0.6--78.0 keV are well modeled with a weak disk black-body emission, dominant thermal Comptonization and relativistic reflection fraction. We have detected a fluorescent Iron-K$\alpha$ line relativistically broadened, and a Compton hump at $\sim$ 30 keV. We constrain the inner disk radius as well as the disk inclination angle and their values are found to be 4.1$^{+0.8}_{-0.6}$ R$_{ISCO}$ (where R$_{ISCO}\equiv$ radius of the innermost stable circular orbit) or 5.1$^{+1.0}_{-0.7}$ r$_{g}$ (where r$_{g}\equiv$ gravitational radius) and 29.8$^{+3.0}_{-2.7}$ $^\circ$ respectively. The best fit broad-band spectra suggest that the source was in the hard state and evolving. The source emission is best described by weak thermal emission along with strong thermal Comptonization from a relatively cold, optically thick, geometrically thin and ionized accretion disk. X-ray spectral modeling helps us to understand the accretion and ejection properties in the vicinity of the compact object.


title: SN 2017gmr: An energetic Type II-P supernova with asymmetries

abstract: We present high-cadence ultraviolet (UV), optical, and near-infrared (NIR) data on the luminous Type II-P supernova SN 2017gmr from hours after discovery through the first 180 days. SN 2017gmr does not show signs of narrow, high-ionization emission lines in the early optical spectra, yet the optical lightcurve evolution suggests that an extra energy source from circumstellar medium (CSM) interaction must be present for at least 2 days after explosion. Modeling of the early lightcurve indicates a ~500R$_{\odot}$ progenitor radius, consistent with a rather compact red supergiant, and late-time luminosities indicate up to 0.130 $\pm$ 0.026 M$_{\odot}$ of $^{56}$Ni are present, if the lightcurve is solely powered by radioactive decay, although the $^{56}$Ni mass may be lower if CSM interaction contributes to the post-plateau luminosity. Prominent multi-peaked emission lines of H$\alpha$ and [O I] emerge after day 154, as a result of either an asymmetric explosion or asymmetries in the CSM. The lack of narrow lines within the first two days of explosion in the likely presence of CSM interaction may be an example of close, dense, asymmetric CSM that is quickly enveloped by the spherical supernova ejecta.


title: Probing the strange content of the proton with charm production in charged current at LHeC

abstract: We study charm production in charged-current deep-inelastic scattering (DIS) using the xFitter framework. Recent results from the LHC have focused renewed attention on the determination of the strange-quark parton distribution function (PDF) and the DIS charm process provides important complementary constraints on this quantity. We examine the current PDF uncertainty, and use LHeC pseudodata to estimate the potential improvement from this proposed facility. As xFitter implements both fixed-flavor- and variable-flavor-number schemes, we can compare the impact of these different theoretical choices; this highlights some interesting aspects of multi-scale calculations. We find that the high-statistics LHeC data covering a wide kinematic range could substantially reduce the strange PDF uncertainty.


title: Bulk superconductivity and role of fluctuations in the iron-based superconductor FeSe at high pressures

abstract: The iron-based superconductor FeSe offers a unique possibility to study the interplay of superconductivity with purely nematic as well magnetic-nematic order by pressure (p) tuning. By measuring specific heat under p up to 2.36GPa, we study the multiple phases in FeSe using a thermodynamic probe. We conclude that superconductivity is bulk across the entire p range and competes with magnetism. Our analysis suggests that superconducting and magnetic fluctuations exist over a wide temperature range above the respective bulk transition temperatures, whenever magnetism is present. These observations highlight similarities between FeSe and underdoped cuprate superconductors where fluctuations play a crucial role.


title: Topological edge states of interacting photon pairs emulated in a topolectrical circuit

abstract: Topological physics opens up a plethora of exciting phenomena allowing to engineer disorder-robust unidirectional flows of light. Recent advances in topological protection of electromagnetic waves suggest that even richer functionalities can be achieved by realizing topological states of quantum light. This area, however, remains largely uncharted due to the number of experimental challenges. Here, we take an alternative route and design a classical structure based on topolectrical circuits which serves as a simulator of a quantum-optical one-dimensional system featuring the topological state of two photons induced by the effective photon-photon interaction. Employing the correspondence between the eigenstates of the original problem and circuit modes, we use the designed simulator to extract the frequencies of bulk and edge two-photon bound states and evaluate the topological invariant directly from the measurements. Furthermore, we perform a reconstruction of the two-photon probability distribution for the topological state associated with one of the circuit eigenmodes.


title: Discovery of recombining plasma inside the extended gamma-ray supernova remnant HB9

abstract: We present the results from the Suzaku X-ray Imaging Spectrometer observation of the mixed-morphology supernova remnant (SNR) HB9 (G160.9+2.6). We discovered recombining plasma (RP) in the western Suzaku observation region and the spectra here are well described by a model having collisional ionization equilibrium (CIE) and RP components. On the other hand, the X-ray spectra from the eastern Suzaku observation region are best reproduced by the CIE and non-equilibrium ionization model. We discuss possible scenarios to explain the origin of the RP emission based on the observational properties and concluded that the rarefaction scenario is a possible explanation for the existence of RP. In addition, the gamma-ray emission morphology and spectrum within the energy range of 0.2$-$300 GeV are investigated using $\sim$10 years of data from the Fermi Large Area Telescope (LAT). The gamma-ray morphology of HB9 is best described by the spatial template of radio continuum emission. The spectrum is well-fit to a log-parabola function and its detection significance was found to be $\sim$25$\sigma$. Moreover, a new gamma-ray point source located just outside the south-east region of the SNR's shell was detected with a significance of $\sim$6$\sigma$. We also investigated the archival HI and CO data and detected an expanding shell structure in the velocity range of $-10.5$ and $+1.8$ km s$^{-1}$ that is coinciding with a region of gamma-ray enhancement at the southern rim of the HB9 shell.


title: On the Conditional Complexity of Sets of Strings

abstract: Given a set X of finite strings, one interesting question to ask is whether there exists a member of X which is simple conditional to all other members of X. Conditional simplicity is measured by low conditional Kolmogorov complexity. We prove the affirmative to this question for sets that have low mutual information with the halting sequence. There are two results with respect to this question. One is dependent on the maximum conditional complexity between two elements of X, the other is dependent on the maximum expected value of the conditional complexity of a member of X relative to each member of X.


title: Understanding Fault Scenarios and Impacts through Fault Injection Experiments in Cielo

abstract: We present a set of fault injection experiments performed on the ACES (LANL/SNL) Cray XE supercomputer Cielo. We use this experimental campaign to improve the understanding of failure causes and propagation that we observed in the field failure data analysis of NCSA's Blue Waters. We use the data collected from the logs and from network performance counter data 1) to characterize the fault-error-failure sequence and recovery mechanisms in the Gemini network and in the Cray compute nodes, 2) to understand the impact of failures on the system and the user applications at different scale, and 3) to identify and recreate fault scenarios that induce unrecoverable failures, in order to create new tests for system and application design. The faults were injected through special input commands to bring down network links, directional connections, nodes, and blades. We present extensions that will be needed to apply our methodologies of injection and analysis to the Cray XC (Aries) systems.


title: Chasing star-planet magnetic interactions: the case of Kepler-78

abstract: Observational evidence of star-planet magnetic interactions (SPMI) in compact exo-systems have been looked for in the past decades. Indeed planets in close-in orbit can be magnetically connected to their host star, and channel Alfv\'en waves carrying large amounts of energy towards the central star. The strength and temporal modulation of SPMIs are primarily set by the magnetic topology of the host star and the orbital characteristics of the planet. As a result, SPMI signals can be modulated over the rotational period of the star, the orbital period of the planet, or a complex combination of the two. The detection of SPMI thus have to rely on multiple-epochs and multiple-wavelengths observational campaigns. We present a new method to characterize SPMIs and apply it to Kepler-78, a late G star with a super-Earth on an 8.5 hours orbit. We model the corona of Kepler-78 using the large-scale magnetic topology of the star observed with Zeeman-Doppler-Imaging. We show that the closeness of Kepler-78b allows the interaction to channel energy flux densities up to a few kW m$^{-2}$ towards the central star. We show that this flux is large enough to be detectable in classical activity tracers such as H$\alpha$. It is nonetheless too weak to explain the modulation observed by \citet{Moutou2016}. We furthermore demonstrate how to predict the temporal modulation of SPMI signals in observed systems such as Kepler-78. The methodology presented here thus paves the road towards denser, specific observational campaigns that would allow a proper identification of SPMIs in compact star-planet systems.


title: Like-Charge Attraction between Metal Nanoparticles in a 1:1 Electrolyte Solution

abstract: We calculate the force between two spherical metal nanoparticles of charge Q 1 and Q 2 in a dilute 1:1 electrolyte solution. Numerically solving the non-linear Poisson-Boltzmann equation, we find that metal nanoparticles with the same sign of charge can attract one another. This is fundamentally different from what is found for like-charged, non-polarizable, colloidal particles, the two body interaction potential for which is always repulsive inside a dilute 1:1 electrolyte. Furthermore, existence of like-charge attraction between spherical metal nanoparticles is even more surprising in view of the result that such attraction is impossible between parallel metal slabs, showing the fundamental importance of curvature. To overcome a slow convergence of the numerical solution of the full non-linear Poisson-Boltzmann equation, we developed a modified Derjaguin approximation which allows us to accurate and rapidly calculate the interaction potential between two metal nanoparticles, or between a metal nanoparticle and a phospholipid membrane.


title: Rare Disease Detection by Sequence Modeling with Generative Adversarial Networks

abstract: Rare diseases affecting 350 million individuals are commonly associated with delay in diagnosis or misdiagnosis. To improve those patients' outcome, rare disease detection is an important task for identifying patients with rare conditions based on longitudinal medical claims. In this paper, we present a deep learning method for detecting patients with exocrine pancreatic insufficiency (EPI) (a rare disease). The contribution includes 1) a large longitudinal study using 7 years medical claims from 1.8 million patients including 29,149 EPI patients, 2) a new deep learning model using generative adversarial networks (GANs) to boost rare disease class, and also leveraging recurrent neural networks to model patient sequence data, 3) an accurate prediction with 0.56 PR-AUC which outperformed benchmark models in terms of precision and recall.


title: Diminishing the Effect of Adversarial Perturbations via Refining Feature Representation

abstract: Deep neural networks are highly vulnerable to adversarial examples, which imposes severe security issues for these state-of-the-art models. Many defense methods have been proposed to mitigate this problem. However, a lot of them depend on modification or additional training of the target model. In this work, we analytically investigate each layer's representation of non-perturbed and perturbed images and show the effect of perturbations on each of these representations. Accordingly, a method based on whitening coloring transform is proposed in order to diminish the misrepresentation of any desirable layer caused by adversaries. Our method can be applied to any layer of any arbitrary model without the need of any modification or additional training. Due to the fact that the full whitening of the layer's representation is not easily differentiable, our proposed method is superbly robust against white-box attacks. Furthermore, we demonstrate the strength of our method against some state-of-the-art black-box attacks.


title: Kayotee: A Fault Injection-based System to Assess the Safety and Reliability of Autonomous Vehicles to Faults and Errors

abstract: Fully autonomous vehicles (AVs), i.e., AVs with autonomy level 5, are expected to dominate road transportation in the near-future and contribute trillions of dollars to the global economy. The general public, government organizations, and manufacturers all have significant concern regarding resiliency and safety standards of the autonomous driving system (ADS) of AVs . In this work, we proposed and developed (a) `Kayotee' - a fault injection-based tool to systematically inject faults into software and hardware components of the ADS to assess the safety and reliability of AVs to faults and errors, and (b) an ontology model to characterize errors and safety violations impacting reliability and safety of AVs. Kayotee is capable of characterizing fault propagation and resiliency at different levels - (a) hardware, (b) software, (c) vehicle dynamics, and (d) traffic resilience. We used Kayotee to study a proprietary ADS technology built by Nvidia corporation and are currently applying Kayotee to other open-source ADS systems.


title: Neutral B-meson mixing from full lattice QCD at the physical point

abstract: We calculate the bag parameters for neutral $B$-meson mixing in and beyond the Standard Model, in full four-flavour lattice QCD for the first time. We work on gluon field configurations that include the effect of $u$, $d$, $s$ and $c$ sea quarks with the Highly Improved Staggered Quark (HISQ) action at three values of the lattice spacing and with three $u/d$ quark masses going down to the physical value. The valence $b$ quarks use the improved NRQCD action and the valence light quarks, the HISQ action. Our analysis was blinded. Our results for the bag parameters for all five operators are the most accurate to date. For the Standard Model operator between $B_s$ and $B_d$ mesons we find: $\hat{B}_{B_s}=1.232(53)$, $\hat{B}_{B_d}=1.222(61)$. Combining our results with lattice QCD calculations of the decay constants using HISQ quarks from the Fermilab/MILC collaboration and with experimental values for $B_s$ and $B_d$ oscillation frequencies allows determination of the CKM elements $V_{ts}$ and $V_{td}$. We find $V_{ts} = 0.04189(93)$, $V_{td} = 0.00867(23)$ and $V_{ts}/V_{td} = 0.2071(27)$. Our results agree well (within $2\sigma$) with values determined from CKM unitarity constraints based on tree-level processes (only). Using a ratio to $\Delta M$ in which CKM elements cancel in the Standard Model, we determine the branching fractions ${\text{Br}}(B_s\rightarrow \mu^+\mu^-) = 3.81(18) \times 10^{-9}$ and ${\text{Br}}(B_d\rightarrow \mu^+\mu^-) = 1.031(54) \times 10^{-10}$. We also give results for matrix elements of the operators $R_0$, $R_1$ and $\tilde{R}_1$ that contribute to neutral $B$-meson width differences.


title: Synergies between space telescopes in the photometric characterization of the atmospheres of Hot Jupiters

abstract: Previous generation of instruments have the opportunity to discover thousands of extra-solar planets and more will come with the current and future planet-search missions. In order to go one step further in the characterization of exoplanets, in this paper we describe a way to compare the photometric observation of Hot Jupiters done with space telescopes such as HST, CHEOPS, TESS, PLATO, and JWST and give the first-hand characterization on their atmospheres. We analyze a set of planetary systems hosting a Hot Jupiter for which an atmospheric template is available in the literature. For each system, we simulate the transit light curves observed by different instruments, convolving the incoming spectrum with the corresponding instrumental throughput. For each instrument, we thus measure the expected transit depth and estimate the associated uncertainty. Finally, we compare the transit depths as seen by the selected instruments and we quantify the effect of the planetary atmosphere on multi-band transit photometry. We also analyze a set of simulated scenarios with different stellar magnitudes, activity levels, transit durations and atmospheric templates to find the best cases for this kind of observational approach. We find in general that current and especially future space telescopes provide enough photometric precision to detect significant differences between the transit depths at different wavelengths. In particular, we find that the chromatic effect due to the atmosphere of the Hot Jupiters is maximized at later spectral types and that the effect of stellar activity is smaller than the measurement uncertainties.


title: Low-memory, discrete ordinates, discontinuous Galerkin methods for radiative transport

abstract: The discrete ordinates discontinuous Galerkin ($S_N$-DG) method is a well-established and practical approach for solving the radiative transport equation. In this paper, we study a low-memory variation of the upwind $S_N$-DG method. The proposed method uses a smaller finite element space that is constructed by coupling spatial unknowns across collocation angles, thereby yielding an approximation with fewer degrees of freedom than the standard method. Like the original $S_N$-DG method, the low memory variation still preserves the asymptotic diffusion limit and maintains the characteristic structure needed for mesh sweeping algorithms. While we observe second-order convergence in scattering dominated, diffusive regime, the low-memory method is in general only first-order accurate. To address this issue, we use upwind reconstruction to recover second-order accuracy. For both methods, numerical procedures based on upwind sweeps are proposed to reduce the system dimension in the underlying Krylov solver strategy.


title: Three stable phases and thermodynamic anomaly in a binary mixture of hard particles

abstract: While the realistically modeling of the thermodynamic behavior of fluids usually demands elaborated atomistic models, much have been learned from simplified ones. Here, we investigate a model where point-like particles (with activity $z_0$) are mixed with molecules that exclude their first and second neighbors (i.e., cubes of lateral size $\lambda=\sqrt{3}a$, with activity $z_2$), both placed on the sites of a simple cubic lattice with parameter $a$. Only hard-core interactions exist among the particles, so that the model is athermal. Despite its simplicity, the grand-canonical solution of this model on a Husimi lattice built with cubes revels a fluid-fluid demixing, yielding a phase diagram with two fluid phases (one of them dominated by small particles - $F0$) and a solid-like phase coexisting at a triple-point. Moreover, the fluid-fluid coexistence line ends at a critical point. An anomaly in the total density ($\rho_T$) of particles is also found, which is hallmarked by minima in the isobaric curves of $\rho_T$ versus $z_0$ (or $z_2$). Interestingly, the line of minimum density cross the phase diagram starting inside the region where both fluid phases are stable, passing through the $F0$ one and ending deep inside its metastable region, in a point where the spinodals of both fluid phases cross each other.


title: Influence of atomic motion on the collective effects in dense and cold atomic ensembles

abstract: We show that atomic motion leads not only to noticeable quantitative, but in some cases also to qualitative modification of collective effects in dense and cold atomic ensembles even in the case when the characteristic Doppler shifts are tens of times smaller than the natural linewidth. The observed influence is explained as a result of the suppression of the impact of sub-radiant collective states caused by the displacement of the atoms.


title: LSTM Language Models for LVCSR in First-Pass Decoding and Lattice-Rescoring

abstract: LSTM based language models are an important part of modern LVCSR systems as they significantly improve performance over traditional backoff language models. Incorporating them efficiently into decoding has been notoriously difficult. In this paper we present an approach based on a combination of one-pass decoding and lattice rescoring. We perform decoding with the LSTM-LM in the first pass but recombine hypothesis that share the last two words, afterwards we rescore the resulting lattice. We run our systems on GPGPU equipped machines and are able to produce competitive results on the Hub5'00 and Librispeech evaluation corpora with a runtime better than real-time. In addition we shortly investigate the possibility to carry out the full sum over all state-sequences belonging to a given word-hypothesis during decoding without recombination.


title: Condition-based Maintenance for Multi-component Systems:Modeling, Structural Properties, and Algorithms

abstract: Condition-based maintenance (CBM) is an effective maintenance strategy to improve system performance while lowering operating and maintenance costs. Real-world systems typically consist of a large number of components with various interactions between components. However, existing studies on CBM focus on single-component systems. Multi-component condition-based maintenance, which joins the components' stochastic degradation processes and the combinatorial maintenance grouping problem, remains an open issue in the literature. In this paper, we study the CBM optimization problem for multi-component systems. We first develop a multi-stage stochastic integer model with the objective of minimizing the total maintenance cost over a finite planning horizon. We then investigate the structural properties of a two-stage model. Based on the structural properties, two efficient algorithms are designed to solve the two-stage model. Algorithm 1 solves the problem to its optimality and Algorithm 2 heuristically searches for high-quality solutions based on Algorithm 1. Our computational studies show that Algorithm 1 obtains optimal solutions in a reasonable amount of time and Algorithm 2 can find high-quality solutions quickly. The multi-stage problem is solved using a rolling horizon approach based on the algorithms for the two-stage problem.


title: On Slicing Sorted Integer Sequences

abstract: Representing sorted integer sequences in small space is a central problem for large-scale retrieval systems such as Web search engines. Efficient query resolution, e.g., intersection or random access, is achieved by carefully partitioning the sequences. In this work we describe and compare two different partitioning paradigms: partitioning by cardinality and partitioning by universe. Although the ideas behind such paradigms have been known in the coding and algorithmic community since many years, inverted index compression has extensively adopted the former paradigm, whereas the latter has received only little attention. As a result, an experimental comparison between these two is missing for the setting of inverted index compression. We also propose and implement a solution that recursively slices the universe of representation of a sequence to achieve compact storage and attain to fast query execution. Albeit larger than some state-of-the-art representations, this slicing approach substantially improves the performance of list intersections and unions while operating in compressed space, thus offering an excellent space/time trade-off for the problem.


title: Kinematic Constraints on Spatial Curvature from Supernovae Ia and Cosmic Chronometers

abstract: An approach to estimate the spatial curvature $\Omega_k$ from data independently of dynamical models is suggested, through kinematic parameterizations of the comoving distance ($D_{C}(z)$) with third degree polynomial, of the Hubble parameter ($H(z)$) with a second degree polynomial and of the deceleration parameter ($q(z)$) with first order polynomial. All these parameterizations were done as function of redshift $z$. We used SNe Ia dataset from Pantheon compilation with 1048 distance moduli estimated in the range $0.01<z<2.3$ with systematic and statistical errors and a compilation of 31 $H(z)$ data estimated from cosmic chronometers. The spatial curvature found for $D_C(z)$ parametrization was $\Omega_{k}=-0.03^{+0.24+0.56}_{-0.30-0.53}$. The parametrization for deceleration parameter $q(z)$ resulted in $\Omega_{k}=-0.08^{+0.21+0.54}_{-0.27-0.45}$. The $H(z)$ parametrization has shown incompatibilities between $H(z)$ and SNe Ia data constraints, so these analyses were not combined. The $D_C(z)$ and $q(z)$ parametrizations are compatible with the spatially flat Universe as predicted by many inflation models and data from CMB. This type of analysis is very appealing as it avoids any bias because it does not depend on assumptions about the matter content of the Universe for estimating $\Omega_k$.


title: Learning to aggregate feature representations

abstract: The Algonauts challenge requires to construct a multi-subject encoder of images to brain activity. Deep networks such as ResNet-50 and AlexNet trained for image classification are known to produce feature representations along their intermediate stages which closely mimic the visual hierarchy. However the challenges introduced in the Algonauts project, including combining data from multiple subjects, relying on very few similarity data points, solving for various ROIs, and multi-modality, require devising a flexible framework which can efficiently accommodate them. Here we build upon a recent state-of-the-art classification network (SE-ResNeXt-50) and construct an adaptive combination of its intermediate representations. While the pretrained network serves as a backbone of our model, we learn how to aggregate feature representations along five stages of the network. During learning, our method enables to modulate and screen outputs from each stage along the network as governed by the optimized objective. We applied our method to the Algonauts2019 fMRI and MEG challenges. Using the combined fMRI and MEG data, our approach was rated among the leading five for both challenges. Surprisingly we find that for both the lower and higher order areas (EVC and IT) the adaptive aggregation favors features stemming at later stages of the network.


title: Capacity Analysis for Joint Radar-Communication Capable Coherent MIMO Radars

abstract: Recently, huge attention is attracted to the concept of integrating communication and radar missions within the same platform. Joint radar-communications (JRC) system gives an important opportunity to reduce spectrum usage and product cost while doing concurrent operation, as target sensing via radar processing and establishing communication links. A JRC-capable coherent MIMO radar system have been proposed recently in the literature. Several methods are introduced to reach dual goal as a notable null level towards the direction of interest of the radar and MIMO radar waveform orthogonality. Due to the limitations originated form the JRC operation, communication channel may encounter unwanted amplitude variations. This unwanted modulation normally affects the communication performance by its nature, due to the fades on radiated signal amplitude towards the direction of communication. However, the effect of this unintentional modulation on communication channel is yet to be investigated. In this paper, the communication channel for JRC capable phase-coded coherent MIMO radars is analyzed and investigated under additive white Gaussian noise and Rayleigh/Rician fading conditions. Communication capacity is evaluated for each channel condition. The results reveal that, using the single-side limited null direction fixed waveform generation method displays the best capacity performance under all channel conditions.


title: On the multiloop soft theorem of the dilaton in the bosonic string

abstract: In this note we show that by fixing the multiloop Green function in the closed bosonic string to be Arakelov's Green function, one obtains factorization of scattering amplitudes with a softly emitted dilaton to the same level as with a graviton to all loop order. This extends our previous analysis at one loop to all loop orders and confirms that some high-energy quantum symmetry in the bosonic string protects the factorization of amplitudes with softly emitted dilatons.


title: Tropical geometry over the tropical hyperfield

abstract: In this text, we merge ideas around the tropical hyperfield with the theory of ordered blueprints to give a new formulation of tropical scheme theory. The key insight is that a nonarchimedean absolute value can be considered as a morphism into the tropical hyperfield. In turn, ordered blueprints make it possible to consider the base change of a classical variety to the tropical hyperfield. We call this base change the scheme theoretic tropicalization of the classical variety. Our first main result describes the Berkovich analytification and the tropicalization of a classical variety as sets of rational points of scheme theoretic tropicalizations, including a characterization of the respective topologies. Our second main result shows that the Giansiracusa bend relations can be derived by a natural construction from the scheme theoretic tropicalization.


title: AVFI: Fault Injection for Autonomous Vehicles

abstract: Autonomous vehicle (AV) technology is rapidly becoming a reality on U.S. roads, offering the promise of improvements in traffic management, safety, and the comfort and efficiency of vehicular travel. With this increasing popularity and ubiquitous deployment, resilience has become a critical requirement for public acceptance and adoption. Recent studies into the resilience of AVs have shown that though the AV systems are improving over time, they have not reached human levels of automation. Prior work in this area has studied the safety and resilience of individual components of the AV system (e.g., testing of neural networks powering the perception function). However, methods for holistic end-to-end resilience assessment of AV systems are still non-existent.


title: Revealing the strokes of autonomous quantum heat engines with work and heat fluctuations

abstract: We analyze an autonomous thermoelectric engine composed of two superconducting qubits coupled to separate heat baths and connected by a Josephson junction. Work and heat are process quantities and not observables of the engine quantum system, but their rates can be derived from the steady state expectation value of appropriate system observables and their fluctuations are given by correlation functions determined with the GKLS master equation and quantum regression theorem. The correlation functions also reveal a cyclical, dynamical transfer of energy - the strokes of the engine - hiding underneath the steady state.


title: The Sensitivity of Counterfactual Fairness to Unmeasured Confounding

abstract: Causal approaches to fairness have seen substantial recent interest, both from the machine learning community and from wider parties interested in ethical prediction algorithms. In no small part, this has been due to the fact that causal models allow one to simultaneously leverage data and expert knowledge to remove discriminatory effects from predictions. However, one of the primary assumptions in causal modeling is that you know the causal graph. This introduces a new opportunity for bias, caused by misspecifying the causal model. One common way for misspecification to occur is via unmeasured confounding: the true causal effect between variables is partially described by unobserved quantities. In this work we design tools to assess the sensitivity of fairness measures to this confounding for the popular class of non-linear additive noise models (ANMs). Specifically, we give a procedure for computing the maximum difference between two counterfactually fair predictors, where one has become biased due to confounding. For the case of bivariate confounding our technique can be swiftly computed via a sequence of closed-form updates. For multivariate confounding we give an algorithm that can be efficiently solved via automatic differentiation. We demonstrate our new sensitivity analysis tools in real-world fairness scenarios to assess the bias arising from confounding.


title: Natural Language Understanding with the Quora Question Pairs Dataset

abstract: This paper explores the task Natural Language Understanding (NLU) by looking at duplicate question detection in the Quora dataset. We conducted extensive exploration of the dataset and used various machine learning models, including linear and tree-based models. Our final finding was that a simple Continuous Bag of Words neural network model had the best performance, outdoing more complicated recurrent and attention based models. We also conducted error analysis and found some subjectivity in the labeling of the dataset.


title: OSNR limitations of chip-based optical frequency comb sources for WDM coherent communications

abstract: Optical frequency combs have the potential to become key building blocks of optical communication subsystems. The strictly equidistant, narrow-band spectral lines of a frequency comb can serve both as carriers for massively parallel data transmission and as local oscillator for coherent reception. Recent experiments have demonstrated the viability of various chip-based comb generator concepts for communication applications, offering transmission capacities of tens of Tbit/s. Here, we investigate the influence of the comb line power and of the carrier-to-noise power ratio on the performance of a frequency comb in a WDM system. We distinguish two regimes of operation depending on whether the comb source or the transmission link limits the performance of the system, i.e., defines the link reach, restricts the choice of modulation format and sets the maximum symbol rate. Finally, we investigate the achievable OSNR and channel capacity when using the tones of a soliton Kerr frequency comb as multi-wavelength carriers for WDM systems.


title: Spin transport in an electrically-driven magnon gas near Bose-Einstein condensation: Hartree-Fock-Keldysh theory

abstract: An easy-plane ferromagnetic insulator in a uniform external magnetic field and in contact with a phonon bath and a normal metal bath is studied theoretically in the presence of dc spin current injection via the spin Hall effect in the metal. The Keldysh path integral formalism is used to model the magnon gas driven into a nonequilibrium steady state by mismatched bath temperatures and/or electrical injection, and we analyze the magnon system in the normal (uncondensed) state, but close to the instability to Bose-Einstein condensation (BEC), within the self-consistent Hartree-Fock approximation. We find that the steady state magnon distribution function generally has a non-thermal form that cannot be described by a single effective chemical potential and effective temperature. We also show that the BEC instability in the electrically-driven magnon system is signaled by a sign change in the imaginary part of the poles for long-wavelength magnon modes and by the divergence of the nonequilibrium magnon distribution function. In the presence of two bath temperatures, we find that the correlation length of the superfluid order parameter fluctuations exhibits nontrivial finite temperature crossover behaviors that are richer than the standard crossover behaviors obtained for the vacuum-superfluid transition in an equilibrium dilute Bose gas. We study the consequences of these thermal crossovers on the magnon spin conductivity and obtain an inverse square-root divergence in the spin conductivity in the vicinity of the electrically-induced BEC instability. A spintronics device capable of testing our spin transport predictions is discussed.


title: Dissipative Particle Dynamics for Directed Self-Assembly of Block Copolymers

abstract: The dissipative particle dynamics (DPD) simulation method has been shown to be a promising tool to study self-assembly of soft matter systems. In particular, it has been used to study block copolymer (BCP) self-assembly. However, previous parametrizations of this model are not able to capture most of the rich phase behaviors of block copolymers in thin films nor in directed self-assembly (chemoepitaxy or graphoepitaxy). Here we extend the applicability of the DPD method for BCPs to make it applicable to thin films and directed self-assembly. Our new reparametrization is able to reproduce the bulk phase behavior, but also manages to predict thin film structures obtained experimentally from chemoepitaxy or graphoepitaxy. A number of different complex structures, such as bilayer nanomeshes, 90{\deg} bend structures, circular cylinders/lamellae and Frank-Kasper phases directed by trenches, post arrays or chemically patterned substrate have all been reproduced in this work. This reparametrized DPD model should serves as a powerful tool to predict BCP self-assembly, especially in some complex systems where it is difficult to implement SCFT.


title: Magnon decay theory of Gilbert damping in metallic antiferromagnets

abstract: Gilbert damping is a key property governing magnetization dynamics in ordered magnets. We present a theoretical study of intrinsic Gilbert damping induced by magnon decay in antiferromagnetic metals through $s$-$d$ exchange interaction. Our theory delineates the qualitative features of damping in metallic antiferromagnets owing to their bipartite nature, in addition to providing analytic expressions for the damping parameters. Magnon-induced intraband electron scattering is found to predominantly cause magnetization damping, whereas the N\'eel field is found to be damped via disorder. Depending on the conduction electron band structure, we predict that magnon-induced interband electron scattering around band crossings may be exploited to engineer a strong N\'eel field damping.


title: A Scalable Architecture for Power Consumption Monitoring in Industrial Production Environments

abstract: Detailed knowledge about the electrical power consumption in industrial production environments is a prerequisite to reduce and optimize their power consumption. Today's industrial production sites are equipped with a variety of sensors that, inter alia, monitor electrical power consumption in detail. However, these environments often lack an automated data collation and analysis. We present a system architecture that integrates different sensors and analyzes and visualizes the power consumption of devices, machines, and production plants. It is designed with a focus on scalability to support production environments of various sizes and to handle varying loads. We argue that a scalable architecture in this context must meet requirements for fault tolerance, extensibility, real-time data processing, and resource efficiency. As a solution, we propose a microservice-based architecture augmented by big data and stream processing techniques. Applying the fog computing paradigm, parts of it are deployed in an elastic, central cloud while other parts run directly, decentralized in the production environment. A prototype implementation of this architecture presents solutions how different kinds of sensors can be integrated and their measurements can be continuously aggregated. In order to make analyzed data comprehensible, it features a single-page web application that provides different forms of data visualization. We deploy this pilot implementation in the data center of a medium-sized enterprise, where we successfully monitor the power consumption of 16~servers. Furthermore, we show the scalability of our architecture with 20,000~simulated sensors.


title: Investigating The Piece-Wise Linearity And Benchmark Related To Koczy-Hirota Fuzzy Linear Interpolation

abstract: Fuzzy Rule Interpolation (FRI) reasoning methods have been introduced to address sparse fuzzy rule bases and reduce complexity. The first FRI method was the Koczy and Hirota (KH) proposed "Linear Interpolation". Besides, several conditions and criteria have been suggested for unifying the common requirements FRI methods have to satisfy. One of the most conditions is restricted the fuzzy set of the conclusion must preserve a Piece-Wise Linearity (PWL) if all antecedents and consequents of the fuzzy rules are preserving on PWL sets at {\alpha}-cut levels. The KH FRI is one of FRI methods which cannot satisfy this condition. Therefore, the goal of this paper is to investigate equations and notations related to PWL property, which is aimed to highlight the problematic properties of the KH FRI method to prove its efficiency with PWL condition. In addition, this paper is focusing on constructing benchmark examples to be a baseline for testing other FRI methods against situations that are not satisfied with the linearity condition for KH FRI.


title: Rate Distortion Theorem and the Multicritical Point of Spin Glass

abstract: A spin system can be thought of as an information coding system that transfers information of the interaction configuration into information of the equilibrium state of the spin variables. Hence it can be expected that the relations between the interaction configuration and equilibrium states are consistent with the known laws of information theory. We show that Shannon's rate-distortion theorem can be used to obtain an universal constraint on neighboring spin correlations for a broad range of Ising spin systems with two-body spin interactions. Remarkably, this constraint gives a bound for the multicritical point in the phase diagram, when a mean-field behavior for the neighboring spin pairs can be expected in the paramagnetic phase.


title: Permutation inference with a finite number of heterogeneous clusters

abstract: I introduce a simple permutation procedure to test conventional (non-sharp) hypotheses about the effect of a binary treatment in the presence of a finite number of large, heterogeneous clusters when the treatment effect is identified by comparisons across clusters. The procedure asymptotically controls size by applying a level-adjusted permutation test to a suitable statistic. The adjustments needed for most empirically relevant situations are tabulated in the paper. The adjusted permutation test is easy to implement in practice and performs well at conventional levels of significance with at least four treated clusters and a similar number of control clusters. It is particularly robust to situations where some clusters are much more variable than others. Examples and an empirical application are provided.


title: On the properties of a Static and a Stationary Charged Black Hole in $f(R)$ Gravity

abstract: In a recent paper, Nashed and Capozziello presented a new class of charged, spherically symmetric black hole solutions of $f(R)$ gravity with an asymptotic flat or (anti-)de Sitter behavior. These metrics depend on a dimensional parameter $\alpha$ and are interesting because they cannot reduce to general relativity solutions. In this paper, we present a corrected study of their physical and thermodynamic properties and generalize these solutions to obtain a new set of stationary, axisymmetric black holes in the $f(R)$ scenario. Some of our results show that the entropy is always positive within the allowed values of parameter $\alpha$ and due to the well behaved quantities such as the Gibbs free energy, we conclude that there is no such a phase transition as discussed in the work of Nashed and Capozziello. We also study the geodesics in these spacetimes and particularly, the stability of the circular orbits to obtain the radius of the Innermost Stable Circular Orbit.


title: ML-based Fault Injection for Autonomous Vehicles: A Case for Bayesian Fault Injection

abstract: The safety and resilience of fully autonomous vehicles (AVs) are of significant concern, as exemplified by several headline-making accidents. While AV development today involves verification, validation, and testing, end-to-end assessment of AV systems under accidental faults in realistic driving scenarios has been largely unexplored. This paper presents DriveFI, a machine learning-based fault injection engine, which can mine situations and faults that maximally impact AV safety, as demonstrated on two industry-grade AV technology stacks (from NVIDIA and Baidu). For example, DriveFI found 561 safety-critical faults in less than 4 hours. In comparison, random injection experiments executed over several weeks could not find any safety-critical faults


title: Structural phase diagram and magnetic properties of Sc-substituted rare earth ferrites R1-xScxFeO3 (R=Lu, Yb, Er, and Ho)

abstract: We have studied the structural stability of Sc-substituted rare earth (R) ferrites R1-xScxFeO3, and constructed a structural phase diagram for different R and x. While RFeO3 and ScFeO3 adopt the orthorhombic and the bixbyite structure respectively, the substituted compound R1-xScxFeO3 may be stable in a different structure. Specifically, for R0.5Sc0.5FeO3, the hexagonal structure can be stable for small R, such as Lu and Yb, while the garnet structure is stable for larger R, such as Er and Ho. The formation of garnet structure of the R0.5Sc0.5FeO3 compounds which requires that Sc occupies both the rare earth and the Fe sites, is corroborated by their magnetic properties.


title: Tropical Grassmannians, cluster algebras and scattering amplitudes

abstract: We provide a cluster-algebraic approach to the computation of the recently introduced generalised biadjoint scalar amplitudes related to Grassmannians ${\rm Gr}(k,n)$. A finite cluster algebra provides a natural triangulation for the tropical Grassmannian whose volume computes the scattering amplitudes. Using this method one can construct the entire colour-ordered amplitude via mutations starting from a single term.


title: Unweighted Donaldson-Thomas theory of the banana 3-fold with section classes

abstract: We further the study of the Donaldson-Thomas theory of the banana threefolds which were recently discovered and studied in [Bryan'19]. These are smooth proper Calabi-Yau threefolds which are fibred by Abelian surfaces such that the singular locus of a singular fibre is a non-normal toric curve known as a "banana configuration". In [Bryan'19] the Donaldson-Thomas partition function for the rank 3 sub-lattice generated by the banana configurations is calculated. In this article we provide calculations with a view towards the rank 4 sub-lattice generated by a section and the banana configurations. We relate the findings to the Pandharipande-Thomas theory for a rational elliptic surface and present new Gopakumar-Vafa invariants for the banana threefold.


title: Is artificial data useful for biomedical Natural Language Processing algorithms?

abstract: A major obstacle to the development of Natural Language Processing (NLP) methods in the biomedical domain is data accessibility. This problem can be addressed by generating medical data artificially. Most previous studies have focused on the generation of short clinical text, and evaluation of the data utility has been limited. We propose a generic methodology to guide the generation of clinical text with key phrases. We use the artificial data as additional training data in two key biomedical NLP tasks: text classification and temporal relation extraction. We show that artificially generated training data used in conjunction with real training data can lead to performance boosts for data-greedy neural network algorithms. We also demonstrate the usefulness of the generated data for NLP setups where it fully replaces real training data.


title: Optimal Bookmaking

abstract: We introduce a general framework for continuous-time betting markets, in which a bookmaker can dynamically control the prices of bets on outcomes of random events. In turn, the prices set by the bookmaker affect the rate or intensity of bets placed by gamblers. The bookmaker seeks a price process that maximizes his expected (utility of) terminal wealth. We obtain explicit solutions or characterizations to the bookmaker's optimal bookmaking problem in various interesting models.


title: Computing an order complete basis for $M^{\infty}(N)$ and Applications

abstract: This paper gives a quick way to construct all modular functions for the group $\Gamma_0(N)$ having only a pole at $\tau = i \infty$. We assume that we are given two modular functions $f,g$ for $\Gamma_0(N)$ with poles only at $i \infty$ and coprime pole orders. As an application we obtain two new identities from which one can derive that $p(11n+6)\equiv 0\pmod{11}$, here $p(n)$ is the usual partition function.


title: Associative Embedding for Game-Agnostic Team Discrimination

abstract: Assigning team labels to players in a sport game is not a trivial task when no prior is known about the visual appearance of each team. Our work builds on a Convolutional Neural Network (CNN) to learn a descriptor, namely a pixel-wise embedding vector, that is similar for pixels depicting players from the same team, and dissimilar when pixels correspond to distinct teams. The advantage of this idea is that no per-game learning is needed, allowing efficient team discrimination as soon as the game starts. In principle, the approach follows the associative embedding framework introduced in arXiv:1611.05424 to differentiate instances of objects. Our work is however different in that it derives the embeddings from a lightweight segmentation network and, more fundamentally, because it considers the assignment of the same embedding to unconnected pixels, as required by pixels of distinct players from the same team. Excellent results, both in terms of team labelling accuracy and generalization to new games/arenas, have been achieved on panoramic views of a large variety of basketball games involving players interactions and occlusions. This makes our method a good candidate to integrate team separation in many CNN-based sport analytics pipelines.


title: Probing Non-Standard Neutrino Interactions with Supernova Neutrinos at Hyper-K

abstract: Non-standard neutrino self interactions (NSSI) could be stronger than Fermi interactions. We investigate the ability to constrain these four-neutrino interactions by their effect on the flux of neutrinos originating from a galactic supernova. In the dense medium of a core collapse supernova, these new self interactions can have a significant impact on neutrino oscillations, leading to changes at the flavor evolution and spectra level. We use simulations of the neutrino flux from a 13 solar mass, core collapse supernova at 10 kpc away, and numerically propagate these neutrinos through the stellar medium taking into account vacuum/MSW oscillations, SM $\nu-\nu$ scattering as well as $\nu-\nu$ interactions that arise from NSSI. We pass the resulting neutrino flux to a simulation of the future Hyper-Kamiokande detector to see what constraints on NSSI parameters are possible when the next galactic supernova becomes visible. We find that these constraints depend strongly on the neutrino mass hierarchy and if the NSSI is flavor-violating or preserving. Sensitivity to NSSI in the normal hierarchy (NH) at Hyper-K is limited by the experiment's ability to efficiently detect $\nu_{e}$, but deviations from no NSSI could be seen if the NSSI is particularly strong. In the inverted hierarchy (IH) scenario, Hyper-K can significantly improve constraints on flavor-violating NSSI down to $\mathcal{O}(10^{-1})G_{F}$.


title: Lecture Notes on Stochastic Processes

abstract: This is lecture notes on the course "Stochastic Processes". In this format, the course was taught in the spring semesters 2017 and 2018 for third-year bachelor students of the Department of Control and Applied Mathematics, School of Applied Mathematics and Informatics at Moscow Institute of Physics and Technology. The base of this course was formed and taught for decades by professors from the Department of Mathematical Foundations of Control A.A. Natan, S.A. Guz, and O.G. Gorbachev. Besides standard chapters of stochastic processes theory (correlation theory, Markov processes) in this book (and lectures) the following chapters are included: von Neumann-Birkhoff-Khinchin ergodic theorem, macrosystem equilibrium concept, Markov Chain Monte Carlo, Markov decision processes and the secretary problem.


title: Thermoacoustic Tomography with Circular Integrating Detectors and Variable Wave Speed

abstract: We explore Thermoacoustic Tomography with circular integrating detectors assuming variable, smooth wave speed. We show that the measurement operator in this case is a Fourier Integral Operator and examine how the singularities in initial data and measured data are related through the canonical relation of this operator. We prove which of those singularities in the initial data are visible from a fixed open subset of the set on which measurements are taken. In addition, numerical results are shown for both full and partial data.


title: DeepTEGINN: Deep Learning Based Tools to Extract Graphs from Images of Neural Networks

abstract: In the brain, the structure of a network of neurons defines how these neurons implement the computations that underlie the mind and the behavior of animals and humans. Provided that we can describe the network of neurons as a graph, we can employ methods from graph theory to investigate its structure or use cellular automata to mathematically assess its function. Although, software for the analysis of graphs and cellular automata are widely available. Graph extraction from the image of networks of brain cells remains difficult. Nervous tissue is heterogeneous, and differences in anatomy may reflect relevant differences in function. Here we introduce a deep learning based toolbox to extracts graphs from images of brain tissue. This toolbox provides an easy-to-use framework allowing system neuroscientists to generate graphs based on images of brain tissue by combining methods from image processing, deep learning, and graph theory. The goals are to simplify the training and usage of deep learning methods for computer vision and facilitate its integration into graph extraction pipelines. In this way, the toolbox provides an alternative to the required laborious manual process of tracing, sorting and classifying. We expect to democratize the machine learning methods to a wider community of users beyond the computer vision experts and improve the time-efficiency of graph extraction from large brain image datasets, which may lead to further understanding of the human mind.


title: GPU-based Parallel Computation Support for Stan

abstract: This paper details an extensible OpenCL framework that allows Stan to utilize heterogeneous compute devices. It includes GPU-optimized routines for the Cholesky decomposition, its derivative, other matrix algebra primitives and some commonly used likelihoods, with more additions planned for the near future. Stan users can now benefit from large speedups offered by GPUs with little effort and without changes to their existing Stan code. We demonstrate the practical utility of our work with two examples - logistic regression and Gaussian Process regression.


title: Channel-Correlation-Enabled Transmit Optimization for MISO Wiretap Channels

abstract: An artificial-noise (AN)-aided beamformer specific to correlated main and wiretap channels is designed in this paper. We consider slow-fading multiple-input-single-output (MISO) wiretap channels with a passive single-antenna eavesdropper, in which independent transmitter-side and correlated receiver-side are assumed. Additionally, the source has accurate main channel information and statistical wiretap channel information. To reduce the secrecy loss due to receiver-side correlation, this paper proposes the scheme of channel-correlation-enabled transmit optimization. Particularly, the correlation is viewed as a resource to acquire more knowledge about wiretap channel. Based on this, the statistical distribution of wiretap channel is described more precisely. Then, the power of AN in the null space of main channel is placed more reasonably instead of simple uniform distribution and an elaborate beamformer for the information-bearing signal is designed. Finally, an efficient algorithm for power allocation between the information-bearing signal and the AN is developed. Simulation results show that the secrecy rate under transmit power and secrecy outage constraint is improved.


title: Global existence of the nonisentropic compressible Euler equations with vacuum boundary surrounding a variable entropy state

abstract: Global existence for the nonisentropic compressible Euler equations with vacuum boundary for all adiabatic constants $\gamma > 1$ is shown through perturbations around a rich class of background nonisentropic affine motions. The notable feature of the nonisentropic motion lies in the presence of non-constant entropies, and it brings a new mathematical challenge to the stability analysis of nonisentropic affine motions. In particular, the estimation of the curl terms requires a careful use of algebraic, nonlinear structure of the pressure. With suitable regularity of the underlying affine entropy, we are able to adapt the weighted energy method developed for the isentropic Euler by Had\v{z}i\'c and Jang to the nonisentropic problem. For large $\gamma$ values, inspired by Shkoller and Sideris, we use time-dependent weights that allow some of the top-order norms to potentially grow as the time variable tends to infinity. We also exploit coercivity estimates here via the fundamental theorem of calculus in time variable for norms which are not top-order.


title: On two-to-one mappings over finite fields

abstract: Two-to-one ($2$-to-$1$) mappings over finite fields play an important role in symmetric cryptography. In particular they allow to design APN functions, bent functions and semi-bent functions. In this paper we provide a systematic study of two-to-one mappings that are defined over finite fields. We characterize such mappings by means of the Walsh transforms. We also present several constructions, including an AGW-like criterion, constructions with the form of $x^rh(x^{(q-1)/d})$, those from permutation polynomials, from linear translators and from APN functions. Then we present $2$-to-$1$ polynomial mappings in classical classes of polynomials: linearized polynomials and monomials, low degree polynomials, Dickson polynomials and Muller-Cohen-Matthews polynomials, etc. Lastly, we show applications of $2$-to-$1$ mappings over finite fields for constructions of bent Boolean and vectorial bent functions, semi-bent functions, planar functions and permutation polynomials. In all those respects, we shall review what is known and provide several new results.


title: Spinon Confinement and a Sharp Longitudinal Mode in Yb$_2$Pt$_2$Pb in Magnetic Fields

abstract: The fundamental excitations in an antiferromagnetic chain of spins-1/2 are spinons, de-confined fractional quasiparticles that when combined in pairs, form a triplet excitation continuum. In an Ising-like spin chain the continuum is gapped and the ground state is N{\'e}el ordered. Here, we report high resolution neutron scattering experiments, which reveal how a magnetic field closes this gap and drives the spin chains in \YPP\ to a critical, disordered Luttinger-liquid state. In \YPP\ the effective spins-1/2 describe the dynamics of large, Ising-like Yb magnetic moments, ensuring that the measured excitations are exclusively longitudinal, which we find to be well described by time-dependent density matrix renormalization group calculations. The inter-chain coupling leads to the confinement of spinons, a condensed matter analog of quark confinement in quantum chromodynamics. Insensitive to transverse fluctuations, our measurements show how a gapless, dispersive longitudinal mode arises from confinement and evolves with magnetic order.


title: Augmenting and Tuning Knowledge Graph Embeddings

abstract: Knowledge graph embeddings rank among the most successful methods for link prediction in knowledge graphs, i.e., the task of completing an incomplete collection of relational facts. A downside of these models is their strong sensitivity to model hyperparameters, in particular regularizers, which have to be extensively tuned to reach good performance [Kadlec et al., 2017]. We propose an efficient method for large scale hyperparameter tuning by interpreting these models in a probabilistic framework. After a model augmentation that introduces per-entity hyperparameters, we use a variational expectation-maximization approach to tune thousands of such hyperparameters with minimal additional cost. Our approach is agnostic to details of the model and results in a new state of the art in link prediction on standard benchmark data.


title: Development and high-power testing of an X-band dielectric-loaded power extractor

abstract: Dielectric loaded structures are promising candidates for use in the structure wakefield acceleration (SWFA) technique, for both the collinear wakefield and the two-beam acceleration (CWA and TBA respectively) approaches, due to their low fabrication cost, low rf losses, and the potential to withstand high gradient. A short pulse (<=20 ns) TBA program is under development at the Argonne Wakefield Accelerator (AWA) facility where dielectric loaded structures are being used for both the power extractor/transfer structure (PETS) and the accelerator. In this study, an X-band 11.7 GHz dielectric PETS was developed and tested at the AWA facility to demonstrate high power wakefield generation. The PETS was driven by a train of eight electron bunches separated by 769.2 ps (9 times of the X-band rf period) in order to achieve coherent wakefield superposition. A total train charge of 360 nC was passed through the PETS structure to generate ~200 MW, ~3 ns flat-top rf pulses without rf breakdown. A future experiment is being planned to increase the generated rf power to approximately ~1 GW by optimizing the structure design and improving the drive beam quality.


title: A Semi-Supervised Self-Organizing Map for Clustering and Classification

abstract: There has been an increasing interest in semi-supervised learning in the recent years because of the great number of datasets with a large number of unlabeled data but only a few labeled samples. Semi-supervised learning algorithms can work with both types of data, combining them to obtain better performance for both clustering and classification. Also, these datasets commonly have a high number of dimensions. This article presents a new semi-supervised method based on self-organizing maps (SOMs) for clustering and classification, called Semi-Supervised Self-Organizing Map (SS-SOM). The method can dynamically switch between supervised and unsupervised learning during the training according to the availability of the class labels for each pattern. Our results show that the SS-SOM outperforms other semi-supervised methods in conditions in which there is a low amount of labeled samples, also achieving good results when all samples are labeled.


title: Online Charge Scheduling for Electric Vehicles in Autonomous Mobility on Demand Fleets

abstract: In this paper, we study an online charge scheduling strategy for fleets of autonomous-mobility-on-demand electric vechicles (AMoD EVs). We consider the case where vehicles complete trips and then enter a between-ride state throughout the day, with their information becoming available to the fleet operator in an online fashion. In the between-ride state, the vehicles must be scheduled for charging and then routed to their next passenger pick-up locations. Additionally, due to the unknown daily sequences of ride requests, the problem cannot be solved by any offline approach. As such, we study an online welfare maximization heuristic based on primal-dual methods that allocates limited fleet charging resources and rebalances the vehicles while avoiding congestion at charging facilities and pick-up locations. We discuss a competitive ratio result comparing the performance of our online solution to the clairvoyant offline solution and provide numerical results highlighting the performance of our heuristic.


title: $\omega$-Lyndon words

abstract: Let $\A$ be a finite non-empty set and $\preceq $ a total order on $\A^\nats$ verifying the following lexicographic like condition: For each $n\in \nats$ and $u, v\in \A^n,$ if $u^\omega \prec v^\omega$ then $ux\prec vy$ for all $x, y \in \A^\nats.$ A word $x\in \A^\nats$ is called $\omega$-Lyndon if $x\prec y$ for each proper suffix $y$ of $x.$ A finite word $w\in \A^+$ is called $\omega$-Lyndon if $w^\omega \prec v^\omega$ for each proper suffix $v$ of $w.$ In this note we prove that every infinite word may be written uniquely as a non-increasing product of $\omega$-Lyndon words.


title: On the generation of rank 3 simple matroids with an application to Terao's freeness conjecture

abstract: In this paper we describe a parallel algorithm for generating all non-isomorphic rank $3$ simple matroids with a given multiplicity vector. We apply our implementation in the HPC version of GAP to generate all rank $3$ simple matroids with at most $14$ atoms and a splitting characteristic polynomial. We have stored the resulting matroids alongside with various useful invariants in a publicly available, ArangoDB-powered database. As a byproduct we show that the smallest divisionally free rank $3$ arrangement which is not inductively free has $14$ hyperplanes and exists in all characteristics distinct from $2$ and $5$. Another database query proves that Terao's freeness conjecture is true for rank $3$ arrangements with $14$ hyperplanes in any characteristic.


title: Ion survival in grazing collisions of H$^-$ with vicinal nanosurfaces probes subband electronic structures

abstract: We study the electron dynamics at a monocrystalline Pd(111) surface with stepped vicinal nanostructures modeled in a simple Kronig-Penney scheme. The unoccupied bands of the surface are resonantly excited \textit{via} the resonant charge transfer (RCT) interaction of the surface with a hydrogen anion reflected at grazing angles. The interaction dynamics is simulated numerically in a quantum mechanical wave packet propagation approach. Visualization of the wave packet density shows that, when the electron is transferred to the metal, the surface and image subband states are the most likely locations of the electron as it evolves through the superlattice. The survival probability of the interacting ion exhibits strong modulations as a function of the vicinal-terrace size and shows peaks at those energies that access the image state subband dispersions. A simple square well model producing standing waves between the steps on the surface suggests the application of such ion-scattering at shallow angles to map electronic substructures in vicinal surfaces. The work also serves as the first proof-of-principle in the utility of our computational method to address, via RCT, surfaces with nanometric patterns.


title: Simulation smoothing for nowcasting with large mixed-frequency VARs

abstract: There is currently an increasing interest in large vector autoregressive (VAR) models. VARs are popular tools for macroeconomic forecasting and use of larger models has been demonstrated to often improve the forecasting ability compared to more traditional small-scale models. Mixed-frequency VARs deal with data sampled at different frequencies while remaining within the realms of VARs. Estimation of mixed-frequency VARs makes use of simulation smoothing, but using the standard procedure these models quickly become prohibitive in nowcasting situations as the size of the model grows. We propose two algorithms that alleviate the computational efficiency of the simulation smoothing algorithm. Our preferred choice is an adaptive algorithm, which augments the state vector as necessary to sample also monthly variables that are missing at the end of the sample. For large VARs, we find considerable improvements in speed using our adaptive algorithm. The algorithm therefore provides a crucial building block for bringing the mixed-frequency VARs to the high-dimensional regime.


title: The Polynomial Complexity of Vector Addition Systems with States

abstract: Vector addition systems are an important model in theoretical computer science and have been used in a variety of areas. In this paper, we consider vector addition systems with states over a parameterized initial configuration. For these systems, we are interested in the standard notion of computational complexity, i.e., we want to understand the length of the longest trace for a fixed vector addition system with states depending on the size of the initial configuration. We show that the asymptotic complexity of a given vector addition system with states is either $\Theta(N^k)$ for some computable integer $k$, where $N$ is the size of the initial configuration, or at least exponential. We further show that $k$ can be computed in polynomial time in the size of the considered vector addition system. Finally, we show that $1 \le k \le 2^n$, where $n$ is the dimension of the considered vector addition system.


title: 5G NR CA-Polar Maximum Likelihood Decoding by GRAND

abstract: CA-Polar codes have been selected for all control channel communications in 5G NR, but accurate, computationally feasible decoders are still subject to development. Here we report the performance of a recently proposed class of optimally precise Maximum Likelihood (ML) decoders, GRAND, that can be used with any block-code. As published theoretical results indicate that GRAND is computationally efficient for short-length, high-rate codes and 5G CA-Polar codes are in that class, here we consider GRAND's utility for decoding them. Simulation results indicate that decoding of 5G CA-Polar codes by GRAND, and a simple soft detection variant, is a practical possibility.


title: Quantization in Compressive Sensing: A Signal Processing Approach

abstract: Influence of the finite-length registers and quantization effects on the reconstruction of sparse and approximately sparse signals is analyzed in this paper. For the nonquantized measurements, the compressive sensing (CS) framework provides highly accurate reconstruction algorithms that produce negligible errors when the reconstruction conditions are met. However, hardware implementations of signal processing algorithms involve the finite-length registers and quantization of the measurements. An analysis of the effects related to the measurements quantization with an arbitrary number of bits is the topic of this paper. A unified mathematical model for the analysis of the quantization noise and the signal nonsparsity on the CS reconstruction is presented. An exact formula for the expected energy of error in the CS-based reconstructed signal is derived. The theory is validated through various numerical examples with quantized measurements, including the cases of approximately sparse signals, noise folding, and floating-point arithmetics.


title: Control of dopant crystallinity in electrochemically treated cuprate thin films

abstract: We present a methodology based on \textit{ex-situ} (post-growth) electrochemistry to control the oxygen concentration in thin films of the superconducting oxide La$_2$CuO$_{4+y}$ grown epitaxially on substrates of isostructural LaSrAlO$_4$. The superconducting transition temperature, which depends on the oxygen concentration, can be tuned by adjusting the pH level of the base solution used for the electrochemical reaction. As our main finding, we demonstrate that the dopant oxygens can either occupy the interstitial layer in an orientationally disordered state or organize into a crystalline phase via a mechanism in which dopant oxygens are inserted into the substrate, changing the lattice symmetry of both the substrate and the epitaxial film. We discuss this mechanism, and describe the resulting methodology as a platform to be explored in thin films of other transition metal oxides.


title: The Number of Gr\"obner Bases in Finite Fields

abstract: In the field of algebraic systems biology, the number of minimal polynomial models constructed using discretized data from an underlying system is related to the number of distinct reduced Gr\"obner bases for the ideal of the data points. While the theory of Gr\"obner bases is extensive, what is missing is a closed form for their number for a given ideal. This work contributes connections between the geometry of data points and the number of Gr\"obner bases associated to small data sets. Furthermore we improve an existing upper bound for the number of Gr\"obner bases specialized for data over a finite field.


title: Private Authentication with Physical Identifiers Through Broadcast Channel Measurements

abstract: A basic model for key agreement with biometric or physical identifiers is extended to include measurements of a hidden source through a general broadcast channel (BC). An inner bound for strong secrecy, maximum key rate, and minimum privacy-leakage and database-storage rates is proposed. The inner bound is shown to be tight for physically-degraded and less-noisy BCs.

categories: cs.IT cs.CR cs.MM eess.SP math.IT math.PR


title: BAT AGN Spectroscopic Survey: XVI. General Physical Characteristics of BAT Blazars

abstract: The recently released 105-month {\it Swift}-Burst Alert Telescope (BAT) all-sky hard X-ray survey catalog presents an opportunity to study astrophysical objects detected in the deepest look at the entire hard X-ray (14$-$195 keV) sky. Here we report the results of a multifrequency study of 146 blazars from this catalog, quadrupling the number compared to past studies, by utilizing recent data from the {\it Fermi}-Large Area Telescope (LAT), Swift-BAT, and archival measurements. In our $\gamma$-ray analysis of $\sim$10 years of the LAT data, 101 are found as $\gamma$-ray emitters, whereas, 45 remains LAT undetected. We model the broadband spectral energy distributions with a synchrotron-inverse Compton radiative model. On average, BAT detected sources host massive black holes ($M_{\rm bh}\sim10^9$ M$_{\odot}$) and luminous accretion disks ($L_{\rm d}\sim10^{46}$ erg s$^{-1}$). At high-redshifts ($z>2$), BAT blazars host more powerful jets with luminous accretion disks compared to those detected only with the {\it Fermi}-LAT. We find good agreement in the black hole masses derived from the single-epoch optical spectroscopic measurements and standard accretion disk modeling approaches. Other physical properties of BAT blazars are similar to those known for {\it Fermi}-LAT detected objects.


title: The independent set problem is FPT for even-hole-free graphs

abstract: The class of even-hole-free graphs is very similar to the class of perfect graphs, and was indeed a cornerstone in the tools leading to the proof of the Strong Perfect Graph Theorem. However, the complexity of computing a maximum independent set (MIS) is a long-standing open question in even-hole-free graphs. From the hardness point of view, MIS is W[1]-hard in the class of graphs without induced 4-cycle (when parameterized by the solution size). Halfway of these, we show in this paper that MIS is FPT when parameterized by the solution size in the class of even-hole-free graphs. The main idea is to apply twice the well-known technique of augmenting graphs to extend some initial independent set.


title: Regularity of linear and polynomial images of Skorohod differentiable measures

abstract: In this paper we study the regularity properties of linear and polynomial images of Skorohod differentiable measures. Firstly, we obtain estimates for the Skorohod derivative norm of a projection of a product of Scorohod differentiable measures. In the second part of the paper we prove Nikolskii--Besov regularity of a polynomial image of a Skorohod differentiable measure on $\mathbb{R}^n$.


title: Multiview Aggregation for Learning Category-Specific Shape Reconstruction

abstract: We investigate the problem of learning category-specific 3D shape reconstruction from a variable number of RGB views of previously unobserved object instances. Most approaches for multiview shape reconstruction operate on sparse shape representations, or assume a fixed number of views. We present a method that can estimate dense 3D shape, and aggregate shape across multiple and varying number of input views. Given a single input view of an object instance, we propose a representation that encodes the dense shape of the visible object surface as well as the surface behind line of sight occluded by the visible surface. When multiple input views are available, the shape representation is designed to be aggregated into a single 3D shape using an inexpensive union operation. We train a 2D CNN to learn to predict this representation from a variable number of views (1 or more). We further aggregate multiview information by using permutation equivariant layers that promote order-agnostic view information exchange at the feature level. Experiments show that our approach is able to produce dense 3D reconstructions of objects that improve in quality as more views are added.


title: A Semi-Supervised Self-Organizing Map with Adaptive Local Thresholds

abstract: In the recent years, there is a growing interest in semi-supervised learning, since, in many learning tasks, there is a plentiful supply of unlabeled data, but insufficient labeled ones. Hence, Semi-Supervised learning models can benefit from both types of data to improve the obtained performance. Also, it is important to develop methods that are easy to parameterize in a way that is robust to the different characteristics of the data at hand. This article presents a new method based on Self-Organizing Map (SOM) for clustering and classification, called Adaptive Local Thresholds Semi-Supervised Self-Organizing Map (ALTSS-SOM). It can dynamically switch between two forms of learning at training time, according to the availability of labels, as in previous models, and can automatically adjust itself to the local variance observed in each data cluster. The results show that the ALTSS-SOM surpass the performance of other semi-supervised methods in terms of classification, and other pure clustering methods when there are no labels available, being also less sensitive than previous methods to the parameters values.


title: On simple $Z_2$-invariant and corner function germs

abstract: V.I.Arnold has classified simple (i.e. having no modules for the classification) singularities (function germs), and also simple boundary singularities (function germs invariant with respect to the action $\sigma(x_1; y_1, \ldots, y_n)=(-x_1; y_1, \ldots, y_n)$ of the group $Z_2$. In particular, it was shown that a function germ (respectively a boundary singularity germ) is simple if and only if the intersection form (respectively the restriction of the intersection form to the subspace to anti-invariant cycles) of a germ in $3+4s$ variables stable equivalent to the one under consideration is negative definite and if and only if the (equivariant) monodromy group on the corresponding subspace is finite. We formulate and prove analogues of these statements for function germs invariant with respect to an arbitrary action of the group $Z_2$ and also for corner singularities.


title: The importance of thermal gradients on the vortex dynamics and magnetic behavior of mesoscopic superconducting samples

abstract: Usually, the measurements of electronic and magnetic properties of superconducting samples are carried out under a constant temperature bath. On the other hand, thermal gradients induce local variation of the superconducting order parameter, and the vortex dynamics can present interesting behaviors. In this work, we solved the time-dependent Ginzburg-Landau equations simulating samples under two different thermal gradients, and considering two values of the Ginzburg-Landau parameter, \k{appa}. We find out that both parameters, i.e., \k{appa} and thermal gradients, play an important role on the vortex dynamics and on the magnetization behavior of the samples.


title: Robust degenerate unfoldings of cycles and tangencies

abstract: We construct open sets of degenerate unfoldings of heterodimensional cycles of any co-index $c>0$ and homoclinic tangencies of arbitrary codimension $c>0$. These sets are known to be the support of unexpected phenomena in families of diffeomorphisms, such as the Kolmogorov typical co-existence of infinitely many attractors. As a prerequisite we also construct robust homoclinic tangencies of large codimension which cannot be inside a strong partially hyperbolic set.


title: A Genetic Algorithm for Astroparticle Physics Studies

abstract: Precision measurements of charged cosmic rays have recently been carried out by space-born (e.g. AMS-02), or ground experiments (e.g. HESS). These measured data are important for the studies of astro-physical phenomena, including supernova remnants, cosmic ray propagation, solar physics and dark matter. Those scenarios usually contain a number of free parameters that need to be adjusted by observed data. Some techniques, such as Markov Chain Monte Carlo and MultiNest, are developed in order to solve the above problem. However, it is usually required a computing farm to apply those tools. In this paper, a genetic algorithm for finding the optimum parameters for cosmic ray injection and propagation is presented. We find that this algorithm gives us the same best fit results as the Markov Chain Monte Carlo but consuming less computing power by nearly 2 orders of magnitudes.


title: Equivariant instanton homology

abstract: We define four versions of equivariant instanton Floer homology ($I^+, I^-, I^\infty$ and $\widetilde I$) for a class of 3-manifolds and $SO(3)$-bundles over them including all rational homology spheres. These versions are analogous to the four flavors of monopole and Heegaard Floer homology theories. This construction is functorial for a large class of 4-manifold cobordisms, and agrees with Donaldson's definition of equivariant instanton homology for integer homology spheres. Furthermore, one of our invariants is isomorphic to Floer's instanton homology for admissible bundles, and we calculate $I^\infty$ in all cases it is defined, away from characteristic 2. The appendix, possibly of independent interest, defines an algebraic construction of three equivariant homology theories for dg-modules over a dg-algebra, the equivariant homology $H^+(A,M)$, the coBorel homology $H^-(A,M)$, and the Tate homology $H^\infty(A,M)$. The constructions of the appendix are used to define our invariants.


title: Multiscale dynamics of a heterotypic cancer cell population within a fibrous extracellular matrix

abstract: Local cancer cell invasion is a complex process involving many cellular and tissue interactions and is an important prerequisite for metastatic spread, the main cause of cancer related deaths. Occurring over many different temporal and spatial scales, the first stage of local invasion is the secretion of matrix-degrading enzymes (MDEs) and the resulting degradation of the extra-cellular matrix (ECM). This process creates space in which the cells can invade and thus enlarge the tumour. As a tumour increases in malignancy, the cancer cells adopt the ability to mutate into secondary cell subpopulations giving rise to a heterogeneous tumour. This new cell subpopulation often carries higher invasive qualities and permits a quicker spread of the tumour. Building upon the recent multiscale modelling framework for cancer invasion within a fibrous ECM introduced in Shuttleworth and Trucu (2019), in this paper we consider the process of local invasion by a heterotypic tumour consisting of two cancer cell populations mixed with a two-phase ECM. To that end, we address the double feedback link between the tissue-scale cancer dynamics and the cell-scale molecular processes through the development of a two-part modelling framework that crucially incorporates the multiscale dynamic redistribution of oriented fibres occurring within a two-phase extra-cellular matrix and combines this with the multiscale leading edge dynamics exploring key matrix-degrading enzymes molecular processes along the tumour interface that drive the movement of the cancer boundary. The modelling framework will be accompanied by computational results that explore the effects of the underlying fibre network on the overall pattern of cancer invasion.


title: Electroweak pion-production on nuclei within the extended factorization scheme

abstract: We have applied the extended factorization scheme to investigate the electroweak pion production on nuclei. The ANL-Osaka model, which was obtained by analyzing the data of $\pi N$, $\gamma N$, $N(e,e'\pi)$ and $N(\nu,\mu\,\pi) $ reactions up to invariant mass $W=$ 2 GeV, is used to generate the matrix elements of current operators relevant to pion-production off the nucleon. Medium effects on the $\Delta$ (1232) component of meson-exchange current are included by using a $\Delta$-nucleus potential determined from the previous $\Delta$-hole model studies of pion-nucleus reactions. Nuclear correlations in the initial target state and in the spectator system(s) are modeled using realistic hole spectral functions. As a first step, we show that the data of $^{12}$C$(e,e')$ up to the $\Delta$ (1232) region can be described reasonably well. The interplay between the pion production and two-body meson-exchange mechanisms is shown to be essential in improving the agreement with the data in the ``dip'' region, between the quasielastic and the $\Delta$ (1232) peaks. Predictions for $^{12}$C$(\nu,\mu\,\pi)$ have also been made. They can be used to estimate pion-emission rates in neutrino-nucleus cross section, which constitutes an important systematic uncertainty to the reconstructed neutrino energy. With further improvements of the Metropolis Monte-Carlo techniques to account for final states comprised of more than two particles, our approach can be employed up to $W=$ 2 GeV, where two-pion production and higher mass nucleon resonances must be included for analyzing the data from accelerator-based neutrino-oscillation experiments.


title: A multiresolution algorithm to generate images of generalized fuzzy fractal attractors

abstract: We provide a new algorithm to generate images of the generalized fuzzy fractal attractors described in Oliveira-2017. We also provide some important results on the approximation of fractal operators to discrete subspaces with application to discrete versions of the deterministic algorithm for fractal image generation in the cases of IFS recovering the classical images from Barnsley et al., Fuzzy IFS from Cabrelli-1992 and GIFS's from Jaros-2016.


title: Advanced Cauchy Mutation for Differential Evolution in Numerical Optimization

abstract: Among many evolutionary algorithms, differential evolution (DE) has received much attention over the last two decades. DE is a simple yet powerful evolutionary algorithm that has been used successfully to optimize various real-world problems. Since it was introduced, many researchers have developed new methods for DE, and one of them makes use of a mutation based on the Cauchy distribution to increase the convergence speed of DE. The method monitors the results of each individual in the selection operator and performs the Cauchy mutation on consecutively failed individuals, which generates mutant vectors by perturbing the best individual with the Cauchy distribution. Therefore, the method can locate the consecutively failed individuals to new positions close to the best individual. Although this approach is interesting, it fails to take into account establishing a balance between exploration and exploitation. In this paper, we propose a sigmoid based parameter control that alters the failure threshold for performing the Cauchy mutation in a time-varying schedule, which can establish a good ratio between exploration and exploitation. Experiments and comparisons have been done with six conventional and six advanced DE variants on a set of 30 benchmark problems, which indicate that the DE variants assisted by the proposed algorithm are highly competitive, especially for multimodal functions.


title: The Pinpoint Comets: 133P/Elst-Pizarro, 249P/LINEAR, 331P/Gibbs, 62412 and 6478 Gault

abstract: From two Active Asteroid (AA) known in 1979 we have advanced to 37 members at the beginning of 2019. More surprisingly, in the first three months of 2019 five new members were added to the list, one of them, 6478 Gault, with a curious, out of the ordinary tail. This prompted a visual search for similar objects in images depositories. We have identified five cometary objects that have the same morphological appearance: A thin long tail coming out of a star-like nucleus with no gas. The members of this new class (Pinpoint Comet Group, PCG) are, 133P/Elst-Pizarro, 249P/LINEAR, 331P/Gibbs, 62412 and 6478 Gault. We look for differences and similarities among them. Four of them belong to the main belt and one to the Jupiter Family of comets. One of them has left the group but will return soon. Two of them are double. All exhibit densities well above the Jupiter Family of comets values. This tells us something about their origin and evolution. In this work we study the secular and rotational light curves of these objects and their phase and SLC plots. We will investigate the meaning of these characteristics in terms of their physical evolution.


title: Variable thrust and high efficiency propulsion with oscillating foils at high Reynolds numbers

abstract: Bio-inspired oscillatory foil propulsion has the ability to traverse various propulsive modes by dynamically changing the foil's heave and pitch kinematics. This research characterizes the propulsion properties and wake dynamics of a symmetric oscillating foil, specifically targeting the high Reynolds number operation of small to medium surface vessels whose propulsive specifications have a broad range of loads and speeds. An unsteady Reynolds-averaged Navier-Stokes (URANS) solver with a k-$\omega$ SST turbulence model is used to sweep through pitch amplitude and frequency at two heave amplitudes of $h_0/c=1$ and $h_0/c=2$ at $Re=10^6$. At $h_0/c=2$, the maximum thrust coefficient is $C_T=8.2$ due to the large intercepted flow area of the foil, whereas at a decreased Strouhal number the thrust coefficient decreases and the maximum propulsive efficiency reaches 75%. Results illustrate the kinematics required to transition between the high-efficiency and high-thrust regimes at high Reynolds number and the resulting changes to the vortex wake structure. The unsteady vortex dynamics throughout the heave-pitch cycle strongly influence the characterization of thrust and propulsive efficiency, and are classified into flow regimes based on performance and vortex structure.