title: The Environment of Lyman Break Analogues (ELBA) survey: Star-forming galaxies in small groups

abstract: The Environment of Lyman Break Analogues (ELBA) survey is an imaging survey of 33 $deg^{2}$ of the southern sky. The survey was observed in {\it u}, {\it g}, {\it r}, and {\it i} bands with the Dark Energy Camera (DECam) on the Blanco telescope. The main goal of this project is to investigate the environment of Lyman break analogues (LBAs), low-redshift (z $\sim $0.2) galaxies that are remarkably similar to typical star-forming galaxies at z $\sim$ 3. We explore whether the environment has any influence on the observed properties of these galaxies, providing valuable insight on the formation and evolution of galaxies over cosmic time. Using the Nearest Neighbour method, we measure the local density of each object ranging from small to large scales (clusters of galaxies). Comparing the environment around LBAs with that of the general galaxy population in the field, we conclude that LBAs, on average, populate denser regions at small scales ($\sim$ $1.5Mpc$), but are located in similar environment to other star-forming galaxies at larger scales ($\sim$ $3.0 Mpc$). This offers evidence that nearby encounters such as mergers may influence the star formation activity in LBAs, before infall onto larger galaxy clusters. We interpret this an indication of galaxy preprocessing, in agreement with theoretical expectations for galaxies at z $\sim$ 2 -3 where the gravitational interactions are more intense in early formation processes of this objects


title: Arguments against using $h^{-1}{\rm Mpc}$ units in observational cosmology

abstract: It is common to express cosmological measurements in units of $h^{-1}{\rm Mpc}$. Here, we review some of the complications that originate from this practice. A crucial problem caused by these units is related to the normalization of the matter power spectrum, which is commonly characterized in terms of the linear-theory rms mass fluctuation in spheres of radius $8\,h^{-1}{\rm Mpc}$, $\sigma_8$. This parameter does not correctly capture the impact of $h$ on the amplitude of density fluctuations. We show that the use of $\sigma_8$ has caused critical misconceptions for both the so-called $\sigma_8$ tension regarding the consistency between low-redshift probes and cosmic microwave background data, and the way in which growth-rate estimates inferred from redshift-space distortions are commonly expressed. We propose to abandon the use of $h^{-1}{\rm Mpc}$ units in cosmology and to characterize the amplitude of the matter power spectrum in terms of $\sigma_{12}$, defined as the mass fluctuation in spheres of radius $12\,{\rm Mpc}$, whose value is similar to the standard $\sigma_8$ for $h\sim 0.67$.


title: Linearization of transition functions of a semi-positive line bundle along a certain submanifold

abstract: Let $X$ be a complex manifold and $L$ be a holomorphic line bundle on $X$. Assume that $L$ is semi-positive, namely $L$ admits a smooth Hermitian metric with semi-positive Chern curvature. Let $Y$ be a compact K\"ahler submanifold of $X$ such that the restriction of $L$ to $Y$ is topologically trivial. We investigate the obstruction for $L$ to be unitary flat on a neighborhood of $Y$ in $X$. As an application, for example, we show the existence of nef, big, and non semi-positive line bundle on a non-singular projective surface.


title: Newly Improved Ionization Corrections for the Neutral Interstellar Medium: Enabling Accurate Abundance Determinations in Star-forming Galaxies throughout the Universe

abstract: Studies measuring the chemical abundances of the neutral gas in star-forming galaxies (SFGs) require ionization correction factors (ICFs) to accurately measure their metal contents. In the work presented here we calculate newly improved ICFs for a sample of SFGs. These new corrections include both the contaminating ionized gas along the line of sight (ICF$_{\rm ionized}$) and unaccounted higher ionization stages in the neutral gas (ICF$_{\rm neutral}$). We make use of recently acquired spectroscopic observations taken with the Cosmic Origins Spectrograph (COS) on board Hubble to measure column densities for Fe II and Fe III. Using the Fe III/Fe II ratios as well as other physical properties (i.e. $\log$[L$_{\rm UV}$], $N$(H I), T, and $Z$) we generate ad-hoc photoionization models with CLOUDY to quantify the corrections required for each of the targets. We identify a luminosity threshold of $\log$[L$_{\rm UV}$]$\sim$ 40.75 erg s$^{-1}$ above which the ICF$_{\rm neutral}$ values for nitrogen are relatively higher (ICF$_{\rm neutral}=0.05$-0.7) than those for the rest of the elements (ICF$_{\rm neutral}\sim 0.01$). This behavior indicates that for the high UV luminosity objects, N II is found in non-negligible quantities in the neutral gas, making these ICF$_{\rm neutral}$ corrections critical for determining the true abundances in the interstellar medium. In addition, we calculate ICFs from a uniform grid of models covering a wide range of physical properties typically observed in studies of SFGs and extragalactic H II regions. We provide the community with tabulated ICF values for the neutral gas abundances measured from a variety of environments and applicable to chemical studies of the high redshift universe.


title: Closed-Form Optimal Impulsive Control of Spacecraft Relative Motion Using Reachable Set Theory

abstract: This paper addresses the spacecraft relative orbit reconfiguration problem of minimizing the delta-v cost of impulsive control actions while achieving a desired state in fixed time. The problem is posed in relative orbit element (ROE) space, which yields insight into relative motion geometry and allows for the straightforward inclusion of perturbations in linear time-variant form. Reachable set theory is used to translate the cost-minimization problem into a geometric path-planning problem and formulate the reachable delta-v minimum, a new metric to assess optimality and quantify reachability of a maneuver scheme. Next, this paper presents a methodology to compute maneuver schemes that meet this new optimality criteria and achieve a prescribed reconfiguration. Though the methodology is applicable to any linear time-variant system, this paper leverages a state representation in ROE to derive new globally optimal maneuver schemes in orbits of arbitrary eccentricity. The methodology is also used to generate quantifiably sub-optimal solutions when the optimal solutions are unreachable. Further, this paper determines the mathematical impact of uncertainties on achieving the desired end state and provides a geometric visualization of those effects on the reachable set. The proposed algorithms are tested in realistic reconfiguration scenarios and validated in a high-fidelity simulation environment.


title: Higher-Order Label Homogeneity and Spreading in Graphs

abstract: Do higher-order network structures aid graph semi-supervised learning? Given a graph and a few labeled vertices, labeling the remaining vertices is a high-impact problem with applications in several tasks, such as recommender systems, fraud detection and protein identification. However, traditional methods rely on edges for spreading labels, which is limited as all edges are not equal. Vertices with stronger connections participate in higher-order structures in graphs, which calls for methods that can leverage these structures in the semi-supervised learning tasks. To this end, we propose Higher-Order Label Spreading (HOLS) to spread labels using higher-order structures. HOLS has strong theoretical guarantees and reduces to standard label spreading in the base case. Via extensive experiments, we show that higher-order label spreading using triangles in addition to edges is up to 4.7% better than label spreading using edges alone. Compared to prior traditional and state-of-the-art methods, the proposed method leads to statistically significant accuracy gains in all-but-one cases, while remaining fast and scalable to large graphs.


title: Quantum dynamics in strongly driven random dipolar magnets

abstract: The random dipolar magnet LiHo$_x$Y$_{1-x}$F$_4$ enters a strongly frustrated regime for small Ho$^{3+}$ concentrations with $x<0.05$. In this regime, the magnetic moments of the Ho$^{3+}$ ions experience small quantum corrections to the common Ising approximation of LiHo$_x$Y$_{1-x}$F$_4$, which lead to a $Z_2$-symmetry breaking and small, degeneracy breaking energy shifts between different eigenstates. Here we show that destructive interference between two almost degenerate excitation pathways burns spectral holes in the magnetic susceptibility of strongly driven magnetic moments in LiHo$_x$Y$_{1-x}$F$_4$. Such spectral holes in the susceptibility, microscopically described in terms of Fano resonances, can already occur in setups of only two or three frustrated moments, for which the driven level scheme has the paradigmatic $\Lambda$-shape. For larger clusters of magnetic moments, the corresponding level schemes separate into almost isolated many-body $\Lambda$-schemes, in the sense that either the transition matrix elements between them are negligibly small or the energy difference of the transitions is strongly off-resonant to the drive. This enables the observation of Fano resonances, caused by many-body quantum corrections to the common Ising approximation also in the thermodynamic limit. We discuss its dependence on the driving strength and frequency as well as the crucial role that is played by lattice dissipation.


title: Sublimation Kinetics for Individual Graphite and Graphene Nano-particles (NPs): NP-to-NP Variations and Evolving Structure-Kinetics and Structure-Emissivity Relationships

abstract: A single nanoparticle (NP) mass spectrometry method was used to measure sublimation rates as a function of nanoparticle temperature (TNP) for a number of individual graphite and graphene NPs. Initially, the NP sublimation rates were ca. 400 times faster than that for bulk graphite, and there were large NP-to-NP variations. Over time, the rate slowed substantially, though remaining well above the bulk rate. The initial activation energies (Eas) were correspondingly low and doubled as a few monolayer's worth of material were sublimed from the surfaces. The high initial rates and low Eas are attributed to large numbers of edge and other low coordination sites on the NP surfaces, and the changes are attributed to atomic-scale "smoothing" of the surface by preferential sublimation of the less stable sites. The emissivity of the NPs also changed after heating, most frequently increasing. The emissivity and sublimation rates were anti-correlated, leading to the conclusion that high densities of low-coordination sites on the NP surfaces enhances sublimation but suppresses emissivity.


title: Theoretical Convergence of Multi-Step Model-Agnostic Meta-Learning

