title: CAR T cell therapy in B-cell acute lymphoblastic leukaemia: Insights from mathematical models

abstract: Immunotherapies use components of the patient immune system to selectively target cancer cells. The use of CAR T cells to treat B-cell malignancies --leukaemias and lymphomas-- is one of the most successful examples, with many patients experiencing long-lasting complete responses to this therapy. This treatment works by extracting the patient's T cells and adding them the CAR group, which enables them to recognize and target cells carrying the antigen CD19+, that is expressed in these haematological tumors. Here we put forward a mathematical model describing the time response of leukaemias to the injection of CAR T-cells. The model accounts for mature and progenitor B-cells, tumor cells, CAR T cells and side effects by incorporating the main biological processes involved. The model explains the early post-injection dynamics of the different compartments and the fact that the number of CAR T cells injected does not critically affect the treatment outcome. An explicit formula is found that provides the maximum CAR T cell expansion in-vivo and the severity of side effects. Our mathematical model captures other known features of the response to this immunotherapy. It also predicts that CD19+ tumor relapses could be the result of the competition between tumor and CAR T cells analogous to predator-prey dynamics. We discuss this fact on the light of available evidences and the possibility of controlling relapses by early re-challenging of the tumor with stored CAR T cells.


title: A Gridless Compressive Sensing Based Channel Estimation for Millimeter Wave Massive MIMO Systems from 1-Bit Measurements

abstract: This paper considers the problem of estimating the sparse millimeter wave (mmWave) massive multiple-input multiple-output (MIMO) OFDM channel from 1-bit quantized measurements. Unlike previous quantized on-grid approaches to compressive sensing, we introduce an off-grid convex scheme which uses binary atomic norm minimization (BiANM) to estimate sparse channel form 1-bit measurements. Moreover, in this paper, we propose a decoupled angle-delay version of BiANM called decoupled binary atomic norm minimization (DeBiANM) in order to reduce complexity. Further, for improving the accuracy of the estimation, we develop reweighted version of BiANM named reweighted binary atomic norm minimization (ReBiANM) that is a trade-off between atomic $\ell_0$ norm and atomic $\ell_1$ norm. Also, for improvement the performance of the DeBiANM, we introduce the reweighted version of DeBiANM, termed as reweighted decoupled binary atomic minimization (ReDeBiANM). Finally, through the simulation results, the performance of the proposed methods are studied.


title: Multi-Person Pose Estimation with Enhanced Feature Aggregation and Selection

abstract: We propose a novel Enhanced Feature Aggregation and Selection network (EFASNet) for multi-person 2D human pose estimation. Due to enhanced feature representation, our method can well handle crowded, cluttered and occluded scenes. More specifically, a Feature Aggregation and Selection Module (FASM), which constructs hierarchical multi-scale feature aggregation and makes the aggregated features discriminative, is proposed to get more accurate fine-grained representation, leading to more precise joint locations. Then, we perform a simple Feature Fusion (FF) strategy which effectively fuses high-resolution spatial features and low-resolution semantic features to obtain more reliable context information for well-estimated joints. Finally, we build a Dense Upsampling Convolution (DUC) module to generate more precise prediction, which can recover missing joint details that are usually unavailable in common upsampling process. As a result, the predicted keypoint heatmaps are more accurate. Comprehensive experiments demonstrate that the proposed approach outperforms the state-of-the-art methods and achieves the superior performance over three benchmark datasets: the recent big dataset CrowdPose, the COCO keypoint detection dataset and the MPII Human Pose dataset. Our code will be released upon acceptance.


title: Finite metric and $k$-metric bases on ultrametric spaces

abstract: Given a metric space $(X,d)$, a set $S\subseteq X$ is called a $k$-\emph{metric generator} for $X$ if any pair of different points of $X$ is distinguished by at least $k$ elements of $S$. A $k$-\emph{metric basis} is a $k$-metric generator of the minimum cardinality in $X$. We prove that ultrametric spaces do not have finite $k$-metric bases for $k>2$. We also characterize when the metric and 2-metric bases of an ultrametric space are finite and, when they are finite, we characterize them. Finally, we prove that an ultrametric space can be easily recovered knowing only the metric basis and the coordinates of the points in it.


title: Carbon on the nanoscale: Ultrastiffness and unambiguous definition of incompressibility

abstract: In the presented work, the features of mechanical stiffness of carbon nanoparticles (nanodiamonds and fullerenes) in a wide range of sizes are considered. The enhancement of nanodiamonds stiffness (comparing to bulk diamond) is studied and explained in the terms of average bond stiffness $k_0$. It is shown that $k_0$ can be useful in the description of various carbon nanostructures and gives reliable estimates of their incompressibility. Moreover, we found that $k_0$ can be well estimated based only on relaxed atomic geometry.


title: Suppression of radiation loss in high kinetic inductance superconducting co-planar waveguides

abstract: We present a novel lab-on-chip technique to measure the very low losses in superconducting transmission lines at (sub-) mm wavelengths. The chips consist of a 100 nm thick NbTiN Coplanar Waveguide (CPW) Fabry Perot (FP) resonator, coupled on one side to an antenna and on the other side to a Microwave Kinetic Inductance detector. Using a single frequency radiation source allows us to measure the frequency response of the FP around 350 GHz and deduce its losses. We show that the loss is dominated by radiation loss inside the CPW line that forms the FP and that it decreases with decreasing line width and increasing kinetic inductance as expected. The results can be quantitatively understood using SONNET simulations. The lowest loss is observed for a CPW with a total width of $6\ \mathrm{\mu m}$ and corresponds to a Q-factor of $\approx15,000$.


title: Eigenrays in 3D heterogeneous anisotropic media: Part IV -- Geometric spreading from traveltime Hessian

abstract: Considering 3D heterogeneous anisotropic media, in Part III of this study we obtain the finite-element solution for the stationary ray path between two fixed endpoints using the finite element approach. It involves computation of the global (all-node) traveltime gradient vector and Hessian matrix, with respect to the nodal location and direction components. The global traveltime Hessian of the resolved stationary ray also plays an important role for obtaining the dynamic characteristics along the ray. It is a band matrix that includes the spatial, directional and mixed second derivatives of the traveltime at all nodes and nodal pairs. In Parts V, VI and VII of this study we explicitly use the global traveltime Hessian for solving the dynamic ray tracing equation, to obtain the geometric spreading at each point along the ray. Moreover, the dynamic ray tracing makes it possible to identify/classify caustics along the ray. In this part (Part IV) we propose an original two-stage approach for computing the relative geometric spreading of the entire ray path without explicitly performing the dynamic ray tracing. The first stage consists of an efficient algorithm for reducing the already computed global traveltime Hessian into a endpoint spatial traveltime Hessian. The second stage involves application of a known workflow for computing the geometrical spreading using the endpoint traveltime Hessian. Note that the method proposed in this part (Part IV) does not deliver information about caustics, nor the geometric spreading at intermediate points of the stationary ray path. Throughout the examples presented in Part VII we demonstrate the accuracy of the method proposed in this part over a set of isotropic and anisotropic examples by comparing the geometric spreading computed from this method based on the traveltime Hessian reduction with the one computed from the dynamic ray tracing.


title: Eigenrays in 3D heterogeneous anisotropic media: Part V -- Dynamics, Variational formulation

