title: Surprising identities for the greedy independent set on Cayley trees

abstract: We prove a surprising symmetry between the law of the size $G_n$ of the greedy independent set on a uniform Cayley tree $ \mathcal{T}_n$ of size $n$ and that of its complement. We show that $G_n$ has the same law as the number of vertices at even height in $ \mathcal{T}_n$ rooted at a uniform vertex. This enables us to compute the exact law of the $G_n$. We also give a Markovian construction of the greedy independent set, which highlights the symmetry of $G_n$ and whose proof uses a new Markovian exploration of rooted Cayley trees which is of independent interest.


title: Error-Correction for Sparse Support Recovery Algorithms

abstract: Consider the compressed sensing setup where the support $s^*$ of an $m$-sparse $d$-dimensional signal $x$ is to be recovered from $n$ linear measurements with a given algorithm. Suppose that the measurements are such that the algorithm does not guarantee perfect support recovery and that true features may be missed. Can they efficiently be retrieved? This paper addresses this question through a simple error-correction module referred to as LiRE. LiRE takes as input an estimate $s_{in}$ of the true support $s^*$, and outputs a refined support estimate $s_{out}$. In the noiseless measurement setup, sufficient conditions are established under which LiRE is guaranteed to recover the entire support, that is $s_{out}$ contains $s^*$. These conditions imply, for instance, that in the high-dimensional regime LiRE can correct a sublinear in $m$ number of errors made by Orthogonal Matching Pursuit (OMP). The computational complexity of LiRE is $O(mnd)$. Experimental results with random Gaussian design matrices show that LiRE substantially reduces the number of measurements needed for perfect support recovery via Compressive Sampling Matching Pursuit, Basis Pursuit (BP), and OMP. Interestingly, adding LiRE to OMP yields a support recovery procedure that is more accurate and significantly faster than BP. This observation carries over in the noisy measurement setup. Finally, as a standalone support recovery algorithm with a random initialization, experiments show that LiRE's reconstruction performance lies between OMP and BP. These results suggest that LiRE may be used generically, on top of any suboptimal baseline support recovery algorithm, to improve support recovery or to operate with a smaller number of measurements, at the cost of a relatively small computational overhead. Alternatively, LiRE may be used as a standalone support recovery algorithm that is competitive with respect to OMP.


title: The Response of Laser Interferometric Gravitational Wave Detectors Beyond the Eikonal Equation

abstract: The response of Michelson interferometers to weak plane gravitational waves is computed at one order of accuracy beyond the eikonal equation. The modulation of the electromagnetic field amplitude and polarisation are taken into account by solving the transport equations of geometrical optics with boundary conditions adapted to laser interferometry. Considering both DC and balanced homodyne readout schemes, explicit formulae for the interferometer output signals are derived. These signals comprise perturbations of the optical path length, frequency and amplitude, and are shown to be insensitive to polarisation perturbations.


title: Director Deformations, Geometric Frustration, and Modulated Phases in Liquid Crystals

abstract: This article analyzes modulated phases in liquid crystals, from the long-established cholesteric and blue phases to the recently discovered twist-bend, splay-bend, and splay nematic phases, as well as the twist-grain-boundary (TGB) and helical nanofilament variations on smectic phases. The analysis uses the concept of four fundamental modes of director deformation: twist, bend, splay, and a fourth mode related to saddle-splay. Each mode is coupled to a specific type of molecular order: chirality, polarization perpendicular and parallel to the director, and octupolar order. When the liquid crystal develops one type of spontaneous order, the ideal local structure becomes nonuniform, with the corresponding director deformation. In general, the ideal local structure is frustrated; it cannot fill space. As a result, the liquid crystal must form a complex global phase, which may have a combination of deformation modes, and may have a periodic array of defects. Thus, the concept of an ideal local structure under geometric frustration provides a unified framework to understand the wide variety of modulated phases.


title: Quantum algorithm for the Navier Stokes equations by using the streamfunction vorticity formulation and the lattice Boltzmann method

abstract: A new algorithm for solving the Navier-Stokes equations (NSE) on a quantum device is presented. For the fluid flow equations the stream function-vorticity formulation is adopted, while the lattice Boltzmann method (LBM) is utilized for solving the corresponding system of equations numerically for one time step. Following the nature of the lattice Boltzmann method, the proposed quantum algorithm consists of five major sections: initialization, collision, propagation, boundary condition implementation, and calculation of macroscopic quantities. The collision and boundary condition step is quantumly implemented by applying the standard-form encoding approach, while the quantum walk procedure is applied for the propagation step. The algorithm is implemented by using IBM's quantum computing software development framework Qiskit, while for the verification purposes two-dimensional (2D) cavity flow is simulated and compared with classical code.


title: Topological Linear System Identification via Moderate Deviations Theory

abstract: Two dynamical systems are topologically equivalent when their phase-portraits can be morphed into each other by a homeomorphic coordinate transformation on the state space. The induced equivalence classes capture qualitative properties such as stability or the oscillatory nature of the state trajectories, for example. In this paper we develop a method to learn the topological class of an unknown stable system from a single trajectory of finitely many state observations. Using a moderate deviations principle for the least squares estimator of the unknown system matrix $\theta$, we prove that the probability of misclassification decays exponentially with the number of observations at a rate that is proportional to the square of the smallest singular value of $\theta$.


title: MalBERT: Using Transformers for Cybersecurity and Malicious Software Detection

abstract: In recent years we have witnessed an increase in cyber threats and malicious software attacks on different platforms with important consequences to persons and businesses. It has become critical to find automated machine learning techniques to proactively defend against malware. Transformers, a category of attention-based deep learning techniques, have recently shown impressive results in solving different tasks mainly related to the field of Natural Language Processing (NLP). In this paper, we propose the use of a Transformers' architecture to automatically detect malicious software. We propose a model based on BERT (Bidirectional Encoder Representations from Transformers) which performs a static analysis on the source code of Android applications using preprocessed features to characterize existing malware and classify it into different representative malware categories. The obtained results are promising and show the high performance obtained by Transformer-based models for malicious software detection.


title: Discovery of the propargyl radical (CH2CCH) in TMC-1: one of the most abundant radicals ever found and a key species for cyclization to benzene in cold dark clouds

abstract: We present the first identification in interstellar space of the propargyl radical (CH2CCH). This species was observed in the cold dark cloud TMC-1 using the Yebes 40m telescope. The six strongest hyperfine components of the 2,0,2-1,0,1 rotational transition, lying at 37.46~GHz, were detected with signal-to-noise ratios in the range 4.6-12.3 sigma. We derive a column density of 8.7e13 cm-2 for CH2CCH, which translates to a fractional abundance relative to H2 of 8.7e-9. This radical has a similar abundance to methyl acetylene, with an abundance ratio CH2CCH/CH3CCH close to one. The propargyl radical is thus one of the most abundant radicals detected in TMC-1, and it is probably the most abundant organic radical with a certain chemical complexity ever found in a cold dark cloud. We constructed a gas-phase chemical model and find calculated abundances that agree with, or fall two orders of magnitude below, the observed value depending on the poorly constrained low-temperature reactivity of CH2CCH with neutral atoms. According to the chemical model, the propargyl radical is essentially formed by the C + C2H4 reaction and by the dissociative recombination of C3Hn+ ions with n = 4-6. The propargyl radical is believed to control the synthesis of the first aromatic ring in combustion processes, and it probably plays a key role in the synthesis of large organic molecules and cyclization processes to benzene in cold dark clouds.


title: Two-step reinforcement learning for model-free redesign of nonlinear optimal regulator

abstract: In many practical control applications, the performance level of a closed-loop system degrades over time due to the change of plant characteristics. Thus, there is a strong need for redesigning a controller without going through the system modeling process, which is often difficult for closed-loop systems. Reinforcement learning (RL) is one of the promising approaches that enable model-free redesign of optimal controllers for nonlinear dynamical systems based only on the measurement of the closed-loop system. However, the learning process of RL usually requires a considerable number of trial-and-error experiments using the poorly controlled system that may accumulate wear on the plant. To overcome this limitation, we propose a model-free two-step design approach that improves the transient learning performance of RL in an optimal regulator redesign problem for unknown nonlinear systems. Specifically, we first design a linear control law that attains some degree of control performance in a model-free manner, and then, train the nonlinear optimal control law with online RL by using the designed linear control law in parallel. We introduce an offline RL algorithm for the design of the linear control law and theoretically guarantee its convergence to the LQR controller under mild assumptions. Numerical simulations show that the proposed approach improves the transient learning performance and efficiency in hyperparameter tuning of RL.


title: PalmTree: Learning an Assembly Language Model for Instruction Embedding

abstract: Deep learning has demonstrated its strengths in numerous binary analysis tasks, including function boundary detection, binary code search, function prototype inference, value set analysis, etc. When applying deep learning to binary analysis tasks, we need to decide what input should be fed into the neural network model. More specifically, we need to answer how to represent an instruction in a fixed-length vector. The idea of automatically learning instruction representations is intriguing, however the existing schemes fail to capture the unique characteristics of disassembly. These schemes ignore the complex intra-instruction structures and mainly rely on control flow in which the contextual information is noisy and can be influenced by compiler optimizations. In this paper, we propose to pre-train an assembly language model called PalmTree for generating general-purpose instruction embeddings by conducting self-supervised training on large-scale unlabeled binary corpora. PalmTree utilizes three pre-training tasks to capture various characteristics of assembly language. These training tasks overcome the problems in existing schemes, thus can help to generate high-quality representations. We conduct both intrinsic and extrinsic evaluations, and compare PalmTree with other instruction embedding schemes. PalmTree has the best performance for intrinsic metrics, and outperforms the other instruction embedding schemes for all downstream tasks.


title: Quantum study of reaction O(3P) + H2(v,j) $\rightarrow$ OH + H: OH formation in strongly UV-irradiated gas