abstract: As a popular meta-learning approach, the model-agnostic meta-learning (MAML) algorithm has been widely used due to its simplicity and effectiveness. However, the convergence of the general multi-step MAML still remains unexplored. In this paper, we develop a new theoretical framework to provide such convergence guarantee for two types of objective functions that are of interest in practice: (a) resampling case (e.g., reinforcement learning), where loss functions take the form in expectation and new data are sampled as the algorithm runs; and (b) finite-sum case (e.g., supervised learning), where loss functions take the finite-sum form with given samples. For both cases, we characterize the convergence rate and the computational complexity to attain an $\epsilon$-accurate solution for multi-step MAML in the general nonconvex setting. In particular, our results suggest that an inner-stage stepsize needs to be chosen inversely proportional to the number $N$ of inner-stage steps in order for $N$-step MAML to have guaranteed convergence. From the technical perspective, we develop novel techniques to deal with the nested structure of the meta gradient for multi-step MAML, which can be of independent interest.


title: Incremental Nonlinear Fault-Tolerant Control of a Quadrotor with Complete Loss of Two Opposing Rotors

abstract: In order to further expand the flight envelope of quadrotors under actuator failures, we design a nonlinear sensor-based fault-tolerant controller to stabilize a quadrotor with failure of two opposing rotors in the high-speed flight condition (> 8m/s). The incremental nonlinear dynamic inversion (INDI) approach which excels in handling model uncertainties is adopted to compensate for the significant unknown aerodynamic effects. The internal dynamics of such an underactuated system have been analyzed, and subsequently stabilized by re-defining the control output. The proposed method can be generalized to control a quadrotor under single-rotor-failure and nominal conditions. For validation, flight tests have been carried out in a large-scale open jet wind tunnel. The position of a damaged quadrotor can be controlled in the presence of significant wind disturbances. A linear quadratic regulator (LQR) approach from the literature has been compared to demonstrate the advantages of the proposed nonlinear method in the windy and high-speed flight condition.


title: Cyberattack Action-Intent-Framework for Mapping Intrusion Observables

abstract: The techniques and tactics used by cyber adversaries are becoming more sophisticated, ironically, as defense getting stronger and the cost of a breach continuing to rise. Understanding the thought processes and behaviors of adversaries is extremely challenging as high profile or even amateur attackers have no incentive to share the trades associated with their illegal activities. One opportunity to observe the actions the adversaries perform is through the use of Intrusion Detection Systems (IDS) which generate alerts in the event that suspicious behavior was detected. The alerts raised by these systems typically describe the suspicious actions via the form of attack 'signature', which do not necessarily reveal the true intent of the attacker performing the action. Meanwhile, several high level frameworks exist to describe the sequence or chain of action types an adversary might perform. These frameworks, however, do not connect the action types to observables of standard intrusion detection systems, nor describing the plausible intents of the adversarial actions. To address these gaps, this work proposes the Action-Intent Framework (AIF) to complement existing Cyber Attack Kill Chains and Attack Taxonomies. The AIF defines a set of Action-Intent States (AIS) at two levels of description: the Macro-AIS describes 'what' the attacker is trying to achieve and the Micro-AIS describes "how" the intended goal is achieved. A full description of both the Macro is provided along with a set of guiding principals of how the AIS is derived and added to the framework.


title: Is Local SGD Better than Minibatch SGD?

abstract: We study local SGD (also known as parallel SGD and federated averaging), a natural and frequently used stochastic distributed optimization method. Its theoretical foundations are currently lacking and we highlight how all existing error guarantees in the convex setting are dominated by a simple baseline, minibatch SGD. (1) For quadratic objectives we prove that local SGD strictly dominates minibatch SGD and that accelerated local SGD is minimax optimal for quadratics; (2) For general convex objectives we provide the first guarantee that at least sometimes improves over minibatch SGD; (3) We show that indeed local SGD does not dominate minibatch SGD by presenting a lower bound on the performance of local SGD that is worse than the minibatch SGD guarantee.


title: Sparse Hop Spanners for Unit Disk Graphs

abstract: A unit disk graph $G$ on a given set $P$ of points in the plane is a geometric graph where an edge exists between two points $p,q \in P$ if and only if $|pq| \leq 1$. A spanning subgraph $G'$ of $G$ is a $k$-hop spanner if and only if for every edge $pq\in G$, there is a path between $p,q$ in $G'$ with at most $k$ edges. We obtain the following results for unit disk graphs in the plane. (I) Every $n$-vertex unit disk graph has a $5$-hop spanner with at most $5.5n$ edges. We analyze the family of spanners constructed by Biniaz (2020) and improve the upper bound on the number of edges from $9n$ to $5.5n$. (II) Using a new construction, we show that every $n$-vertex unit disk graph has a $3$-hop spanner with at most $11n$ edges. (III) Every $n$-vertex unit disk graph has a $2$-hop spanner with $O(n\log n)$ edges. This is the first nontrivial construction of $2$-hop spanners. (IV) For every sufficiently large positive integer $n$, there exists a set $P$ of $n$ points on a circle, such that every plane hop spanner on $P$ has hop stretch factor at least $4$. Previously, no lower bound greater than $2$ was known. (V) For every finite point set on a circle, there exists a plane (i.e., crossing-free) $4$-hop spanner. As such, this provides a tight bound for points on a circle. (VI) The maximum degree of $k$-hop spanners cannot be bounded from above by a function of $k$ for any positive integer $k$.


title: Performance of a reconciliation method operating on a discreet quantum key distribution system

abstract: Reconciliation is a mechanism allowing to weed out the discrepancies between two correlated variables. It has great role in every Quantum Key Distribution protocol where the key has to be transmitted through a noisy channel or as in our case of study in presence of an eavesdropping. In this paper, we show that for discrete-variable QKD protocols, this problem can be advantageously solved with Turbo codes. In particular, we demonstrate that our method leads to a significant improvement of Bit Error Rate, may divide it by three in presence of a eavesdropper even with great eavesdropping capability.


title: Analyzing Grant-Free Access for URLLC Service

abstract: 5G New Radio (NR) is expected to support new ultra-reliable low-latency communication (URLLC) service targeting at supporting the small packets transmissions with very stringent latency and reliability requirements. Current Long Term Evolution (LTE) system has been designed based on grantbased (GB) (i.e., dynamic grant) random access, which can hardly support the URLLC requirements. Grant-free (GF) (i.e., configured grant) access is proposed as a feasible and promising technology to meet such requirements, especially for uplink transmissions, which effectively saves the time of requesting/waiting for a grant. While some basic GF access features have been proposed and standardized in NR Release-15, there is still much space to improve. Being proposed as 3GPP study items, three GF access schemes with Hybrid Automatic Repeat reQuest (HARQ) retransmissions including Reactive, K-repetition, and Proactive, are analyzed in this paper. Specifically, we present a spatiotemporal analytical framework for the contention-based GF access analysis. Based on this framework, we define the latent access failure probability to characterize URLLC reliability and latency performances. We propose a tractable approach to derive and analyze the latent access failure probability of the typical UE under three GF HARQ schemes. Our results show that under shorter latency constraints, the Proactive scheme provides the lowest latent access failure probability, whereas, under longer latency constraints, the K-repetition scheme achieves the lowest latent access failure probability, which depends on K. If K is overestimated, the Proactive scheme provides lower latent access failure probability than the K-repetition scheme.


title: Polarisation-insensitive generation of vector modes using a digital micromirror device

abstract: In recent time there has been an increasing amount of interest in developing novel techniques for the generation of complex vector light beams. Amongst these, digital holography stands out as one of the most flexible and versatile with almost unlimited freedom to generate scalar and vector light beams with arbitrary polarisation distributions and spatial transverse profile. Recently, we put forward a novel method to quantify the non-separability of vector modes in which we reported first measurements of a compact and robust device to generate such vector modes that fully exploits the polarisation-independence of Digital Micromirror Devices (DMDs). In this manuscript we fully characterise this device and provide qualitative and quantitative analysis of the generated modes. First by reconstructing their transverse polarisation distribution, using stokes polarimetry, followed by a measure of their degree of non-separability, determined through the concurrence.


title: Capacity-achieving Spatially Coupled Sparse Superposition Codes with AMP Decoding

abstract: Sparse superposition codes, also called sparse regression codes (SPARCs), are a class of codes for efficient communication over the AWGN channel at rates approaching the channel capacity. In a standard SPARC, codewords are sparse linear combinations of columns of an i.i.d. Gaussian design matrix, while in a spatially coupled SPARC the design matrix has a block-wise structure, where the variance of the Gaussian entries can be varied across blocks. A well-designed spatial coupling structure can significantly enhance the error performance of iterative decoding algorithms such as Approximate Message Passing (AMP). In this paper, we obtain a non-asymptotic bound on the probability of error of spatially coupled SPARCs with AMP decoding. Applying this bound to a simple band-diagonal design matrix, we prove that spatially coupled SPARCs with AMP decoding achieve the capacity of the AWGN channel. The bound also highlights how the decay of error probability depends on each design parameter of the spatially coupled SPARC. An attractive feature of AMP decoding is that its asymptotic mean squared error (MSE) can be predicted via a deterministic recursion called state evolution. Our result provides the first proof that the MSE concentrates on the state evolution prediction for spatially coupled designs. Combined with the state evolution prediction, this result implies that spatially coupled SPARCs with the proposed band-diagonal design are capacity-achieving. Using the proof technique used to establish the main result, we also obtain a concentration inequality for the MSE of AMP applied to compressed sensing with spatially coupled design matrices. Finally we provide numerical simulation results that demonstrate the finite length error performance of spatially coupled SPARCs. The performance is compared with coded modulation schemes that use LDPC codes from the DVB-S2 standard.


title: Interpretable Multi-Headed Attention for Abstractive Summarization at Controllable Lengths

abstract: Abstractive summarization at controllable lengths is a challenging task in natural language processing. It is even more challenging for domains where limited training data is available or scenarios in which the length of the summary is not known beforehand. At the same time, when it comes to trusting machine-generated summaries, explaining how a summary was constructed in human-understandable terms may be critical. We propose Multi-level Summarizer (MLS), a supervised method to construct abstractive summaries of a text document at controllable lengths. The key enabler of our method is an interpretable multi-headed attention mechanism that computes attention distribution over an input document using an array of timestep independent semantic kernels. Each kernel optimizes a human-interpretable syntactic or semantic property. Exhaustive experiments on two low-resource datasets in the English language show that MLS outperforms strong baselines by up to 14.70% in the METEOR score. Human evaluation of the summaries also suggests that they capture the key concepts of the document at various length-budgets.


