Search

Noncollinear magnets are notoriously difficult to describe within first-principles approaches based on density-functional theory (DFT) because of the presence of low-lying spin excitations. At the level of ground-state calculations, several methods exist to constrain the magnetic moments to a predetermined configuration, and thereby accelerate convergence towards self-consistency. Their use in a perturbative context, however, remains very limited. Here we present a general methodological framework to achieve parametric control over the local spin moments at the linear-response level. Our strategy builds on the concept of Legendre transform to switch between various flavors of magnetic functionals, and to relate their second derivatives via simple linear-algebra operations. Thereby, we can address an arbitrary response function at the time-dependent DFT level of theory with optimal accuracy and minimal computational effort. In the low frequency limit, we identify the leading correction to the existing adiabatic formulation of the problem [S. Ren \emph{et al.}, Phys. Rev. X {\bf 14}, 011041 (2024)], consisting in a renormalization of the phonon and magnon masses due to electron inertia. As a demonstration, we apply our methodology to the study of the THz optical response in bulk CrI$_3$, where we identify a hybrid electromagnon with mixed spin-lattice character.

The attainment of accurate numerical solutions of ill-conditioned linear algebraic problems involving totally positive matrices has been gathering considerable attention among researchers over the last years. In parallel, the interest of $q$-calculus has been steadily growing in the literature. In this work the $q$-analogue of the Abel polynomial basis is studied. The total positivity of the matrix of change of basis between monomial and $q$-Abel bases is characterized, providing its bidiagonal factorization. Moreover, well-known high relative accuracy results of Vandermonde matrices corresponding to increasing positive nodes are extended to the decreasing negative case. This further allows to solve with high relative accuracy several algebraic problems concerning collocation, Wronskian and Gramian matrices of $q$-Abel polynomials. Finally, a series of numerical tests support the presented theoretical results and illustrate the goodness of the method where standard approaches fail to deliver accurate solutions.

Ultrasound follicle tracking is an important part of cycle monitoring. OSIS Ovary (Online System for Image Segmentation for the Ovary) has been conceived aiming to aid the management of the workflow in follicle tracking, one of the most iterative procedures in cycle monitoring during ovarian stimulation. In the present study, we compared OSIS Ovary (as three-dimensional ultrasound-based automated system) with the two-dimensional manual standard measurement method, in order to assess the reliability of the main measurements obtained to track follicle growth during ovarian stimulation cycles, the follicle size and count. Based on the mean follicle diameter and follicle count values obtained, the Pearson/intraclass correlation coefficients were 0.976/0.987 and 0.804/0.889 in >=10mm follicles, 0.989/0.994 and 0.809/0.867 in >=13mm follicles and 0.995/0.997 and 0.791/0.840 in >=16mm follicles. The mean difference (MnD) for the mean diameter and follicle count was, respectively, 0.759/0.161 in >=10mm follicles, 0.486/1.033 in >=13mm follicles and 0.784/0.486 in >=16mm follicles. The upper and lower limits of agreement (ULA and LLA) were 3.641/2.123 and 5.392/3.070 in >=10mm follicles, 3.496/2.522 and 4.285/2.218 in >=13mm follicles, and 3.723/2.153 and 2.432/1.459 in >=16mm follicles. The limits of agreement range (LoAR) were 5.764/8.462 in >=10mm follicles, 6.048/6.503 in >=13mm follicles and 5.876/3.891 in >=16mm follicles. P<0.05 was considered for all calculations. As three-dimensional ultrasound-based automated system in comparison with two-dimensional manual method standard, we found OSIS Ovary as a reliable tool to track follicle growth during ovarian stimulation cycles, Comment: 21 pages, 4 figures, 2 tables

In this paper, we propose and analyze two different stream ciphers based on a Skew Tent Map and a Modified Logistic Map respectively. In order to improve the randomness of these systems, a single method for increasing the period length of the generated sequences has been applied. The results prove that the randomness of these systems can be severally increased by using this method, making these systems suitable for secure communications., Comment: Proceedings of 2016 12th Conference on Ph.D. Research in Microelectronics and Electronics (PRIME)

In recent years many efforts have been devoted to finding bidiagonal factorizations of nonsingular totally positive matrices, since their accurate computation allows to numerically solve several important algebraic problems with great precision, even for large ill-conditioned matrices. In this framework, the present work provides the factorization of the collocation matrices of Newton bases -- of relevance when considering the Lagrange interpolation problem -- together with an algorithm that allows to numerically compute it to high relative accuracy. This further allows to determine the coefficients of the interpolating polynomial and to compute the singular values and the inverse of the collocation matrix. Conditions that guarantee high relative accuracy for these methods and, in the former case, for the classical recursion formula of divided differences, are determined. Numerical errors due to imprecise computer arithmetic or perturbed input data in the computation of the factorization are analyzed. Finally, numerical experiments illustrate the accuracy and effectiveness of the proposed methods with several algebraic problems, in stark contrast with traditional approaches., Comment: Numer Algor (2024)

Enabling on-demand control of heat flow is key for the development of next-generation electronic devices, solid-state heat pumps, and thermal logic. However, precise and agile tuning of the relevant microscopic material parameters for adjusting thermal conductivities remains elusive. Here, we study several single crystals of lanthanum aluminate (LaAlO$_{3}$) with different domain structures and show that ferroelastic domain walls behave as boundaries that act like efficient controllers to govern thermal conductivity. At low temperature (3 K), we demonstrate a fivefold reduction in thermal conductivity induced by domain walls orthogonal to the heat flow and a twofold reduction when they are parallel to the heat flow. Atomistic calculations fully support this experimental observation. By breaking down phonon scattering mechanisms, we also analyze the temperature dependence of the thermal conductivity to derive a quantitative relation between thermal conductivity variations and domain wall organization and density.