abstract: The reaction between atomic oxygen and molecular hydrogen is an important one in astrochemistry as it regulates the abundance of the hydroxyl radical and serves to open the chemistry of oxygen in diverse astronomical environments. However, the existence of a high activation barrier in the reaction with ground state oxygen atoms limits its efficiency in cold gas. In this study we calculate the dependence of the reaction rate coefficient on the rotational and vibrational state of H$_2$ and evaluate the impact on the abundance of OH in interstellar regions strongly irradiated by far-UV photons, where H2 can be efficiently pumped to excited vibrational states. We use a recently calculated potential energy surface and carry out time-independent quantum mechanical scattering calculations to compute rate coefficients for the reaction O(3P) + H2(v,j) -> OH + H, with H2 in vibrational states v = 0-7 and rotational states j = 0-10. We find that the reaction becomes significantly faster with increasing vibrational quantum number of H2, although even for high vibrational states of H2 (v = 4-5) for which the reaction is barrierless, the rate coefficient does not strictly attain the collision limit and still maintains a positive dependence with temperature. We implemented the calculated state-specific rate coefficients in the Meudon PDR code to model the Orion Bar PDR and evaluate the impact on the abundance of the OH radical. We find the fractional abundance of OH is enhanced by up to one order of magnitude in regions of the cloud corresponding to Av = 1.3-2.3, compared to the use of a thermal rate coefficient for O + H2, although the impact on the column density of OH is modest, of about 60 %. The calculated rate coefficients will be useful to model and interpret JWST observations of OH in strongly UV-illuminated environments.


title: Low Stress Ion Conductance Microscopy of Sub-Cellular Stiffness

abstract: Directly examining subcellular mechanics whilst avoiding excessive strain of a live cell requires the precise control of light stress on very small areas, which is fundamentally difficult. Here we use a glass nanopipet out of contact with the plasma membrane to both exert the stress on the cell and also accurately monitor cellular compression. This allows the mapping of cell stiffness at a lateral resolution finer than 100 nm. We calculate the stress a nanopipet exerts on a cell as the sum of the intrinsic pressure between the tip face and the plasma membrane plus its direct pressure on any glycocalyx, both evaluated from the gap size in terms of the ion current decrease. A survey of cell types confirms that an intracellular pressure of approximately 120 Pa begins to detach the plasma membrane from the cytoskeleton and reveals that the first 660 +/- 90 nm of compression of a neuron cell body is much softer than previous methods have been able to detect.


title: Absolute continuity of the solution to the stochastic Burgers equation

abstract: We prove the existence and the Besov regularity of the density of the solution to a general parabolic SPDE which includes the stochastic Burgers equation on an unbounded domain. We use an elementary approach based on the fractional integration by parts.


title: A lower bound on the spectral gap of Schr\"odinger operators with weak potentials of compact support

abstract: In this paper we continue the study of the spectral gap of Schr\"odinger operators on large intervals and subject to Neumann boundary conditions. The main goal is to derive a lower bound on the spectral gap which is polynomial in the interval length. This bound is derived for a class of bounded potentials of compact support which are weak enough in a suitable sense.


title: Rubbing and Drawing: Generic Ways to Improve the Thermoelectric Power Factor of Organic Semiconductors?

abstract: Highly oriented polymer films can show considerable anisotropy in the thermoelectric properties leading to power factors beyond those predicted by the widely obeyed power law linking the thermopower $S$ and the electrical conductivity $\sigma$ as $S\propto\sigma^{-1/4}$. This has led to encouraging practical results with respect to the electrical conductivity, notwithstanding that the conditions necessary to enhance $\sigma$ and $S$ simultaneously are less clear. Here, kinetic Monte Carlo simulations are used to study the impact of structural anisotropy on the thermoelectric properties of disordered organic semiconductors. We find that stretching is a suitable strategy to improve the conductivity along the direction of strain, while the effect on the power factor depends on the morphology the polymer crystallizes. In general, crystalline polymers show a simultaneous increase in $\sigma$ and $S$ which is not the case for amorphous polymers. Moreover, we show that the trends resulting from simulations based on variable-range hopping are in good agreement with experiments and can describe the different functional dependencies in the $S$ versus $\sigma$ behaviour of different directions.


title: Generalized Emergent Dark Energy Model and the Hubble Constant Tension

abstract: We investigate a generalized form of the phenomenologically emergent dark energy model, known as generalized emergent dark energy (GEDE), introduced by Li and Shafieloo [Astrophys. J. {\bf 902}, 58 (2020)] in light of a series of cosmological probes and considering the evolution of the model at the level of linear perturbations. This model introduces a free parameter $\Delta$ that can discriminate between the $\Lambda$CDM (corresponds to $\Delta=0$) or the phenomenologically emergent dark energy (PEDE) (corresponds to $\Delta=1$) models, allowing us to determine which model is preferred most by the fit of the observational datasets. We find evidence in favor of the GEDE model for Planck alone and in combination with R19, while the Bayesian model comparison is inconclusive when Supernovae Type Ia or BAO data are included. In particular, we find that $\Lambda$CDM model is disfavored at more than $2\sigma$ CL for most of the observational datasets considered in this work and PEDE is in agreement with Planck 2018+BAO+R19 combination within $1\sigma$ CL.


title: Bioluminescence modeling for deep sea experiments

abstract: We develop a modeling framework for bioluminescence light found in the deep sea near neutrino telescopes by combining a hydrodynamic model with a stochastic one. The bioluminescence is caused by organisms when exposed to a non-constant water flow, such as past the neutrino telescopes. We model the flow using the incompressible Navier-Stokes equations for Reynolds numbers between 4000 and 23000. The discretization relies on a finite element method which includes upwind-stabilization for the velocity field. On top of the flow model, we simulate a population of random microscopic organisms. Their movement and emission are stochastic processes which we model using Monte Carlo methods. We observe unique time-series for the photon counts depending on the flow velocity and detector specifications. This opens up the possibility of categorizing organisms using neutrino detectors. We show that the average light-yield and pulse shapes require precise flow modeling, while the emission timing is chaotic. From this we construct a fast modeling scheme, requiring only a subset of computationally expensive flow and population modeling.


title: Proactive and AoI-aware Failure Recovery for Stateful NFV-enabled Zero-Touch 6G Networks: Model-Free DRL Approach

abstract: In this paper, we propose a Zero-Touch, deep reinforcement learning (DRL)-based Proactive Failure Recovery framework called ZT-PFR for stateful network function virtualization (NFV)-enabled networks. To this end, we formulate a resource-efficient optimization problem minimizing the network cost function including resource cost and wrong decision penalty. As a solution, we propose state-of-the-art DRL-based methods such as soft-actor-critic (SAC) and proximal-policy-optimization (PPO). In addition, to train and test our DRL agents, we propose a novel impending-failure model. Moreover, to keep network status information at an acceptable freshness level for appropriate decision-making, we apply the concept of age of information to strike a balance between the event and scheduling based monitoring. Several key systems and DRL algorithm design insights for ZT-PFR are drawn from our analysis and simulation results. For example, we use a hybrid neural network, consisting long short-term memory layers in the DRL agents structure, to capture impending-failures time dependency.


title: Time-varying $\ell_0$ optimization for Spike Inference from Multi-Trial Calcium Recordings

abstract: Optical imaging of genetically encoded calcium indicators is a powerful tool to record the activity of a large number of neurons simultaneously over a long period of time from freely behaving animals. However, determining the exact time at which a neuron spikes and estimating the underlying firing rate from calcium fluorescence data remains challenging, especially for calcium imaging data obtained from a longitudinal study. We propose a multi-trial time-varying $\ell_0$ penalized method to jointly detect spikes and estimate firing rates by robustly integrating evolving neural dynamics across trials. Our simulation study shows that the proposed method performs well in both spike detection and firing rate estimation. We demonstrate the usefulness of our method on calcium fluorescence trace data from two studies, with the first study showing differential firing rate functions between two behaviors and the second study showing evolving firing rate function across trials due to learning.


title: Hybrid Point Cloud Semantic Compression for Automotive Sensors: A Performance Evaluation

abstract: In a fully autonomous driving framework, where vehicles operate without human intervention, information sharing plays a fundamental role. In this context, new network solutions have to be designed to handle the large volumes of data generated by the rich sensor suite of the cars in a reliable and efficient way. Among all the possible sensors, Light Detection and Ranging (LiDAR) can produce an accurate 3D point cloud representation of the surrounding environment, which in turn generates high data rates. For this reason, efficient point cloud compression is paramount to alleviate the burden of data transmission over bandwidth-constrained channels and to facilitate real-time communications. In this paper, we propose a pipeline to efficiently compress LiDAR observations in an automotive scenario. First, we leverage the capabilities of RangeNet++, a Deep Neural Network (DNN) used to semantically infer point labels, to reduce the channel load by selecting the most valuable environmental data to be disseminated. Second, we compress the selected points using Draco, a 3D compression algorithm which is able to obtain compression up to the quantization error. Our experiments, validated on the Semantic KITTI dataset, demonstrate that it is possible to compress and send the information at the frame rate of the LiDAR, thus achieving real-time performance.


title: AnswerQuest: A System for Generating Question-Answer Items from Multi-Paragraph Documents

abstract: One strategy for facilitating reading comprehension is to present information in a question-and-answer format. We demo a system that integrates the tasks of question answering (QA) and question generation (QG) in order to produce Q&A items that convey the content of multi-paragraph documents. We report some experiments for QA and QG that yield improvements on both tasks, and assess how they interact to produce a list of Q&A items for a text. The demo is accessible at qna.sdl.com.


title: Fast Interactive Video Object Segmentation with Graph Neural Networks

abstract: Pixelwise annotation of image sequences can be very tedious for humans. Interactive video object segmentation aims to utilize automatic methods to speed up the process and reduce the workload of the annotators. Most contemporary approaches rely on deep convolutional networks to collect and process information from human annotations throughout the video. However, such networks contain millions of parameters and need huge amounts of labeled training data to avoid overfitting. Beyond that, label propagation is usually executed as a series of frame-by-frame inference steps, which is difficult to be parallelized and is thus time consuming. In this paper we present a graph neural network based approach for tackling the problem of interactive video object segmentation. Our network operates on superpixel-graphs which allow us to reduce the dimensionality of the problem by several magnitudes. We show, that our network possessing only a few thousand parameters is able to achieve state-of-the-art performance, while inference remains fast and can be trained quickly with very little data.