title: Andrews-Bressoud series and Wronskians

abstract: The Andrews-Bressoud identities are one of many families of $q$-series identities relating an infinite sum to an infinite product. While the original motivation for studying these series relates to partitions, they can also be viewed in relation to irreducible characters of minimal models in the theory of vertex operator algebras. Furthermore, considering certain Wronskians of the Andrews-Bressoud series for a given modulus produces additional $q$-series, which themselves exhibit interesting and predictable modularity properties. In this paper, we connect the Andrews-Bressoud series to modular forms and prove results about the modularity of their associated Wronskians.


title: The Solar wind prevents re-accretion of debris after Mercury's giant impact

abstract: The planet Mercury possesses an anomalously large iron core, and a correspondingly high bulk density. Numerous hypotheses have been proposed in order to explain such a large iron content. A long-standing idea holds that Mercury once possessed a larger silicate mantle which was removed by a giant impact early in the the Solar system's history. A central problem with this idea has been that material ejected from Mercury is typically re-accreted onto the planet after a short (~Myr) timescale. Here, we show that the primordial Solar wind would have provided sufficient drag upon ejected debris to remove them from Mercury-crossing trajectories before re-impacting the planet's surface. Specifically, the young Sun likely possessed a stronger wind, fast rotation and strong magnetic field. Depending upon the time of the giant impact, the ram pressure associated with this wind would push particles outward into the Solar system, or inward toward the Sun, on sub-Myr timescales, depending upon the size of ejected debris. Accordingly, the giant impact hypothesis remains a viable pathway toward the removal of planetary mantles, both on Mercury and extrasolar planets, particularly those close to young stars with strong winds.


title: Kinetic roughening of the urban skyline

abstract: \begin{abstract} In this Letter we follow the asymptotic spatial correlation of buildings' heights $G_{\infty}(r)$ in the whole of the Netherlands \cite{bag3d}, which comprises $\approx 10,000,000$ buildings, for the purpose of recovering its scaling with respect to space and time given respectively by the exponents $r^{2 \alpha}$ and $t^{2\beta}$. This allows us to identify the universality class of the evolution of the urban skyline seen as a dynamically evolving interface. Two major classes of cities were identified based on the recovered value of $\alpha =0.4$ and $\alpha = 0$, which correspond respectively to the KPZ and the EW universality classes. Picking a discrete model from each of these classes and mapping it to physical rules for constructions in cities we conclude that imposed restrictions on buildings' heights are reflected in the exponent $\alpha$ and thus have implications on how the skyline evolves.


title: Design of Millimeter-Wave Single-Shot Beam Training for True-Time-Delay Array

abstract: Beam training is one of the most important and challenging tasks in millimeter-wave and sub-terahertz communications. Novel transceiver architectures and signal processing techniques are required to avoid prohibitive training overhead when large antenna arrays with narrow beams are used. In this work, we leverage recent developments in wide range true-time-delay (TTD) analog arrays and frequency dependent probing beams to accelerate beam training. We propose an algorithm that achieves high-accuracy angle of arrival estimation with a single training symbol. Further, the impact of TTD front-end impairments on beam training accuracy is investigated, including the impact of gain, phase, and delay errors. Lastly, the study on impairments and required specifications of resolution and range of analog delay taps are used to provide a design insight of energy efficient TTD array, which employs a novel architecture with discrete-time sampling based TTD elements.


title: No Massive Companion to the Coherent Radio-Emitting M Dwarf GJ 1151

abstract: The recent detection of circularly polarized, long-duration (>8 hr) low-frequency (~150 MHz) radio emission from the M4.5 dwarf GJ 1151 has been interpreted as arising from a star-planet interaction via the electron cyclotron maser instability. The existence or parameters of the proposed planets have not been determined. Using 20 new HARPS-N observations, we put 99th-percentile upper limits on the mass of any close companion to GJ 1151 at Msini < 5.6 M earth. With no stellar, brown dwarf, or giant planet companion likely in a close orbit, our data are consistent with detected radio emission emerging from a magnetic interaction between a short-period terrestrial-mass planet and GJ 1151.


title: 1324- and 2143-avoiding Kazhdan-Lusztig immanants and k-positivity

abstract: Immanants are functions on square matrices generalizing the determinant and permanent. Kazhdan-Lusztig immanants, which are indexed by permutations, involve $q=1$ specializations of Type A Kazhdan-Lusztig polynomials, and were defined in (Rhoades-Skandera, 2006). Using results of (Haiman, 1993) and (Stembridge, 1991), Rhoades and Skandera showed that Kazhdan-Lusztig immanants are nonnegative on matrices whose minors are nonnegative. We investigate which Kazhdan-Lusztig immanants are positive on $k$-positive matrices (matrices whose minors of size $k \times k$ and smaller are positive). We show that the Kazhdan-Lusztig immanant indexed by $v$ is positive on $k$-positive matrices when $v$ avoids 1324 and 2143 and for all non-inversions $i<j$ of $v$, either $j-i \leq k$ or $v_j-v_i\leq k$. Our main tool is Lewis Carroll's identity.


title: A Mystery in Chamaeleon: Serendipitous Discovery of a Galactic Symbiotic Nova

abstract: We present the serendipitous discovery of a low luminosity nova occurring in a symbiotic binary star system in the Milky Way. We lay out the extensive archival data alongside new follow-up observations related to the stellar object V$^*$ CN Cha in the constellation of Chamaeleon. The object had long period ($\sim\! 250\,$day), high amplitude ($\sim\! 3\,$mag) optical variability in its recent past, preceding an increase in optical brightness by $\sim\! 8\,$magnitudes and a persistence at this luminosity for about 3 years, followed by a period of $\sim\! 1.4\,{\rm mag}\,{\rm yr}^{-1}$ dimming. The object's current optical luminosity seems to be dominated by H$\alpha$ emission, which also exhibits blue-shifted absorption (a P-Cygni-like profile). After consideration of a number of theories to explain these myriad observations, we determine that V$^*$ CN Cha is most likely a symbiotic (an evolved star-white dwarf binary) system which has undergone a long-duration, low luminosity, nova. Interpreted in this way, the outburst in V$^*$ CN Cha is among the lowest luminosity novae ever observed.


title: The Capacity and Optimal Signaling for Gaussian MIMO Channels Under Interference Constraints (full version)

abstract: Gaussian MIMO channel under total transmit and interference power constraints (TPC and IPC) is considered. A closed-form solution for the optimal transmit covariance matrix in the general case is obtained using the KKT-based approach (up to dual variables). While closed-from solutions for optimal dual variables are possible in special cases, an iterative bisection algorithm (IBA) is proposed to find the optimal dual variables in the general case and its convergence is proved for some special cases. Numerical experiments illustrate its efficient performance. Bounds for the optimal dual variables are given, which facilitate numerical solutions. An interplay between the TPC and IPC is studied, including the transition from power-limited to interference-limited regimes as the total transmit power increases. Sufficient and necessary conditions for each constraint to be redundant are given. A number of explicit closed-form solutions are obtained, including full-rank and rank-1 (beamforming) cases as well as the case of identical eigenvectors (typical for massive MIMO settings). A bound on the rank of optimal covariance is established. A number of unusual properties of optimal covariance matrix are pointed out.


title: Detecting Pair-Instability Supernovae at z<5 with the James Webb Space Telescope

abstract: Pair-instability supernovae (PISNe) are the ultimate cosmic lighthouses, capable of being observed at z<25 and revealing the properties of primordial stars at cosmic dawn. But it is now understood that the spectra and light curves of these events evolved with redshift as the universe became polluted with heavy elements because chemically enriched stars in this mass range typically lose most of their hydrogen envelopes and explode as bare helium cores. The light curves of such transients can be considerably dimmer in the near infrared (NIR) today than those of primordial PISNe of equal energy and progenitor mass. Here, we calculate detection rates for PISNe whose progenitors lost their outer layers to either line-driven winds or rotation at z<10, their detection limit in redshift for the James Webb Space Telescope (JWST). We find that JWST may be able to detect only Pop II (metal-poor) PISNe over the redshift range of z<4, but not their Pop III (metal-free) counterparts.


title: Logarithmic Holder continuous mappings and Beltrami equation

abstract: The paper is devoted to the study of mappings satisfying the inverse Poletsky inequality. We study the local behavior of these mappings, moreover, we are most interested in the case when the corresponding majorant is integrable on some set of spheres of positive linear measure. The most important result is a logarithmic H\"{o}lder continuity of such mappings at inner points. As a corollary, we obtained the existence of a continuous $ACL$-solution of the Beltrami equation, which is logarithmic H\"{o}lder continuous.


title: Irreversible hydrodynamic trapping by surface rollers

abstract: A colloidal particle driven by externally actuated rotation can self-propel parallel to a rigid boundary by exploiting the hydrodynamic coupling that surfaces induce between translation and rotation. As such a roller moves along the boundary it generates local vortical flows, which can be used to trap and transport passive cargo particles. However, the details and conditions for this trapping mechanism have not yet been fully understood. Here, we show that the trapping of cargo is accomplished through time-irreversible interactions between the cargo and the boundary, leading to its migration across streamlines into a steady flow vortex next to the roller. The trapping mechanism is explained analytically with a two dimensional model, investigated numerically in three dimensions for a wide range of parameters and is shown to be analogous to the deterministic lateral displacement (DLD) technique used in microfluidics for the separation of differently sized particles. The several geometrical parameters of the problem are analysed and we predict that thin, disc-like rollers offer the most favourable trapping conditions.


title: DFSSD: Deep Faults and Shallow State Duality, A Provably Strong Obfuscation Solution for Circuits with Restricted Access to Scan Chain