Starting with a manifold $M$ and a semi-free action of $S^3$ on it, we have the Smith-Gysin sequence: $$ \cdots \to H^{*}( M) \to H^{*-3}(M/S^3, M^{S^3}) \oplus H^{*} (M^{S^3}) \to H^{*+1}(M/S^3, M^{S^3}) \to H^{*+1}(M) \to \cdots $$ In this paper, we construct a Smith-Gysin sequence that does not require the semi-free condition. This sequence includes a new term, referred to as the "exotic term," which depends on the subset $M^{S^1}$: $$ \cdots \to H^{*}(M) \to H^{*-3} (M/S^3, \Sigma/S^3) \oplus H^{*}(M^{S^3}) \oplus \left( H^{*-2}(M^{S^1})\right)^{-\mathbb{Z}_2} \to H^{*+1}(M/S^3,M^{S^3}) \to H^{*+1}(M) \to \cdots $$ Here, $\Sigma \subset M$ is the subset of points in $M$ whose isotropy groups are infinite. The group $\mathbb{Z}_2$ acts on $M^{S^1}$ by $j \in S^3$., Comment: We have changed the abstrast, which was unreadable

The Hard Lefschetz Property (HLP) has recently been formulated in the context of isometric flows without singularities on manifolds. In this category, two versions of the HLP (transverse and not) have been proven to be equivalent, thus generalizing what happens in the important cases of both K-contact and Sasakian manifolds. In this work we define both versions of the HLP for almost-free S3 -actions, and prove that they agree for actions satisfying a cohomological condition, which includes the important category of 3-Sasakian manifolds, where those two versions of the HLP are shown to be held. We also provide a family of examples of free actions of the 3-sphere which are not 3-Sasakian manifolds, but satisfy the HLP., Comment: 8 pages

Time-resolved illumination provides rich spatio-temporal information for applications such as accurate depth sensing or hidden geometry reconstruction, becoming a useful asset for prototyping and as input for data-driven approaches. However, time-resolved illumination measurements are high-dimensional and have a low signal-to-noise ratio, hampering their applicability in real scenarios. We propose a novel method to compactly represent time-resolved illumination using mixtures of exponentially-modified Gaussians that are robust to noise and preserve structural information. Our method yields representations two orders of magnitude smaller than discretized data, providing consistent results in applications such as hidden scene reconstruction and depth estimation, and quantitative improvements over previous approaches.

Non-line-of-sight (NLOS) imaging methods are capable of reconstructing complex scenes that are not visible to an observer using indirect illumination. However, they assume only third-bounce illumination, so they are currently limited to single-corner configurations, and present limited visibility when imaging surfaces at certain orientations. To reason about and tackle these limitations, we make the key observation that planar diffuse surfaces behave specularly at wavelengths used in the computational wave-based NLOS imaging domain. We call such surfaces virtual mirrors. We leverage this observation to expand the capabilities of NLOS imaging using illumination beyond the third bounce, addressing two problems: imaging single-corner objects at limited visibility angles, and imaging objects hidden behind two corners. To image objects at limited visibility angles, we first analyze the reflections of the known illuminated point on surfaces of the scene as an estimator of the position and orientation of objects with limited visibility. We then image those limited visibility objects by computationally building secondary apertures at other surfaces that observe the target object from a direct visibility perspective. Beyond single-corner NLOS imaging, we exploit the specular behavior of virtual mirrors to image objects hidden behind a second corner by imaging the space behind such virtual mirrors, where the mirror image of objects hidden around two corners is formed. No specular surfaces were involved in the making of this paper.

We predict a large in-plane polarization response to bending in a broad class of trigonal two-dimensional crystals. We define and compute the relevant flexoelectric coefficients from first principles as linear-response properties of the undistorted layer, by using the primitive crystal cell. The ensuing response (evaluated for SnS$_{2}$, silicene, phosphorene and RhI$_{3}$ monolayers and for a hexagonal BN bilayer) is up to one order of magnitude larger than the out-of-plane components in the same material. We illustrate the topological implications of our findings by calculating the polarization textures that are associated with a variety of rippled and bent structures. We also determine the longitudinal electric fields induced by a flexural phonon at leading order in amplitude and momentum.

Robots deployed to the real world must be able to interact with other agents in their environment. Dynamic game theory provides a powerful mathematical framework for modeling scenarios in which agents have individual objectives and interactions evolve over time. However, a key limitation of such techniques is that they require a-priori knowledge of all players' objectives. In this work, we address this issue by proposing a novel method for learning players' objectives in continuous dynamic games from noise-corrupted, partial state observations. Our approach learns objectives by coupling the estimation of unknown cost parameters of each player with inference of unobserved states and inputs through Nash equilibrium constraints. By coupling past state estimates with future state predictions, our approach is amenable to simultaneous online learning and prediction in receding horizon fashion. We demonstrate our method in several simulated traffic scenarios in which we recover players' preferences for, e.g., desired travel speed and collision-avoidance behavior. Results show that our method reliably estimates game-theoretic models from noise-corrupted data that closely matches ground-truth objectives, consistently outperforming state-of-the-art approaches., Comment: arXiv admin note: text overlap with arXiv:2106.03611