title: Optical constants of a solar system organic analog and the Allende meteorite in the near and mid-infrared (1.5-13 {\mu}m)

abstract: Measurements of visible and near-infrared reflection (0.38-5 {\mu}m) and mid to far infrared emission (5-200 {\mu}m) from telescope and satellite remote sensing instruments make it possible to investigate the composition of planetary surfaces via electronic transitions and vibrational modes of chemical bonds. Red spectral slopes at visible and near infrared wavelengths and absorption features at 3.3 and 3.4 {\mu}m observed in circumstellar disks, the interstellar medium, and on the surfaces of solar-system bodies are interpreted to be due to the presence of organic material and other carbon compounds. Identifying the origin of these features requires measurements of the optical properties of a variety of relevant analog and planetary materials. Spectroscopic models of dust within circumstellar disks and the interstellar medium as well as planetary regoliths often incorporate just one such laboratory measurement despite the wide variation in absorption and extinction properties of organic and other carbon-bearing materials. Here we present laboratory measurements of transmission spectra in the 1.5-13 {\mu}m region and use these to derive real and imaginary indices of refraction for two samples: 1) an analog to meteoritic insoluble organic matter and 2) a powdered Allende meteorite sample. We also test our refractive index retrieval method on a previously published transmission spectrum of an Mg-rich olivine. We compare optical measurements of the insoluble organic-matter analog to those of other solar-system and extrasolar organic analogs, such as amorphous carbon and tholins, and find that the indices of refraction of the newly characterized material differ significantly from other carbonaceous samples.


title: Local models for the trianguline variety and partially classical families

abstract: We generalize Breuil-Hellmann-Schraen's local model for the trianguline variety to certain points with non-regular Hodge-Tate weights. With the local models we are able to prove, under the Taylor-Wiles hypothesis, the existence of certain companion points on the global eigenvariety and the appearance of related companion constituents in the completed cohomology for non-regular crystalline Galois representations. The new ingredients in the proof of the global applications are results relating the partial classicality of locally analytic representations (the existence of non-zero locally algebraic vectors in the non-Borel parabolic Emerton's Jacquet modules), the partially de Rham properties of Galois representations (the de Rhamness of graded pieces along the paraboline filtrations of the associated $(\varphi,\Gamma)$-modules over the Robba rings) and the relevant properties of cycles on the generalized Steinberg varieties. We prove that partial classicality implies partial de Rhamness in finite slope cases using Ding's partial eigenvarieties.


title: On single-point inversions of magnetic dipole lines in the corona

abstract: Prompted by a recent paper by Dima and Schad, we re-consider the problem of inferring magnetic properties of the corona using polarimetric observations of magnetic dipole (M1) lines. Dima and Schad point to a potential source of degeneracy in a formalism developed by Plowman, which under some circumstances can lead to the solution being under-determined. Here we clarify the nature of the problem. Its resolution lies in solving for the scattering geometry using the elongation of the observed region of the corona. We discuss some conceptual problems that arise when casting the problem for inversion in the observer's reference frame, and satisfactorily resolve difficulties identified by Plowman, Dima and Schad.


title: Real-Time Forecasting of Driver-Vehicle Dynamics on 3D Roads: a Deep-Learning Framework Leveraging Bayesian Optimisation

abstract: Most state-of-the-art works in trajectory forecasting for automotive target predicting the pose and orientation of the agents in the scene. This represents a particularly useful problem, for instance in autonomous driving, but it does not cover a spectrum of applications in control and simulation that require information on vehicle dynamics features other than pose and orientation. Also, multi-step dynamic simulation of complex multibody models does not seem to be a viable solution for real-time long-term prediction, due to the high computational time required. To bridge this gap, we present a deep-learning framework to model and predict the evolution of the coupled driver-vehicle system dynamics jointly on a complex road geometry. It consists of two components. The first, a neural network predictor, is based on Long Short-Term Memory autoencoders and fuses the information on the road geometry and the past driver-vehicle system dynamics to produce context-aware predictions. The second, a Bayesian optimiser, is proposed to tune some significant hyperparameters of the network. These govern the network complexity, as well as the features importance. The result is a self-tunable framework with real-time applicability, which allows the user to specify the features of interest. The approach has been validated with a case study centered on motion cueing algorithms, using a dataset collected during test sessions of a non-professional driver on a dynamic driving simulator. A 3D track with complex geometry has been employed as driving environment to render the prediction task challenging. Finally, the robustness of the neural network to changes in the driver and track was investigated to set guidelines for future works.


title: Entropy-stable discontinuous Galerkin difference methods for hyperbolic conservation laws

abstract: The paper describes the construction of entropy-stable discontinuous Galerkin difference (DGD) discretizations for hyperbolic conservation laws on unstructured grids. The construction takes advantage of existing theory for entropy-stable summation-by-parts (SBP) discretizations. In particular, the paper shows how DGD discretizations -- both linear and nonlinear -- can be constructed by defining the SBP trial and test functions in terms of interpolated DGD degrees of freedom. In the case of entropy-stable discretizations, the entropy variables rather than the conservative variables must be interpolated to the SBP nodes. A fully-discrete entropy-stable scheme is obtained by adopting the relaxation Runge-Kutta version of the midpoint method. In addition, DGD matrix operators for the first derivative are shown to be dense-norm SBP operators. Numerical results are presented to verify the accuracy and entropy-stability of the DGD discretization in the context of the Euler equations. The results suggest that DGD and SBP solution errors are similar for the same number of degrees of freedom. Finally, an investigation of the DGD spectra shows that spectral radius is relatively insensitive to discretization order; however, the high-order methods do suffer from the linear instability reported for other entropy-stable discretizations.


title: Multi-Session Visual SLAM for Illumination Invariant Re-Localization in Indoor Environments

abstract: For robots navigating using only a camera, illumination changes in indoor environments can cause re-localization failures during autonomous navigation. In this paper, we present a multi-session visual SLAM approach to create a map made of multiple variations of the same locations in different illumination conditions. The multi-session map can then be used at any hour of the day for improved re-localization capability. The approach presented is independent of the visual features used, and this is demonstrated by comparing re-localization performance between multi-session maps created using the RTAB-Map library with SURF, SIFT, BRIEF, BRISK, KAZE, DAISY and SuperPoint visual features. The approach is tested on six mapping and six localization sessions recorded at 30 minute intervals during sunset using a Google Tango phone in a real apartment.


title: Ricci curvature of Bruhat orders

abstract: We study the Ricci curvature of the Hasse diagrams of the Bruhat order of finite irreducible Coxeter groups. For this purpose we compute the maximum degree of these graphs for types $B_n$ and $D_n$. The proof uses a new graph $\Gamma(\pi)$ defined for any element $\pi$ in the corresponding group.


title: Onboarding in Software Product Lines: ConceptMaps as Welcome Guides

abstract: With a volatile labour and technological market, onboarding is becoming increasingly important. The process of incorporating a new developer, a.k.a. the newcomer, into a software development team is reckoned to be lengthy, frustrating and expensive. Newcomers face personal, interpersonal, process and technical barriers during their incorporation, which, in turn, affects the overall productivity of the whole team. This problem exacerbates for Software Product Lines (SPLs), where their size and variability combine to make onboarding even more challenging, even more so for developers that are transferred from the Application Engineering team into the Domain Engineering team, who will be our target newcomers. This work presents concept maps on the role of sensemaking scaffolds to help to introduce these newcomers into the SPL domain. Concept maps, used as knowledge visualisation tools, have been proven to be helpful for meaningful learning. Our main insight is to capture concepts of the SPL domain and their interrelationships in a concept map, and then, present them incrementally, helping newcomers grasp the SPL and aiding them in exploring it in a guided manner while avoiding information overload. This work's contributions are four-fold. First, concept maps are proposed as a representation to introduce newcomers into the SPL domain. Second, concept maps are presented as the means for a guided exploration of the SPL core assets. Third, a feature-driven concept map construction process is introduced. Last, the usefulness of concept maps as guides for SPL onboarding is tested through a formative evaluation. Link to the online demo: url="https://rebrand.ly/wacline-cmap"


title: Certificates of quantum many-body properties assisted by machine learning

abstract: Computationally intractable tasks are often encountered in physics and optimization. Such tasks often comprise a cost function to be optimized over a so-called feasible set, which is specified by a set of constraints. This may yield, in general, to difficult and non-convex optimization tasks. A number of standard methods are used to tackle such problems: variational approaches focus on parameterizing a subclass of solutions within the feasible set; in contrast, relaxation techniques have been proposed to approximate it from outside, thus complementing the variational approach by providing ultimate bounds to the global optimal solution. In this work, we propose a novel approach combining the power of relaxation techniques with deep reinforcement learning in order to find the best possible bounds within a limited computational budget. We illustrate the viability of the method in the context of finding the ground state energy of many-body quantum systems, a paradigmatic problem in quantum physics. We benchmark our approach against other classical optimization algorithms such as breadth-first search or Monte-Carlo, and we characterize the effect of transfer learning. We find the latter may be indicative of phase transitions, with a completely autonomous approach. Finally, we provide tools to generalize the approach to other common applications in the field of quantum information processing.


title: Tor circuit fingerprinting defenses using adaptive padding

abstract: Online anonymity and privacy has been based on confusing the adversary by creating indistinguishable network elements. Tor is the largest and most widely deployed anonymity system, designed against realistic modern adversaries. Recently, researchers have managed to fingerprint Tor's circuits -- and hence the type of underlying traffic -- simply by capturing and analyzing traffic traces. In this work, we study the circuit fingerprinting problem, isolating it from website fingerprinting, and revisit previous findings in this model, showing that accurate attacks are possible even when the application-layer traffic is identical. We then proceed to incrementally create defenses against circuit fingerprinting, using a generic adaptive padding framework for Tor based on WTF-PAD. We present a simple defense which delays a fraction of the traffic, as well as a more advanced one which can effectively hide onion service circuits with zero delays. We thoroughly evaluate both defenses, both analytically and experimentally, discovering new subtle fingerprints, but also showing the effectiveness of our defenses.