abstract: In this paper, we introduce DFSSD, a novel logic locking solution for sequential and FSM circuits with a restricted (locked) access to the scan chain. DFSSD combines two techniques for obfuscation: (1) Deep Faults, and (2) Shallow State Duality. Both techniques are specifically designed to resist against sequential SAT attacks based on bounded model checking. The shallow state duality prevents a sequential SAT attack from taking a shortcut for early termination without running an exhaustive unbounded model checker to assess if the attack could be terminated. The deep fault, on the other hand, provides a designer with a technique for building deep, yet key recoverable faults that could not be discovered by sequential SAT (and bounded model checker based) attacks in a reasonable time.


title: Grid-based minimization at scale: Feldman-Cousins corrections for SBN

abstract: We present a computational model for the construction of Feldman-Cousins (FC) corrections frequently used in High Energy Physics (HEP) analysis. The program contains a grid-based minimization and is written in C++. Our algorithms exploit vectorization through Eigen3, yielding a single-core speed-up of 350 compared to the original implementation, and achieve MPI data parallelism by using DIY. We demonstrate the application to scale very well at High Performance Computing (HPC) sites. We use HDF5 in conjunction with HighFive to write results of the calculation to file.


title: A strong law of large numbers for scrambled net integration

abstract: This article provides a strong law of large numbers for integration on digital nets randomized by a nested uniform scramble. The motivating problem is optimization over some variables of an integral over others, arising in Bayesian optimization. This strong law requires that the integrand have a finite moment of order $p$ for some $p>1$. Previously known results implied a strong law only for Riemann integrable functions. Previous general weak laws of large numbers for scrambled nets require a square integrable integrand. We generalize from $L^2$ to $L^p$ for $p>1$ via the Riesz-Thorin interpolation theorem


title: Work statistics across a quantum phase transition

abstract: We investigate the statistics of the work performed during a quench across a quantum phase transition using the adiabatic perturbation theory. It is shown that all the cumulants of work exhibit universal scaling behavior analogous to the Kibble-Zurek scaling for the average density of defects. Two kinds of transformations are considered: quenches between two gapped phases in which a critical point is traversed, and quenches that end near the critical point. In contrast to the scaling behavior of the density of defects, the scaling behavior of the work cumulants are shown to be qualitatively different for these two kinds of quenches. However, in both cases the corresponding exponents are fully determined by the dimension of the system and the critical exponents of the transition, as in the traditional Kibble-Zurek mechanism (KZM). Thus, our study deepens our understanding about the nonequilibrium dynamics of a quantum phase transition by revealing the imprint of the KZM on the work statistics.


title: Homogeneous Einstein metrics on non-K\"ahler C-spaces

abstract: We study homogeneous Einstein metrics on indecomposable non-K\"ahlerian C-spaces, i.e. even-dimensional torus bundles $M=G/H$ with $\mathsf{rank} G>\mathsf{rank} H$ over flag manifolds $F=G/K$ of a compact simple Lie group $G$. Based on the theory of painted Dynkin diagrams we present the classification of such spaces. Next we focus on the family \[ M_{\ell, m, n}:=\mathsf{SU}(\ell+m+n)/\mathsf{SU}(\ell)\times\mathsf{SU}(m)\times\mathsf{SU}(n)\,,\quad \ell, m, n\in\mathbb{Z}_{+} \] and examine several of its geometric properties. We show that invariant metrics on $M_{\ell, m, n}$ are not diagonal and beyond certain exceptions their parametrization depends on six real parameters. By using such an invariant Riemannian metric, we compute the diagonal and the non-diagonal part of the Ricci tensor and present explicitly the algebraic system of the homogeneous Einstein equation. For general positive integers $\ell, m, n$, by applying mapping degree theory we provide the existence of at least one $\mathsf{SU}(\ell+m+n)$-invariant Einstein metric on $M_{\ell, m, n}$. For $\ell=m$ we show the existence of two $\mathsf{SU}(2m+n)$ invariant Einstein metrics on $M_{m, m, n}$, and for $\ell=m=n$ we obtain four $\mathsf{SU}(3n)$-invariant Einstein metrics on $M_{n, n, n}$. We also examine the isometry problem for these metrics, while for a plethora of cases induced by fixed $\ell, m, n$, we provide the numerical form of all non-isometric invariant Einstein metrics.


title: VAT Compliance Incentives

abstract: In this work I clarify VAT evasion incentives through a game theoretical approach. Traditionally, evasion has been linked to the decreasing risk aversion in higher revenues (Allingham and Sandmo (1972), Cowell (1985) (1990)). I claim tax evasion to be a rational choice when compliance is stochastically more expensive than evading, even in absence of controls and sanctions. I create a framework able to measure the incentives for taxpayers to comply. The incentives here are deductions of specific VAT documented expenses from the income tax. The issue is very well known and deduction policies at work in many countries. The aim is to compute the right parameters for each precise class of taxpayers. VAT evasion is a collusive conduct between the two counterparts of the transaction. I therefore first explore the convenience for the two private counterparts to agree on the joint evasion and to form a coalition. Crucial is that compliance incentives break the agreement among the transaction participants' coalition about evading. The game solution leads to boundaries for marginal tax rates or deduction percentages, depending on parameters, able to create incentives to comply The stylized example presented here for VAT policies, already in use in many countries, is an attempt to establish a more general method for tax design, able to make compliance the "dominant strategy", satisfying the "outside option" constraint represented by evasion, even in absence of audit and sanctions. The theoretical results derived here can be easily applied to real data for precise tax design engineering.


title: Learning Similarity Metrics for Numerical Simulations

abstract: We propose a neural network-based approach that computes a stable and generalizing metric (LSiM) to compare data from a variety of numerical simulation sources. We focus on scalar time-dependent 2D data that commonly arises from motion and transport-based partial differential equations (PDEs). Our method employs a Siamese network architecture that is motivated by the mathematical properties of a metric. We leverage a controllable data generation setup with PDE solvers to create increasingly different outputs from a reference simulation in a controlled environment. A central component of our learned metric is a specialized loss function that introduces knowledge about the correlation between single data samples into the training process. To demonstrate that the proposed approach outperforms existing metrics for vector spaces and other learned, image-based metrics, we evaluate the different methods on a large range of test data. Additionally, we analyze generalization benefits of an adjustable training data difficulty and demonstrate the robustness of LSiM via an evaluation on three real-world data sets.

categories: cs.LG physics.data-an physics.flu-dyn stat.ML


title: Statistical study of magnetosheath jet-driven bow waves

abstract: When a magnetosheath jet (localized dynamic pressure enhancements) compresses ambient magnetosheath at a (relative) speed faster than the local magnetosonic speed, a bow wave or shock can form ahead of the jet. Such bow waves or shocks were recently observed to accelerate particles, thus contributing to magnetosheath heating and particle acceleration in the extended environment of Earth bow shock. To further understand the characteristics of jet-driven bow waves, we perform a statistical study to examine which solar wind conditions favor their formation and whether it is common for them to accelerate particles. We identified 364 out of 2859 (13%) magnetosheath jets to have a bow wave or shock ahead of them with Mach number typically larger than 1.1. We show that large solar wind plasma beta, weak interplanetary magnetic field (IMF) strength, large solar wind Alfven Mach number, and strong solar wind dynamic pressure present favorable conditions for their formation. We also show that magnetosheath jets with bow waves or shocks are more frequently associated with higher maximum ion and electron energies than those without them, confirming that it is common for these structures to accelerate particles. In particular, magnetosheath jets with bow waves have electron energy flux enhanced on average by a factor of 2 compared to both those without bow waves and the ambient magnetosheath. Our study implies that magnetosheath jets can contribute to shock acceleration of particles especially for high Mach number shocks. Therefore, shock models should be generalized to include magnetosheath jets and concomitant particle acceleration.


title: Electron acceleration by magnetosheath jet-driven bow waves

abstract: Magnetosheath jets are localized fast flows with enhanced dynamic pressure. When they supermagnetosonically compress the ambient magnetosheath plasma, a bow wave or shock can form ahead of them. Such a bow wave was recently observed to accelerate ions and possibly electrons. The ion acceleration process was previously analyzed, but the electron acceleration process remains largely unexplored. Here we use multi-point observations by Time History of Events and Macroscale during Substorms from three events to determine whether and how magnetosheath jet-driven bow waves can accelerate electrons. We show that when suprathermal electrons in the ambient magnetosheath convect towards a bow wave, some electrons are shock-drift accelerated and reflected towards the ambient magnetosheath and others continue moving downstream of the bow wave resulting in bi-directional motion. Our study indicates that magnetosheath jet-driven bow waves can result in additional energization of suprathermal electrons in the magnetosheath. It implies that magnetosheath jets can increase the efficiency of electron acceleration at planetary bow shocks or other similar astrophysical environments.


title: Generalised Igusa-Todorov functions and Lat-Igusa-Todorov algebras

abstract: In this paper we study a generalisation of the Igusa-Todorov functions which gives rise to a vast class of algebras satisfying the finitistic dimension conjecture. This class of algebras is called Lat-Igusa-Todorov and includes, among others, the Igusa-Todorov algebras (defined by J. Wei) and the self-injective algebras which in general are not Igusa-Todorov algebras. Finally, some applications of the developed theory are given in order to relate the different homological dimensions which have been discussed through the paper.


title: Global Convergence of Deep Networks with One Wide Layer Followed by Pyramidal Topology

abstract: Recent works have shown that gradient descent can find a global minimum for over-parameterized neural networks where the widths of all the hidden layers scale polynomially with $N$ ($N$ being the number of training samples). In this paper, we prove that, for deep networks, a single layer of width $N$ following the input layer suffices to ensure a similar guarantee. In particular, all the remaining layers are allowed to have constant widths, and form a pyramidal topology. We show an application of our result to the widely used LeCun's initialization and obtain an over-parameterization requirement for the single wide layer of order $N^2.$


title: High-precision quantum algorithms for partial differential equations