Given a smooth action of the sphere $\mathbb S^3$ on a manifold $M$, we have previously constructed a Gysin sequence relating the cohomology of the manifold $M$ and that of the orbit space $M/\mathbb S^3$. This sequence involves an exotic term depending on the subset $M^{\mathbb S^1}$. Notice that the orbit space is a stratified pseudomanifold, a kind of singular spaces where intersection cohomology applies. In the case where the the action is semi-free, the first author has already constructed a Gysin sequence relating the cohomology of $M$ and the intersection cohomology of $M/\mathbb S^3$. What happens if the action is not semi-free? This is the goal of this work. The situation is more complicated and we do not find a Gysin sequence but a Gysin braid relating the cohomology of $M$ and the intersection cohomology of $M/\mathbb S^3$. This braid also contains an exotic term depending this time on the intersection cohomology of $M^{\mathbb S^1}$.

The Observatorio Astrof\'isico de Javalambre (OAJ{\dag}1) in Spain is a young astronomical facility, conceived and developed from the beginning as a fully automated observatory with the main goal of optimizing the processes in the scientific and general operation of the Observatory. The OAJ has been particularly conceived for carrying out large sky surveys with two unprecedented telescopes of unusually large fields of view (FoV): the JST/T250, a 2.55m telescope of 3deg field of view, and the JAST/T80, an 83cm telescope of 2deg field of view. The most immediate objective of the two telescopes for the next years is carrying out two unique photometric surveys of several thousands square degrees, J-PAS{\dag}2 and J-PLUS{\dag}3, each of them with a wide range of scientific applications, like e.g. large structure cosmology and Dark Energy, galaxy evolution, supernovae, Milky Way structure, exoplanets, among many others. To do that, JST and JAST are equipped with panoramic cameras under development within the J-PAS collaboration, JPCam and T80Cam respectively, which make use of large format (~ 10k x 10k) CCDs covering the entire focal plane. This paper describes in detail, from operations point of view, a comparison between the detailed cost of the global automation of the Observatory and the standard automation cost for astronomical facilities, in reference to the total investment and highlighting all benefits obtained from this approach and difficulties encountered. The paper also describes the engineering development of the overall facilities and infrastructures for the fully automated observatory and a global overview of current status, pinpointing lessons learned in order to boost observatory operations performance, achieving scientific targets, maintaining quality requirements, but also minimizing operation cost and human resources., Comment: Global Observatory Control System GOCS

Ancient solutions of the Ricci flow arise naturally as models for singularity formation. There has been significant progress towards the classification of such solutions under natural geometric assumptions. Nonnegatively curved solutions in dimensions 2 and 3, and uniformly PIC solutions in higher dimensions are now well understood. We consider ancient solutions of arbitrary dimension which are complete and have Type~I curvature growth. We show that a $\kappa$-noncollapsed Type~I ancient solution which is noncompact and has nonnegative sectional curvature necessarily splits at least one Euclidean factor. It follows that a $\kappa$-noncollapsed Type~I ancient solution which is weakly PIC2 is a locally symmetric space.

The sensitivity of the rare decays $\eta^{(\prime)}\to\pi^{0}\gamma\gamma$ and $\eta^{\prime}\to\eta\gamma\gamma$ to signatures of a leptophobic $B$ boson in the MeV-GeV mass range is analyzed in this work. By adding an explicit $B$-boson resonance exchange, $\eta\to B\gamma\to\pi^{0}\gamma\gamma$, to the Standard Model contributions from vector and scalar meson exchanges, and employing experimental data for the associated branching ratios, it allows us to improve the current constraints on the $B$-boson mass $m_{B}$ and coupling to Standard Model particles $\alpha_{B}$. From these constraints and the analysis of the available experimental $\gamma\gamma$ invariant mass distribution, we show that a $B$-boson signature in the resonant mass range $m_{\pi^{0}}\lesssim m_{B}\lesssim m_{\eta}$ is strongly suppressed and would be very difficult to experimentally identify, assuming that the leptophobic $B$ boson only decays to Standard Model particles. In contrast, the limits outside this mass window are less stringent and the corresponding $t$- and $u$-channel signatures may still be observable in the data, as it occurs with the nonresonant Standard Model $\rho$, $\omega$ and $\phi$ meson exchanges. In addition, we make use of experimental data from the $\eta^{\prime}\to\pi^{0}\gamma\gamma$ and $\eta^{\prime}\to\eta\gamma\gamma$ decays to explore larger $B$-boson masses. Our results are relevant for the $B$-boson search programs at existing and forthcoming light-meson facilities, such as KLOE(-II) and Jefferson Lab Eta Factory experiments., Comment: 15 pages, 12 figures, 1 table; v2 includes a discussion of the uncertainty of the Standard Model estimate and an improvement in the overall presentation

Charge transport plays a crucial role in manifold potential applications of two-dimensional materials, including field effect transistors, solar cells, and transparent conductors. At most operating temperatures, charge transport is hindered by scattering of carriers by lattice vibrations. Assessing the intrinsic phonon-limited carrier mobility is thus of paramount importance to identify promising candidates for next-generation devices. Here we provide a framework to efficiently compute the drift and Hall carrier mobility of two-dimensional materials through the Boltzmann transport equation by relying on a Fourier-Wannier interpolation. Building on a recent formulation of long-range contributions to dynamical matrices and phonon dispersions [Phys. Rev. X 11, 041027 (2021)], we extend the approach to electron-phonon coupling including the effect of dynamical dipoles and quadrupoles. We identify an unprecedented contribution associated with the Berry connection that is crucial to preserve the Wannier-gauge covariance of the theory. This contribution is not specific to 2D crystals, but also concerns the 3D case, as we demonstrate via an application to bulk SrO. We showcase our method on a wide selection of relevant monolayers ranging from SnS2 to MoS2, graphene, BN, InSe, and phosphorene. We also discover a non-trivial temperature evolution of the Hall hole mobility in InSe whereby the mobility increases with temperature above 150 K due to the mexican-hat electronic structure of the InSe valence bands. Overall, we find that dynamical quadrupoles are essential and can impact the carrier mobility in excess of 75%., Comment: 31 pages and 16 figures