title: Performance and Complexity Analysis of bi-directional Recurrent Neural Network Models vs. Volterra Nonlinear Equalizers in Digital Coherent Systems

abstract: We investigate the complexity and performance of recurrent neural network (RNN) models as post-processing units for the compensation of fibre nonlinearities in digital coherent systems carrying polarization multiplexed 16-QAM and 32-QAM signals. We evaluate three bi-directional RNN models, namely the bi-LSTM, bi-GRU and bi-Vanilla-RNN and show that all of them are promising nonlinearity compensators especially in dispersion unmanaged systems. Our simulations show that during inference the three models provide similar compensation performance, therefore in real-life systems the simplest scheme based on Vanilla-RNN units should be preferred. We compare bi-Vanilla-RNN with Volterra nonlinear equalizers and exhibit its superiority both in terms of performance and complexity, thus highlighting that RNN processing is a very promising pathway for the upgrade of long-haul optical communication systems utilizing coherent detection.


title: Improving the Utility of Poisson-Distributed, Differentially Private Synthetic Data via Prior Predictive Truncation with an Application to CDC WONDER

abstract: CDC WONDER is a web-based tool for the dissemination of epidemiologic data collected by the National Vital Statistics System. While CDC WONDER has built-in privacy protections, they do not satisfy formal privacy protections such as differential privacy and thus are susceptible to targeted attacks. Given the importance of making high-quality public health data publicly available while preserving the privacy of the underlying data subjects, we aim to improve the utility of a recently developed approach for generating Poisson-distributed, differentially private synthetic data by using publicly available information to truncate the range of the synthetic data. Specifically, we utilize county-level population information from the U.S. Census Bureau and national death reports produced by the CDC to inform prior distributions on county-level death rates and infer reasonable ranges for Poisson-distributed, county-level death counts. In doing so, the requirements for satisfying differential privacy for a given privacy budget can be reduced by several orders of magnitude, thereby leading to substantial improvements in utility. To illustrate our proposed approach, we consider a dataset comprised of over 26,000 cancer-related deaths from the Commonwealth of Pennsylvania belonging to over 47,000 combinations of cause-of-death and demographic variables such as age, race, sex, and county-of-residence and demonstrate the proposed framework's ability to preserve features such as geographic, urban/rural, and racial disparities present in the true data.


title: Intercensal updating using structure-preserving methods and satellite imagery

abstract: Censuses are fundamental building blocks of most modern-day societies, yet collected every ten years at best. We propose an extension of the widely popular census updating technique Structure Preserving Estimation by incorporating auxiliary information in order to take ongoing subnational population shifts into account. We apply our method by incorporating satellite imagery as additional source to derive annual small-area updates of multidimensional poverty indicators from 2013 to 2020 for a population at risk: female-headed households in Senegal. We evaluate the performance of our proposal using data from two different census periods.


title: Digital Interference Mitigation in Space Division Multiplexing Self-Homodyne Coherent Detection

abstract: We propose a digital interference mitigation scheme to reduce the impact of mode coupling in space division multiplexing self-homodyne coherent detection and experimentally verify its effectiveness in 240-Gbps mode-multiplexed transmission over 3-mode multimode fiber.


title: Human Activity Recognition using Deep Learning Models on Smartphones and Smartwatches Sensor Data

abstract: In recent years, human activity recognition has garnered considerable attention both in industrial and academic research because of the wide deployment of sensors, such as accelerometers and gyroscopes, in products such as smartphones and smartwatches. Activity recognition is currently applied in various fields where valuable information about an individual's functional ability and lifestyle is needed. In this study, we used the popular WISDM dataset for activity recognition. Using multivariate analysis of covariance (MANCOVA), we established a statistically significant difference (p<0.05) between the data generated from the sensors embedded in smartphones and smartwatches. By doing this, we show that smartphones and smartwatches don't capture data in the same way due to the location where they are worn. We deployed several neural network architectures to classify 15 different hand and non-hand-oriented activities. These models include Long short-term memory (LSTM), Bi-directional Long short-term memory (BiLSTM), Convolutional Neural Network (CNN), and Convolutional LSTM (ConvLSTM). The developed models performed best with watch accelerometer data. Also, we saw that the classification precision obtained with the convolutional input classifiers (CNN and ConvLSTM) was higher than the end-to-end LSTM classifier in 12 of the 15 activities. Additionally, the CNN model for the watch accelerometer was better able to classify non-hand oriented activities when compared to hand-oriented activities.


title: Inverse design of Raman amplifier in frequency and distance domain using Convolutional Neural Networks

abstract: We present a Convolutional Neural Network (CNN) architecture for inverse Raman amplifier design. This model aims at finding the pump powers and wavelengths required for a target signal power evolution, both in distance along the fiber and in frequency. Using the proposed framework, the prediction of the pump configuration required to achieve a target power profile is demonstrated numerically with high accuracy in C-band considering both counter-propagating and bidirectional pumping schemes. For a distributed Raman amplifier based on a 100 km single-mode fiber, a low mean set (0.51, 0.54 and 0.64 dB) and standard deviation set (0.62, 0.43 and 0.38 dB) of the maximum test error are obtained numerically employing 2 and 3 counter, and 4 bidirectional propagating pumps, respectively.

categories: eess.SP eess.IV physics.app-ph physics.optics


title: On the invariance properties of Vaidya-Bonner geodesics via symmetry operators

abstract: In the present paper, we try to investigate the Noether symmetries and Lie point symmetries of the Vaidya-Bonner geodesics. Classification of one-dimensional subalgebras of Lie point symmetries are considered. In fact, the collection of pairwise non-conjugate one-dimensional subalgebras that are called the optimal system of subalgebras is determined. Moreover, as illustrative examples, the symmetry analysis is implemented on two special cases of the system.


title: Hot-carrier dynamics in InAs/AlAsSb multiple-quantum wells

abstract: A type-II InAs/AlAs$_{0.16}$Sb$_{0.84}$ multiple-quantum well sample is investigated for the photoexcited carrier dynamics as a function of excitation photon energy and lattice temperature. Time-resolved measurements are performed using a near-infrared pump pulse, with photon energies near to and above the band gap, probed with a terahertz probe pulse. The transient terahertz absorption is characterized by a multi-rise, multi-decay function that captures long-lived decay times and a metastable state of for an excess-photon energy of $>100$ meV. For sufficient excess-photon energy, excitation of the metastable state is followed by a transition to the long-lived states. Excitation dependence of the long-lived states map onto a near-direct band gap ($E{_g}$) density of states with an Urbach tail below $E{_g}$. As temperature increases, the long-lived decay times increase $<E{_g}$, due to the increased phonon interaction of the unintentional defect states, and by phonon stabilization of the hot carriers $>E{_g}$. Additionally, Auger (and/or trap-assisted Auger) scattering above the onset of the plateau may also contribute to longer hot-carrier lifetimes. Meanwhile, the initial decay component shows strong dependence on excitation energy and temperature, reflecting the complicated initial transfer of energy between valence-band and defect states, indicating methods to further prolong hot carriers for technological applications.


title: Self-Supervised Longitudinal Neighbourhood Embedding

abstract: Longitudinal MRIs are often used to capture the gradual deterioration of brain structure and function caused by aging or neurological diseases. Analyzing this data via machine learning generally requires a large number of ground-truth labels, which are often missing or expensive to obtain. Reducing the need for labels, we propose a self-supervised strategy for representation learning named Longitudinal Neighborhood Embedding (LNE). Motivated by concepts in contrastive learning, LNE explicitly models the similarity between trajectory vectors across different subjects. We do so by building a graph in each training iteration defining neighborhoods in the latent space so that the progression direction of a subject follows the direction of its neighbors. This results in a smooth trajectory field that captures the global morphological change of the brain while maintaining the local continuity. We apply LNE to longitudinal T1w MRIs of two neuroimaging studies: a dataset composed of 274 healthy subjects, and Alzheimer's Disease Neuroimaging Initiative (ADNI, N=632). The visualization of the smooth trajectory vector field and superior performance on downstream tasks demonstrate the strength of the proposed method over existing self-supervised methods in extracting information associated with normal aging and in revealing the impact of neurodegenerative disorders. The code is available at \url{https://github.com/ouyangjiahong/longitudinal-neighbourhood-embedding.git}.


title: Generating Images with Sparse Representations

abstract: The high dimensionality of images presents architecture and sampling-efficiency challenges for likelihood-based generative models. Previous approaches such as VQ-VAE use deep autoencoders to obtain compact representations, which are more practical as inputs for likelihood-based models. We present an alternative approach, inspired by common image compression methods like JPEG, and convert images to quantized discrete cosine transform (DCT) blocks, which are represented sparsely as a sequence of DCT channel, spatial location, and DCT coefficient triples. We propose a Transformer-based autoregressive architecture, which is trained to sequentially predict the conditional distribution of the next element in such sequences, and which scales effectively to high resolution images. On a range of image datasets, we demonstrate that our approach can generate high quality, diverse images, with sample metric scores competitive with state of the art methods. We additionally show that simple modifications to our method yield effective image colorization and super-resolution models.


title: Overcoming Poor Word Embeddings with Word Definitions

abstract: Modern natural language understanding models depend on pretrained subword embeddings, but applications may need to reason about words that were never or rarely seen during pretraining. We show that examples that depend critically on a rarer word are more challenging for natural language inference models. Then we explore how a model could learn to use definitions, provided in natural text, to overcome this handicap. Our model's understanding of a definition is usually weaker than a well-modeled word embedding, but it recovers most of the performance gap from using a completely untrained word.


title: Finite element discretization methods for velocity-pressure and stream function formulations of surface Stokes equations