abstract: Parts V, VI and VII of this study are dedicated to the computation of paraxial rays and dynamic characteristics along the stationary rays obtained numerically in Part III. In this part (Part V), we formulate the linear, second-order, Jacobi dynamic ray tracing equation. In Part VI, we compare the Lagrangian and Hamiltonian approaches to the dynamic ray tracing in heterogeneous isotropic and anisotropic media; we demonstrate that the two approaches are compatible and derive the relationships between the Lagrangian's and Hamiltonian's Hessian matrices. In Part VII, we apply a similar finite-element solver, as used for the kinematic ray tracing, to compute the dynamic characteristics between the source and any point along the ray. The dynamic characteristics in our study include the relative geometric spreading and the phase correction factor due to caustics (i.e., the amplitude and the phase of the Green's function of waves propagating in 3D heterogeneous anisotropic media). The basic solution of the Jacobi equation is a shift vector of a paraxial ray in the plane normal to the ray direction at each point along the central ray. A general paraxial ray is defined by a linear combination of up to four basic vector solutions, each corresponds to specific initial conditions related to the ray coordinates at the source. We define the four basic solutions with two pairs of initial condition sets: point-source and plane-wave. For the proposed point-source ray coordinates and initial conditions, we derive the ray Jacobian and relate it to the relative geometric spreading for general anisotropy. Finally, we introduce a new dynamic parameter, referred as the normalized geometrical spreading, considered a measure of complexity of the propagated wave/ray phenomena, and hence, we propose using a criterion based on this parameter as a qualifying factor associated with the given ray solution.


title: Eigenrays in 3D heterogeneous anisotropic media: Part VI -- Dynamics, Lagrangian vs. Hamiltonian approaches

abstract: In Part V of this study, we presented an original Lagrangian approach for computing the dynamic characteristics along stationary rays, by solving the linear, second-order Jacobi differential equation, considering four sets of initial conditions as the basic solutions. We then focused on the computation of the geometric spreading and identification of caustics, where only the two point-source basic solutions with their corresponding initial conditions are required. Solutions of the Jacobi equation represent the normal shift vectors of the paraxial rays and define the geometry of the ray tube with respect to the stationary central ray. Rather than the Lagrangian approach, the dynamic characteristics are traditionally computed with the Hamiltonian approach, formulated normally in terms of two first-order differential equations, where the solution variables are the paraxial shifts and paraxial slowness changes along the ray. In this part (Part VI), we compare and relate the two approaches. We first combine the two first-order Hamiltonian dynamic equations, eliminating the paraxial variations of the slowness vector. This leads to a second-order differential equation in terms of the Hamiltonian shift alone, whose ray-normal counterpart coincides with the normal shift of the corresponding Lagrangian solution, while the ray-tangent component does not affect the Jacobian and the geometric spreading. Comparing the proposed Lagrangian approach to the dynamic ray tracing with the "classical" Hamiltonian approach, we demonstrate that they are fully compatible for general anisotropy. We then derive the two-way relationships between the Hamiltonian's and Lagrangian's Hessians, which are the core computational elements of dynamic ray theory. Finally, we demonstrate the relationships between these two types of the Hessians numerically, for a general triclinic medium.


title: Dichotomy between Deterministic and Probabilistic Models in Countably Additive Effectus Theory

abstract: Effectus theory is a relatively new approach to categorical logic that can be seen as an abstract form of generalized probabilistic theories (GPTs). While the scalars of a GPT are always the real unit interval [0,1], in an effectus they can form any effect monoid. Hence, there are quite exotic effectuses resulting from more pathological effect monoids. In this paper we introduce sigma-effectuses, where certain countable sums of morphisms are defined. We study in particular sigma-effectuses where unnormalized states can be normalized. We show that a non-trivial sigma-effectus with normalization has as scalars either the two-element effect monoid 0,1 or the real unit interval [0,1]. When states and/or predicates separate the morphisms we find that in the 0,1 case the category must embed into the category of sets and partial functions (and hence the category of Boolean algebras), showing that it implements a deterministic model, while in the [0,1] case we find it embeds into the category of Banach order-unit spaces and of Banach pre-base-norm spaces (satisfying additional properties), recovering the structure present in GPTs. Hence, from abstract categorical and operational considerations we find a dichotomy between deterministic and convex probabilistic models of physical theories.


title: Nanofocusing optics for an X-ray free-electron laser generating an extreme intensity of 100 EW/cm$^2$ using total reflection mirrors

abstract: A nanofocusing optical system referred to as $\textit{100 exa}$ for an X-ray free-electron laser (XFEL) was developed to generate an extremely high intensity of 100 EW/cm$^2$ (10$^2$$^0$ W/cm$^2$) using total reflection mirrors. The system is based on Kirkpatrick-Baez geometry, with 250 mm long elliptically figured mirrors optimized for the SPring-8 Angstrom Compact Free-Electron Laser (SACLA) XFEL facility. The nano-precision surface employed is coated with rhodium and offers a high reflectivity of 80%, with a photon energy of up to 12 keV, under total reflection conditions. Incident X-rays on the optics are reflected with a large spatial acceptance of over 900 $\mu$m. The focused beam is 210 nm $\times$ 120 nm (full width at half maximum) and was evaluated at a photon energy of 10 keV. The optics developed for $\textit{100 exa}$ efficiently achieved an intensity of 1 $\times$ 10$^2$$^0$ W/cm$^2$ with a pulse duration of 7 fs and a pulse energy of 150 $\mu$J (25% of the pulse energy generated at the light source). The experimental chamber, which can provide varied stage arrangements and sample conditions, including vacuum environments and atmospheric pressure helium, was set up with the focusing optics to meet the experimental requirements.


title: Linear Time-Varying MPC for Nonprehensile Object Manipulation with a Nonholonomic Mobile Robot

abstract: This paper proposes a technique to manipulate an object with a nonholonomic mobile robot by pushing, which is a nonprehensile manipulation motion primitive. Such a primitive involves unilateral constraints associated with the friction between the robot and the manipulated object. Violating this constraint produces the slippage of the object during the manipulation, preventing the correct achievement of the task. A linear time-varying model predictive control is designed to include the unilateral constraint within the control action properly. The approach is verified in a dynamic simulation environment through a Pioneer 3-DX wheeled robot executing the pushing manipulation of a package.


title: Wise Sliding Window Segmentation: A classification-aided approach for trajectory segmentation

abstract: Large amounts of mobility data are being generated from many different sources, and several data mining methods have been proposed for this data. One of the most critical steps for trajectory data mining is segmentation. This task can be seen as a pre-processing step in which a trajectory is divided into several meaningful consecutive sub-sequences. This process is necessary because trajectory patterns may not hold in the entire trajectory but on trajectory parts. In this work, we propose a supervised trajectory segmentation algorithm, called Wise Sliding Window Segmentation (WS-II). It processes the trajectory coordinates to find behavioral changes in space and time, generating an error signal that is further used to train a binary classifier for segmenting trajectory data. This algorithm is flexible and can be used in different domains. We evaluate our method over three real datasets from different domains (meteorology, fishing, and individuals movements), and compare it with four other trajectory segmentation algorithms: OWS, GRASP-UTS, CB-SMoT, and SPD. We observed that the proposed algorithm achieves the highest performance for all datasets with statistically significant differences in terms of the harmonic mean of purity and coverage.


title: Using Deep Reinforcement Learning Methods for Autonomous Vessels in 2D Environments

abstract: Unmanned Surface Vehicles technology (USVs) is an exciting topic that essentially deploys an algorithm to safely and efficiently performs a mission. Although reinforcement learning is a well-known approach to modeling such a task, instability and divergence may occur when combining off-policy and function approximation. In this work, we used deep reinforcement learning combining Q-learning with a neural representation to avoid instability. Our methodology uses deep q-learning and combines it with a rolling wave planning approach on agile methodology. Our method contains two critical parts in order to perform missions in an unknown environment. The first is a path planner that is responsible for generating a potential effective path to a destination without considering the details of the root. The latter is a decision-making module that is responsible for short-term decisions on avoiding obstacles during the near future steps of USV exploitation within the context of the value function. Simulations were performed using two algorithms: a basic vanilla vessel navigator (VVN) as a baseline and an improved one for the vessel navigator with a planner and local view (VNPLV). Experimental results show that the proposed method enhanced the performance of VVN by 55.31 on average for long-distance missions. Our model successfully demonstrated obstacle avoidance by means of deep reinforcement learning using planning adaptive paths in unknown environments.


title: Plasma Surface Metallurgy of Materials Based on Double Glow Discharge Phenomenon

abstract: Plasma Surface Metallurgy/Alloying is a kind of surface metallurgy/alloying to employ low temperature plasma produced by glow discharge to diffuse alloying elements into the surface of substrate material to form an alloy layer. The first plasma surface metallurgy technology is plasma nitriding invented by German scientist Dr. Bernard Berghuas in 1930. He was the first person to apply glow discharge to realize the surface alloying. In order to break the limitation of plasma nitriding technology, which can only be applied to a few non-metallic gaseous elements such as nitrogen, carbon, sulfur, the "Double Glow Discharge Phenomenon"was found in 1978. Based on this phenomenon the "Double Glow Plasma Surface Metallurgy Technology", also known as the "Xu-Tec Process" was invented in 1980. It can utilize any chemical elements in the periodic table including solid metallic, gas non-metallic elements and their combination to realize plasma surface alloying, hence greatly expanded the field of surface alloying. Countless surface alloys with high hardness, wear resistance and corrosion resistance, such as high speed steels, nickel base alloys and burn resistant alloys have been produced on the surfaces of a variety of materials. This technology may greatly improve the surface properties of metal materials, comprehensively improve the quality of mechanical products, save a lot of precious alloy elements for human beings. Based on the plasma nitriding technology, the Xu-Tec Process has opened up a new material engineering field of "Plasma Surface Metallurgy". This Review Article briefly presents the history of glow discharge and surface alloying, double glow discharge phenomenon, basic principle and current status of Double Glow Plasma Surface Metallurgy/Alloying. Industrial applications, advantages and future potential of the Xu-Tec process are also presented.


title: Coweight lattice $A^*_n$ and lattice simplices

abstract: There exist as many index-$k$ sublattices of the hexagonal lattice up to isometry as there exist lattice triangles with normalized volume $k$ up to unimodular equivalence, which can be explained using orbifolds. In dimension 3, it was noted that the number of sublattices of the fcc and the bcc lattices and the number of lattice tetrahedra all seem to be the same. We provide a bijection between the sublattices of the coweight lattice $A^*_n$ and the $n$-dimensional lattice simplices. It explains, proves, and generalizes the observed coincidences to arbitrary dimension.


title: A note on the ternary purely exponential Diophantine equation $A^x+B^y=C^z$ with $A+B=C^2$

abstract: Let $\ell, m, r$ be fixed positive integers such that $2\nmid \ell$, $3\nmid \ell m$, $\ell>r$ and $3\mid r$. In this paper, using the BHV theorem on the existence of primitive divisors of Lehmer numbers, we prove that if $\min\{r\ell m^2-1,(\ell-r)\ell m^2+1\}>30$, then the equation $(r\ell m^2-1)^x+((\ell -r)\ell m^2+1)^y=(\ell m)^z$ has only the positive integer solution $(x,y,z)=(1,1,2)$.


title: Non-linear dynamical response of interacting bosons to synthetic electric field

abstract: We theoretically study the non-linear response of interacting neutral bosonic gas in a synthetically driven one-dimensional optical lattice. In particular, we examine the bosonic analogue of electronic higher harmonic generation in a strong time-dependent synthetic vector potential manifesting itself as the synthetic electric field. We show that the vector potential can generate reasonably high harmonics in the insulating regime, while the superfluid regime exhibits only a few harmonics. In the insulating regime, the number of harmonics increases with the increase in the strength of the vector potential. This originates primarily due to the field-driven resonant and non-resonant excitations in the neutral Mott state and their recombination with the ground state. If the repulsive interaction between two atoms ($U$) is close to the strength of the gauge potential ($A_0$), the resonant quasiparticle-quasihole pairs on nearest-neighbor sites, namely dipole states are found to a play a dominant role in the generating higher harmonics. However, in the strong-field limit $A_0\gg U$, the nonresonant states where quasiparticle-quasihole pairs are not on nearest-neighbor sites give rise to higher harmonics.


title: Fast Cross-domain Data Augmentation through Neural Sentence Editing

abstract: Data augmentation promises to alleviate data scarcity. This is most important in cases where the initial data is in short supply. This is, for existing methods, also where augmenting is the most difficult, as learning the full data distribution is impossible. For natural language, sentence editing offers a solution - relying on small but meaningful changes to the original ones. Learning which changes are meaningful also requires large amounts of training data. We thus aim to learn this in a source domain where data is abundant and apply it in a different, target domain, where data is scarce - cross-domain augmentation. We create the Edit-transformer, a Transformer-based sentence editor that is significantly faster than the state of the art and also works cross-domain. We argue that, due to its structure, the Edit-transformer is better suited for cross-domain environments than its edit-based predecessors. We show this performance gap on the Yelp-Wikipedia domain pairs. Finally, we show that due to this cross-domain performance advantage, the Edit-transformer leads to meaningful performance gains in several downstream tasks.


title: Characterizing uninorms on bounded lattices

abstract: This paper establishes some equivalent conditions of a uninorm, extending an arbitrary triangular norm on [0, e] or an arbitrary triangular conorm on [e, 1] to the whole lattice.


title: Self-similar solutions of the spherically symmetric Euler equations for general equations of state

abstract: The study of spherically symmetric motion is important for the theory of explosion waves. In this paper, we construct rigorously self-similar solutions to the Riemann problem of the spherically symmetric Euler equations for general equations of state. We used the assumption of self-similarity to reduce the spherically symmetric Euler equations to a system of nonlinear ordinary differential equations, from which we obtain detailed structures of solutions besides their existence.


title: Grant-Free Non-Orthogonal Multiple Access: A Key Enabler for 6G-IoT

abstract: The proliferating number of devices with short payloads as well as low power budget has already driven researchers away from classical grant-based access schemes that are notorious for their large signalling overhead as well as power-consuming retransmissions. Instead, light-weight random access protocols have been re-investigated and their throughput has been improved in orders of magnitude with sophisticated yet still low-complex transceiver algorithms. In fact, grant-free access has been identified as a key medium access control technique for providing massive connectivity in machine type communications in cellular networks. In this paper, we show that grant-free access combined with non-orthogonal transmission schemes is a promising solution for 6G Internet of Things (IoT). We present novel and promising results for deep learning (DL)-based techniques for joint user detection and decoding. Then, we propose a multi-layered model for GF-NOMA for power-efficient communications. We also discuss resource allocation issues to enable the co-existence of GF-NOMA with other orthogonal or even grant-based schemes. Finally, we conclude with proposed research directions for medium access towards enabling 6G-IoT.


title: Sample-Specific Output Constraints for Neural Networks

abstract: Neural networks reach state-of-the-art performance in a variety of learning tasks. However, a lack of understanding the decision making process yields to an appearance as black box. We address this and propose ConstraintNet, a neural network with the capability to constrain the output space in each forward pass via an additional input. The prediction of ConstraintNet is proven within the specified domain. This enables ConstraintNet to exclude unintended or even hazardous outputs explicitly whereas the final prediction is still learned from data. We focus on constraints in form of convex polytopes and show the generalization to further classes of constraints. ConstraintNet can be constructed easily by modifying existing neural network architectures. We highlight that ConstraintNet is end-to-end trainable with no overhead in the forward and backward pass. For illustration purposes, we model ConstraintNet by modifying a CNN and construct constraints for facial landmark prediction tasks. Furthermore, we demonstrate the application to a follow object controller for vehicles as a safety-critical application. We submitted an approach and system for the generation of safety-critical outputs of an entity based on ConstraintNet at the German Patent and Trademark Office with the official registration mark DE10 2019 119 739.


title: Spatio-temporal filtering in laser Doppler holography for retinal blood flow imaging

abstract: Laser Doppler holography (LDH) is a full-field interferometric imaging technique recently applied in ophthalmology to measure blood flow, a parameter of high clinical interest. From the temporal fluctuations of digital holograms acquired at ultrafast frame rates, LDH reveals retinal and choroidal blood flow with a few milliseconds of temporal resolution. However LDH experiences difficulties to detect slower blood flow as it requires to work with low Doppler frequency shifts which are corrupted by eye motion. We here demonstrate the use of a spatio-temporal decomposition adapted from Doppler ultrasound that provides a basis appropriate to the discrimination of blood flow from eye motion. A singular value decomposition (SVD) can be used as a simple, robust, and efficient way to separate the Doppler fluctuations of blood flow from those of strong spatial coherence such as eye motion. We show that the SVD outperforms the conventional Fourier based filter to reveal slower blood flow, and dramatically improves the ability of LDH to reveal vessels of smaller size or with a pathologically reduced blood flow.


title: Radial Star Formation Histories in 32 Nearby Galaxies

abstract: The spatially resolved star formation histories are studied for 32 normal star-forming galaxies drawn from the the Spitzer Extended Disk Galaxy Exploration Science survey. At surface brightness sensitivities fainter than 28 mag arcsec$^{-2}$, the new optical photometry is deep enough to complement archival ultraviolet and infrared imaging and to explore the properties of the emission well beyond the traditional optical extents of these nearby galaxies. Fits to the spectral energy distributions using a delayed star formation history model indicate a subtle but interesting average radial trend for the spiral galaxies: the inner stellar systems decrease in age with increasing radius, consistent with inside-out disk formation, but the trend reverses in the outermost regions with the stellar age nearly as old as the innermost stars. These results suggest an old stellar outer disk population formed through radial migration and/or the cumulative history of minor mergers and accretions of satellite dwarf galaxies. The subset of S0 galaxies studied here show the opposite trend compared to what is inferred for spirals: characteristic stellar ages that are increasingly older with radius for the inner portions of the galaxies, and increasingly younger stellar ages for the outer portions. This result suggests that either S0 galaxies are not well modeled by a delayed-$\tau$ model, and/or that S0 galaxies have a more complicated formation history than spiral galaxies.


title: Distributed projected-reflected-gradient algorithms for stochastic generalized Nash equilibrium problems

abstract: We consider the stochastic generalized Nash equilibrium problem (SGNEP) with joint feasibility constraints and expected-value cost functions. We propose a distributed stochastic projected reflected gradient algorithm and show its almost sure convergence when the pseudogradient mapping is monotone and the solution is unique. The algorithm is based on monotone operator splitting methods tailored for SGNEPs when the expected-value pseudogradient mapping is approximated at each iteration via an increasing number of samples of the random variable. Finally, we show that a preconditioned variant of our proposed algorithm has convergence guarantees when the pseudogradient mapping is cocoercive.


title: Search for non-thermal X-ray emission in the colliding wind binary Cyg OB2 #8A

abstract: Cyg OB2 #8a is a massive O-type binary displaying strong non-thermal radio emission. Owing to the compactness of this binary, emission of non-thermal X-ray photons via inverse Compton scattering is expected. We first revised the orbital solution for Cyg OB2 #8a using new optical spectra. We then reduced and analysed X-ray spectra obtained with XMM-Newton, Swift, INTEGRAL, and NuSTAR. The analysis of the XMM-Newton and Swift data allows us to better characterise the X-ray emission from the stellar winds and colliding winds region at energies below 10 keV. We confirm the variation of the broad-band light curve of Cyg OB2 #8a along the orbit with, for the first time, the observation of the maximum emission around phase 0.8. The minimum ratio of the X-ray to bolometric flux of Cyg OB2 #8a remains well above the level expected for single O-type stars, indicating that the colliding wind region is not disrupted during the periastron passage. The analysis of the full set of publicly available INTEGRAL observations allows us to refine the upper limit on the non-thermal X-ray flux of the Cyg OB2 region between 20 and 200 keV. Two NuSTAR observations (phases 0.028 and 0.085) allow us to study the Cyg OB2 #8a spectrum up to 30 keV. These data do not provide evidence of the presence of non-thermal X-rays, but bring more stringent constraints on the flux of a putative non-thermal component. Finally, we computed, thanks to a new dedicated model, the anisotropic inverse Compton emission generated in the wind shock region. The theoretical non-thermal emission appears to be compatible with observational limits and the kinetic luminosity computed from these models is in good agreement with the unabsorbed flux observed below 10 keV.


title: Anti M-Weierstrass function sequences

abstract: Large algebraic structures are found inside the space of sequences of continuous functions on a compact interval having the property that, the series defined by each sequence converges absolutely and uniformly on the interval but the series of the upper bounds diverges. So showing that there exist many examples satisfying the conclusion but not the hypothesis of the Weierstrass M-test.


title: A numerical approach for fluid deformable surfaces

abstract: Fluid deformable surfaces show a solid-fluid duality which establishes a tight interplay between tangential flow and surface deformation. We derive the governing equations as a thin film limit and provide a general numerical approach for their solution. The simulation results demonstrate the rich dynamics resulting from this interplay, where in the presence of curvature any shape change is accompanied by a tangential flow and, vice versa, the surface deforms due to tangential flow. However, they also show that the only possible stable stationary state in the considered setting is a sphere with zero velocity.


title: Phase transitions in perturbative walking dynamics

abstract: In this paper, we investigate the dynamics of the confinement-deconfinement phase transition in a toy model where the walking dynamics is realized perturbatively. We study the properties of the phase transition focusing on the possible cosmological signatures it can provide. Interestingly the model is well under perturbative control only when the mass of the lightest field - the dilaton/scalon is much lighter than the rest of the fields and the phase transition proceeds slowly leading to strong signals in the stochastic gravitational wave spectrum.


title: A nearly-neutral biallelic Moran model with biased mutation and linear and quadratic selection

abstract: In this article, a biallelic reversible mutation model with linear and quadratic selection is analyzed. The approach reconnects to one proposed by Kimura ( Possibility of extensive neutral evolution under stabilizing selection with special reference to nonrandom use of codons (PNAS,1981)), who starts from a diffusion model and derives its equilibrium distribution up to a constant. We use a boundary-mutation Moran model, which approximates a general mutation model for small effective mutation rates, and derive its equilibrium distribution for polymorphic and monomorphic variants in small to moderately sized populations. Using this model, we show that biased mutation rates and linear selection alone can cause patterns of polymorphism rates within and substitution rates between populations that are usually ascribed to balancing or overdominant selection. We illustrate this using a data set of short introns and fourfold degenerate sites from Drosophila simulans and Drosophila melanogaster.


title: Associated Derived Invariants for Geometric Mappings of Non-Symmetric Affine Connection Spaces

abstract: The invariants of the Thomas and the Weyl type for a mapping between non-symmetric affine connection spaces are obtained with respect to the factored deformation tensor in this paper. Motivated by two invariants of the Weyl type obtained in (N. O. Vesi\'c, Basic Invariants of Geometric Mappings, [17]), we founded novel invariants of the Weyl type. Invariants for almost geodesic mappings of the third type are searched at the end of this paper.


title: A method to identify geochemical mineralization on linear transect

abstract: Mineral exploration in biogeochemistry is related to the detection of anomalies in soil, which is driven by many factors and thus a complex problem. Mik\v{s}ov\'a, Rieser, and Filzmoser (2019) have introduced a method for the identification of spatial patterns with increased element concentrations in samples along a linear sampling transect. This procedure is based on fitting Generalized Additive Models (GAMs) to the concentration data, and computing a curvature measure from the pairwise log-ratios of these fits. The higher the curvature, the more likely one or both elements of the pair indicate local mineralization. This method is applied on two geochemical data sets which have been collected specifically for the purpose of mineral exploration. The aim is to test the technique for its ability to identify pathfinder elements to detect mineralized zones, and to verify whether the method can indicate which sampling material is best suited for this purpose. Reference: Mik\v{s}ov\'a D., Rieser C., Filzmoser P. (2019). "Identification of mineralization in geochemistry along a transect based on the spatial curvature of log-ratios." arXiv, (1912.02867).


title: A Block Coordinate Descent-based Projected Gradient Algorithm for Orthogonal Non-negative Matrix Factorization

abstract: This article utilizes the projected gradient method (PG) for a non-negative matrix factorization problem (NMF), where one or both matrix factors must have orthonormal columns or rows. We penalise the orthonormality constraints and apply the PG method via a block coordinate descent approach. This means that at a certain time one matrix factor is fixed and the other is updated by moving along the steepest descent direction computed from the penalised objective function and projecting onto the space of non-negative matrices. Our method is tested on two sets of synthetic data for various values of penalty parameters. The performance is compared to the well-known multiplicative update (MU) method from Ding (2006), and with a modified global convergent variant of the MU algorithm recently proposed by Mirzal (2014). We provide extensive numerical results coupled with appropriate visualizations, which demonstrate that our method is very competitive and usually outperforms the other two methods.


title: Mobility-aware Beam Steering in Metasurface-based Programmable Wireless Environments

abstract: Programmable wireless environments (PWEs) utilize electromagnetic metasurfaces to transform wireless propagation into a software-controlled resource. In this work we study the effects of user device mobility on the efficiency of PWEs. An analytical model is proposed, which describes the potential misalignment between user-emitted waves and the active PWE configuration, and can constitute the basis for studying queuing problems in PWEs. Subsequently, a novel, beam steering approach is proposed which can effectively mitigate the misalignment effects. Ray-tracing-based simulations evaluate the proposed scheme.


title: A Nonconvex Low-Rank Tensor Completion Model for Spatiotemporal Traffic Data Imputation

abstract: Sparsity and missing data problems are very common in spatiotemporal traffic data collected from various sensing systems. Making accurate imputation is critical to many applications in intelligent transportation systems. In this paper, we formulate the missing data imputation problem in spatiotemporal traffic data in a low-rank tensor completion (LRTC) framework and define a novel truncated nuclear norm (TNN) on traffic tensors of location$\times$day$\times$time of day. In particular, we introduce an universal rate parameter to control the degree of truncation on all tensor modes in the proposed LRTC-TNN model, and this allows us to better characterize the hidden patterns in spatiotemporal traffic data. Based on the framework of the Alternating Direction Method of Multipliers (ADMM), we present an efficient algorithm to obtain the optimal solution for each variable. We conduct numerical experiments on four spatiotemporal traffic data sets, and our results show that the proposed LRTC-TNN model outperforms many state-of-the-art imputation models with missing rates/patterns. Moreover, the proposed model also outperforms other baseline models in extreme missing scenarios.


title: Stellar wind effects on the atmospheres of close-in giants: a possible reduction in escape instead of increased erosion

abstract: The atmospheres of highly irradiated exoplanets are observed to undergo hydrodynamic escape. However, due to strong pressures, stellar winds can confine planetary atmospheres, reducing their escape. Here, we investigate under which conditions atmospheric escape of close-in giants could be confined by the large pressure of their host star's winds. For that, we simulate escape in planets at a range of orbital distances ([0.04, 0.14] au), planetary gravities ([36%, 87%] of Jupiter's gravity), and ages ([1, 6.9] Gyr). For each of these simulations, we calculate the ram pressure of these escaping atmospheres and compare them to the expected stellar wind external pressure to determine whether a given atmosphere is confined or not. We show that, although younger close-in giants should experience higher levels of atmospheric escape, due to higher stellar irradiation, stellar winds are also stronger at young ages, potentially reducing escape of young exoplanets. Regardless of the age, we also find that there is always a region in our parameter space where atmospheric escape is confined, preferably occurring at higher planetary gravities and orbital distances. We investigate confinement of some known exoplanets and find that the atmosphere of several of them, including pi Men c, should be confined by the winds of their host stars, thus potentially preventing escape in highly irradiated planets. Thus, the lack of hydrogen escape recently reported for pi Men c could be caused by the stellar wind.


title: On the minimal set of controllers and sensors for linear power flow

abstract: We consider a linear power flow model with interval-bounded nodal power injections and limited line power flows. We determine the minimal number of power injections to control based on a minimal set of measurements, such that the overall system is feasible for all assignments of the non-controlled power injections. For the important case where the possible measurements are the nodal power injections, we show that the problem can be solved efficiently as a mixed-integer linear program (MILP). When also line power flows are considered as potential measurements, we derive an iterative, greedy algorithm that provides a feasible, but potentially conservative solution. We apply the developed algorithms to both a small microgrid and a modified version of the IEEE 118 bus test power system. We show that in both cases a sparse solution in terms of the number of required controllers and measurements can be obtained. Moreover, the number of required measurements can be reduced significantly if line flow measurements are considered additionally to nodal power injections.


title: Second-order Post-Minkowskian scattering in arbitrary dimensions

abstract: We extract the long-range gravitational potential between two scalar particles with arbitrary masses from the two-to-two elastic scattering amplitude at 2nd Post-Minkowskian order in arbitrary dimensions. In contrast to the four-dimensional case, in higher dimensions the classical potential receives contributions from box topologies. Moreover, the kinematical relation between momentum and position on the classical trajectory contains a new term which is quadratic in the tree-level amplitude. A precise interplay between this new relation and the formula for the scattering angle ensures that the latter is still linear in the classical part of the scattering amplitude, to this order, matching an earlier calculation in the eikonal approach. We point out that both the eikonal exponentiation and the reality of the potential to 2nd post-Minkowskian order can be seen as a consequence of unitarity. We finally present closed-form expressions for the scattering angle given by leading-order gravitational potentials for dimensions ranging from four to ten.


title: Cross-domain Object Detection through Coarse-to-Fine Feature Adaptation

abstract: Recent years have witnessed great progress in deep learning based object detection. However, due to the domain shift problem, applying off-the-shelf detectors to an unseen domain leads to significant performance drop. To address such an issue, this paper proposes a novel coarse-to-fine feature adaptation approach to cross-domain object detection. At the coarse-grained stage, different from the rough image-level or instance-level feature alignment used in the literature, foreground regions are extracted by adopting the attention mechanism, and aligned according to their marginal distributions via multi-layer adversarial learning in the common feature space. At the fine-grained stage, we conduct conditional distribution alignment of foregrounds by minimizing the distance of global prototypes with the same category but from different domains. Thanks to this coarse-to-fine feature adaptation, domain knowledge in foreground regions can be effectively transferred. Extensive experiments are carried out in various cross-domain detection scenarios. The results are state-of-the-art, which demonstrate the broad applicability and effectiveness of the proposed approach.


title: Double-EIT Ground-State Cooling of Stationary Two-Dimensional Ion Lattices

abstract: We theoretically and experimentally investigate double electromagnetically induced transparency (double-EIT) cooling of two-dimensional ion crystals confined in a Paul trap. The double-EIT ground-state cooling is observed for \Yb ions with clock state, for which EIT cooling has not been realized like many other ions with a simple $\Lambda$-scheme. A cooling rate of $\dot{\bar n}=34~(\pm1.8)~\rm{ms}^{-1}$ and a cooling limit of $\bar n=0.06~(\pm 0.059)$ are observed for a single ion. The measured cooling rate and limit are consistent with theoretical predictions. We apply double-EIT cooling to the transverse modes of two-dimensional (2D) crystals with up to 12 ions. In our 2D crystals, the micromotion and the transverse mode directions are perpendicular, which makes them decoupled. Therefore, the cooling on transverse modes is not disturbed by micromotion, which is confirmed in our experiment. For the center of mass mode of a 12 ions crystal, we observe a cooling rate and a cooling limit that are consistent with those of a single ion, including heating rates proportional to the number of ions. This method can be extended to other hyperfine qubits, and near ground-state cooling of stationary 2D crystals with large numbers of ions may advance the field of quantum information sciences.


title: On the diameter of the polytope of the stable marriage with ties

abstract: The stable marriage problem with ties is a well-studied and interesting problem in game theory. We are given a set of men and a set of women. Each individual has a preference ordering on the opposite group, which can possibly contain ties. A stable marriage is given by a matching between men and women for which there is no blocking pair, i.e., a men and a women who strictly prefer each other to their current partner in the matching. In this paper, we study the diameter of the polytope given by the convex hull of characteristic vectors of stable marriages, in the setting with ties. We prove an upper bound of $\lfloor \frac{n}{3}\rfloor$ on the diameter, where $n$ is the total number of men and women, and give a family of instances for which the bound holds tight. Our result generalizes the bound on the diameter of the standard stable marriage polytope (i.e., the well-known polytope that describes the setting without ties), developed previously in the literature.


title: Using variability and VLBI to measure cosmological distances

abstract: In this paper, we propose a new approach to determining cosmological distances to active galactic nuclei (AGN) via light travel-time arguments, which can be extended from nearby sources to very high redshift sources. The key assumption is that the variability seen in AGN is constrained by the speed of light and therefore provides an estimate of the linear size of an emitting region. This can then be compared with the angular size measured with very long baseline interferometry (VLBI) in order to derive a distance. We demonstrate this approach on a specific well studied low redshift (z = 0.0178) source 3C84 (NGC 1275), which is the bright radio core of the Perseus Cluster. We derive an angular diameter distance including statistical errors of $D_{A} = 72^{+5}_{-6}$ Mpc for this source, which is consistent with other distance measurements at this redshift. Possible sources of systematic errors and ways to correct for them are discussed.


title: IETF protocol suite for the Internet of Things: Overview and Recent Advancements

abstract: Internet of Things (IoT) is a rapidly growing technological domain and the effort that many Standards Developing Organizations (SDOs) and alliances are dedicating to it is constantly increasing. The Internet Engineering Task Force (IETF) is certainly one of the most active SDOs operating in this broad domain. With the aim of building a comprehensive, interoperable and streamlined IoT stack, several IETF working groups are developing protocols, in multiple technological domains, which are already nowadays very relevant to the IoT. This article gives a concise but comprehensive survey of the IETF efforts in the IoT. We discuss mainstream standardization and research activities, as well as other IETF-connected supporting initiatives provided still in the context of IoT.


title: Graph Neural Networks for Decentralized Controllers

abstract: Dynamical systems comprised of autonomous agents arise in many relevant problems such as multi-agent robotics, smart grids, or smart cities. Controlling these systems is of paramount importance to guarantee a successful deployment. Optimal centralized controllers are readily available but face limitations in terms of scalability and practical implementation. Optimal decentralized controllers, on the other hand, are difficult to find. In this paper, we propose a framework using graph neural networks (GNNs) to learn decentralized controllers from data. While GNNs are naturally distributed architectures, making them perfectly suited for the task, we adapt them to handle delayed communications as well. Furthermore, they are equivariant and stable, leading to good scalability and transferability properties. The problem of flocking is explored to illustrate the potential of GNNs in learning decentralized controllers.


title: Accurate Optimization of Weighted Nuclear Norm for Non-Rigid Structure from Motion

abstract: Fitting a matrix of a given rank to data in a least squares sense can be done very effectively using 2nd order methods such as Levenberg-Marquardt by explicitly optimizing over a bilinear parameterization of the matrix. In contrast, when applying more general singular value penalties, such as weighted nuclear norm priors, direct optimization over the elements of the matrix is typically used. Due to non-differentiability of the resulting objective function, first order sub-gradient or splitting methods are predominantly used. While these offer rapid iterations it is well known that they become inefficent near the minimum due to zig-zagging and in practice one is therefore often forced to settle for an approximate solution. In this paper we show that more accurate results can in many cases be achieved with 2nd order methods. Our main result shows how to construct bilinear formulations, for a general class of regularizers including weighted nuclear norm penalties, that are provably equivalent to the original problems. With these formulations the regularizing function becomes twice differentiable and 2nd order methods can be applied. We show experimentally, on a number of structure from motion problems, that our approach outperforms state-of-the-art methods.


title: BVI-CC: A Dataset for Research on Video Compression and Quality Assessment

abstract: The video technology scenery has been very vivid over the past years, with novel video coding technologies introduced that promise improved compression performance over state-of-the-art technologies. Despite the fact that a lot of video datasets are available, representative content of the wide parameter space along with subjective evaluations of variations of encoded content from an unpartial end is required. In response to this requirement, this paper features a dataset, the BVI-CC. Three video codecs were deployed to create the variations of the encoded sequences: High Efficiency Video Coding (HEVC) Test Model (HM), AOMedia Video 1 (AV1), and Versatile Video Coding (VVC) Test Model (VTM). Nine source video sequences were carefully selected to offer both diversity and representativeness in the spatio-temporal domain. Different spatial resolution versions of the sequences were created and encoded by all three codecs at pre-defined target bit rates. The compression efficiency of the codecs was evaluated with commonly used objective quality metrics, and the subjective quality of their reconstructed content was also evaluated through psychophysical experiments. Additionally, an adaptive bit rate (convex hull rate-distortion optimization across spatial resolutions) test case was assessed using both objective and subjective evaluations. Finally, the computational complexities of the tested codecs were examined. All data have been made publicly available as part of the dataset, which can be used for coding performance evaluation and video quality metric development.


title: Breaking BEC

abstract: In this work quantum corrections to the classical evolution of a relativistic scalar condensate are studied. The problem is approached by means of two different perturbative approaches: the 2-particle-irreducible (2PI) effective action and the expansion in the self-coupling. In the weak coupling regime, the decoherence of the classical state is observed. The corresponding timescale is identified with the quantum break-time.


title: Relativistic stars in 4D Einstein-Gauss-Bonnet gravity

abstract: In the present paper we investigate the structure of relativistic stars in 4D Einstein-Gauss-Bonnet gravity. The mass-radius relations are obtained for realistic hadronic and for strange quark star equations of state, and for a wide range of the Gauss-Bonnet coupling parameter $\alpha$. Even though the deviations from general relativity for nonzero values of $\alpha$ can be large enough, they are still comparable with the variations due to different modern realistic equations of state if we restrict ourselves to moderate value of $\alpha$. That is why the current observations of the neutron star masses and radii alone can not impose stringent constraints on the value of the parameter $\alpha$. Nevertheless some rough constraints on $\alpha$ can be put. The existence of stable stellar mass black holes imposes $\sqrt{\alpha}\lesssim 2.6 {\rm km}$ for $\alpha>0$ while the requirement that the maximum neutron star mass should be greater than two solar masses gives $\sqrt{|\alpha|}\lesssim 3.9 {\rm km}$ for $\alpha<0$. We also present an exact solution describing the structure of relativistic stars with uniform energy density in 4D Einstein-Gauss-Bonnet gravity.


title: Projected Estimation for Large-dimensional Matrix Factor Models

abstract: In this study, we propose a projection estimation method for large-dimensional matrix factor models with cross-sectionally spiked eigenvalues. By projecting the observation matrix onto the row or column factor space, we simplify factor analysis for matrix series to that for a lower-dimensional tensor. This method also reduces the magnitudes of the idiosyncratic error components, thereby increasing the signal-to-noise ratio, because the projection matrix linearly filters the idiosyncratic error matrix. We theoretically prove that the projected estimators of the factor loading matrices achieve faster convergence rates than existing estimators under similar conditions. Asymptotic distributions of the projected estimators are also presented. A novel iterative procedure is given to specify the pair of row and column factor numbers. Extensive numerical studies verify the empirical performance of the projection method. Two real examples in finance and macroeconomics reveal factor patterns across rows and columns, which coincides with financial, economic, or geographical interpretations.


title: PathVQA: 30000+ Questions for Medical Visual Question Answering

abstract: Is it possible to develop an "AI Pathologist" to pass the board-certified examination of the American Board of Pathology? To achieve this goal, the first step is to create a visual question answering (VQA) dataset where the AI agent is presented with a pathology image together with a question and is asked to give the correct answer. Our work makes the first attempt to build such a dataset. Different from creating general-domain VQA datasets where the images are widely accessible and there are many crowdsourcing workers available and capable of generating question-answer pairs, developing a medical VQA dataset is much more challenging. First, due to privacy concerns, pathology images are usually not publicly available. Second, only well-trained pathologists can understand pathology images, but they barely have time to help create datasets for AI research. To address these challenges, we resort to pathology textbooks and online digital libraries. We develop a semi-automated pipeline to extract pathology images and captions from textbooks and generate question-answer pairs from captions using natural language processing. We collect 32,799 open-ended questions from 4,998 pathology images where each question is manually checked to ensure correctness. To our best knowledge, this is the first dataset for pathology VQA. Our dataset will be released publicly to promote research in medical VQA.


title: A Time Series Approach To Player Churn and Conversion in Videogames

abstract: Players of a free-to-play game are divided into three main groups: non-paying active users, paying active users and inactive users. A State Space time series approach is then used to model the daily conversion rates between the different groups, i.e., the probability of transitioning from one group to another. This allows, not only for predictions on how these rates are to evolve, but also for a deeper understanding of the impact that in-game planning and calendar effects have. It is also used in this work for the detection of marketing and promotion campaigns about which no information is available. In particular, two different State Space formulations are considered and compared: an Autoregressive Integrated Moving Average process and an Unobserved Components approach, in both cases with a linear regression to explanatory variables. Both yield very close estimations for covariate parameters, producing forecasts with similar performances for most transition rates. While the Unobserved Components approach is more robust and needs less human intervention in regards to model definition, it produces significantly worse forecasts for non-paying user abandonment probability. More critically, it also fails to detect a plausible marketing and promotion campaign scenario.


title: Energy Optimization in Binary Star Systems: Explanation for Equal Mass Members in Close Orbits

abstract: Observations indicate that members of close stellar binaries often have mass ratios close to unity, while longer-period systems exhibit a more uniform mass-ratio distribution. This paper provides a theoretical explanation for this finding by determining the tidal equilibrium states for binary star systems --- subject to the constraints of conservation of angular momentum and constant total mass. This work generalizes previous treatments by including the mass fraction as a variable in the optimization problem. The results show that the lowest energy state accessible to the system corresponds to equal mass stars on a circular orbit, where the stellar spin angular velocities are both synchronized and aligned with the orbit. These features are roughly consistent with observed properties of close binary systems. We also find the conditions required for this minimum energy state to exist: [1] The total angular momentum must exceed a critical value, [2] the orbital angular momentum must be three times greater than the total spin angular momentum, and [3] the semimajor axis is bounded from above. The last condition implies that sufficiently wide binaries are not optimized with equal mass stars, where the limiting binary separation occurs near $a_0\approx16R_\ast$.


title: Embedding technique and network analysis of scientific innovations emergence in an arXiv-based concept network

abstract: Novelty is an inherent part of innovations and discoveries. Such processes may be considered as an appearance of new ideas or as an emergence of atypical connections between the existing ones. The importance of such connections hints for investigation of innovations through network or graph representation in the space of ideas. In such representation, a graph node corresponds to the relevant concept (idea), whereas an edge between two nodes means that the corresponding concepts have been used in a common context. In this study we address the question about a possibility to identify the edges between existing concepts where the innovations may emerge. To this end, we use a well-documented scientific knowledge landscape of 1.2M arXiv.org manuscripts dated starting from April 2007 and until September 2019. We extract relevant concepts for them using the ScienceWISE.info platform. Combining approaches developed in complex networks science and graph embedding, we discuss the predictability of edges (links) on the scientific knowledge landscape where the innovations may appear.


title: Performance Analysis of Millimeter Wave Wireless Power Transfer With Imperfect Beam Alignment

abstract: In this paper, the impact of imperfect beam alignment (IBA) on millimeter wave (mmWave) wireless power transfer (WPT) is investigated. We consider a mmWave WPT network, where the location of the energy transmitters follows a Poisson point process. Instead of the mostly used flat-top antenna model, we adopt the Gaussian antenna model suggested by the 3rd Generation Partnership Project (3GPP) for better accuracy. Two beam alignment error (BAE) models, i.e., truncated Gaussian and uniform models, are introduced to represent different BAE sources. We derive the probability density functions (PDFs) of the cascaded antenna gain with both BAE models and then provide the approximated PDFs for tractability. With the help of Fox's H function, the analytic expression for the energy coverage probability with nonlinear energy harvesting model is derived. Besides, we deduce a closed-form expression of the average harvested radio frequency (RF) energy. The simulation results verify our theoretical results and demonstrate the performance degradation incurred by BAE. It also shows that the flat-top antenna model cannot always provide accurate performance evaluation in the presence of BAE.


title: Diffusive dynamics of elongated particles in active colloidal suspensions of motile algae

abstract: Swimming microorganisms can influence the diffusion of passive particles. The effect of this swimmer-particle interaction depends on different properties, such as the hydrodynamic field of the swimmer and the relative sizes of microorganisms and particles. We investigated an enhancement of the diffusion of silica doublets in a suspension of microalgae Chlamydomonas reinhardtii in a flat capillary. Depending on the concentration of microswimmers, the translational and rotational diffusion constants increase by several orders of magnitude in the presence of the swimming algae. For low concentrations of algae, the doublets exhibit Brownian motion in a fluctuating flow field generated by multiple swimmers. One can observe strong, diffusive transport caused by occasional large displacements. At high swimmer concentration, the algae form dense clusters, where the rotational motion of the doublets shows a subdiffusive behaviour while the translational motion remains diffusive.


title: Observation of Time-Reversal Invariant Helical Edge-Modes in Bilayer Graphene/WSe$_2$ Heterostructure

abstract: Topological insulators, along with Chern insulators and Quantum Hall insulator phases, are considered as paradigms for symmetry protected topological phases of matter. This article reports the experimental realization of the time-reversal invariant helical edge-modes in bilayer graphene/monolayer WSe$_2$-based heterostructures -- a phase generally considered as a precursor to the field of generic topological insulators. Our observation of this elusive phase depended crucially on our ability to create mesoscopic devices comprising both a moir\'e superlattice potential and strong spin-orbit coupling; this resulted in materials whose electronic band structure could be tuned from trivial to topological by an external displacement field. We find that the topological phase is characterized by a bulk bandgap and by helical edge-modes with electrical conductance quantized exactly to $2e^2/h$ in zero external magnetic field. We put the helical edge-modes on firm grounds through supporting experiments, including the verification of predictions of the Landauer-B$\mathrm{\ddot{u}}$ttiker model for quantum transport in multi-terminal mesoscopic devices. Our non-local transport properties measurements show that the helical edge-modes are dissipationless and equilibrate at the contact probes. We achieved the tunability of the different topological phases with electric and magnetic fields, which allowed us to achieve topological phase transitions between trivial and multiple, distinct topological phases. We also present results of a theoretical study of a realistic model which, in addition to replicating our experimental results, explains the origin of the topological insulating bulk and helical edge-modes. Our experimental and theoretical results establish a viable route to realizing the time-reversal invariant $\mathbb{Z}_2$ topological phase of matter.


title: Dissecting the $\Delta I= 1/2$ rule at large $N_c$

abstract: We study the scaling of kaon decay amplitudes with the number of colours, $N_c$, in a theory with four degenerate flavours, $N_f=4$. In this scenario, two current-current operators, $Q^\pm$, mediate $\Delta S=1$ transitions, such as the two isospin amplitudes of non-leptonic kaon decays for $K\to (\pi\pi)_{I=0,2}$, $A_0$ and $A_2$. In particular, we concentrate on the simpler $K\to\pi$ amplitudes, $A^\pm$, mediated by these two operators. A diagrammatic analysis of the large-$N_c$ scaling of these observables is presented, which demonstrates the anticorrelation of the leading ${\mathcal O}(1/N_c)$ and ${\mathcal O}(N_f/N_c^2)$ corrections in both amplitudes. Using our new $N_f=4$ and previous quenched data, we confirm this expectation and show that these corrections are $naturally$ large and may be at the origin of the $\Delta I=1/2$ rule. The evidence for the latter is indirect, based on the matching of the amplitudes to their prediction in Chiral Perturbation Theory, from which the LO low-energy couplings of the chiral weak Hamiltonian, $g^\pm$, can be determined. A NLO estimate of the $K \to (\pi\pi)_{I=0,2}$ isospin amplitudes can then be derived, which is in good agreement with the experimental value.


title: Optimising Game Tactics for Football

abstract: In this paper we present a novel approach to optimise tactical and strategic decision making in football (soccer). We model the game of football as a multi-stage game which is made up from a Bayesian game to model the pre-match decisions and a stochastic game to model the in-match state transitions and decisions. Using this formulation, we propose a method to predict the probability of game outcomes and the payoffs of team actions. Building upon this, we develop algorithms to optimise team formation and in-game tactics with different objectives. Empirical evaluation of our approach on real-world datasets from 760 matches shows that by using optimised tactics from our Bayesian and stochastic games, we can increase a team chances of winning by up to 16.1\% and 3.4\% respectively.


title: Interdisciplinarity metric based on the co-citation network

abstract: Quantifying the interdisciplinarity of a research is a relevant problem in the evaluative bibliometrics. The concept of interdisciplinarity is ambiguous and multidimensional. Thus, different measures of interdisciplinarity have been propose in the literature. However, few studies have proposed interdisciplinary metrics without previously defining classification sets, and no one use the co-citation network for this purpose. In this study we propose an interdisciplinary metric based on the co-citation network. This is a way to define the publication's field without resorting to pre-defined classification sets. We present a characterization of a publication's field and then we use this definition to propose a new metric of the interdisciplinarity degree for publications (papers) and journals as units of analysis. The proposed measure has an aggregative property that makes it scalable from a paper individually to a set of them (journal) without more than adding the numerators and denominators in the proportions that define this new indicator. Moreover, the aggregated value of two or more units is strictly among all the individual values.


title: Adaptive Name Entity Recognition under Highly Unbalanced Data

abstract: For several purposes in Natural Language Processing (NLP), such as Information Extraction, Sentiment Analysis or Chatbot, Named Entity Recognition (NER) holds an important role as it helps to determine and categorize entities in text into predefined groups such as the names of persons, locations, quantities, organizations or percentages, etc. In this report, we present our experiments on a neural architecture composed of a Conditional Random Field (CRF) layer stacked on top of a Bi-directional LSTM (BI-LSTM) layer for solving NER tasks. Besides, we also employ a fusion input of embedding vectors (Glove, BERT), which are pre-trained on the huge corpus to boost the generalization capacity of the model. Unfortunately, due to the heavy unbalanced distribution cross-training data, both approaches just attained a bad performance on less training samples classes. To overcome this challenge, we introduce an add-on classification model to split sentences into two different sets: Weak and Strong classes and then designing a couple of Bi-LSTM-CRF models properly to optimize performance on each set. We evaluated our models on the test set and discovered that our method can improve performance for Weak classes significantly by using a very small data set (approximately 0.45\%) compared to the rest classes.


title: An approach to parameter identifiability for a class of nonlinear models represented in LPV form

abstract: In several model-based system maintenance problems, parameters are used to represent unknown characteristics of a component, equipment degradation, etc. This allows for modelling constant, slow-varying terms. The identifiability of these parameters is an important condition to estimate them. Linear Parameter Varying (LPV) models are being increasingly used in the industries as a bridge between linear and nonlinear models. Techniques exist that can rewrite some nonlinear models in LPV form. However, the problem of identifiability of these models is still at a nascent stage. In this paper, we propose an approach to verify identifiability of unknown parameters for LPV state-space models. It makes use of a parity-space like formulation to eliminate the states of the model. The resulting input-output-parameter equation is analysed to verify the identifiability of the original model or a subset of unknown parameters. This approach provides a framework for both continuous-time and discrete-time models and we illustrate it using examples.


title: Analysis of a semi-implicit structure-preserving finite element method for the nonstationary incompressible Magnetohydrodynamics equations

abstract: We revise the structure-preserving finite element method in [K. Hu, Y. MA and J. Xu. (2017) Stable finite element methods preserving $\nabla \cdot \mathbf{B}=0$ exactly for MHD models. Numer. Math.,135, 371-396]. The revised method is semi-implicit in time-discretization. We prove the linearized scheme preserves the divergence free property for the magnetic field exactly at each time step. Further, we showed the linearized scheme is unconditionally stable and we obtain optimal convergence in the energy norm of the revised method even for solutions with low regularity.


title: Robust Medical Instrument Segmentation Challenge 2019

abstract: Intraoperative tracking of laparoscopic instruments is often a prerequisite for computer and robotic-assisted interventions. While numerous methods for detecting, segmenting and tracking of medical instruments based on endoscopic video images have been proposed in the literature, key limitations remain to be addressed: Firstly, robustness, that is, the reliable performance of state-of-the-art methods when run on challenging images (e.g. in the presence of blood, smoke or motion artifacts). Secondly, generalization; algorithms trained for a specific intervention in a specific hospital should generalize to other interventions or institutions. In an effort to promote solutions for these limitations, we organized the Robust Medical Instrument Segmentation (ROBUST-MIS) challenge as an international benchmarking competition with a specific focus on the robustness and generalization capabilities of algorithms. For the first time in the field of endoscopic image processing, our challenge included a task on binary segmentation and also addressed multi-instance detection and segmentation. The challenge was based on a surgical data set comprising 10,040 annotated images acquired from a total of 30 surgical procedures from three different types of surgery. The validation of the competing methods for the three tasks (binary segmentation, multi-instance detection and multi-instance segmentation) was performed in three different stages with an increasing domain gap between the training and the test data. The results confirm the initial hypothesis, namely that algorithm performance degrades with an increasing domain gap. While the average detection and segmentation quality of the best-performing algorithms is high, future research should concentrate on detection and segmentation of small, crossing, moving and transparent instrument(s) (parts).


title: A 75kb SRAM in 65nm CMOS for In-Memory Computing Based Neuromorphic Image Denoising

abstract: This paper presents an in-memory computing (IMC) architecture for image denoising. The proposed SRAM based in-memory processing framework works in tandem with approximate computing on a binary image generated from neuromorphic vision sensors. Implemented in TSMC 65nm process, the proposed architecture enables approximately 2000X energy savings (approximately 222X from IMC) compared to a digital implementation when tested with the video recordings from a DAVIS sensor and achieves a peak throughput of 1.25-1.66 frames/us.


title: $\Lambda$-submodules of finite index of anticyclotomic plus and minus Selmer groups of elliptic curves

abstract: Let $p$ be an odd prime and $K$ an imaginary quadratic field where $p$ splits. Under appropriate hypotheses, Bertolini showed that the Selmer group of a $p$-ordinary elliptic curve over the anticyclotomic $\mathbb Z_p$-extension of $K$ does not admit any proper $\Lambda$-submodule of finite index, where $\Lambda$ is a suitable Iwasawa algebra. We generalize this result to the plus and minus Selmer groups (in the sense of Kobayashi) of $p$-supersingular elliptic curves. In particular, in our setting the plus/minus Selmer groups have $\Lambda$-corank one, so they are not $\Lambda$-cotorsion. As an application of our main theorem, we prove results in the vein of Greenberg-Vatsal on Iwasawa invariants of $p$-congruent elliptic curves, extending to the supersingular case results for $p$-ordinary elliptic curves due to Hatley-Lei.


title: The collision frequency of electron-neutral-particle in the weakly ionized plasma with the power-law velocity distribution

abstract: We study the collision frequency of electron-neutral-particle in the weakly ionized plasma with the power-law velocity q-distribution and derive the formulation of the average collision frequency. We find that the average collision frequency in the q-distributed plasma also depends strongly on the q-parameter and thus is generally different from that in the Maxwell-distributed plasma, which therefore modifies the transport coefficients in the previous studies of the weakly ionized plasmas with the power-law velocity distributions.


title: From Bit To Bedside: A Practical Framework For Artificial Intelligence Product Development In Healthcare

abstract: Artificial Intelligence (AI) in healthcare holds great potential to expand access to high-quality medical care, whilst reducing overall systemic costs. Despite hitting the headlines regularly and many publications of proofs-of-concept, certified products are failing to breakthrough to the clinic. AI in healthcare is a multi-party process with deep knowledge required in multiple individual domains. The lack of understanding of the specific challenges in the domain is, therefore, the major contributor to the failure to deliver on the big promises. Thus, we present a decision perspective framework, for the development of AI-driven biomedical products, from conception to market launch. Our framework highlights the risks, objectives and key results which are typically required to proceed through a three-phase process to the market launch of a validated medical AI product. We focus on issues related to Clinical validation, Regulatory affairs, Data strategy and Algorithmic development. The development process we propose for AI in healthcare software strongly diverges from modern consumer software development processes. We highlight the key time points to guide founders, investors and key stakeholders throughout their relevant part of the process. Our framework should be seen as a template for innovation frameworks, which can be used to coordinate team communications and responsibilities towards a reasonable product development roadmap, thus unlocking the potential of AI in medicine.


title: Attention U-Net Based Adversarial Architectures for Chest X-ray Lung Segmentation

abstract: Chest X-ray is the most common test among medical imaging modalities. It is applied for detection and differentiation of, among others, lung cancer, tuberculosis, and pneumonia, the last with importance due to the COVID-19 disease. Integrating computer-aided detection methods into the radiologist diagnostic pipeline, greatly reduces the doctors' workload, increasing reliability and quantitative analysis. Here we present a novel deep learning approach for lung segmentation, a basic, but arduous task in the diagnostic pipeline. Our method uses state-of-the-art fully convolutional neural networks in conjunction with an adversarial critic model. It generalized well to CXR images of unseen datasets with different patient profiles, achieving a final DSC of 97.5% on the JSRT dataset.


title: Twisted Hochschild Homology of Quantum Flag Manifolds and K\"ahler Forms

abstract: We study the twisted Hochschild homology of quantum flag manifolds, the twist being the modular automorphism of the Haar state. We prove that every quantum flag manifold admits a non-trivial class in degree two, with an explicit representative defined in terms of a certain projection. The corresponding classical two-form, via the Hochschild-Kostant-Rosenberg theorem, is identified with a K\"ahler form on the flag manifold.


title: Safe Crossover of Neural Networks Through Neuron Alignment

abstract: One of the main and largely unexplored challenges in evolving the weights of neural networks using genetic algorithms is to find a sensible crossover operation between parent networks. Indeed, naive crossover leads to functionally damaged offspring that do not retain information from the parents. This is because neural networks are invariant to permutations of neurons, giving rise to multiple ways of representing the same solution. This is often referred to as the competing conventions problem. In this paper, we propose a two-step safe crossover(SC) operator. First, the neurons of the parents are functionally aligned by computing how well they correlate, and only then are the parents recombined. We compare two ways of measuring relationships between neurons: Pairwise Correlation (PwC) and Canonical Correlation Analysis (CCA). We test our safe crossover operators (SC-PwC and SC-CCA) on MNIST and CIFAR-10 by performing arithmetic crossover on the weights of feed-forward neural network pairs. We show that it effectively transmits information from parents to offspring and significantly improves upon naive crossover. Our method is computationally fast,can serve as a way to explore the fitness landscape more efficiently and makes safe crossover a potentially promising operator in future neuroevolution research and applications.


title: Partitions of non-complete Baire metric spaces

abstract: We investigate the properties of ideals associated with Kuratowski partitions of non-complete Baire metric spaces. We show that such an ideal can be precipitous.


title: A Developmental Neuro-Robotics Approach for Boosting the Recognition of Handwritten Digits

abstract: Developmental psychology and neuroimaging research identified a close link between numbers and fingers, which can boost the initial number knowledge in children. Recent evidence shows that a simulation of the children's embodied strategies can improve the machine intelligence too. This article explores the application of embodied strategies to convolutional neural network models in the context of developmental neuro-robotics, where the training information is likely to be gradually acquired while operating rather than being abundant and fully available as the classical machine learning scenarios. The experimental analyses show that the proprioceptive information from the robot fingers can improve network accuracy in the recognition of handwritten Arabic digits when training examples and epochs are few. This result is comparable to brain imaging and longitudinal studies with young children. In conclusion, these findings also support the relevance of the embodiment in the case of artificial agents' training and show a possible way for the humanization of the learning process, where the robotic body can express the internal processes of artificial intelligence making it more understandable for humans.


title: Distributed Gradient Methods for Nonconvex Optimization: Local and Global Convergence Guarantees

abstract: The article discusses distributed gradient-descent algorithms for computing local and global minima in nonconvex optimization. For local optimization, we focus on distributed stochastic gradient descent (D-SGD)--a simple network-based variant of classical SGD. We discuss local minima convergence guarantees and explore the simple but critical role of the stable-manifold theorem in analyzing saddle-point avoidance. For global optimization, we discuss annealing-based methods in which slowly decaying noise is added to D-SGD. Conditions are discussed under which convergence to global minima is guaranteed. Numerical examples illustrate the key concepts in the paper.


title: Dynamical Environment and Surface Characteristics of Asteroid (16) Psyche

abstract: Radar observations show that (16) Psyche is one of the largest and most massive asteroids of the M-class located in the main belt, with a diameter of approximately 230 km. This fact makes Psyche a unique object since observations indicated an iron-nickel composition. It is believed that this body may be what was left of a metal core of an early planet that would have been fragmented over millions of years due to violent collisions. In this work we study a variety of dynamical aspects related to the surface, as well as, the environment around this asteroid. We use computational tools to explore the gravitational field generated by this body, assuming constant values for its density and rotation period. We then determine a set of physical and dynamical characteristics over its entire surface. The results include the geometric altitude, geopotential altitude, tilt, slope, among others. We also explore the neighborhood around the asteroid (16) Psyche, so that the location and linear stability of the equilibrium points were found. We found four external equilibrium points, two of them linearly stable. We confirmed the stability of these points by performing numerical simulations of massless particles around the asteroid, which also showed an asymmetry in the size of the stable regions. In addition, we integrate a cloud of particles in the vicinity of (16) Psyche in order to verify in which regions of its surface the particles are most likely to collide.


title: Galactic Center IRS13E: Colliding Stellar Winds or an Intermediate Mass Black Hole?

abstract: A small cluster of massive stars residing in the Galactic center, collectively known as IRS13E, is of special interest due to its close proximity to Sgr A* and the possibility that an embedded intermediate-mass black hole (IMBH) binds its member stars. It has been suggested that colliding winds from two member stars, both classified as Wolf-Rayet type, are responsible for the observed X-ray, infrared and radio emission from IRS13E. We have conducted an in-depth study of the X-ray spatial, temporal and spectral properties of IRS13E, based on 5.6 Ms of ultra-deep Chandra observations obtained over 20 years. These X-ray observations show no significant evidence for source variability. We have also explored the kinematics of the cluster members, using Keck near-infrared imaging and spectroscopic data on a 14-yr baseline that considerably improve the accuracy of stars' proper motions. The observations are interpreted using 3-dimensional hydrodynamical simulations of colliding winds tailored to match the physical conditions of IRS13E, leading us to conclude that the observed X-ray spectrum and morphology can be well explained by the colliding wind scenario, in the meantime offering no support for the presence of a putative IMBH. An IMBH more massive than a few $10^3{\rm~M_\odot}$ is also strongly disfavored by the stellar kinematics.


title: A termination criterion for stochastic gradient descent for binary classification

abstract: We propose a new, simple, and computationally inexpensive termination test for constant step-size stochastic gradient descent (SGD) applied to binary classification on the logistic and hinge loss with homogeneous linear predictors. Our theoretical results support the effectiveness of our stopping criterion when the data is Gaussian distributed. This presence of noise allows for the possibility of non-separable data. We show that our test terminates in a finite number of iterations and when the noise in the data is not too large, the expected classifier at termination nearly minimizes the probability of misclassification. Finally, numerical experiments indicate for both real and synthetic data sets that our termination test exhibits a good degree of predictability on accuracy and running time.


title: Asymmetric Partisan Voter Turnout Games

abstract: Since Downs proposed that the act of voting is irrational in 1957, myriad models have been proposed to explain voting and account for observed turnout patterns. We propose a model in which partisans consider both the instrumental and expressive benefits of their vote when deciding whether or not to abstain in an election, introducing an asymmetry that most other models do not consider. Allowing learning processes within our electorate, we analyze what turnout states are rationalizable under various conditions. Our model predicts comparative statics that are consistent with voter behavior. Furthermore, relaxing some of our preliminary assumptions eliminates some of the discrepancies between our model and empirical voter behavior.


title: A remnant planetary core in the hot-Neptune desert

abstract: The interiors of giant planets remain poorly understood. Even for the planets in the Solar System, difficulties in observation lead to large uncertainties in the properties of planetary cores. Exoplanets that have undergone rare evolutionary processes provide a route to understanding planetary interiors. Planets found in and near the typically barren hot-Neptune 'desert' (a region in mass-radius space that contains few planets) have proved to be particularly valuable in this regard. These planets include HD149026b, which is thought to have an unusually massive core, and recent discoveries such as LTT9779b and NGTS-4b, on which photoevaporation has removed a substantial part of their outer atmospheres. Here we report observations of the planet TOI-849b, which has a radius smaller than Neptune's but an anomalously large mass of $39.1^{+2.7}_{-2.6}$ Earth masses and a density of $5.2^{+0.7}_{-0.8}$ grams per cubic centimetre, similar to Earth's. Interior structure models suggest that any gaseous envelope of pure hydrogen and helium consists of no more than $3.9^{+0.8}_{-0.9}$ per cent of the total planetary mass. The planet could have been a gas giant before undergoing extreme mass loss via thermal self-disruption or giant planet collisions, or it could have avoided substantial gas accretion, perhaps through gap opening or late formation. Although photoevaporation rates cannot account for the mass loss required to reduce a Jupiter-like gas giant, they can remove a small (a few Earth masses) hydrogen and helium envelope on timescales of several billion years, implying that any remaining atmosphere on TOI-849b is likely to be enriched by water or other volatiles from the planetary interior. We conclude that TOI-849b is the remnant core of a giant planet.


title: Adversarial Attacks on Monocular Depth Estimation

abstract: Recent advances of deep learning have brought exceptional performance on many computer vision tasks such as semantic segmentation and depth estimation. However, the vulnerability of deep neural networks towards adversarial examples have caused grave concerns for real-world deployment. In this paper, we present to the best of our knowledge the first systematic study of adversarial attacks on monocular depth estimation, an important task of 3D scene understanding in scenarios such as autonomous driving and robot navigation. In order to understand the impact of adversarial attacks on depth estimation, we first define a taxonomy of different attack scenarios for depth estimation, including non-targeted attacks, targeted attacks and universal attacks. We then adapt several state-of-the-art attack methods for classification on the field of depth estimation. Besides, multi-task attacks are introduced to further improve the attack performance for universal attacks. Experimental results show that it is possible to generate significant errors on depth estimation. In particular, we demonstrate that our methods can conduct targeted attacks on given objects (such as a car), resulting in depth estimation 3-4x away from the ground truth (e.g., from 20m to 80m).


title: Extended Existence Probability Using Digital Maps for Object Verification

abstract: A main task for automated vehicles is an accurate and robust environment perception. Especially, an error-free detection and modeling of other traffic participants is of great importance to drive safely in any situation. For this purpose, multi-object tracking algorithms, based on object detections from raw sensor measurements, are commonly used. However, false object hypotheses can occur due to a high density of different traffic participants in complex, arbitrary scenarios. For this reason, the presented approach introduces a probabilistic model to verify the existence of a tracked object. Therefore, an object verification module is introduced, where the influences of multiple digital map elements on a track's existence are evaluated. Finally, a probabilistic model fuses the various influences and estimates an extended existence probability for every track. In addition, a Bayes Net is implemented as directed graphical model to highlight this work's expandability. The presented approach, reduces the number of false positives, while retaining true positives. Real world data is used to evaluate and to highlight the benefits of the presented approach, especially in urban scenarios.


title: Path Integral Molecular Dynamics for Fermions: Alleviating the Sign Problem with the Bogoliubov Inequality

abstract: We present a method for performing path integral molecular dynamics (PIMD) simulations for fermions and address its sign problem. PIMD simulations are widely used for studying many-body quantum systems at thermal equilibrium. However, they assume that the particles are distinguishable and neglect bosonic and fermionic exchange effects. Interacting fermions play a key role in many chemical and physical systems, such as electrons in quantum dots and ultracold trapped atoms. A direct sampling of the fermionic partition function is impossible using PIMD since its integrand is not positive definite. We show that PIMD simulations for fermions are feasible by employing our recently developed method for bosonic PIMD and reweighting the results to obtain fermionic expectation values. The approach is tested against path integral Monte Carlo (PIMC) simulations for up to 7 electrons in a two-dimensional quantum dot for a range of interaction strengths. However, like PIMC, the method suffers from the sign problem at low temperatures. We propose a simple approach for alleviating it by simulating an auxiliary system with a larger average sign and obtaining an upper bound to the energy of the original system using the Bogoliubov inequality. This allows fermions to be studied at temperatures lower than would otherwise have been feasible using PIMD, as demonstrated in the case of a three-electron quantum dot. Our results extend the boundaries of PIMD simulations of fermions and will hopefully stimulate the development of new approaches for tackling the sign problem.


title: Microlensing of radially pulsating stars

abstract: Here, we study the microlensing of radially pulsating stars. Discerning and characterizing the properties of distant, faint pulsating stars is achievable through high-cadence microlensing observations. Combining stellar variability period with microlensing gives the source distance, type, and radius and helps better determine the lens parameters. Considering periodically variations in their radius and surface temperature, their microlensing light curves are resulted from multiplication of the magnification factor with variable finite size effect by the intrinsic brightness curves of pulsing source. The variable finite source size due to pulsation can be significant for transit and single microlensing and while caustic-crossing features. This kind of deviation in the magnification factor is considerable when the ratio of the source radius to the projected lens-source distance is in the range of $\rho_{\star}/u \in[0.4,10]$ and its duration is short and in the same order of the time of crossing the source radius. Other deviations due to variable source intensity and its area make colored and periodic deviations which are asymmetric with respect to the signs of pulsation phase. The positive phases makes deviations with larger amplitude that negative phase. These deviations dominate in filters with short wave lengths (e.g., $B-$band). The position of magnification peaks in microlensing of variable stars varies and this displacement differs in different filters.


title: The Magellan/PFS Exoplanet Search: A 55-day period dense Neptune transiting the bright ($V=8.6$) star HD 95338

abstract: We report the detection of a transiting, dense Neptune planet candidate orbiting the bright ($V=8.6$) K0.5V star HD 95338. Detection of the 55-day periodic signal comes from the analysis of precision radial velocities from the Planet Finder Spectrograph on the Magellan II Telescope. Follow-up observations with HARPS also confirm the presence of the periodic signal in the combined data. HD 95338 was also observed by the Transiting Exoplanet Survey Satellite ({\it TESS}) where we identify a clear single transit in the photometry. A Markov Chain Monte Carlo period search on the velocities allows strong constraints on the expected transit time, matching well the epoch calculated from \tess{} data, confirming both signals describe the same companion. A joint fit model yields an absolute mass of 42.44$^{+2.22}_{-2.08} M_{\oplus}$ and a radius of 3.89$^{+0.19}_{-0.20}$ $R_{\oplus}$ which translates to a density of 3.98$^{+0.62}_{-0.64}$ \gcm\, for the planet. Given the planet mass and radius, structure models suggest it is composed of a mixture of ammonia, water, and methane. HD 95338\,b is one of the most dense Neptune planets yet detected, indicating a heavy element enrichment of $\sim$90\% ($\sim38\, M_{\oplus}$). This system presents a unique opportunity for future follow-up observations that can further constrain structure models of cool gas giant planets.