Despite considerable efforts, accurate computations of electron-phonon and carrier transport properties of low-dimensional materials from first principles have remained elusive. By building on recent advances in the description of long-range electrostatics, we develop a general approach to the calculation of electron-phonon couplings in two-dimensional materials. We show that the nonanalytic behavior of the electron-phonon matrix elements depends on the Wannier gauge, but that a missing Berry connection restores invariance to quadrupolar order. We showcase these contributions in a MoS$_2$ monolayer, calculating intrinsic drift and Hall mobilities with precise Wannier interpolations. We also find that the contributions of dynamical quadrupoles to the scattering potential are essential, and that their neglect leads to errors of 23% and 76% in the room temperature electron and hole Hall mobilities, respectively., Comment: 6 pages and 2 figures

The study of two-dimensional (2D) van der Waals materials has been an active field of research in the development of new optoelectronics and photonic applications over the last decade. Organic-inorganic layered perovskites are currently some of the most promising 2D van der Waals materials, due to their exceptional optical brightness and enhanced excitonic effects. However, low crystal quality and spectral diffusion usually broaden the exciton linewidth, obscuring the fine structure of the exciton in conventional photoluminescence experiments. Here, we propose a mechanical approach for reducing the effect of spectral diffusion by means of hBN-capping on layered perovskites with different thicknesses, revealing the exciton fine structure. We used a stochastic model to link the reduction of the spectral linewidth with the population of active charge fluctuation centres present in the organic spacer taking part in the dynamical Stark shift. Active fluctuation centres are reduced by a factor of 3.7 to 7.1 when we include hBN-capping according to our direct spectral measurements. This rate is in good agreement with the analysis of the overlap between the squared perovskite lattice and the hexagonal hBN lattice. Van der Waals forces between both lattices cause the partial clamping of the perovskite organic spacer molecules, and hence, the amplitude of the dynamical Stark shift characteristic of the spectral diffusion effect is reduced. Our work provides an easy and low-cost solution to the problem of accessing important fine-structure excitonic state information, along with an explanation of the important carrier dynamics present in the organic spacer that affect the quality of the optical emission.

Two-dimensional (2D) van der Waals nanomaterials have attracted considerable attention for potential use in photonic and optoelectronic applications in the nanoscale, due to their outstanding electrical and optical properties, differing from their bulk state. Currently, 2D perovskite belonging to this group of nanomaterials is widely studied for a wide range of optoelectronic applications. Thanks to their excitonic properties, 2D perovskites are also promising materials for photonics and nonlinear devices working at room temperature. Nevertheless, strong excitonic effects can reduce the photocurrent characteristics when using thinner perovskites phases. In this work, we present solid experimental evidence for the presence of single donor-acceptor pair optical transitions in 2D Lead Halide Perovskites, characterized by sub meV linewidths ($\simeq 120 \mu eV$) and long decay times (5-8 ns). Micro-photoluminescence evidence is supported by detailed Photoemission measurements, and a model simulation. The association of Phenethylammonium with Methylammonium cations, the latter molecule being only present in 2D Halide Perovskites with thicker phases $n \geq 2$, has been identified as the source of the donor-acceptor pair formation, corresponding to the displacement of lead atoms and their replacement by methylamonium. Our seminal study of discrete Donor-Acceptor Pair (DAP) sharp and bright optical transitions in 2D Lead Halide Perovskites opens new routes to implement DAP as the carrier sources for novel designs of optoelectronic devices with 2D perovskites, and will foster the development of future outstanding properties in non-linear quantum technologies.

In this work we investigate the prospect of observing new-physics signatures via $C\!P$ violation in $\eta^{(\prime)}\to\pi^0\mu^+\mu^-$ and $\eta^\prime\to\eta\mu^+\mu^-$ decays at the REDTOP experiment. We make use of the SMEFT to parametrise the new-physics $C\!P$-violating effects and find that the projected REDTOP statistics are not competitive with respect to nEDM experiments. This reasserts the $\eta\to\mu^+\mu^-$ process as the most promising channel to find $C\!P$-violation at this experimental facility., Comment: 12 pages, 1 figure, 1 table; v2 includes a comment about the neutron EDM constraint and possible accidental Wilson coefficient cancellations

Reach-avoid optimal control problems, in which the system must reach certain goal conditions while staying clear of unacceptable failure modes, are central to safety and liveness assurance for autonomous robotic systems, but their exact solutions are intractable for complex dynamics and environments. Recent successes in reinforcement learning methods to approximately solve optimal control problems with performance objectives make their application to certification problems attractive; however, the Lagrange-type objective used in reinforcement learning is not suitable to encode temporal logic requirements. Recent work has shown promise in extending the reinforcement learning machinery to safety-type problems, whose objective is not a sum, but a minimum (or maximum) over time. In this work, we generalize the reinforcement learning formulation to handle all optimal control problems in the reach-avoid category. We derive a time-discounted reach-avoid Bellman backup with contraction mapping properties and prove that the resulting reach-avoid Q-learning algorithm converges under analogous conditions to the traditional Lagrange-type problem, yielding an arbitrarily tight conservative approximation to the reach-avoid set. We further demonstrate the use of this formulation with deep reinforcement learning methods, retaining zero-violation guarantees by treating the approximate solutions as untrusted oracles in a model-predictive supervisory control framework. We evaluate our proposed framework on a range of nonlinear systems, validating the results against analytic and numerical solutions, and through Monte Carlo simulation in previously intractable problems. Our results open the door to a range of learning-based methods for safe-and-live autonomous behavior, with applications across robotics and automation. See https://github.com/SafeRoboticsLab/safety_rl for code and supplementary material., Comment: Accepted in Robotics: Science and Systems (RSS), 2021