abstract: Quantum computers can produce a quantum encoding of the solution of a system of differential equations exponentially faster than a classical algorithm can produce an explicit description. However, while high-precision quantum algorithms for linear ordinary differential equations are well established, the best previous quantum algorithms for linear partial differential equations (PDEs) have complexity $\mathrm{poly}(1/\epsilon)$, where $\epsilon$ is the error tolerance. By developing quantum algorithms based on adaptive-order finite difference methods and spectral methods, we improve the complexity of quantum algorithms for linear PDEs to be $\mathrm{poly}(d, \log(1/\epsilon))$, where $d$ is the spatial dimension. Our algorithms apply high-precision quantum linear system algorithms to systems whose condition numbers and approximation errors we bound. We develop a finite difference algorithm for the Poisson equation and a spectral algorithm for more general second-order elliptic equations.


title: Retention of conceptual learning after an interactive introductory physics course

abstract: The cyclic format of the undergraduate physics curriculum depends on students' ability to recall and utilize material covered in prior courses in order to reliably build on that knowledge in later courses. However, there is evidence to suggest that people often do not retain all, or even most, of what they learned previously. How much information is retained appears to be dependent both on the individuals' approach to learning as well as the style of instruction. In particular, there is evidence to suggest that active engagement techniques in the classroom can improve students' retention of the material over time. Here, we report the findings of a longitudinal investigation of students' retention of conceptual understanding as measured by the Force and Motion Conceptual Evaluation (FMCE) following a first-semester, calculus-based introductory physics course, which features significant active engagement in both lecture and recitation. By administering the FMCE at the end of a first-semester physics course and again at the beginning of the subsequent second-semester physics course, we examine students' knowledge retention over time periods ranging from 1-15 months. We find that the shift in students' FMCE scores between these two courses is positive but corresponds to a small effect size, indicating that students retained effectively all of their conceptual learning (as measured by the FMCE). This finding largely persists even as the length of the gap between the two courses increases. We also find that, when breaking out students' performance on individual questions, the majority of students maintain their score on individual questions. Averaged over all questions, roughly a fifth of the students switched their answers from right to wrong or wrong to right on any given item.


title: Online Parameter Estimation for Safety-Critical Systems with Gaussian Processes

abstract: Parameter estimation is crucial for modeling, tracking, and control of complex dynamical systems. However, parameter uncertainties can compromise system performance under a controller relying on nominal parameter values. Typically, parameters are estimated using numerical regression approaches framed as inverse problems. However, they suffer from non-uniqueness due to existence of multiple local optima, reliance on gradients, numerous experimental data, or stability issues. Addressing these drawbacks, we present a Bayesian optimization framework based on Gaussian processes (GPs) for online parameter estimation. It uses an efficient search strategy over a response surface in the parameter space for finding the global optima with minimal function evaluations. The response surface is modeled as correlated surrogates using GPs on noisy data. The GP posterior predictive variance is exploited for smart adaptive sampling. This balances the exploration versus exploitation trade-off which is key in reaching the global optima under limited budget. We demonstrate our technique on an actuated planar pendulum and safety-critical quadrotor in simulation with changing parameters. We also benchmark our results against solvers using interior point method and sequential quadratic program. By reconfiguring the controller with new optimized parameters iteratively, we drastically improve trajectory tracking of the system versus the nominal case and other solvers.


title: Winding homology of knotoids

abstract: Knotoids were introduced by V. Turaev as open-ended knot-type diagrams that generalize knots. Turaev defined a two-variable polynomial invariant of knotoids which encompasses a generalization of the Jones knot polynomial to knotoids. We define a triply-graded homological invariant of knotoids categorifying the Turaev polynomial, called winding homology. Forgetting one of the three gradings gives a generalization of the Khovanov knot homology to knotoids.


title: Elastica catastrophe machine: theory, design and experiments

abstract: The theory, the design and the experimental validation of a catastrophe machine based on a flexible element are addressed for the first time. A general theoretical framework is developed by extending that of the classical catastrophe machines made up of discrete elastic systems. The new formulation, based on the nonlinear solution of the elastica, is enhanced by considering the concept of the universal snap surface. Among the infinite set of elastica catastrophe machines, two families are proposed and investigated to explicitly assess their features. The related catastrophe locus is disclosed in a large variety of shapes, very different from those generated by the classical counterpart. Substantial changes in the catastrophe locus properties, such as convexity and number of bifurcation points, are achievable by tuning the design parameters of the proposed machines towards the design of very efficient snapping devices. Experiments performed on the physical realization of the elastica catastrophe machine fully validate the present theoretical approach. The developed model can find applications in mechanics at different scales, for instance, in the design of new devices involving actuation or hysteresis loop mechanisms to achieve energy harvesting, locomotion, and wave mitigation.


title: A survey of statistical learning techniques as applied to inexpensive pediatric Obstructive Sleep Apnea data

abstract: Pediatric obstructive sleep apnea affects an estimated 1-5% of elementary-school aged children and can lead to other detrimental health problems. Swift diagnosis and treatment are critical to a child's growth and development, but the variability of symptoms and the complexity of the available data make this a challenge. We take a first step in streamlining the process by focusing on inexpensive data from questionnaires and craniofacial measurements. We apply correlation networks, the Mapper algorithm from topological data analysis, and singular value decomposition in a process of exploratory data analysis. We then apply a variety of supervised and unsupervised learning techniques from statistics, machine learning, and topology, ranging from support vector machines to Bayesian classifiers and manifold learning. Finally, we analyze the results of each of these methods and discuss the implications for a multi-data-sourced algorithm moving forward.


title: Ensemble Deep Learning on Large, Mixed-Site fMRI Datasets in Autism and Other Tasks

abstract: Deep learning models for MRI classification face two recurring problems: they are typically limited by low sample size, and are abstracted by their own complexity (the "black box problem"). In this paper, we train a convolutional neural network (CNN) with the largest multi-source, functional MRI (fMRI) connectomic dataset ever compiled, consisting of 43,858 datapoints. We apply this model to a cross-sectional comparison of autism (ASD) vs typically developing (TD) controls that has proved difficult to characterise with inferential statistics. To contextualise these findings, we additionally perform classifications of gender and task vs rest. Employing class-balancing to build a training set, we trained 3$\times$300 modified CNNs in an ensemble model to classify fMRI connectivity matrices with overall AUROCs of 0.6774, 0.7680, and 0.9222 for ASD vs TD, gender, and task vs rest, respectively. Additionally, we aim to address the black box problem in this context using two visualization methods. First, class activation maps show which functional connections of the brain our models focus on when performing classification. Second, by analyzing maximal activations of the hidden layers, we were also able to explore how the model organizes a large and mixed-centre dataset, finding that it dedicates specific areas of its hidden layers to processing different covariates of data (depending on the independent variable analyzed), and other areas to mix data from different sources. Our study finds that deep learning models that distinguish ASD from TD controls focus broadly on temporal and cerebellar connections, with a particularly high focus on the right caudate nucleus and paracentral sulcus.


title: Lake Ice Monitoring with Webcams and Crowd-Sourced Images

abstract: Lake ice is a strong climate indicator and has been recognised as part of the Essential Climate Variables (ECV) by the Global Climate Observing System (GCOS). The dynamics of freezing and thawing, and possible shifts of freezing patterns over time, can help in understanding the local and global climate systems. One way to acquire the spatio-temporal information about lake ice formation, independent of clouds, is to analyse webcam images. This paper intends to move towards a universal model for monitoring lake ice with freely available webcam data. We demonstrate good performance, including the ability to generalise across different winters and different lakes, with a state-of-the-art Convolutional Neural Network (CNN) model for semantic image segmentation, Deeplab v3+. Moreover, we design a variant of that model, termed Deep-U-Lab, which predicts sharper, more correct segmentation boundaries. We have tested the model's ability to generalise with data from multiple camera views and two different winters. On average, it achieves intersection-over-union (IoU) values of ~71% across different cameras and ~69% across different winters, greatly outperforming prior work. Going even further, we show that the model even achieves 60% IoU on arbitrary images scraped from photo-sharing web sites. As part of the work, we introduce a new benchmark dataset of webcam images, Photi-LakeIce, from multiple cameras and two different winters, along with pixel-wise ground truth annotations.


title: An Entropy-Maximization Approach to Automated Training Set Generation for Interatomic Potentials

abstract: Machine learning (ML)-based interatomic potentials are currently garnering a lot of attention as they strive to achieve the accuracy of electronic structure methods at the computational cost of empirical potentials. Given their generic functional forms, the transferability of these potentials is highly dependent on the quality of the training set, the generation of which is a highly labor-intensive activity. Good training sets should at once contain a very diverse set of configurations while avoiding redundancies that incur cost without providing benefits. We formalize these requirements in a local entropy maximization framework and propose an automated sampling scheme to sample from this objective function. We show that this approach generates much more diverse training sets than unbiased sampling and is competitive with hand-crafted training sets.


title: CBIR using features derived by Deep Learning

abstract: In a Content Based Image Retrieval (CBIR) System, the task is to retrieve similar images from a large database given a query image. The usual procedure is to extract some useful features from the query image, and retrieve images which have similar set of features. For this purpose, a suitable similarity measure is chosen, and images with high similarity scores are retrieved. Naturally the choice of these features play a very important role in the success of this system, and high level features are required to reduce the semantic gap. In this paper, we propose to use features derived from pre-trained network models from a deep-learning convolution network trained for a large image classification problem. This approach appears to produce vastly superior results for a variety of databases, and it outperforms many contemporary CBIR systems. We analyse the retrieval time of the method, and also propose a pre-clustering of the database based on the above-mentioned features which yields comparable results in a much shorter time in most of the cases.


title: A relative-error inertial-relaxed inexact projective splitting algorithm