abstract: In this paper we study parametric TraceFEM and parametric SurfaceFEM (SFEM) discretizations of a surface Stokes problem. These methods are applied both to the Stokes problem in velocity-pressure formulation and in stream function formulation. A class of higher order methods is presented in a unified framework. Numerical efficiency aspects of the two formulations are discussed and a systematic comparison of TraceFEM and SFEM is given. A benchmark problem is introduced in which a scalar reference quantity is defined and numerically determined.


title: The Post-impact Evolution of the X-ray Emitting Gas in SNR 1987A Viewed by XMM-Newton

abstract: Since 1996 the blast wave driven by SN 1987A has been interacting with the dense circumstellar material, which provides us with a unique opportunity to study the early evolution of a newborn supernova remnant (SNR). Based on the XMM-Newton RGS and EPIC-pn X-ray observations from 2007 to 2019, we investigated the post-impact evolution of the X-ray emitting gas in SNR 1987A. The hot plasma is represented by two non-equilibrium ionization components with temperature of $\sim0.6$ keV and $\sim2.5$ keV. The low-temperature plasma has a density $\sim2400$ cm$^{-3}$, which is likely dominated by the lower density gas inside the equatorial ring (ER). The high-temperature plasma with a density $\sim550$ cm$^{-3}$ could be dominated by the H II region and the high-latitude material beyond the ring. In the last few years, the emission measure of the low-temperature plasma has been decreasing, indicating that the blast wave has left the main ER. But the blast wave is still propagating into the high-latitude gas, resulting in the steadily increase of the high-temperature emission measure. In the meantime, the average abundances of N, O, Ne, and Mg are found to be declining, which may reflect the different chemical compositions between two plasma components. We also detected the Fe K lines in most of the observations, showing increasing flux and centroid energy. We interpret the Fe K lines as from a third hot component, which may come from the reflected shock-heated gas or originate from Fe-rich ejecta clumps, shocked by the reverse shock.


title: The gravitational bending angle by static and spherically symmetric black holes in bumblebee gravity

abstract: This work investigates the influence of the Lorentz symmetry breaking in the bending angle of massive particles and light for bumblebee black hole solutions. The solutions analyzed break the Lorentz symmetry due to a non-zero vacuum expectation value of the bumblebee field. We use the Ishihara method, which allows us to study the bending angle of light for finite distances, and it is applicable to non-asymptotically flat spacetimes when considering the receiver viewpoint. In order to analyze the deflection of massive particles, we systematize the Ishihara method for its application in the Jacobi metric. This systematization allows the study of the deflection angle of massive particles using the Gauss-Bonnet theorem. We consider two backgrounds: the first was found by Bertolami et al. and is asymptotically flat. The second was found recently by Maluf et al. and is not asymptotically flat due to an effective cosmological constant.


title: Anomalicious: Automated Detection of Anomalous and Potentially Malicious Commits on GitHub

abstract: Security is critical to the adoption of open source software (OSS), yet few automated solutions currently exist to help detect and prevent malicious contributions from infecting open source repositories. On GitHub, a primary host of OSS, repositories contain not only code but also a wealth of commit-related and contextual metadata - what if this metadata could be used to automatically identify malicious OSS contributions? In this work, we show how to use only commit logs and repository metadata to automatically detect anomalous and potentially malicious commits. We identify and evaluate several relevant factors which can be automatically computed from this data, such as the modification of sensitive files, outlier change properties, or a lack of trust in the commit's author. Our tool, Anomalicious, automatically computes these factors and considers them holistically using a rule-based decision model. In an evaluation on a data set of 15 malware-infected repositories, Anomalicious showed promising results and identified 53.33% of malicious commits, while flagging less than 1% of commits for most repositories. Additionally, the tool found other interesting anomalies that are not related to malicious commits in an analysis of repositories with no known malicious commits.


title: Analytic genericity of diffusing orbits in a priori unstable Hamiltonian systems

abstract: The genericity of Arnold diffusion in the analytic category is an open problem. In this paper, we study this problem in the following a priori unstable Hamiltonian system with a time-periodic perturbation \[\mathcal{H}_\varepsilon(p,q,I,\varphi,t)=h(I)+\sum_{i=1}^n\pm \left(\frac{1}{2}p_i^2+V_i(q_i)\right)+\varepsilon H_1(p,q,I,\varphi, t), \] where $(p,q)\in \mathbb{R}^n\times\mathbb{T}^n$, $(I,\varphi)\in\mathbb{R}^d\times\mathbb{T}^d$ with $n, d\geq 1$, $V_i$ are Morse potentials, and $\varepsilon$ is a small non-zero parameter. The unperturbed Hamiltonian is not necessarily convex, and the induced inner dynamics does not need to satisfy a twist condition. Using geometric methods we prove that Arnold diffusion occurs for generic analytic perturbations $H_1$. Indeed, the set of admissible $H_1$ is $C^\omega$ dense and $C^3$ open (a fortiori, $C^\omega$ open). Our perturbative technique for the genericity is valid in the $C^k$ topology for all $k\in [3,\infty)\cup\{\infty, \omega\}$.


title: Isometries of the quaternionic hyperbolic line

abstract: We give a classification of the matrices in the unitary group U(1,1;H),where H is the division ring of the real quaternions. To this end, we consider the complex representation phi(P) for P in U(1,1;H). Next, we compute the characteristic polynomial f of the 4x4 complex matrix phi(P) and then study the sign of the resultant of f and its derivative f'.


title: Optimal Path Planning using CAMIS: a Continuous Anisotropic Model for Inclined Surfaces

abstract: The optimal traverse of irregular terrains made by ground mobile robots heavily depends on the adequacy of the cost models used to plan the path they follow. The criteria to define optimality may be based on minimizing energy consumption and/or preserving the robot stability. This entails the proper assessment of anisotropy to account for the robot driving on top of slopes with different directions. To fulfill this demand, this paper presents the Continuous Anisotropic Model for Inclined Surfaces, a cost model compatible with anisotropic path planners like the bi-directional Ordered Upwind Method. This model acknowledges how the orientation of the robot with respect to any slope determines its energetic cost, considering the action of gravity and terramechanic effects such as the slippage. Moreover, the proposed model can be tuned to define a trade-off between energy minimization and Roll angle reduction. The results from two simulation tests demonstrate how, to find the optimal path in scenarios containing slopes, in certain situations the use of this model can be more advantageous than relying on isotropic cost functions. Finally, the outcome of a field experiment involving a skid-steering robot that drives on top of a real slope is also discussed.


title: Evolution of Splashback Boundaries and Gaseous Outskirts: Insights from Mergers of Self-similar Galaxy Clusters

abstract: A self-similar spherical collapse model predicts a dark matter (DM) splashback and accretion shock in the outskirts of galaxy clusters while misses a key ingredient of structure formation - processes associated with mergers. To fill this gap, we perform simulations of merging self-similar clusters and investigate their DM and gas evolution in an idealized cosmological context. Our simulations show that the cluster rapidly contracts during the major merger and the splashback radius $r_{\rm sp}$ decreases, approaching the virial radius $r_{\rm vir}$. While $r_{\rm sp}$ correlates with a smooth mass accretion rate (MAR) parameter $\Gamma_{\rm s}$ in the self-similar model, our simulations show a similar trend with the total MAR $\Gamma_{\rm vir}$ (includes both mergers and $\Gamma_{\rm s}$). The scatter of the $\Gamma_{\rm vir}-r_{\rm sp}/r_{\rm vir}$ relation indicates a generally low $\Gamma_{\rm s}\sim1$ in clusters in cosmological simulations. In contrast to the DM, the hot gaseous atmospheres significantly expand by the merger-accelerated (MA-) shocks formed when the runaway merger shocks overtake the outer accretion shock. After a major merger, the MA-shock radius is larger than $r_{\rm sp}$ by a factor of up to $\sim1.7$ for $\Gamma_{\rm s}\lesssim1$ and is $\sim r_{\rm sp}$ for $\Gamma_{\rm s}\gtrsim3$. This implies that (1) mergers could easily generate the MA-shock-splashback offset measured in cosmological simulations, and (2) the smooth MAR is small in regions away from filaments where MA-shocks reside. We further discuss various shocks and contact discontinuities formed at different epochs of the merger, the ram pressure stripping in cluster outskirts, and the dependence of member galaxies' splashback feature on their orbital parameters.


title: SoK: Cryptojacking Malware

abstract: Emerging blockchain and cryptocurrency-based technologies are redefining the way we conduct business in cyberspace. Today, a myriad of blockchain and cryptocurrency systems, applications, and technologies are widely available to companies, end-users, and even malicious actors who want to exploit the computational resources of regular users through \textit{cryptojacking} malware. Especially with ready-to-use mining scripts easily provided by service providers (e.g., Coinhive) and untraceable cryptocurrencies (e.g., Monero), cryptojacking malware has become an indispensable tool for attackers. Indeed, the banking industry, major commercial websites, government and military servers (e.g., US Dept. of Defense), online video sharing platforms (e.g., Youtube), gaming platforms (e.g., Nintendo), critical infrastructure resources (e.g., routers), and even recently widely popular remote video conferencing/meeting programs (e.g., Zoom during the Covid-19 pandemic) have all been the victims of powerful cryptojacking malware campaigns. Nonetheless, existing detection methods such as browser extensions that protect users with blacklist methods or antivirus programs with different analysis methods can only provide a partial panacea to this emerging cryptojacking issue as the attackers can easily bypass them by using obfuscation techniques or changing their domains or scripts frequently. Therefore, many studies in the literature proposed cryptojacking malware detection methods using various dynamic/behavioral features.


title: High thermoelectric power factor of poly(3-hexylthiophene) through in-plane alignment and doping with a molybdenum dithiolene complex