We present the derivation and code implementation of a first-principles methodology to calculate the lattice-mediated contributions to the bulk flexoelectric tensor. The approach is based on our recent analytical long-wavelength extension of density-functional perturbation theory [Royo and Stengel, Phys. Rev. X 9, 021050 (2019)], and avoids the cumbersome numerical derivatives with respect to the wave vector that were adopted in previous implementations. To substantiate our results, we revisit and numerically validate the sum rules that relate flexoelectricity and uniform elasticity by generalizing them to regimes where finite forces and stresses are present. We also revisit the definition of the elastic tensor under stress, especially in regards to the existing linear-response implementation. We demonstrate the performance of our method by applying it to representative cubic crystals and to the tetragonal low-temperature polymorph of SrTiO$_3$, obtaining excellent agreement with the available literature data.

Two-dimensional (2D) insulators are a key element in the design and fabrication of van der Waals heterostructures. They are vital as transparent dielectric spacers whose thickness can influence both the photonic, electronic, and optoelectronic properties of 2D devices. Simultaneously, they provide protection of the active layers in the heterostructure. For these critical roles, hexagonal Boron Nitride (hBN) is the dominant choice due to its large bandgap, atomic flatness, low defect density, and encapsulation properties. However, the broad catalogue of 2D insulators offers exciting opportunities to replace hBN in certain applications that require transparent thin layers with additional optical degrees of freedom. Here we investigate the potential of single-crystalline Molybdenum Oxide (MoO$_3$) as an alternative 2D insulator for the design of nanodevices that require precise adjustment of the light polarization at the nanometer scale. First, we measure the wavelength-dependent refractive indices of MoO$_3$ along its three main crystal axes and determine the in-plane and out-of-plane anisotropy of its optical properties. We find the birefringence in MoO$_3$ nanosheets compares favorably with other 2D materials that exhibit strong birefringence, such as black phosphorus, ReS$_2$, or ReSe$_2$, in particular in the visible spectral range where MoO$_3$ has the unique advantage of transparency. Finally, we demonstrate the suitability of MoO$_3$ for dielectric encapsulation by reporting linewidth narrowing and reduced inhomogeneous broadening of 2D excitons and optically active quantum emitters, respectively, in a prototypical monolayer transition-metal dichalcogenide semiconductor. These results show the potential of MoO$_3$ as a 2D dielectric layer for manipulation of the light polarization in vertical 2D heterostructures., Comment: 7 pages, 5 figures

The paper proposes a method to convert a deep learning object detector into an equivalent spiking neural network. The aim is to provide a conversion framework that is not constrained to shallow network structures and classification problems as in state-of-the-art conversion libraries. The results show that models of higher complexity, such as the RetinaNet object detector, can be converted with limited loss in performance., Comment: 5 pages, submitted to ISCMI 2021 conference

Robots and autonomous systems must interact with one another and their environment to provide high-quality services to their users. Dynamic game theory provides an expressive theoretical framework for modeling scenarios involving multiple agents with differing objectives interacting over time. A core challenge when formulating a dynamic game is designing objectives for each agent that capture desired behavior. In this paper, we propose a method for inferring parametric objective models of multiple agents based on observed interactions. Our inverse game solver jointly optimizes player objectives and continuous-state estimates by coupling them through Nash equilibrium constraints. Hence, our method is able to directly maximize the observation likelihood rather than other non-probabilistic surrogate criteria. Our method does not require full observations of game states or player strategies to identify player objectives. Instead, it robustly recovers this information from noisy, partial state observations. As a byproduct of estimating player objectives, our method computes a Nash equilibrium trajectory corresponding to those objectives. Thus, it is suitable for downstream trajectory forecasting tasks. We demonstrate our method in several simulated traffic scenarios. Results show that it reliably estimates player objectives from a short sequence of noise-corrupted partial state observations. Furthermore, using the estimated objectives, our method makes accurate predictions of each player's trajectory., Comment: Submitted to RSS2021

Detecting pedestrians is a crucial task in autonomous driving systems to ensure the safety of drivers and pedestrians. The technologies involved in these algorithms must be precise and reliable, regardless of environment conditions. Relying solely on RGB cameras may not be enough to recognize road environments in situations where cameras cannot capture scenes properly. Some approaches aim to compensate for these limitations by combining RGB cameras with TOF sensors, such as LIDARs. However, there are few works that address this problem using exclusively the 3D geometric information provided by LIDARs. In this paper, we propose a PointNet++ based architecture to detect pedestrians in dense 3D point clouds. The aim is to explore the potential contribution of geometric information alone in pedestrian detection systems. We also present a semi-automatic labeling system that transfers pedestrian and non-pedestrian labels from RGB images onto the 3D domain. The fact that our datasets have RGB registered with point clouds enables label transferring by back projection from 2D bounding boxes to point clouds, with only a light manual supervision to validate results. We train PointNet++ with the geometry of the resulting 3D labelled clusters. The evaluation confirms the effectiveness of the proposed method, yielding precision and recall values around 98%.

The Hard Lefschetz Property (HLP) is an important property which has been studied in several categories of the symplectic world. For Sasakian manifolds, this duality is satisfied by the basic cohomology (so, it is a transverse property), but a new version of the HLP has been recently given in terms of duality of the cohomology of the manifold itself in arXiv:1306.2896. Both properties were proved to be equivalent (see arXiv:1311.1431) in the case of K-contact flows. In this paper we extend both versions of the HLP (transverse and not) to the more general category of isometric flows, and show that they are equivalent. We also give some explicit examples which illustrate the categories where the HLP could be considered., Comment: Accepted for publication in Transformation Groups. Minor changes in the exposition of the proof of Lemma 2.6