abstract: For solving structured monotone inclusion problems involving the sum of finitely many maximal monotone operators, we propose and study a relative-error inertial-relaxed inexact projective splitting algorithm. The proposed algorithm benefits from a combination of inertial and relaxation effects, which are both controlled by parameters within a certain range. We propose sufficient conditions on these parameters and study the interplay between them in order to guarantee weak convergence of sequences generated by our algorithm. Additionally, the proposed algorithm also benefits from inexact subproblem solution within a relative-error criterion. Illustrative numerical experiments on LASSO problems indicate some improvement when compared with previous (noninertial and exact) versions of projective splitting.


title: A note on devising HDG+ projections on polyhedral elements

abstract: In this note, we propose a simple way of constructing HDG+ projections on polyhedral elements. The projections enable us to analyze the Lehrenfeld-Sch\"{o}berl HDG (HDG+) methods in a very concise manner, and make many existing analysis techniques of standard HDG methods reusable for HDG+. The novelty here is an alternative way of constructing the projections without using $M$-decompositions as a middle step. This extends our previous results [S. Du and F.-J. Sayas, SpringerBriefs in Mathematics (2019)] (elliptic problems) and [S. Du and F.-J. Sayas, Math.~Comp.~89~(2020), 1745-1782] (elasticity) to polyhedral meshes.


title: The interconnectedness of the economic content in the speeches of the US Presidents

abstract: The speeches stated by influential politicians can have a decisive impact on the future of a country. In particular, the economic content of such speeches affects the economy of countries and their financial markets. For this reason, we examine a novel dataset containing the economic content of 951 speeches stated by 45 US Presidents from George Washington (April 1789) to Donald Trump (February 2017). In doing so, we use an economic glossary carried out by means of text mining techniques. The goal of our study is to examine the structure of significant interconnections within a network obtained from the economic content of presidential speeches. In such a network, nodes are represented by talks and links by values of cosine similarity, the latter computed using the occurrences of the economic terms in the speeches. The resulting network displays a peculiar structure made up of a core (i.e. a set of highly central and densely connected nodes) and a periphery (i.e. a set of non-central and sparsely connected nodes). The presence of different economic dictionaries employed by the Presidents characterize the core-periphery structure. The Presidents' talks belonging to the network's core share the usage of generic (non-technical) economic locutions like "interest" or "trade". While the use of more technical and less frequent terms characterizes the periphery (e.g. "yield" ). Furthermore, the speeches close in time share a common economic dictionary. These results together with the economics glossary usages during the US periods of boom and crisis provide unique insights on the economic content relationships among Presidents' speeches.


title: Decoration of growth sector boundaries with single nitrogen vacancy centres in as-grown single crystal HPHT synthetic diamond

abstract: Large (> 100 mm$^3$), relatively pure (type II) and low birefringence single crystal diamond can be produced by high pressure high temperature (HPHT) synthesis. In this study we examine a HPHT sample of good crystalline perfection, containing less than 1 ppb (part per billion carbon atoms) of boron impurity atoms in the {001} growth sector and only tens of ppb of nitrogen impurity atoms. It is shown that the boundaries between {111} and {113} growth sectors are decorated by negatively charged nitrogen vacancy centres (NV$^-$): no decoration is observed at any other type of growth sector interface. This decoration can be used to calculated the relative {111} and {113} growth rates. The bulk (001) sector contains concentrations of luminescent point defects (excited with 488 and 532 nm wavelengths) below 10$^{11}$ cm$^{-3}$ (10$^{-3}$ ppb). We observe the negatively charged silicon-vacancy (SiV$^-$) defect in the bulk {111} sectors along with a zero phonon line emission associated with a nickel defect at 884 nm (1.40 eV). No preferential orientation is seen for either NV$^-$ or SiV$^-$ defects, but the nickel related defect is oriented with its trigonal axis along the <111> sector growth direction. Since the NV$^-$ defect is expected to readily re-orientate at HPHT diamond growth temperatures, no preferential orientation is expected for this defect but the lack of preferential orientation of SiV$^-$ in {111} sectors is not explained.


title: Microstructure Design of Low-Melting-Point Alloy (LMPA)/ Polymer Composites for Dynamic Dry Adhesion Tuning in Soft Gripping

abstract: Tunable dry adhesion is a crucial mechanism in compliant manipulation. The gripping force, mainly originated from the van der Waals force between the adhesive composite and the object to be gripped, can be controlled by reversibly varying the physical properties (e.g., stiffness) of the composite via external stimuli. The maximal gripping force Fmax and its tunability depend on, among other factors, the stress distribution on the gripping interface and its fracture dynamics (during detaching), which in turn are determined by the composite microstructure. Here, we present a computational framework for the modeling and design of a class of binary smart composites containing a porous low-melting-point alloy (LMPA) phase and a polymer phase, in order to achieve desirable dynamically tunable dry adhesion. In particular, we employ spatial correlation functions to quantify, model and represent the complex bi-continuous microstructure of the composites, from which a wide spectrum of realistic virtual 3D composite microstructures can be generated using stochastic optimization. A recently developed volume-compensated lattice-particle (VCLP) method is then employed to model the dynamic interfacial fracture process to compute Fmax for different composite microstructures. We focus on the interface defect tuning (IDT) mechanism for dry adhesion tuning enabled by the composite, in which the thermal expansion of the LMPA phase due to Joule heating initializes small cracks on the adhesion interface, subsequently causing the detachment of the gripper from the object due to interfacial fracture. We find that for an optimal microstructure among the ones studied here, a 10-fold dynamic tuning of Fmax before and after the thermal expansion of the LMPA phase can be achieved. Our computational results can provide valuable guidance for experimental fabrication of the LMPA-polymer composites.


title: QoS Evaluation and Prediction for C-V2X Communication in Commercially-Deployed LTE and Mobile Edge Networks

abstract: Cellular vehicle-to-everything (C-V2X) communication is a key enabler for future cooperative automated driving and safety-related applications. The requirements they demand in terms of Quality of Service (QoS) performance vary according to the use case. For instance, Day-1 applications such as Emergency Brake Light warning have less strict requirements than remote driving. In this paper, we seek to answer two questions: Are current LTE networks ready to support Day-1 applications at all times? And, can underperforming situations be reliably predicted based on GPS and network-related information? To address these questions, we first implement a system that collects positioning data and LTE key performance indicators (KPIs) with a higher time resolution than commercial off-the-shelf LTE modems, while simultaneously measuring the end-to-end (E2E) delay of an LTE network. We then use this system to assess the readiness of multiple mobile network operators (MNOs) and a live Mobile Edge Computing (MEC) deployment in an urban scenario. For evaluating whether an adaptable operation is possible in adverse circumstances, e.g., by performing hybrid networking or graceful degradation, we finally use Machine Learning to generate a client-based QoS predictor and forecast the achievable QoS levels.


title: Observational nonidentifiability, generalized likelihood and free energy

abstract: We study the parameter estimation problem in mixture models with observational nonidentifiability: the full model (also containing hidden variables) is identifiable, but the marginal (observed) model is not. Hence global maxima of the marginal likelihood are (infinitely) degenerate and predictions of the marginal likelihood are not unique. We show how to generalize the marginal likelihood by introducing an effective temperature, and making it similar to the free energy. This generalization resolves the observational nonidentifiability, since its maximization leads to unique results that are better than a random selection of one degenerate maximum of the marginal likelihood or the averaging over many such maxima. The generalized likelihood inherits many features from the usual likelihood, e.g. it holds the conditionality principle, and its local maximum can be searched for via suitably modified expectation-maximization method. The maximization of the generalized likelihood relates to entropy optimization.

categories: stat.ML cond-mat.stat-mech cs.LG physics.data-an


title: Photon directional profile from stimulated decay of axion clouds with arbitrary momentum distributions

abstract: We model clusters of axions with spherically symmetric spacial but arbitrary momentum distributions and study the directional profile of photos produced in their evolution through spontaneous and stimulated decay of axions via the process a --> gamma + gamma. Several specific examples are presented.


title: Numerical study of Zakharov-Kuznetsov equations in two dimensions

abstract: We present a detailed numerical study of solutions to the (generalized) Zakharov-Kuznetsov equation in two spatial dimensions with various power nonlinearities. In the $L^{2}$-subcritical case, numerical evidence is presented for the stability of solitons and the soliton resolution for generic initial data. In the $L^2$-critical and supercritical cases, solitons appear to be unstable against both dispersion and blow-up. It is conjectured that blow-up happens in finite time and that blow-up solutions have some resemblance of being self-similar, i.e., the blow-up core forms a rightward moving self-similar type rescaled profile with the blow-up happening at infinity in the critical case and at a finite location in the supercritical case. In the $L^{2}$-critical case, the blow-up appears to be similar to the one in the $L^{2}$-critical generalized Korteweg-de Vries equation with the profile being a dynamically rescaled soliton.


title: Singular radial solutions for the Lin-Ni-Takagi equation

abstract: We study singular radially symmetric solution to the Lin-Ni-Takagi equation for a supercritical power non-linearity in dimension $N\geq 3$. It is shown that for any ball and any $k \geq 0$, there is a singular solution that satisfies Neumann boundary condition and oscillates at least $k$ times around the constant equilibrium. Moreover, we show that the Morse index of the singular solution is finite or infinite if the exponent is respectively larger or smaller than the Joseph--Lundgren exponent.


title: DLITE: The Discounted Least Information Theory of Entropy

abstract: We propose an entropy-based information measure, namely the Discounted Least Information Theory of Entropy (DLITE), which not only exhibits important characteristics expected as an information measure but also satisfies conditions of a metric. Classic information measures such as Shannon Entropy, KL Divergence, and Jessen-Shannon Divergence have manifested some of these properties while missing others. This work fills an important gap in the advancement of information theory and its application, where related properties are desirable.


title: Strategyproof Facility Location Mechanisms with Richer Action Spaces

abstract: We study facility location problems where agents control multiple locations and when reporting their locations can choose to hide some locations (hiding), report some locations more than once (replication) and lie about their locations (manipulation). We fully characterize all facility location mechanisms that are anonymous, efficient, and strategyproof with respect to the richer strategic behavior for this setting. We also provide a characterization with respect to manipulation only. This is the first, to the best of our knowledge, characterization result for the strategyproof facility location mechanisms where each agent controls multiple locations.