abstract: Here we report a record thermoelectric power factor of up to 160 $\mu$ W m-1 K-2 for the conjugated polymer poly(3-hexylthiophene) (P3HT). This result is achieved through the combination of high-temperature rubbing of thin films together with the use of a large molybdenum dithiolene p-dopant with a high electron affinity. Comparison of the UV-vis-NIR spectra of the chemically doped samples to electrochemically oxidized material reveals an oxidation level of 10%, i.e. one polaron for every 10 repeat units. The high power factor arises due to an increase in the charge-carrier mobility and hence electrical conductivity along the rubbing direction. We conclude that P3HT, with its facile synthesis and outstanding processability, should not be ruled out as a potential thermoelectric material.


title: Quantum control of a nanoparticle optically levitated in cryogenic free space

abstract: Tests of quantum mechanics on a macroscopic scale require extreme control over mechanical motion and its decoherence. Quantum control of mechanical motion has been achieved by engineering the radiation-pressure coupling between a micromechanical oscillator and the electromagnetic field in a resonator. Furthermore, measurement-based feedback control relying on cavity-enhanced detection schemes has been used to cool micromechanical oscillators to their quantum ground states. In contrast to mechanically tethered systems, optically levitated nanoparticles are particularly promising candidates for matter-wave experiments with massive objects, since their trapping potential is fully controllable. In this work, we optically levitate a femto-gram dielectric particle in cryogenic free space, which suppresses thermal effects sufficiently to make the measurement backaction the dominant decoherence mechanism. With an efficient quantum measurement, we exert quantum control over the dynamics of the particle. We cool its center-of-mass motion by measurement-based feedback to an average occupancy of 0.65 motional quanta, corresponding to a state purity of 43%. The absence of an optical resonator and its bandwidth limitations holds promise to transfer the full quantum control available for electromagnetic fields to a mechanical system. Together with the fact that the optical trapping potential is highly controllable, our experimental platform offers a route to investigating quantum mechanics at macroscopic scales.


title: A Pilot Study on Visually Stimulated Cognitive Tasks for EEG-Based Dementia Recognition

abstract: In the status quo, dementia is yet to be cured. Precise diagnosis prior to the onset of the symptoms can prevent the rapid progression of the emerging cognitive impairment. Recent progress has shown that Electroencephalography (EEG) is the promising and cost-effective test to facilitate the detection of neurocognitive disorders. However, most of the existing works have been using only resting-state EEG. The efficiencies of EEG signals from various cognitive tasks, for dementia classification, have yet to be thoroughly investigated. In this study, we designed four cognitive tasks that engage different cognitive performances: attention, working memory, and executive function. We investigated these tasks by using statistical analysis on both time and frequency domains of EEG signals from three classes of human subjects: Dementia (DEM), Mild Cognitive Impairment (MCI), and Normal Control (NC). We also further evaluated the classification performances of two features extraction methods: Principal Component Analysis (PCA) and Filter Bank Common Spatial Pattern (FBCSP). We found that the working memory related tasks yielded good performances for dementia recognition in both cases using PCA and FBCSP. Moreover, FBCSP with features combination from four tasks revealed the best sensitivity of 0.87 and the specificity of 0.80. To our best knowledge, this is the first work that concurrently investigated several cognitive tasks for dementia recognition using both statistical analysis and classification scores. Our results yielded essential information to design and aid in conducting further experimental tasks to early diagnose dementia patients.


title: Shadows of parametrized axially symmetric black holes allowing for separation of variables

abstract: Metric of axially symmetric asymptotically flat black holes in an arbitrary metric theory of gravity can be represented in the general form which depends on infinite number of parameters. We constrain this general class of metrics by requiring the existence of additional symmetries, which lead to the separation of variables in the Hamilton-Jacobi and Klein-Gordon equations, and show that once the metric functions change sufficiently moderately in some region near the black hole, the black-hole shadow depends on a few deformation parameters only. We analyze the influence of these parameters on the black-hole shadow. We also show that the shadow of the rotating black hole in the Einstein-dilaton-Gauss-Bonnet theory is well approximated if the terms violating the separation of variables are neglected in the metric.


title: Experimentally Validated Hopping-Transport Model for Energetically Disordered Organic Semiconductors

abstract: Charge transport in disordered organic semiconductors occurs by hopping of charge carriers between localized sites that are randomly distributed in a strongly energy dependent density of states. Extracting disorder and hopping parameters from experimental data like temperature dependent current-voltage characteristics typically relies on parametrized mobility functionals that are integrated in a drift-diffusion solver. Surprisingly, the functional based on the extended Gaussian disorder model (eGDM) has been extremely successful at this, despite it being based on the assumption of nearest neighbor hopping (nnH) on a regular lattice. We here propose a variable range hopping (VRH) model that has been integrated in a freeware drift-diffusion solver. The mobility model has been calibrated using kinetic Monte Carlo calculations and shows good agreement with the Monte Carlo calculations over the experimentally relevant part of the parameter space. The model is applied to temperature-dependent space charge limited current (SCLC) measurements of different systems. In contrast to the eGDM, the VRH model provides a consistent description of both p-type and n-type devices. We find a critical ratio of aNN/$\alpha$ (mean inter-site distance / localization radius) of ~3 below which hopping to non-nearest neighbors becomes important around room temperature and the eGDM cannot be used for parameter extraction. Typical (Gaussian) disorder values in the range 45-120 meV are found, without any clear correlation with photovoltaic performance when the same active layer is used in an organic solar cell.


title: Revisiting the g-null paradox

abstract: The parametric g-formula is an approach to estimating causal effects of sustained treatment strategies from observational data. An often cited limitation of the parametric g-formula is the g-null paradox: a phenomenon in which model misspecification in the parametric g-formula is guaranteed under the conditions that motivate its use (i.e., when identifiability conditions hold and measured time-varying confounders are affected by past treatment). Many users of the parametric g-formula know they must acknowledge the g-null paradox as a limitation when reporting results but still require clarity on its meaning and implications. Here we revisit the g-null paradox to clarify its role in causal inference studies. In doing so, we present analytic examples and a simulation-based illustration of the bias of parametric g-formula estimates under the conditions associated with this paradox. Our results highlight the importance of avoiding overly parsimonious models for the components of the g-formula when using this method.


title: Three-dimensional realizations of flood flow in large-scale rivers using the neural fuzzy-based machine-learning algorithm

abstract: Machine learning methods have been extensively used to study the dynamics of complex fluid flows. One such algorithm, known as adaptive neural fuzzy inference system (ANFIS), can generate data-driven predictions for flow fields but has not been applied to natural geophysical flows in large-scale rivers. Herein, we demonstrate the potential of ANFIS to produce three-dimensional (3D) realizations of the instantaneous flood flow field in several large-scale, virtual meandering rivers. The 3D dynamics of flood flow in large-scale rivers were obtained using large-eddy simulation (LES). The LES results, i.e., the 3D velocity components, were employed to train the learnable coefficients of an ANFIS. The trained ANFIS, along with a few time-steps of LES results (precursor data) were then used to produce 3D realizations of flood flow fields in large-scale rivers with geometries other than the one the ANFIS was trained with. We also used the trained ANFIS to generate 3D realizations of river flow at a discharge other than that the ANFIS was trained with. The flow field results obtained from ANFIS were validated using separate LES runs to assess the accuracy of the 3D instantaneous realizations of the machine learning algorithm. An error analysis was conducted to quantify the discrepancies among the ANFIS and LES results for various flood flow predictions in large-scale rivers.


title: Acoustic Su-Schrieffer-Heeger lattice: Direct mapping of acoustic waveguides to Su-Schrieffer-Heeger model

abstract: Topological physics strongly relies on prototypical lattice model with particular symmetries. We report here on a theoretical and experimental work on acoustic waveguides that is directly mapped to the one-dimensional Su-Schrieffer-Heeger chiral model. Starting from the continuous two dimensional wave equation we use a combination of monomadal approximation and the condition of equal length tube segments to arrive at the wanted discrete equations. It is shown that open or closed boundary conditions topological leads automatically to the existence of edge modes. We illustrate by graphical construction how the edge modes appear naturally owing to a quarter-wavelength condition and the conservation of flux. Furthermore, the transparent chirality of our system, which is ensured by the geometrical constraints allows us to study chiral disorder numerically and experimentally. Our experimental results in the audible regime demonstrate the predicted robustness of the topological edge modes.


title: A Learning-Based Approach to Address Complexity-Reliability Tradeoff in OS Decoders

abstract: In this paper, we study the tradeoffs between complexity and reliability for decoding large linear block codes. We show that using artificial neural networks to predict the required order of an ordered statistics based decoder helps in reducing the average complexity and hence the latency of the decoder. We numerically validate the approach through Monte Carlo simulations.


title: Topology-driven ordering of flocking matter

abstract: When interacting motile units self-organize into flocks, they realize one of the most robust ordered state found in nature. However, after twenty five years of intense research, the very mechanism controlling the ordering dynamics of both living and artificial flocks has remained unsettled. Here, combining active-colloid experiments, numerical simulations and analytical work, we explain how flocking liquids heal their spontaneous flows initially plagued by collections of topological defects to achieve long-ranged polar order even in two dimensions. We demonstrate that the self-similar ordering of flocking matter is ruled by a living network of domain walls linking all $\pm 1$ vortices, and guiding their annihilation dynamics. Crucially, this singular orientational structure echoes the formation of extended density patterns in the shape of interconnected bow ties. We establish that this double structure emerges from the interplay between self-advection and density gradients dressing each $-1$ topological charges with four orientation walls. We then explain how active Magnus forces link all topological charges with extended domain walls, while elastic interactions drive their attraction along the resulting filamentous network of polarization singularities. Taken together our experimental, numerical and analytical results illuminate the suppression of all flow singularities, and the emergence of pristine unidirectional order in flocking matter.


title: Harnessing Geometric Constraints from Emotion Labels to improve Face Verification

abstract: For the task of face verification, we explore the utility of harnessing auxiliary facial emotion labels to impose explicit geometric constraints on the embedding space when training deep embedding models. We introduce several novel loss functions that, in conjunction with a standard Triplet Loss [43], or ArcFace loss [10], provide geometric constraints on the embedding space; the labels for our loss functions can be provided using either manually annotated or automatically detected auxiliary emotion labels. Our method is implemented purely in terms of the loss function and does not require any changes to the neural network backbone of the embedding function.