The manipulation of light emitted by two-dimensional semiconductors grounds forthcoming technologies in the field of on-chip communications. However, these technologies require from the so elusive out-of-plane photon sources to achieve an efficient coupling of radiated light into planar devices. Here we propose a versatile spectroscopic method that enables the identification of the out-of-plane component of dipoles. The method is based on the selective coupling of light emitted by in-plane and out-of-plane dipoles to the whispering gallery modes of spherical dielectric microresonators, in close contact to them. We have applied this method to demonstrate the existence of dipoles with an out-of-plane orientation in monolayer WSe2 at room temperature. Micro-photoluminescent measurements, numerical simulations based on finite element methods, and ab-initio calculations have identified trions as the source responsible for this out-of-plane emission, opening new routes for realizing on-chip integrated systems with applications in information processing and quantum communications., Comment: 21 pages, 4 figures

In theories with a warped extra dimension, composite fermions, as e.g. the right-handed top quark, can be very strongly coupled to Kaluza-Klein (KK) fields. In particular, the KK gluons in the presence of such composite fields become very broad resonances, thus remarkably modifying their experimental signatures. We have computed the pole mass and the pole width of the KK gluon, triggered by its interaction with quarks, as well as the prediction for proton-proton cross-sections using the full propagator and compared it with that obtained from the usual Breit-Wigner approximation. We compare both approaches, along with the existing experimental data from ATLAS and CMS, for the $t\bar t$, $t\bar t W$, $t\bar t Z$, $t\bar t H$, and $t\bar t t\bar t$ channels. We have found differences between the two approaches of up to about 100%, highlighting that the effect of broad resonances can be dramatic on present, and mainly future, experimental searches. The channel $t\bar t t\bar t$ is particularly promising because the size of the cross-section signal is of the same order of magnitude as the Standard Model prediction, and future experimental analyses in this channel, especially for broad resonances, can shed light on the nature of possible physics beyond the Standard Model., Comment: 32 pages, 14 figures; v2 includes additional references and correct units for cross-section on vertical axis of Fig. 14

We develop a fundamental theory of the long-range electrostatic interactions in two-dimensional crystals by performing a rigorous study of the nonanalyticities of the Coulomb kernel. We find that the dielectric functions are best represented by $2\times 2$ matrices, with nonuniform macroscopic potentials that are two-component hyperbolic functions of the out-of-plane coordinate, $z$. We demonstrate our arguments by deriving the long-range interatomic forces in the adiabatic regime, where we identify a formerly overlooked dipolar coupling involving the out-of-plane components of the dynamical charges. The resulting formula is exact up to an arbitrary multipolar order, which we illustrate in practice via the explicit inclusion of dynamical quadrupoles. By performing numerical tests on monolayer BN, SnS$_2$ and BaTiO$_3$ membranes, we show that our method allows for a drastic improvement in the description of the long-range electrostatic interactions, with comparable benefits to the quality of the interpolated phonon band structure.

In this work we present a multi-armed bandit framework for online expert selection in Markov decision processes and demonstrate its use in high-dimensional settings. Our method takes a set of candidate expert policies and switches between them to rapidly identify the best performing expert using a variant of the classical upper confidence bound algorithm, thus ensuring low regret in the overall performance of the system. This is useful in applications where several expert policies may be available, and one needs to be selected at run-time for the underlying environment., Comment: In proceedings of the 59th IEEE Conference on Decision and Control 2020. arXiv admin note: text overlap with arXiv:1707.05714

Building on recent developments in electronic-structure methods, we define and calculate the flexoelectric response of two-dimensional (2D) materials fully from first principles. In particular, we show that the open-circuit voltage response to a flexural deformation is a fundamental linear-response property of the crystal that can be calculated within the primitive unit cell of the flat configuration. Applications to graphene, silicene, phosphorene, BN and transition-metal dichalcogenide monolayers reveal that two distinct contributions exist, respectively of purely electronic and lattice-mediated nature. Within the former, we identify a key $metric$ term, consisting in the quadrupolar moment of the unperturbed charge density. We propose a simple continuum model to connect our findings with the available experimental measurements of the converse flexoelectric effect., Comment: To appear in Phys. Rev. Lett

A complete theoretical analysis of the $C$-conserving semileptonic decays $\eta^{(\prime)}\to\pi^0l^+l^-$ and $\eta^\prime\to\eta l^+l^-$ ($l=e$ or $\mu$) is carried out within the framework of the Vector Meson Dominance (VMD) model. An existing phenomenological model is used to parametrise the VMD coupling constants and the associated numerical values are obtained from an optimisation fit to $V\to P\gamma$ and $P\to V\gamma$ radiative decays ($V=\rho^0$, $\omega$, $\phi$ and $P=\pi^0$, $\eta$, $\eta^{\prime}$). The decay widths and dilepton energy spectra for the two $\eta\to\pi^0l^+l^-$ processes obtained using this approach are compared and found to be in good agreement with other results available in the published literature. Theoretical predictions for the four $\eta^{\prime}\to\pi^0l^+l^-$ and $\eta^\prime\to\eta l^+l^-$ decay widths and dilepton energy spectra are calculated and presented for the first time in this work., Comment: 11 pages, 2 figures, 1 table; v2 includes estimated errors coming from the dependence on the model and a detailed explanation of the structure of the $VP\gamma$ vertex; v3 corrects eqs. (6) and (15), includes discussions on the Feynman integral numerical computations and the use of the energy-dependent width for the $\rho$ propagator, as well as corrected numerical results in Table I and Fig. 2

The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged), Comment: The miniJPAS data and associated value added catalogues are publicly accessible via this url: http://archive.cefca.es/catalogues/minijpas-pdr201912