title: Informative Path Planning for Mobile Sensing with Reinforcement Learning

abstract: Large-scale spatial data such as air quality, thermal conditions and location signatures play a vital role in a variety of applications. Collecting such data manually can be tedious and labour intensive. With the advancement of robotic technologies, it is feasible to automate such tasks using mobile robots with sensing and navigation capabilities. However, due to limited battery lifetime and scarcity of charging stations, it is important to plan paths for the robots that maximize the utility of data collection, also known as the informative path planning (IPP) problem. In this paper, we propose a novel IPP algorithm using reinforcement learning (RL). A constrained exploration and exploitation strategy is designed to address the unique challenges of IPP, and is shown to have fast convergence and better optimality than a classical reinforcement learning approach. Extensive experiments using real-world measurement data demonstrate that the proposed algorithm outperforms state-of-the-art algorithms in most test cases. Interestingly, unlike existing solutions that have to be re-executed when any input parameter changes, our RL-based solution allows a degree of transferability across different problem instances.


title: Towards Query-Efficient Black-Box Adversary with Zeroth-Order Natural Gradient Descent

abstract: Despite the great achievements of the modern deep neural networks (DNNs), the vulnerability/robustness of state-of-the-art DNNs raises security concerns in many application domains requiring high reliability. Various adversarial attacks are proposed to sabotage the learning performance of DNN models. Among those, the black-box adversarial attack methods have received special attentions owing to their practicality and simplicity. Black-box attacks usually prefer less queries in order to maintain stealthy and low costs. However, most of the current black-box attack methods adopt the first-order gradient descent method, which may come with certain deficiencies such as relatively slow convergence and high sensitivity to hyper-parameter settings. In this paper, we propose a zeroth-order natural gradient descent (ZO-NGD) method to design the adversarial attacks, which incorporates the zeroth-order gradient estimation technique catering to the black-box attack scenario and the second-order natural gradient descent to achieve higher query efficiency. The empirical evaluations on image classification datasets demonstrate that ZO-NGD can obtain significantly lower model query complexities compared with state-of-the-art attack methods.


title: Fair Clustering with Multiple Colors

abstract: A fair clustering instance is given a data set $A$ in which every point is assigned some color. Colors correspond to various protected attributes such as sex, ethnicity, or age. A fair clustering is an instance where membership of points in a cluster is uncorrelated with the coloring of the points. Of particular interest is the case where all colors are equally represented. If we have exactly two colors, Chierrichetti, Kumar, Lattanzi and Vassilvitskii (NIPS 2017) showed that various $k$-clustering objectives admit a constant factor approximation. Since then, a number of follow up work has attempted to extend this result to a multi-color case, though so far, the only known results either result in no-constant factor approximation, apply only to special clustering objectives such as $k$-center, yield bicrititeria approximations, or require $k$ to be constant. In this paper, we present a simple reduction from unconstrained $k$-clustering to fair $k$-clustering for a large range of clustering objectives including $k$-median, $k$-means, and $k$-center. The reduction loses only a constant factor in the approximation guarantee, marking the first true constant factor approximation for many of these problems.


title: ReaxFF reactive force field study of polymerization of polymer matrix in carbon nanotube-composite system

abstract: Human transport to Mars and deep space explorations demand the development of new materials with extraordinary high performance-to-mass ratios. Promising candidates to fulfill these requirements are ultrahigh strength lightweight (UHSL) materials, which consist of polymer matrices fortified by pristine carbon nanotubes (CNTs). Previous investigations have showed that with an increase in CNT diameter, its preferred configuration changes from a circular form to a flattened shape that can be obtained in high pressure or tension conditions. The ReaxFF reactive force field can reveal detailed chemical interactions at the atomistic scale. To enable ReaxFF simulations on CNT/polymer interfaces, we trained force field parameters to capture the proper structure of flattened carbon nanotubes (flCNTs), i.e. dumbbell-like shape CNTs, against available polymer consistent force field -- interface force field (PCFF-IFF) data which had good proximity to density functional theory (DFT) data. In this study we used accelerated ReaxFF molecular dynamics simulation using the optimized force field to study the polymerization of diglycidyl ether of bisphenol F (Bis F) and diethyltoluenediamine (DEDTA) molecules in vicinity of circular and flattened CNTs. Our results indicate that the flat regions of flCNT are more favorable spots for the polymers to settle compared to curved regions due to higher binding energies. Moreover, higher dimer generation around flCNT results in more effective coating of the CNT which leads to higher load transfer in compared to circular CNT. According to our results there is a high alignment between polymers and CNT surface which is due to strong pi-pi interactions of aromatic carbon rings in the polymers and CNTs. These ReaxFF simulations indicate the capability of this method to simultaneously observe the polymerization of monomers along with their interactions with CNTs.


title: Fragmentation modelling of the August 2019 impact on Jupiter

abstract: On 7th August 2019, an impact flash lasting $\sim1$s was observed on Jupiter. The video of this event was analysed to obtain the lightcurve and determine the energy release and initial mass. We find that the impactor released a total energy of $96-151$ kilotons of TNT, corresponding to an initial mass between $190-260$ metric tonnes with a diameter between $4-10$m. We developed a fragmentation model to simulate the atmospheric breakup of the object and reproduce the lightcurve. We model three different materials: cometary, stony and metallic at speeds of $60$, $65 $ and $70$ km/s to determine the material makeup of the impacting object. The slower cases are best fit by a strong, metallic object while the faster cases require a weaker material.


title: Bivariate Continuous q-Hermite Polynomials and Deformed Quantum Serre Relations

abstract: We introduce bivariate versions of the continuous q-Hermite polynomials. We obtain algebraic properties for them (generating function, explicit expressions in terms of the univariate ones, backward difference equations and recurrence relations) and analytic properties (determining the orthogonality measure). We find a direct link between bivariate continuous q-Hermite polynomials and the star product method of [Kolb and Yakimov, Adv. Math. 2020] for quantum symmetric pairs to establish deformed quantum Serre relations for quasi-split quantum symmetric pairs of Kac-Moody type. We prove that these defining relations are obtained from the usual quantum Serre relations by replacing all monomials by multivariate orthogonal polynomials.


title: Concurrent probing of electron-lattice dephasing induced by photoexcitation in 1T-TaSeTe using ultrafast electron diffraction

abstract: It has been technically challenging to concurrently probe the electrons and the lattices in materials during non-equilibrium processes, allowing their correlations to be determined. Here, in a single set of ultrafast electron diffraction patterns taken on the charge-density-wave (CDW) material 1T-TaSeTe, we discover a temporal shift in the diffraction intensity measurements as a function of scattering angle. With the help of dynamic models and theoretical calculations, we show that the ultrafast electrons probe both the valence-electron and lattice dynamic processes, resulting in the temporal shift measurements. Our results demonstrate unambiguously that the CDW is not merely a result of the periodic lattice deformation ever-present in 1T-TaSeTe but has significant electronic origin. This method demonstrates a novel approach for studying many quantum effects that arise from electron-lattice dephasing in molecules and crystals for next-generation devices.


title: LocoGAN -- Locally Convolutional GAN

abstract: In the paper we construct a fully convolutional GAN model: LocoGAN, which latent space is given by noise-like images of possibly different resolutions. The learning is local, i.e. we process not the whole noise-like image, but the sub-images of a fixed size. As a consequence LocoGAN can produce images of arbitrary dimensions e.g. LSUN bedroom data set. Another advantage of our approach comes from the fact that we use the position channels, which allows the generation of fully periodic (e.g. cylindrical panoramic images) or almost periodic ,,infinitely long" images (e.g. wall-papers).


title: Deep Transform and Metric Learning Network: Wedding Deep Dictionary Learning and Neural Networks

abstract: On account of its many successes in inference tasks and denoising applications, Dictionary Learning (DL) and its related sparse optimization problems have garnered a lot of research interest. While most solutions have focused on single layer dictionaries, the improved recently proposed Deep DL (DDL) methods have also fallen short on a number of issues. We propose herein, a novel DDL approach where each DL layer can be formulated as a combination of one linear layer and a Recurrent Neural Network (RNN). The RNN is shown to flexibly account for the layer-associated and learned metric. Our proposed work unveils new insights into Neural Networks and DDL and provides a new, efficient and competitive approach to jointly learn a deep transform and a metric for inference applications. Extensive experiments are carried out to demonstrate that the proposed method can not only outperform existing DDL but also state-of-the-art generic CNNs.


title: Transporting Experimental Results with Entropy Balancing

abstract: We show how entropy balancing can be used for transporting experimental treatment effects from a trial population onto a target population. This method is doubly-robust in the sense that if either the outcome model or the probability of trial participation is correctly specified, then the estimate of the target population average treatment effect is consistent. Furthermore, we only require the sample moments of the effect modifiers drawn from the target population to consistently estimate the target population average treatment effect. We compared the finite-sample performance of entropy balancing with several alternative methods for transporting treatment effects between populations. Entropy balancing techniques are efficient and robust to violations of model misspecification. We also examine the results of our proposed method in an applied analysis of the Action to Control Cardiovascular Risk in Diabetes Blood Pressure (ACCORD-BP) trial transported to a sample of US adults with diabetes taken from the National Health and Nutrition Examination Survey (NHANES) cohort.


title: Trust Repairing for Human-Swarm Cooperation inDynamic Task Response