title: KMOS study of the mass accretion rate from Class I to Class II in NGC 1333

abstract: The mass accretion rate is the fundamental parameter to understand the process of mass assembly that results in the formation of a low-mass star. This parameter has been largely studied in Classical TTauri stars in star-forming regions with ages of 1-10Myr. However, little is known about the accretion properties of young stellar objects (YSOs) in younger regions and early stages of star formation, such as in the Class0/I phases. We present new NIR spectra of 17 ClassI/Flat and 35 ClassII sources located in the young (<1Myr) NGC1333 cluster, acquired with the KMOS instrument at the VLT. Our goal is to study whether the mass accretion rate evolves with age, as suggested by the widely adopted viscous evolution model, by comparing the properties of the NGC1333 members with samples of older regions. We measured the stellar parameters and accretion rates of our sample, finding a correlation between accretion and stellar luminosity, and between mass accretion rate and stellar mass. Both correlations are compatible within the errors with the older Lupus star-forming region, while only the latter is consistent with results from ChamaeleonI. The ClassI sample shows larger accretion luminosities with respect to the ClassII stars of the same cloud. However, the derived accretion rates are not sufficiently high to build up the inferred stellar masses, assuming steady accretion during the ClassI lifetime. This suggests that the sources are not in their main accretion phase and that most of their mass has already been accumulated during a previous stage and/or that the accretion is an episodic phenomenon. We show that some of the targets originally classified as Class I through Spitzer photometry are in fact evolved or low accreting objects. This evidence can have implications for the estimated protostellar phase lifetimes. Further observations are needed to determine if this is a general result.


title: Learning to Extend Molecular Scaffolds with Structural Motifs

abstract: Recent advancements in deep learning-based modeling of molecules promise to accelerate in silico drug discovery. A plethora of generative models is available, building molecules either atom-by-atom and bond-by-bond or fragment-by-fragment. However, many drug discovery projects require a fixed scaffold to be present in the generated molecule, and incorporating that constraint has only recently been explored. Here, we propose MoLeR, a graph-based model that naturally supports scaffolds as initial seed of the generative procedure, which is possible because it is not conditioned on the generation history. Our experiments show that MoLeR performs comparably to state-of-the-art methods on unconstrained molecular optimization tasks, and outperforms them on scaffold-based tasks, while being an order of magnitude faster to train and sample from than existing approaches. Furthermore, we show the influence of a number of seemingly minor design choices on the overall performance.


title: A combinatorial statistic for labeled threshold graphs

abstract: Consider the collection of hyperplanes in $\mathbb{R}^n$ whose defining equations are given by $\{x_i + x_j = 0\mid 1\leq i<j\leq n\}$. This arrangement is called the threshold arrangement since its regions are in bijection with labeled threshold graphs on $n$ vertices. Zaslavsky's theorem implies that the number of regions of this arrangement is the sum of coefficients of the characteristic polynomial of the arrangement. In the present article we give a combinatorial meaning to these coefficients as the number of labeled threshold graphs with a certain property, thus answering a question posed by Stanley.


title: Towards Automated Benchmark Support for Multi-Blockchain Interoperability-Facilitating Platforms

abstract: Since the introduction of the first Bitcoin blockchain in 2008, different decentralized blockchain systems such as Ethereum, Hyperledger Fabric, and Corda, have emerged with public and private accessibility. It has been widely acknowledged that no single blockchain network will fit all use cases. As a result, we have observed the increasing popularity of multi-blockchain ecosystem in which customers will move toward different blockchains based on their particular requirements. Hence, the efficiency and security requirements of interactions among these heterogeneous blockchains become critical. In realization of this multi-blockchain paradigm, initiatives in building Interoperability-Facilitating Platforms (IFPs) that aim at bridging different blockchains (a.k.a. blockchain interoperability) have come to the fore. Despite current efforts, it is extremely difficult for blockchain customers (organizations, governments, companies) to understand the trade-offs between different IFPs and their suitability for different application domains before adoption. A key reason is due to a lack of fundamental and systematic approaches to assess the variables among different IFPs. To fill this gap, developing new IFP requirements specification and open-source benchmark tools to advance research in distributed, multi-blockchain interoperability, with emphasis on IFP performance and security challenges are required. In this document, we outline a research proposal study to the community to realize this gap.


title: Effects of the shape of curvature peaks on the size of primordial black holes

abstract: We simulate numerically the formation of spherically symmetric primordial black holes (PBHs) seeded by different families of primordial curvature perturbations profiles on a Friedman-Robertson-Walker (FRW) Universe filled by radiation fluid. We have studied the dependency on the curvature profile of the initial mass $M_{\rm BH,i}$ of the PBHs at the time of apparent horizon formation $t_{AH}$, and the final mass $M_{\rm BH,f}$ after the accretion process, using an excision technique, comparing $M_{\rm BH,i}$ to previous analytical estimations obtained using compensated PBHs model approach. The analytical estimations are in agreement with numerical results, except for large values of the initial perturbation amplitude, when the compensated model is less accurate. The masses $M_{\rm BH,f}$ and $M_{\rm BH,i}$ do not depend only on the shape around the compaction function peak, but on the full profile of the initial curvature perturbation. We also estimate the accretion effects, and for those PBHs with masses relevant for the dark matter abundance, with a final mass equal to the horizon crossing mass, we find $M_{\rm BH,f}\approx 3 M_{\rm BH,i} $.


title: Decomposable Submodular Function Minimization via Maximum Flow

abstract: This paper bridges discrete and continuous optimization approaches for decomposable submodular function minimization, in both the standard and parametric settings. We provide improved running times for this problem by reducing it to a number of calls to a maximum flow oracle. When each function in the decomposition acts on $O(1)$ elements of the ground set $V$ and is polynomially bounded, our running time is up to polylogarithmic factors equal to that of solving maximum flow in a sparse graph with $O(\vert V \vert)$ vertices and polynomial integral capacities. We achieve this by providing a simple iterative method which can optimize to high precision any convex function defined on the submodular base polytope, provided we can efficiently minimize it on the base polytope corresponding to the cut function of a certain graph that we construct. We solve this minimization problem by lifting the solutions of a parametric cut problem, which we obtain via a new efficient combinatorial reduction to maximum flow. This reduction is of independent interest and implies some previously unknown bounds for the parametric minimum $s,t$-cut problem in multiple settings.


title: Lyapunov-Regularized Reinforcement Learning for Power System Transient Stability

abstract: Transient stability of power systems is becoming increasingly important because of the growing integration of renewable resources. These resources lead to a reduction in mechanical inertia but also provide increased flexibility in frequency responses. Namely, their power electronic interfaces can implement almost arbitrary control laws. To design these controllers, reinforcement learning (RL) has emerged as a powerful method in searching for optimal non-linear control policy parameterized by neural networks. A key challenge is to enforce that a learned controller must be stabilizing. This paper proposes a Lyapunov regularized RL approach for optimal frequency control for transient stability in lossy networks. Because the lack of an analytical Lyapunov function, we learn a Lyapunov function parameterized by a neural network. The losses are specially designed with respect to the physical power system. The learned neural Lyapunov function is then utilized as a regularization to train the neural network controller by penalizing actions that violate the Lyapunov conditions. Case study shows that introducing the Lyapunov regularization enables the controller to be stabilizing and achieve smaller losses.


title: Low pseudomoments of Euler products

abstract: In this paper, we determine the order of magnitude of the $2q$-th pseudomoment of powers of the Riemann zeta function $\zeta(s)^{\alpha}$ for $0<q\le 1/2$ and $0< \alpha<1$, completing the results of Bondarenko, Heap and Seip, and of Gerspach. Our results also apply to more general Euler products satisfying certain conditions.


title: Modal Reasoning = Metric Reasoning, via Lawvere

abstract: Graded modal types systems and coeffects are becoming a standard formalism to deal with context-dependent computations where code usage plays a central role. The theory of program equivalence for modal and coeffectful languages, however, is considerably underdeveloped if compared to the denotational and operational semantics of such languages. This raises the question of how much of the theory of ordinary program equivalence can be given in a modal scenario. In this work, we show that coinductive equivalences can be extended to a modal setting, and we do so by generalising Abramsky's applicative bisimilarity to coeffectful behaviours. To achieve this goal, we develop a general theory of ternary program relations based on the novel notion of a comonadic lax extension, on top of which we define a modal extension of Abramsky's applicative bisimilarity (which we dub modal applicative bisimilarity). We prove such a relation to be a congruence, this way obtaining a compositional technique for reasoning about modal and coeffectful behaviours. But this is not the end of the story: we also establish a correspondence between modal program relations and program distances. This correspondence shows that modal applicative bisimilarity and (a properly extended) applicative bisimilarity distance coincide, this way revealing that modal program equivalences and program distances are just two sides of the same coin.


title: Rissanen Data Analysis: Examining Dataset Characteristics via Description Length

abstract: We introduce a method to determine if a certain capability helps to achieve an accurate model of given data. We view labels as being generated from the inputs by a program composed of subroutines with different capabilities, and we posit that a subroutine is useful if and only if the minimal program that invokes it is shorter than the one that does not. Since minimum program length is uncomputable, we instead estimate the labels' minimum description length (MDL) as a proxy, giving us a theoretically-grounded method for analyzing dataset characteristics. We call the method Rissanen Data Analysis (RDA) after the father of MDL, and we showcase its applicability on a wide variety of settings in NLP, ranging from evaluating the utility of generating subquestions before answering a question, to analyzing the value of rationales and explanations, to investigating the importance of different parts of speech, and uncovering dataset gender bias.


title: OperA: Attention-Regularized Transformers for Surgical Phase Recognition

abstract: In this paper we introduce OperA, a transformer-based model that accurately predicts surgical phases from long video sequences. A novel attention regularization loss encourages the model to focus on high-quality frames during training. Moreover, the attention weights are utilized to identify characteristic high attention frames for each surgical phase, which could further be used for surgery summarization. OperA is thoroughly evaluated on two datasets of laparoscopic cholecystectomy videos, outperforming various state-of-the-art temporal refinement approaches.