Information over the phonon band structure is crucial to predicting many thermodynamic properties of materials, such as thermal transport coefficients. Highly accurate phonon dispersion curves can be, in principle, calculated in the framework of density-functional perturbation theory (DFPT). However, well-established techniques can run into trouble (or even catastrophically fail) in the case of piezoelectric materials, where the acoustic branches hardly reproduce the physically correct sound velocity. Here we identify the culprit in the higher-order multipolar interactions between atoms, and demonstrate an effective procedure that fixes the aforementioned issue. Our strategy drastically improves the predictive power of perturbative lattice-dynamical calculations in piezoelectric crystals, and is directly implementable for high-throughput generation of materials databases.

An enhanced phenomenological model that includes isospin-symmetry breaking is presented in this letter. The model is then used in a number of statistical fits to the most recent experimental data for the radiative transitions $V\!P\gamma$ ($V=\rho$, $K^*$, $\omega$, $\phi$ and $P=\pi$, $K$, $\eta$, $\eta^{\prime}$) and estimations for the mixing angles amongst the three pseudoscalar states with vanishing third-component of isospin are obtained. The quality of the performed fits is good, e.g. $\chi^2_{\textrm{min}}/\textrm{d.o.f} = 1.9$. The current experimental uncertainties allow for isospin symmetry violations with a confidence level of approximately $2.5\sigma$., Comment: 10 pages, 3 tables

We include the treatment of quadrupolar fields beyond the Fr\"ohlich interaction in the first-principles electron-phonon vertex in semiconductors. Such quadrupolar fields induce long-range interactions that have to be taken into account for accurate physical results. We apply our formalism to Si (nonpolar), GaAs, and GaP (polar) and demonstrate that electron mobilities show large errors if dynamical quadrupoles are not properly treated.

We describe a new approach to compute the electron-phonon self-energy and carrier mobilities in semiconductors. Our implementation does not require a localized basis set to interpolate the electron-phonon matrix elements, with the advantage that computations can be easily automated. Scattering potentials are interpolated on dense $\mathbf{q}$ meshes using Fourier transforms and ab initio models to describe the long-range potentials generated by dipoles and quadrupoles. To reduce significantly the computational cost, we take advantage of crystal symmetries and employ the linear tetrahedron method and double-grid integration schemes, in conjunction with filtering techniques in the Brillouin zone. We report results for the electron mobility in Si, GaAs, and GaP obtained with this new methodology.

Iterative linear-quadratic (ILQ) methods are widely used in the nonlinear optimal control community. Recent work has applied similar methodology in the setting of multiplayer general-sum differential games. Here, ILQ methods are capable of finding local equilibria in interactive motion planning problems in real-time. As in most iterative procedures, however, this approach can be sensitive to initial conditions and hyperparameter choices, which can result in poor computational performance or even unsafe trajectories. In this paper, we focus our attention on a broad class of dynamical systems which are feedback linearizable, and exploit this structure to improve both algorithmic reliability and runtime. We showcase our new algorithm in three distinct traffic scenarios, and observe that in practice our method converges significantly more often and more quickly than was possible without exploiting the feedback linearizable structure., Comment: 7 pages, 5 figures, accepted to IEEE International Conference on Robotics and Automation (2020)

Distance metric learning based on triplet loss has been applied with success in a wide range of applications such as face recognition, image retrieval, speaker change detection and recently recommendation with the CML model. However, as we show in this article, CML requires large batches to work reasonably well because of a too simplistic uniform negative sampling strategy for selecting triplets. Due to memory limitations, this makes it difficult to scale in high-dimensional scenarios. To alleviate this problem, we propose here a 2-stage negative sampling strategy which finds triplets that are highly informative for learning. Our strategy allows CML to work effectively in terms of accuracy and popularity bias, even when the batch size is an order of magnitude smaller than what would be needed with the default uniform sampling. We demonstrate the suitability of the proposed strategy for recommendation and exhibit consistent positive results across various datasets., Comment: SIGIR 2019

In the recent years, singing voice separation systems showed increased performance due to the use of supervised training. The design of training datasets is known as a crucial factor in the performance of such systems. We investigate on how the characteristics of the training dataset impacts the separation performances of state-of-the-art singing voice separation algorithms. We show that the separation quality and diversity are two important and complementary assets of a good training dataset. We also provide insights on possible transforms to perform data augmentation for this task.

This paper describes the Polarimetric and Helioseismic Imager on the Solar Orbiter mission (SO/PHI), the first magnetograph and helioseismology instrument to observe the Sun from outside the Sun-Earth line. It is the key instrument meant to address the top-level science question: How does the solar dynamo work and drive connections between the Sun and the heliosphere? SO/PHI will also play an important role in answering the other top-level science questions of Solar Orbiter, as well as hosting the potential of a rich return in further science. SO/PHI measures the Zeeman effect and the Doppler shift in the FeI 617.3nm spectral line. To this end, the instrument carries out narrow-band imaging spectro-polarimetry using a tunable LiNbO_3 Fabry-Perot etalon, while the polarisation modulation is done with liquid crystal variable retarders (LCVRs). The line and the nearby continuum are sampled at six wavelength points and the data are recorded by a 2kx2k CMOS detector. To save valuable telemetry, the raw data are reduced on board, including being inverted under the assumption of a Milne-Eddington atmosphere, although simpler reduction methods are also available on board. SO/PHI is composed of two telescopes; one, the Full Disc Telescope (FDT), covers the full solar disc at all phases of the orbit, while the other, the High Resolution Telescope (HRT), can resolve structures as small as 200km on the Sun at closest perihelion. The high heat load generated through proximity to the Sun is greatly reduced by the multilayer-coated entrance windows to the two telescopes that allow less than 4% of the total sunlight to enter the instrument, most of it in a narrow wavelength band around the chosen spectral line., Comment: 36 pages, 41 figures