abstract: Emergency happens in human-UAV cooperation, such as criminal activity tracking and urgent needs for ground assistance. Emergency response usually has high requirements on the motion control of the multi-UAV system, by maintaining both the team performance and team behaviors. However When a UAV swarm executes tasks in a real-world environment, because of real-world factors, such as system reliability and environmental disturbances, some robots in the swarm will behave abnormally, such as slow flocking speed, wrong heading direction, or poor spatial relations. In the meanwhile, incorrect trust between human and UAV swarm could map the abnormal behavior of faulty robot to the whole swarm and request a time-consuming intervention from human supervisor, damage the UAV swarm response for a dynamic task, even evolve to a failure of task because of accumulated error. To correct reflect the trust between humans and UAV swarm and rebuild the trust to improve the performance caused by incorrect trust. We propose a dynamic trust repair model. The dynamic trust model focus on human-supervisory UAV system which can help UAV swarm to reduce the negative influence from faulty UAV on the performance of the UAV swarm, get a flexible reaction and stable human-supervisory UAV task performance. Results show that trust model could improve the performance of the swarm for dynamic task response and regain human trust.


title: Scaling Up Electronic Structure Calculations on Quantum Computers: The Frozen Natural Orbital Based Method of Increments

abstract: The method of increments and frozen natural orbital (MI-FNO) framework is introduced to help expedite the application of noisy, intermediate-scale quantum~(NISQ) devices for quantum chemistry simulations. The MI-FNO framework provides a systematic reduction of the occupied and virtual orbital spaces for quantum chemistry simulations. The correlation energies of the resulting increments from the MI-FNO reduction can then be solved by various algorithms, including quantum algorithms such as the phase estimation algorithm and the variational quantum eigensolver (VQE). The unitary coupled-cluster singles and doubles VQE framework is used to obtain correlation energies for the case of small molecules (i.e., BeH$_2$, CH$_4$, NH$_3$, H$_2$O, and HF) using the cc-pVDZ basis set. The quantum resource requirements are estimated for a constrained geometry complex (CGC) catalyst that is utilized in industrial settings for the polymerization of $\alpha$-olefins. We show that the MI-FNO approach provides a significant reduction in the qubit requirements relative to the full system simulations. We propose that the MI-FNO framework can create scalable examples of quantum chemistry problems that are appropriate for assessing the progress of NISQ devices.


title: Mapping the intrinsic photocurrent streamlines through micromagnetic heterostructure devices

abstract: Like air flowing over a wing, optimizing the flow of electronic charge is essential to the operation of nanoscale devices. Unfortunately, the delicate interplay of charge, spin, and heat in complex devices has precluded detailed imaging of charge flow. Here, we report on the visualization of intrinsic charge current streamlines through yttrium iron garnet micromagnetic heterostructures. Scanning photovoltage microscopy of precisely designed devices leads to striking spatial patterns, with prominent photovoltage features emerging in corners and narrow constrictions. These patterns, which evolve continuously with rotation of an external magnetic field, enable rich spatial mapping of fluid-like flow. Taking inspiration from aerodynamic Clark Y airfoils, we engineer micromagnetic wing shaped devices, called electrofoils, which allow us to precisely contort, compress and decompress flowlines of electronic charge.120 (39) e2221815120


title: Temperature-dependent mechanical properties of ZrC and HfC from first principles

abstract: In order to gain insight into the effect of elevated temperature on mechanical performance of ZrC and HfC, the temperature-dependent elastic constants of these two carbides have been systematically studied using two different methods (first-principles qusi-harmonic approximation and first-principles qusi-static approximation). Isoentropic C$_{11}$ gradually decrease and the isoentropic C$_{12}$ slightly increases for ZrC and HfC under temperature. While the isoentropic C$_{44}$ of both ZrC and HfC is insensitive to the temperature. In both ZrC and HfC cases, the decline of temperature-dependent C$_{11}$ calculated from QHA is more pronounced than from QSA. Temperature effects on modulus of elasticity, poisson's ratio, elastic anisotropy, hardness and fracture toughness are further explored and discussed. The results indicate that degradation of B, G and E estimated from QHA is more serious than these from QSA at high temperature. From room temperature to 2500 K, the theoretical decreasing slope (from QHA) of B, G, and E is 0.38 GPa$\cdot$K$^{-1}$, 0.30 GPa$\cdot$K$^{-1}$, and 0.32 GPa$\cdot$K$^{-1}$ for ZrC and 0.33 GPa$\cdot$K$^{-1}$, 0.29 GPa$\cdot$K$^{-1}$, and 0.30 GPa$\cdot$K$^{-1}$ for HfC.


title: Repair rate lower bounds for distributed storage

abstract: One of the primary objectives of a distributed storage system is to reliably store a large amount $dsize$ of source data for a long duration using a large number $N$ of unreliable storage nodes, each with capacity $nsize$. The storage overhead $\beta$ is the fraction of system capacity available beyond $dsize$, i.e., $\beta = 1- \frac{dsize}{N \cdot nsize}$. Storage nodes fail randomly over time and are replaced with initially empty nodes, and thus data is erased from the system at an average rate $erate = \lambda \cdot N \cdot nsize$, where $1/\lambda$ is the average lifetime of a node before failure. To maintain recoverability of the source data, a repairer continually reads data over a network from nodes at some average rate $rrate$, and generates and writes data to nodes based on the read data. The main result is that, for any repairer, if the source data is recoverable at each point in time then it must be the case that $rrate \ge \frac{erate}{2 \cdot \beta}$ asymptotically as $N$ goes to infinity and beta goes to zero. This inequality provides a fundamental lower bound on the average rate that any repairer needs to read data from the system in order to maintain recoverability of the source data.


title: Empirical Policy Evaluation with Supergraphs

abstract: We devise and analyze algorithms for the empirical policy evaluation problem in reinforcement learning. Our algorithms explore backward from high-cost states to find high-value ones, in contrast to forward approaches that work forward from all states. While several papers have demonstrated the utility of backward exploration empirically, we conduct rigorous analyses which show that our algorithms can reduce average-case sample complexity from $O(S \log S)$ to as low as $O(\log S)$.


title: CAUSE: Learning Granger Causality from Event Sequences using Attribution Methods

abstract: We study the problem of learning Granger causality between event types from asynchronous, interdependent, multi-type event sequences. Existing work suffers from either limited model flexibility or poor model explainability and thus fails to uncover Granger causality across a wide variety of event sequences with diverse event interdependency. To address these weaknesses, we propose CAUSE (Causality from AttribUtions on Sequence of Events), a novel framework for the studied task. The key idea of CAUSE is to first implicitly capture the underlying event interdependency by fitting a neural point process, and then extract from the process a Granger causality statistic using an axiomatic attribution method. Across multiple datasets riddled with diverse event interdependency, we demonstrate that CAUSE achieves superior performance on correctly inferring the inter-type Granger causality over a range of state-of-the-art methods.


title: Hawking radiation in optics and beyond

abstract: Hawking radiation was originally proposed in astrophysics, but it has been generalized and extended to other physical systems receiving the name of analogue Hawking radiation. In the last two decades, several attempts have been made to measure it in a laboratory, one of the most successful systems is in optics. Light interacting in a dielectric material causes an analogue Hawking effect, in fact, its stimulated version has already been detected and the search for the spontaneous signal is currently ongoing. We briefly review the general derivation of Hawking radiation, then we focus on the optical analogue and present some novel numerical results. Finally, we call for a generalization of the term Hawking radiation.


title: Framework of goals for writing in physics lab classes

abstract: Writing is an integral part of the process of science. In the undergraduate physics curriculum, the most common place that students engage with scientific writing is in lab classes, typically through lab notebooks, reports, and proposals. There has not been much research on why and how we include writing in physics lab classes, and instructors may incorporate writing for a variety of reasons. Through a broader study of multiweek projects in advanced lab classes, we have developed a framework for thinking about and understanding the role of writing in lab classes. This framework defines and describes the breadth of goals for incorporating writing in lab classes, and is a tool we can use to begin to understand why, and subsequently how, we teach scientific writing in physics.


title: Initial and Boundary Value Problems for the Caputo Fractional Self-Adjoint Difference Equations

abstract: In this paper we develop the theory of initial and boundary value problems for the self-adjoint nabla fractional difference equation containing a Caputo fractional nabla difference that is given by \[ \nabla[p(t+1)\nabla_{a*}^\nu x(t+1)] + q(t)x(t) = h(t), \] where $0 < \nu \leq 1$. We give an introduction to the nabla fractional calculus with Caputo fractional differences. We investigate properties of the specific self-adjoint nabla fractional difference equation given above. We prove existence and uniqueness theorems for both initial and boundary value problems under appropriate conditions. We introduce the definition of a Cauchy function which allows us to give a variation of constants formula for solving initial value problems. We then show that this Cauchy function is important in finding a Green's function for a boundary value problem with Sturm-Liouville type boundary conditions. Several inequalities concerning a certain Green's function are derived. These results are important in using fixed point theorems for proving the existence of solutions to boundary value problems for nonlinear fractional equations related to our linear self-adjoint equation.


title: Surface Susceptibility Synthesis of Metasurface Holograms for creating Electromagnetic Illusions

abstract: A systematic approach is presented to exploit the rich field transformation capabilities of Electromagnetic (EM) metasurfaces for creating a variety of illusions using the concept of metasurface holograms. A system level approach to metasurface hologram synthesis is presented here, in which the hologram is co-designed with the desired object to be projected. A structured approach for the classification of the creation of EM illusions is proposed for better organization and tractability of the overall synthesis problem. The deliniation is in terms of the initial incident (reference) illumination of the object to be recreated (front/back-lit), the position of illusion (posterior/anterior), and the illumination used to create the illusion (front/back). Therefore the classification is based on the specific relationship between the reference object to be recreated, the observer measuring the object, the orientation and placement of the reference and illumination field, and the desired placement of the metasurface hologram creating a virtual image. In the paper a general design procedure to synthesize metasurface holograms is presented based on Integral Equations (IE) and Generalized Sheet Transition Conditions (GSTCs), where the metasurface hologram is described as zero thickness sheet with tensorial surface susceptibility densities. Several selected configurations are chosen to illustrate various aspects of the hologram creation in 2D, along with a novel numerical technique to artificially reverse propagate the scattered fields, required in the synthesis process. Finally, the impact of the metasurface size and the illumination field strength on the quality of the reconstructed scattered fields is also discussed.