title: Measuring Mathematical Problem Solving With the MATH Dataset

abstract: Many intellectual endeavors require mathematical problem solving, but this skill remains beyond the capabilities of computers. To measure this ability in machine learning models, we introduce MATH, a new dataset of 12,500 challenging competition mathematics problems. Each problem in MATH has a full step-by-step solution which can be used to teach models to generate answer derivations and explanations. To facilitate future research and increase accuracy on MATH, we also contribute a large auxiliary pretraining dataset which helps teach models the fundamentals of mathematics. Even though we are able to increase accuracy on MATH, our results show that accuracy remains relatively low, even with enormous Transformer models. Moreover, we find that simply increasing budgets and model parameter counts will be impractical for achieving strong mathematical reasoning if scaling trends continue. While scaling Transformers is automatically solving most other text-based tasks, scaling is not currently solving MATH. To have more traction on mathematical problem solving we will likely need new algorithmic advancements from the broader research community.


title: Genetic Algorithm based hyper-parameters optimization for transfer Convolutional Neural Network

abstract: Hyperparameter optimization is a challenging problem in developing deep neural networks. Decision of transfer layers and trainable layers is a major task for design of the transfer convolutional neural networks (CNN). Conventional transfer CNN models are usually manually designed based on intuition. In this paper, a genetic algorithm is applied to select trainable layers of the transfer model. The filter criterion is constructed by accuracy and the counts of the trainable layers. The results show that the method is competent in this task. The system will converge with a precision of 97% in the classification of Cats and Dogs datasets, in no more than 15 generations. Moreover, backward inference according the results of the genetic algorithm shows that our method can capture the gradient features in network layers, which plays a part on understanding of the transfer AI models.


title: Atomic Carbon [CI]$(^3P_1-^3P_0)$ Mapping of the Nearby Galaxy M83

abstract: Atomic carbon (CI) has been proposed to be a global tracer of the molecular gas as a substitute for CO, however, its utility remains unproven. To evaluate the suitability of CI as the tracer, we performed [CI]$(^3P_1-^3P_0)$ (hereinafter [CI](1-0)) mapping observations of the northern part of the nearby spiral galaxy M83 with the ASTE telescope and compared the distributions of [CI](1-0) with CO lines (CO(1-0), CO(3-2), and $^{13}$CO(1-0)), HI, and infrared (IR) emission (70, 160, and 250$ \mu$m). The [CI](1-0) distribution in the central region is similar to that of the CO lines, whereas [CI](1-0) in the arm region is distributed outside the CO. We examined the dust temperature, $T_{\rm dust}$, and dust mass surface density, $\Sigma_{\rm dust}$, by fitting the IR continuum-spectrum distribution with a single-temperature modified blackbody. The distribution of $\Sigma_{\rm dust}$ shows a much better consistency with the integrated intensity of CO(1-0) than with that of [CI](1-0), indicating that CO(1-0) is a good tracer of the cold molecular gas. The spatial distribution of the [CI] excitation temperature, $T_{\rm ex}$, was examined using the intensity ratio of the two [CI] transitions. An appropriate $T_{\rm ex}$ at the central, bar, arm, and inter-arm regions yields a constant [C]/[H$_2$] abundance ratio of $\sim7 \times 10^{-5}$ within a range of 0.1 dex in all regions. We successfully detected weak [CI](1-0) emission, even in the inter-arm region, in addition to the central, arm, and bar regions, using spectral stacking analysis. The stacked intensity of [CI](1-0) is found to be strongly correlated with $T_{\rm dust}$. Our results indicate that the atomic carbon is a photodissociation product of CO, and consequently, compared to CO(1-0), [CI](1-0) is less reliable in tracing the bulk of "cold" molecular gas in the galactic disk.


title: Neural network-based image reconstruction in swept-source optical coherence tomography using undersampled spectral data

abstract: Optical Coherence Tomography (OCT) is a widely used non-invasive biomedical imaging modality that can rapidly provide volumetric images of samples. Here, we present a deep learning-based image reconstruction framework that can generate swept-source OCT (SS-OCT) images using undersampled spectral data, without any spatial aliasing artifacts. This neural network-based image reconstruction does not require any hardware changes to the optical set-up and can be easily integrated with existing swept-source or spectral domain OCT systems to reduce the amount of raw spectral data to be acquired. To show the efficacy of this framework, we trained and blindly tested a deep neural network using mouse embryo samples imaged by an SS-OCT system. Using 2-fold undersampled spectral data (i.e., 640 spectral points per A-line), the trained neural network can blindly reconstruct 512 A-lines in ~6.73 ms using a desktop computer, removing spatial aliasing artifacts due to spectral undersampling, also presenting a very good match to the images of the same samples, reconstructed using the full spectral OCT data (i.e., 1280 spectral points per A-line). We also successfully demonstrate that this framework can be further extended to process 3x undersampled spectral data per A-line, with some performance degradation in the reconstructed image quality compared to 2x spectral undersampling. This deep learning-enabled image reconstruction approach can be broadly used in various forms of spectral domain OCT systems, helping to increase their imaging speed without sacrificing image resolution and signal-to-noise ratio.


title: Some further q-shift difference results on Hayman conjecture

abstract: In this paper, we investigate the zero distributions of $q$-shift difference-differential polynomials of meromorphic functions with zero-order that extends and generalizes the classical Hayman results of the zeros of differential polynomials to q-shift difference. We also investigate the uniqueness problem of $q$-shift difference-differential polynomials sharing a polynomial value with finite weight.


title: Enhanced X-ray Emission from the Most Radio-Powerful Quasar in the Universe's First Billion Years

abstract: We present deep (265 ks) Chandra X-ray observations of PSO J352.4034$-$15.3373, a quasar at z=5.831 that, with a radio-to-optical flux ratio of R>1000, is one of the radio-loudest quasars in the early universe and is the only quasar with observed extended radio jets of kpc-scale at $z \gtrsim 6$. Modeling the X-ray spectrum of the quasar with a power law, we find a best fit of $\Gamma = 1.99^{+0.29}_{-0.28}$, leading to an X-ray luminosity of $L_{2-10} = 1.26^{+0.45}_{-0.33} \times 10^{45}\ {\rm erg}\ {\rm s}^{-1}$ and an X-ray to UV brightness ratio of $\alpha_{\rm OX} = -1.36 \pm 0.11$. We identify a diffuse structure 50 kpc (${\sim}8^{\prime\prime}$) to the NW of the quasar along the jet axis that corresponds to a $3\sigma$ enhancement in the angular density of emission and can be ruled out as a background fluctuation with a probability of P=0.9985. While with few detected photons the spectral fit of the structure is uncertain, we find that it has a luminosity of $L_{2-10}\sim10^{44}\ {\rm erg}\ {\rm s}^{-1}$. These observations therefore potentially represent the most distant quasar jet yet seen in X-rays. We find no evidence for excess X-ray emission where the previously-reported radio jets are seen (which have an overall linear extent of $0.^{\prime\prime}28$), and a bright X-ray point source located along the jet axis to the SE is revealed by optical and NIR imaging to not be associated with the quasar.


title: Spatially resolved spectroscopy across stellar surfaces. IV. F, G, & K-stars: Synthetic 3D spectra at hyper-high resolution

abstract: High-precision stellar analyses require hydrodynamic 3D modeling. Such models predict changes across stellar disks of spectral line shapes, asymmetries, and wavelength shifts. For testing models in stars other than the Sun, spatially resolved observations are feasible from differential spectroscopy during exoplanet transits, retrieving spectra of stellar surface segments that successively become hidden behind the transiting planet, as shown in Papers I, II, and III. Synthetic high-resolution spectra over extended spectral regions are now available from 3D models. Similar to other ab initio simulations, these data contain patterns that have not been specifically modeled but may be revealed after analyses analogous to those of a large volume of observations. From five 3D models spanning T=3964-6726K (approx. spectral types K8V-F3V), synthetic spectra at hyper-high resolution (R>1,000,000) were analyzed. Selected FeI and FeII lines at various positions across stellar disks were searched for patterns between different lines in the same star and for similar lines between different stars. Such patterns are identified for representative photospheric lines of different strengths, excitation potential, and ionization level, encoding the hydrodynamic 3D structure. Line profiles and bisectors are shown for various stars at different positions across stellar disks. Absolute convective wavelength shifts are obtained as differences to 1D models, where such shifts do not occur. Observable relationships for line properties are retrieved from realistically complex synthetic spectra. Such patterns may also test very detailed 3D modeling, including non-LTE effects. While present results are obtained at hyper-high spectral resolution, the subsequent Paper V examines their practical observability at realistically lower resolutions, and in the presence of noise.


title: Near-IR observations of the young star [BHB2007]-1: A sub-stellar companion opening the gap in the disk

abstract: The presence of planets or sub-stellar objects still embedded in their native protoplanetary disks is indirectly suggested by disk sub-structures like gaps, cavities, and spirals. However, these companions are rarely detected. We present VLT/NACO high-contrast images in $J$, $H$, $K_S$, and $L^{\prime}$ band of the young star [BHB2007]-1 probing the inclined disk in scattered light and revealing the probable presence of a companion. The point source is detected in the $L^{\prime}$ band in spatial correspondence with complementary VLA observations. This object is constrained to have a mass in the range of 37-47 M$_{Jup}$ and is located at 50 au from the central star, inside the 70 au-large disk cavity recently imaged by ALMA, that is absent from our NACO data (down to 20 au). This mass range is compatible with the upper end derived from the size of the ALMA cavity. The NIR disk brightness is highly asymmetric around the minor axis, with the southern side 5.5 times brighter than the northern side. The constant amount of asymmetry across all wavelengths suggests that it is due to a shadow cast by a misaligned inner disk. The massive companion that we detect could, in principle, explain the possible disk misalignment, as well as the different cavity sizes inferred by the NACO and ALMA observations. The confirmation and characterization of the companion is entrusted to future observations.