To use neural networks in safety-critical settings it is paramount to provide assurances on their runtime operation. Recent work on ReLU networks has sought to verify whether inputs belonging to a bounded box can ever yield some undesirable output. Input-splitting procedures, a particular type of verification mechanism, do so by recursively partitioning the input set into smaller sets. The efficiency of these methods is largely determined by the number of splits the box must undergo before the property can be verified. In this work, we propose a new technique based on shadow prices that fully exploits the information of the problem yielding a more efficient generation of splits than the state-of-the-art. Results on the Airborne Collision Avoidance System (ACAS) benchmark verification tasks show a considerable reduction in the partitions generated which substantially reduces computation times. These results open the door to improved verification methods for a wide variety of machine learning applications including vision and control.

Van der Waals materials offer a wide range of atomic layers with unique properties that can be easily combined to engineer novel electronic and photonic devices. A missing ingredient of the van der Waals platform is a two-dimensional crystal with naturally occurring out-of-plane luminescent dipole orientation. Here we measure the far-field photoluminescence intensity distribution of bulk InSe and two-dimensional InSe, WSe$_2$ and MoSe$_2$. We demonstrate, with the support of ab-initio calculations, that layered InSe flakes sustain luminescent excitons with an intrinsic out-of-plane orientation, in contrast with the in-plane orientation of dipoles we find in two-dimensional WSe$_2$ and MoSe$_2$ at room-temperature. These results, combined with the high tunability of the optical response and outstanding transport properties, position layered InSe as a promising semiconductor for novel optoelectronic devices, in particular for hybrid integrated photonic chips which exploit the out-of-plane dipole orientation., Comment: Published in Nature Communications. 12 pages, 5 figures. Link: https://www.nature.com/articles/s41467-019-11920-4

The scalar and vector meson exchange contributions to the doubly radiative decays $\eta^{(\prime)}\to\pi^0\gamma\gamma$ and $\eta^\prime\to\eta\gamma\gamma$ are analysed within the Linear Sigma Model and Vector Meson Dominance frameworks, respectively. Predictions for the diphoton invariant mass spectra and the associated integrated branching ratios are given and compared with current available experimental data. Whilst a satisfactory description of the shape of the $\eta\to\pi^{0}\gamma\gamma$ and $\eta^\prime\to\pi^{0}\gamma\gamma$ decay spectra is obtained, thus supporting the validity of the approach, the corresponding branching ratios cannot be reproduced simultaneously. A first theoretical prediction for the recently measured $\eta^\prime\to\eta\gamma\gamma$ by the BESIII collaboration is also presented., Comment: 16 pages, 4 figures, 2 tables; v2: includes experimental measurement for the $\eta^\prime\to\eta\gamma\gamma$, a better treatment for the $\sigma$ propagator and a new discussion on the values of the $VP\gamma$ couplings

Density-functional perturbation theory (DFPT) is nowadays the method of choice for the accurate computation of linear and non-linear response properties of materials from first principles. A notable advantage of DFPT over alternative approaches is the possibility of treating incommensurate lattice distortions with an arbitrary wavevector, ${\bf q}$, at essentially the same computational cost as the lattice-periodic case. Here we show that ${\bf q}$ can be formally treated as a perturbation parameter, and used in conjunction with established results of perturbation theory (e.g. the "2n+1" theorem) to perform a long-wave expansion of an arbitrary response function in powers of the wavevector components. This provides a powerful, general framework to accessing a wide range of spatial dispersion effects that were formerly difficult to calculate by means of first-principles electronic-structure methods. In particular, the physical response to the spatial gradient of any external field can now be calculated at negligible cost, by using the response functions to $\mathit{uniform}$ perturbations (electric, magnetic or strain fields) as the only input. We demonstrate our method by calculating the flexoelectric and dynamical quadrupole tensors of selected crystalline insulators and model systems.

We address the problem of disambiguating large scale catalogs through the definition of an unknown artist clustering task. We explore the use of metric learning techniques to learn artist embeddings directly from audio, and using a dedicated homonym artists dataset, we compare our method with a recent approach that learn similar embeddings using artist classifiers. While both systems have the ability to disambiguate unknown artists relying exclusively on audio, we show that our system is more suitable in the case when enough audio data is available for each artist in the train dataset. We also propose a new negative sampling method for metric learning that takes advantage of side information such as music genre during the learning phase and shows promising results for the artist clustering task., Comment: published in ISMIR 2018

We consider the task of multimodal music mood prediction based on the audio signal and the lyrics of a track. We reproduce the implementation of traditional feature engineering based approaches and propose a new model based on deep learning. We compare the performance of both approaches on a database containing 18,000 tracks with associated valence and arousal values and show that our approach outperforms classical models on the arousal detection task, and that both approaches perform equally on the valence prediction task. We also compare the a posteriori fusion with fusion of modalities optimized simultaneously with each unimodal model, and observe a significant improvement of valence prediction. We release part of our database for comparison purposes., Comment: Published in ISMIR 2018

In this paper, we propose to infer music genre embeddings from audio datasets carrying semantic information about genres. We show that such embeddings can be used for disambiguating genre tags (identification of different labels for the same genre, tag translation from a tag system to another, inference of hierarchical taxonomies on these genre tags). These embeddings are built by training a deep convolutional neural network genre classifier with large audio datasets annotated with a flat tag system. We show empirically that they makes it possible to retrieve the original taxonomy of a tag system, spot duplicates tags and translate tags from a tag system to another., Comment: published in ISMIR 2018