title: On the Choice of Instruments in Mixed Frequency Specification Tests

abstract: Time averaging has been the traditional approach to handle mixed sampling frequencies. However, it ignores information possibly embedded in high frequency. Mixed data sampling (MIDAS) regression models provide a concise way to utilize the additional information in high-frequency variables. In this paper, we propose a specification test to choose between time averaging and MIDAS models, based on a Durbin-Wu-Hausman test. In particular, a set of instrumental variables is proposed and theoretically validated when the frequency ratio is large. As a result, our method tends to be more powerful than existing methods, as reconfirmed through the simulations.


title: Melding the Data-Decisions Pipeline: Decision-Focused Learning for Combinatorial Optimization

abstract: Creating impact in real-world settings requires artificial intelligence techniques to span the full pipeline from data, to predictive models, to decisions. These components are typically approached separately: a machine learning model is first trained via a measure of predictive accuracy, and then its predictions are used as input into an optimization algorithm which produces a decision. However, the loss function used to train the model may easily be misaligned with the end goal, which is to make the best decisions possible. Hand-tuning the loss function to align with optimization is a difficult and error-prone process (which is often skipped entirely). We focus on combinatorial optimization problems and introduce a general framework for decision-focused learning, where the machine learning model is directly trained in conjunction with the optimization algorithm to produce high-quality decisions. Technically, our contribution is a means of integrating common classes of discrete optimization problems into deep learning or other predictive models, which are typically trained via gradient descent. The main idea is to use a continuous relaxation of the discrete problem to propagate gradients through the optimization procedure. We instantiate this framework for two broad classes of combinatorial problems: linear programs and submodular maximization. Experimental results across a variety of domains show that decision-focused learning often leads to improved optimization performance compared to traditional methods. We find that standard measures of accuracy are not a reliable proxy for a predictive model's utility in optimization, and our method's ability to specify the true goal as the model's training objective yields substantial dividends across a range of decision problems.


title: Modeling Photoionized Turbulent Material in the Circumgalactic Medium

abstract: The circumgalactic medium (CGM) of nearby star-forming galaxies show clear indications of \ion{O}{6} absorption accompanied by little to no \ion{N}{5} absorption. This unusual spectral signature, accompanied by absorption from lower ionization state species whose columns vary by orders of magnitude along \st{difference} \textbf{different} sightlines, indicates that the CGM must be viewed as a dynamic, multiphase medium, such as occurs in the presence of turbulence. To explore this possibility, we carry out a series of chemodynamical simulations of a isotropic turbulent media, using the MAIHEM package. The simulations assume a metallicity of $0.3\ Z_{\odot}$ and a redshift zero metagalatic UV background, and they track ionizations, recombinations, and species-by-species radiative cooling for a wide range of elements. We find that turbulence with a one-dimensional velocity dispersion of $\sigma_{1D} \approx 60$ km/s replicates many of the observed features within the CGM, such as clumping of low ionization-state ions and the existence of \ion{O}{6} at moderate ionization parameters. However, unlike observations, \ion{N}{5} often arises in our simulations with derived column densities of a similar magnitude to those of \ion{O}{6}. While higher values of $\sigma_{1D}$ lead to a thermal runaway in our isotropic simulations, this would not be the case in stratified media, and thus we speculate that more complex models of the turbulence may well match the absence of \ion{N}{5} in the CGM of star-forming galaxies.


title: Discovery of Parabolic Microresonators Produced via Fiber Tapering

abstract: We demonstrate a new method for creation of surface nanoscale axial photonics (SNAP) microresonators with harmonic profiles via fiber tapering in a laser-heated microfurnace. The simple procedure makes microresonators that support hundreds of axial modes with good spacing uniformity, yielding a promising prospective method for fabricating miniature frequency comb generators and dispersionless delay lines.


title: Quintessence Saves Higgs Instability

abstract: We study a model where quintessence potential $e^{-\xi\phi}$ coupled to Higgs potential. We calculate the evolution of the quintessence, and track the running of the effective Higgs self-coupling. We find it slightly larger than that of the standard model in the past. Requiring the electroweak vacuum to be absolutely stable in inflationary era, we find a lower bound $\xi> 0.35\pm 0.05$, where the uncertainty is mainly from the measurement of the top quark mass. This lower bound, together with the upper bound from the observation for dark energy $\xi\lesssim0.6$, narrows down the parameter space and makes it possible to test this model in the near future. Interestingly, the bound on $\xi$, if actually shown to be the case by observation, supports the recently proposed Swampland Conjecture.


title: A non-discrete space $X$ with $C_p(X)$ Menger at infinity

abstract: In a paper by Bella, Tokg\"os and Zdomskyy it is asked whether there exists a Tychonoff space $X$ such that the remainder of $C_p(X)$ in some compactification is Menger but not $\sigma$-compact. In this paper we prove that it is consistent that such space exists and in particular its existence follows from the existence of a Menger ultrafilter.


title: Feasibility and coordination of multiple mobile vehicles with mixed equality and inequality constraints

abstract: We consider the problem of feasible coordination control for multiple homogeneous or heterogeneous mobile vehicles subject to various constraints (nonholonomic motion constraints, holonomic coordination constraints, equality/inequality constraints etc). We develop a general framework involving differential-algebraic equations and viability theory to describe and determine coordination feasibility for a coordinated motion control under heterogeneous vehicle dynamics and different types of coordination constraints. If a solution exists for the derived differential-algebraic equations and/or inequalities, a heuristic algorithm is proposed for generating feasible trajectories for each individual vehicle. In case studies on coordinating two vehicles, we derive analytical solutions to motion generation for two-vehicle groups consisting of car-like vehicles, unicycle vehicles, or vehicles with constant speeds, which serve as benchmark coordination tasks for more complex vehicle groups. We show several simulation experiments on multi-vehicle coordination under various constraints to validate the theory and the effectiveness of the proposed schemes.


title: BPE and computer-extracted parenchymal enhancement for breast cancer risk, response monitoring, and prognosis

abstract: Functional behavior of breast cancer - representing underlying biology - can be analyzed using MRI. The most widely used breast MR imaging protocol is dynamic contrast-enhanced T1-weighted imaging. The cancer enhances on dynamic contrast-enhanced MR imaging because the contrast agent leaks from the leaky vessels into the interstitial space. The contrast agent subsequently leaks back into the vascular space, creating a washout effect. The normal parenchymal tissue of the breast can also enhance after contrast injection. This enhancement generally increases over time. Typically, a radiologist assesses this background parenchymal enhancement (BPE) using the Breast Imaging Reporting and Data System (BI-RADS). According to the BI-RADS, BPE refers to the volume of enhancement and the intensity of enhancement and is divided in four incremental categories: minimal, mild, moderate, and marked. Researchers have developed semi-automatic and automatic methods to extract properties of BPE from MR images. For clarity, in this syllabus the BI-RADS definition will be referred to as BPE, whereas the computer-extracted properties will not. Both BPE and computer-extracted parenchymal enhancement properties have been linked to screening and diagnosis, hormone status and age, risk of development of breast cancer, response monitoring, and prognosis.


title: Dynamic Model of Planar Sliding

abstract: In this paper, we present a principled method to model general planar sliding motion with distributed convex contact patch. The effect of contact patch with indeterminate pressure distribution can be equivalently modeled as the contact wrench at one point contact. We call this point equivalent contact point. Our dynamic model embeds ECP within the equations of slider's motion and friction model which approximates the distributed contact patch, and eventually brings us a system of quadratic equations. This discrete-time dynamic model allows us to solve for the two components of tangential friction impulses, the friction moment and the slip speed. The state of the slider as well as the ECP can be computed by solving a system of linear equations once the contact impulses are computed. In addition, we derive the closed form solutions for the state of slider for quasi-static motion. Furthermore, in pure translation case, based on the discrete-time model, we present the closed form expressions for the friction impulses the slider suffers and the state of it at each time step. Simulation examples are shown to demonstrate the validity of our approach.


title: Topological superconductivity in full shell proximitized nanowires

abstract: We consider a new model system supporting Majorana zero modes based on semiconductor nanowires with a full superconducting shell. We demonstrate that, in the presence of spin-orbit coupling in the semiconductor induced by a radial electric field, the winding of the superconducting order parameter leads to a topological phase supporting Majorana zero modes. The topological phase persists over a large range of chemical potentials and can be induced by a predictable and weak magnetic field piercing the cylinder. The system can be readily realized in semiconductor nanowires covered by a full superconducting shell, opening a pathway for realizing topological quantum computing proposals.


title: Flux-induced Majorana modes in full-shell nanowires

abstract: We demonstrate a novel means of creating Majorana zero modes using magnetic flux applied to a full superconducting shell surrounding a semiconducting nanowire core, unifying approaches based on proximitized nanowires and vortices in topological superconductors. In the destructive Little-Parks regime, reentrant regions of superconductivity are associated with integer number of phase windings in the shell. Tunneling into the core reveals a hard induced gap near zero applied flux, corresponding to zero phase winding, and a gapped region with a discrete zero-energy state for flux around {\Phi}_0 = h/2e, corresponding to 2{\pi} phase winding. Coulomb peak spacing in full-shell islands around one applied flux shows exponentially decreasing deviation from 1e periodicity with device length, consistent with the picture of Majorana modes located at the ends of the wire.


title: Nanoscale design of the local density of optical states

abstract: We propose a design concept for tailoring the local density of optical states (LDOS) in dielectric nanostructures, based on the phase distribution of the scattered optical fields induced by point-like emitters. First we demonstrate that the LDOS can be expressed in terms of a coherent summation of constructive and destructive contributions. By using an iterative approach, dielectric nanostructures can be designed to effectively remove the destructive terms. In this way dielectric Mie resonators, featuring low LDOS for electric dipoles, can be reshaped to enable enhancements of three orders of magnitude. To demonstrate the generality of the method, we also design nanocavities that enhance the radiated power of a circular dipole, a quadrupole and an arbitrary collection of coherent dipoles. Our concept provides a powerful tool for high-performance dielectric resonators, and affords fundamental insights into light-matter coupling at the nanoscale.


title: A Statistical Learning Approach to Ultra-Reliable Low Latency Communication

abstract: Mission-critical applications require Ultra-Reliable Low Latency (URLLC) wireless connections, where the packet error rate (PER) goes down to $10^{-9}$. Fulfillment of the bold reliability figures becomes meaningful only if it can be related to a statistical model in which the URLLC system operates. However, this model is generally not known and needs to be learned by sampling the wireless environment. In this paper we treat this fundamental problem in the simplest possible communication-theoretic setting: selecting a transmission rate over a dynamic wireless channel in order to guarantee high transmission reliability. We introduce a novel statistical framework for design and assessment of URLLC systems, consisting of three key components: (i) channel model selection; (ii) learning the model using training; (3) selecting the transmission rate to satisfy the required reliability. As it is insufficient to specify the URLLC requirements only through PER, two types of statistical constraints are introduced, Averaged Reliability (AR) and Probably Correct Reliability (PCR). The analysis and the evaluations show that adequate model selection and learning are indispensable for designing consistent physical layer that asymptotically behaves as if the channel was known perfectly, while maintaining the reliability requirements in URLLC systems.


title: Spectral data for spin Higgs bundles

abstract: In this paper we determine the spectral data parametrizing Higgs bundles in a generic fiber of the Hitchin map for the case where the structure group is the special Clifford group with fixed Clifford norm. These are spin and "twisted" spin Higgs bundles. The method used relates variations in spectral data with respect to the Hecke transformations for orthogonal bundles introduced by Abe. The explicit description also recovers a result from the geometric Langlands program which states that the fibers of the Hitchin map are the dual abelian varieties to the corresponding fibers of the moduli spaces of projective orthogonal Higgs bundles (in the even case) and projective symplectic Higgs bundles (in the odd case).


title: The hot Jupiter period-mass distribution as a signature of in situ formation

abstract: More than two decades after the widespread detection of Jovian-class planets on short-period orbits around other stars, their dynamical origins remain imperfectly understood. In the traditional narrative, these highly irradiated giant planets, like Jupiter and Saturn, are envisioned to have formed at large stello-centric distances and to have subsequently undergone large-scale orbital decay. Conversely, more recent models propose that a large fraction of hot Jupiters could have formed via rapid gas accretion in their current orbital neighborhood. In this study, we examine the period-mass distribution of close-in giant planets, and demonstrate that the inner boundary of this population conforms to the expectations of the in-situ formation scenario. Specifically, we show that if conglomeration unfolds close to the disk's inner edge, the semi-major axis - mass relation of the emergent planets should follow a power law $a \propto M^{-2/7}$ - a trend clearly reflected in the data. We further discuss corrections to this relationship due to tidal decay of planetary orbits. Although our findings do not discount orbital migration as an active physical process, they suggest that the characteristic range of orbital migration experienced by giant planets is limited.


title: SocialRobot: Towards a Personalized Elderly Care Mobile Robot

abstract: SocialRobot is a collaborative European project, which focuses on providing a practical and interactive solution to improve the quality of life of elderly people. Having this in mind, a state of the art robotic mobile platform has been integrated with virtual social care technology to meet the elderly individual needs and requirements, following a human centered approach. In this short paper, we make an overview of SocialRobot, the developed architecture and the human-robot interactive scenarios being prepared and tested in the framework of the project for dissemination and exploitation purposes.


title: Identifying Quantum Phase Transitions using Artificial Neural Networks on Experimental Data

abstract: Machine learning techniques such as artificial neural networks are currently revolutionizing many technological areas and have also proven successful in quantum physics applications. Here we employ an artificial neural network and deep learning techniques to identify quantum phase transitions from single-shot experimental momentum-space density images of ultracold quantum gases and obtain results, which were not feasible with conventional methods. We map out the complete two-dimensional topological phase diagram of the Haldane model and provide an accurate characterization of the superfluid-to-Mott-insulator transition in an inhomogeneous Bose-Hubbard system. Our work points the way to unravel complex phase diagrams of general experimental systems, where the Hamiltonian and the order parameters might not be known.

categories: cond-mat.quant-gas cond-mat.dis-nn cond-mat.mes-hall quant-ph


title: Studying the morphology of HI isodensity surfaces during reionization using Shapefinders and percolation analysis

abstract: Minkowski functionals and Shapefinders shed light on the connectedness of large-scale structure by determining its topology and morphology. We use a sophisticated code, SURFGEN2, to measure the Minkowski functionals and Shapefinders of individual clusters by modelling cluster surfaces using the 'Marching Cube 33' triangulation algorithm. In this paper, we study the morphology of simulated neutral hydrogen (HI) density fields using Shapefinders at various stages of reionization from the excursion set approach. Accompanying the Shapefinders, we also employ the 'largest cluster statistic' (LCS) to understand the percolation process. Percolation curves demonstrate that the non-Gaussianity in the HI field increases as reionization progresses. The large clusters in both the HI overdense and underdense excursion sets possess similar values of "thickness" ($T$), as well as "breadth" ($B$), but their third Shapefinder - "length" ($L$) - becomes almost proportional to their volume. The large clusters in both HI overdense and underdense segments are overwhelmingly filamentary. The 'cross-section' of a filamentary cluster can be estimated using the product of the first two Shapefinders, $T \times B$. Hence the cross sections of the large clusters at the onset of percolation do not vary much with volume and their sizes only differ in terms of their lengths. This feature appears more vividly in HI overdense regions than in underdense regions and is more pronounced at lower redshifts which correspond to an advanced stage of reionization.


title: Defending Elections Against Malicious Spread of Misinformation

abstract: The integrity of democratic elections depends on voters' access to accurate information. However, modern media environments, which are dominated by social media, provide malicious actors with unprecedented ability to manipulate elections via misinformation, such as fake news. We study a zero-sum game between an attacker, who attempts to subvert an election by propagating a fake new story or other misinformation over a set of advertising channels, and a defender who attempts to limit the attacker's impact. Computing an equilibrium in this game is challenging as even the pure strategy sets of players are exponential. Nevertheless, we give provable polynomial-time approximation algorithms for computing the defender's minimax optimal strategy across a range of settings, encompassing different population structures as well as models of the information available to each player. Experimental results confirm that our algorithms provide near-optimal defender strategies and showcase variations in the difficulty of defending elections depending on the resources and knowledge available to the defender.


title: Deep Compressive Autoencoder for Action Potential Compression in Large-Scale Neural Recording

abstract: Understanding the coordinated activity underlying brain computations requires large-scale, simultaneous recordings from distributed neuronal structures at a cellular-level resolution. One major hurdle to design high-bandwidth, high-precision, large-scale neural interfaces lies in the formidable data streams that are generated by the recorder chip and need to be online transferred to a remote computer. The data rates can require hundreds to thousands of I/O pads on the recorder chip and power consumption on the order of Watts for data streaming alone. We developed a deep learning-based compression model to reduce the data rate of multichannel action potentials. The proposed model is built upon a deep compressive autoencoder (CAE) with discrete latent embeddings. The encoder is equipped with residual transformations to extract representative features from spikes, which are mapped into the latent embedding space and updated via vector quantization (VQ). The decoder network reconstructs spike waveforms from the quantized latent embeddings. Experimental results show that the proposed model consistently outperforms conventional methods by achieving much higher compression ratios (20-500x) and better or comparable reconstruction accuracies. Testing results also indicate that CAE is robust against a diverse range of imperfections, such as waveform variation and spike misalignment, and has minor influence on spike sorting accuracy. Furthermore, we have estimated the hardware cost and real-time performance of CAE and shown that it could support thousands of recording channels simultaneously without excessive power/heat dissipation. The proposed model can reduce the required data transmission bandwidth in large-scale recording experiments and maintain good signal qualities. The code of this work has been made available at https://github.com/tong-wu-umn/spike-compression-autoencoder


title: Postponing the orthogonality catastrophe: efficient state preparation for electronic structure simulations on quantum devices

abstract: Despite significant work on resource estimation for quantum simulation of electronic systems, the challenge of preparing states with sufficient ground state support has so far been largely neglected. In this work we investigate this issue in several systems of interest, including organic molecules, transition metal complexes, the uniform electron gas, Hubbard models, and quantum impurity models arising from embedding formalisms such as dynamical mean-field theory. Our approach uses a state-of-the-art classical technique for high-fidelity ground state approximation. We find that easy-to-prepare single Slater determinants such as the Hartree-Fock state often have surprisingly robust support on the ground state for many applications of interest. For the most difficult systems, single-determinant reference states may be insufficient, but low-complexity reference states may suffice. For this we introduce a method for preparation of multi-determinant states on quantum computers.


title: Extending Neural Generative Conversational Model using External Knowledge Sources

abstract: The use of connectionist approaches in conversational agents has been progressing rapidly due to the availability of large corpora. However current generative dialogue models often lack coherence and are content poor. This work proposes an architecture to incorporate unstructured knowledge sources to enhance the next utterance prediction in chit-chat type of generative dialogue models. We focus on Sequence-to-Sequence (Seq2Seq) conversational agents trained with the Reddit News dataset, and consider incorporating external knowledge from Wikipedia summaries as well as from the NELL knowledge base. Our experiments show faster training time and improved perplexity when leveraging external knowledge.


title: Robustness of Quantum-Enhanced Adaptive Phase Estimation

abstract: As all physical adaptive quantum-enhanced metrology schemes operate under noisy conditions with only partially understood noise characteristics, so a practical control policy must be robust even for unknown noise. We aim to devise a test to evaluate the robustness of AQEM policies and assess the resource used by the policies. The robustness test is performed on QEAPE by simulating the scheme under four phase-noise models corresponding to normal-distribution noise, random-telegraph noise, skew-normal-distribution noise, and log-normal-distribution noise. Control policies are devised either by an evolutionary algorithm under the same noisy conditions, albeit ignorant of its properties, or a Bayesian-based feedback method that assumes no noise. Our robustness test and resource comparison method can be used to determining the efficacy and selecting a suitable policy.


title: Forward trijet production in proton-nucleus collisions

abstract: Using the formalism of the light-cone wave function in perturbative QCD together with the hybrid factorization, we compute the cross-section for three particle production at forward rapidities in proton-nucleus collisions. We focus on the quark channel, in which the three produced partons -- a quark accompanied by a gluon pair, or two quarks plus one antiquark -- are all generated via two successive splittings starting with a quark that was originally collinear with the proton. The three partons are put on-shell by their scattering off the nuclear target, described as a Lorentz-contracted "shockwave". The three-parton component of the quark light-cone wave function that we compute on this occasion is also an ingredient for other interesting calculations, like the next-to-leading order correction to the cross-section for the production of a pair of jets.


title: Gradient descent in higher codimension

abstract: We consider the behavior of gradient flow and of discrete and noisy gradient descent. It is commonly noted that the addition of noise to the process of discrete gradient descent can affect the trajectory of gradient descent. In previous work, we observed such effects. There, we considered the case where the minima had codimension 1. In this note, we do some computer experiments and observe the behavior of noisy gradient descent in the more complex setting of minima of higher codimension.


title: A New Secure Network Architecture to Increase Security Among Virtual Machines in Cloud Computing

abstract: Cloud computing is a new model of computing which provides scalability, flexibility and on-demand service. Virtualization is one of the main components of the cloud, but unfortunately, this technology suffers from many security vulnerabilities. The main purpose of this paper is to present a new secure architecture of Virtual Network machines in order to increase security among virtual machines in a virtualized environment (Xen as a case study). First, we expose the different network modes based on Xen Hypervisor, and then we analyze vulnerabilities and security issues within this kind of environment. Finally, we present in details new secure architecture and demonstrate how it can face the main security network attacks.


title: Entanglement Entropy from Numerical Holography

abstract: In this thesis I present numerical studies of entanglement entropy and the quantum null energy condition in strongly coupled far-from-equilibrium quantum states using holography. I give a careful introduction into the numerical methods used and provide the computer codes to compute entanglement entropy and the quantum null energy condition. These methods are then applied to systems of various degrees of complexity, including homogeneous and isotropic far- from-equilibrium quantum quenches dual to Vaidya spacetimes, to homogeneous and anisotropic finite temperature states dual to anisotropic black brane geometries, and to inhomogeneous and anisotropic states of colliding walls of energy dual to gravitational shock wave collisions in Anti-de Sitter space. For all these scenarios I compute the fully non-linear dynamics of the dual geometry, which requires to numerically solve five-dimensional Einstein equations with negative cosmological constant and asymptotic Anti-de Sitter boundary conditions. The numerical solutions for the geometries allow to extract the time evolution of the holographic energy momentum tensor and provides the background for computing two-point functions, entanglement entropy and the quantum null energy condition. From the anisotropic system one learns that the near-equilibrium dynamic of entanglement entropy has an accurate description in terms of quasinormal modes. In the shock wave system I identify characteristic features of entanglement entropy that allow to discriminate between thick and narrow shocks. All my numerical studies confirm the quantum null energy condition, also the shock wave system, which itself can violate the classical null energy condition for sufficiently narrow shocks. My results also show that the quantum null energy condition can saturate in the far-from-equilibrium regime.


title: Fairness in Online Social Network Timelines: Measurements, Models and Mechanism Design

abstract: Facebook News Feed personalization algorithm has a significant impact, on a daily basis, on the lifestyle, mood and opinion of millions of Internet users. Nonetheless, the behavior of such algorithm lacks transparency, motivating measurements, modeling and analysis in order to understand and improve its properties. In this paper, we propose a reproducible methodology encompassing measurements, an analytical model and a fairness-based News Feed design. The model leverages the versatility and analytical tractability of time-to-live (TTL) counters to capture the visibility and occupancy of publishers over a News Feed. Measurements are used to parameterize and to validate the expressive power of the proposed model. Then, we conduct a what-if analysis to assess the visibility and occupancy bias incurred by users against a baseline derived from the model. Our results indicate that a significant bias exists and it is more prominent at the top position of the News Feed. In addition, we find that the bias is non-negligible even for users that are deliberately set as neutral with respect to their political views, motivating the proposal of a novel and more transparent fairness-based News Feed design.


title: Notes on squeezed states in x-space representation

abstract: In these notes, we discuss squeezed states using the elementary quantum language based on one-dimensional Schr\"odinger equation. No operators are used. The language of quantum optics is mentioned only for a hint to solve a differential equation. Sections II and III (squeezed vacuum and pure squeezed states) can be used in a standard undergraduate course on quantum mechanics after discussion of a harmonic oscillator. Section IV discusses density matrix of mixed squeezed states in $x$-representation. These notes present explicit (therefore somewhat lengthy) step-by-step derivations. These notes are not intended for publication as a journal paper.


title: Continuous and discrete damping reduction for systems with quadratic interaction

abstract: We study the connection between Lagrangian and Hamiltonian descriptions of closed/open dynamics, for a collection of particles with quadratic interaction (closed system) and a sub-collection of particles with linear damping (open system). We consider both continuous and discrete versions of mechanics. We define the Damping Reduction as the mapping from the equations of motion of the closed system to those of the open one. As variational instruments for the obtention of these equations we use the Hamilton's principle (closed dynamics) and Lagrange-d'Alembert principle (open dynamics). We establish the commutativity of the branches Legendre transform + Damping Reduction and Damping Reduction+Legendre transform, where the Legendre transform is the usual mapping between Lagrangian and Hamiltonian mechanics. At a discrete level, this commutativity provides interesting insight about the resulting integrators. More concretely, Discrete Damping Reduction yields particular numerical schemes for linearly damped systems which are not symplectic anymore, but preserve some of the features of their symplectic counterparts from which they proceed (for instance the semi-implicitness in some cases). The theoretical results are illustrated with the examples of the heat bath and transmission lines. In the latter case some simulations are displayed, showing a better performance of the integrators with variational origin.


title: Charged Compact Stars in $f(\mathcal{G})$ Gravity

abstract: This work is devoted to investigate some of the interior configuration of static anisotropic spherical stellar charged structures in the regime of $f(\mathcal{G})$ gravity, where $\mathcal{G}$ is the Gauss Bonnet invariant. The structure of particular charged stars is analyzed with the help of solution obtained by Krori and Barua under different viable models in $f(\mathcal{G})$ gravity theory. The behavior of some physical aspects is investigated with the help of plots and the viability of our modeling is analyzed through different energy conditions. We have also studied some behavior of these realistic charged compact stars and discuss some aspects like density variation, evolution of stresses, different forces, stability of these stars, measure of anisotropy, equation of state parameters and the distribution of charges.


title: Diagnosing numerical Cherenkov instabilities in relativistic plasma simulations based on general meshes

abstract: Numerical Cherenkov radiation (NCR) or instability is a detrimental effect frequently found in electromagnetic particle-in-cell (EM-PIC) simulations involving relativistic plasma beams. NCR is caused by spurious coupling between electromagnetic-field modes and multiple beam resonances. This coupling may result from the slow down of poorly-resolved waves due to numerical (grid) dispersion and from aliasing mechanisms. NCR has been studied in the past for finite-difference-based EM-PIC algorithms on regular (structured) meshes with rectangular elements. In this work, we extend the analysis of NCR to finite-element-based EM-PIC algorithms implemented on unstructured meshes. The influence of different mesh element shapes and mesh layouts on NCR is studied. Analytic predictions are compared against results from finite-element-based EM-PIC simulations of relativistic plasma beams on various mesh types.


title: Riemann-Finsler Geometry and Lorentz-Violating Scalar Fields

abstract: The correspondence between Riemann-Finsler geometries and effective field theories with spin-independent Lorentz violation is explored. We obtain the general quadratic action for effective scalar field theories in any spacetime dimension with Lorentz-violating operators of arbitrary mass dimension. Classical relativistic point-particle lagrangians are derived that reproduce the momentum-velocity and dispersion relations of quantum wave packets. The correspondence to Finsler structures is established, and some properties of the resulting Riemann-Finsler spaces are investigated. The results provide support for open conjectures about Riemann-Finsler geometries associated with Lorentz-violating field theories.


title: Opacity, variability and kinematics of AGN jets

abstract: Synchrotron self-absorption in active galactic nuclei (AGN) jets manifests itself as a time delay between flares observed at high and low radio frequencies. It is also responsible for the observing frequency dependent change in size and position of the apparent base of the jet, aka the core shift effect, detected with very long baseline interferometry (VLBI). We measure the time delays and the core shifts in 11 radio-loud AGN to estimate the speed of their jets without relying on multi-epoch VLBI kinematics analysis. The 15$-$8 GHz total flux density time lags are obtained using Gaussian process regression, the core shift values are measured using VLBI observations and adopted from the literature. A strong correlation is found between the apparent core shift and the observed time delay. Our estimate of the jet speed is higher than the apparent speed of the fastest VLBI components by the median coefficient of 1.4. The coefficient ranges for individual sources from 0.5 to 20. We derive Doppler factors, Lorentz factors and viewing angles of the jets, as well as the corresponding de-projected distance from the jet base to the core. The results support evidence for acceleration of the jets with bulk motion Lorentz factor $\Gamma\propto R^{0.52\pm0.03}$ on de-projected scales $R$ of 0.5$-$500 parsecs.


title: The threshold effects in one-dimensional strongly-interacting systems out of equilibrium

abstract: In this work we investigate the phenomena associated with the new thresholds in the spectrum of excitations arising when different one-dimensional strongly interacting systems are voltage biased and weakly coupled by tunneling. We develop the perturbation theory with respect to tunneling and derive an asymptotic behavior of physical quantities close to threshold energies. We reproduce earlier results for the electron relaxation at the edge of an integer quantum Hall system and for the non-equilibrium Fermi edge singularity phenomenon. In contrast to the previous works, our analysis does not rely on the free-fermionic character of local tunneling, therefore we are able to extend our theory to wider class of systems, without well-defined electron excitations, such as spinless Luttinger liquids and chiral quantum Hall edge states at fractional filling factors.


title: Multi-Conjugate Adaptive Optics Simulator for the Thirty Meter Telescope: Design, Implementation, and Results

abstract: We present a multi-conjugate adaptive optics (MCAO) system simulator bench, HeNOS (Herzberg NFIRAOS Optical Simulator). HeNOS is developed to validate the performance of the MCAO system for the Thirty Meter Telescope, as well as to demonstrate techniques critical for future AO developments. In this paper, we focus on describing the derivations of parameters that scale the 30-m telescope AO system down to a bench experiment and explain how these parameters are practically implemented on an optical bench. While referring other papers for details of AO technique developments using HeNOS, we introduce the functionality of HeNOS, in particular, three different single-conjugate AO modes that HeNOS currently offers: a laser guide star AO with a Shack-Hartmann wavefront sensor, a natural guide star AO with a pyramid wavefront sensor, and a laser guide star AO with a sodium spot elongation on the Shack-Hartmann corrected by a truth wavefront sensing on a natural guide star. Laser tomography AO and ultimate MCAO are being prepared to be implemented in the near future.


title: The Origin of r-process Enhanced Metal-Poor Halo Stars In Now-Destroyed Ultra-Faint Dwarf Galaxies

abstract: The highly r-process enhanced (r-II) metal-poor halo stars we observe today could play a key role in understanding early ultra-faint dwarf galaxies, the smallest building blocks of the Milky Way. If a significant fraction of metal-poor r-II halo stars originated in the ultra-faint dwarf galaxies that merged to help form the Milky Way, observations of r-II stars could help us study these now-destroyed systems and probe the formation history of our Galaxy. To conduct our initial investigation into this possible connection, we use high-resolution cosmological simulations of Milky-Way-mass galaxies from the Caterpillar suite in combination with a simple, empirically motivated treatment of r-process enrichment. We determine the fraction of metal-poor halo stars that could have formed from highly r-process enhanced gas in now-destroyed low-mass ultra-faint dwarf galaxies, the simulated r-II fraction, and compare it to the "as observed" r-II fraction. We find that the simulated fraction, f_{r-II,sim} ~ 1-2%, can account for around half of the "as observed" fraction, f_{r-II,obs} ~ 2-4%. The "as observed" fraction likely overrepresents the fraction of r-II stars due to incomplete sampling, though, meaning f_{r-II,sim} likely accounts for more than half of the true f_{r-II,obs}. Further considering some parameter variations and scatter between individual simulations, the simulated fraction can account for around 20-80% of the "as observed" fraction.


title: Quantum Kibble-Zurek mechanism and critical dynamics on a programmable Rydberg simulator

abstract: Quantum phase transitions (QPTs) involve transformations between different states of matter that are driven by quantum fluctuations. These fluctuations play a dominant role in the quantum critical region surrounding the transition point, where the dynamics are governed by the universal properties associated with the QPT. While time-dependent phenomena associated with classical, thermally driven phase transitions have been extensively studied in systems ranging from the early universe to Bose Einstein Condensates, understanding critical real-time dynamics in isolated, non-equilibrium quantum systems is an outstanding challenge. Here, we use a Rydberg atom quantum simulator with programmable interactions to study the quantum critical dynamics associated with several distinct QPTs. By studying the growth of spatial correlations while crossing the QPT, we experimentally verify the quantum Kibble-Zurek mechanism (QKZM) for an Ising-type QPT, explore scaling universality, and observe corrections beyond QKZM predictions. This approach is subsequently used to measure the critical exponents associated with chiral clock models, providing new insights into exotic systems that have not been understood previously, and opening the door for precision studies of critical phenomena, simulations of lattice gauge theories and applications to quantum optimization.


title: Lampray: Multi-group long characteristics ray tracing for adaptive mesh radiation hydrodynamics

abstract: We present Lampray: a multi-group long characteristics ray tracing method for adaptive mesh radiation hydrodynamics in the Ramses code. It avoids diffusion, captures shadows, and treats colliding beams correctly, and therefore complements existing moment-based ray tracing in Ramses. Lampray includes different options for interpolation between ray and cell domain, and use either integral, Fourier, or an implicit method for hydrogen ionization to solve the radiative transfer. The opacity can either be tabulated or computed through a coupling to the general non-equilibrium astro-chemistry framework Krome. We use an H-He-C-O network with 36 species and 240 reactions to track the photo-chemistry in the interstellar medium across 6 and 10 orders of magnitude in temperature and density. Self-shielding prescriptions for H$_2$ and CO are used together with a new model for the diffuse interstellar UV-field. We also track the dust temperature, formation of H$_2$ on grains, and H$_2$O and CO ices in detail. Lampray is tested against standard benchmarks for molecular cloud and star formation simulations, including the formation of a Str\"omgren sphere, the expansion of an ionization front, the photo-evaporation of a dense clump, and the H-He-C-O chemistry in a static photo-dissociation front. Efficient parallelisation is achieved with a separate domain decomposition for rays where points along a ray reside in the same memory space, and data movement from cell- to ray-domain is done with a direct hash-table lookup algorithm. Point sources are treated without splitting rays, and therefore the method currently only scales to a few point sources, while diffuse radiation has excellent scaling.


title: Modeling the Connection Between Subhalos and Satellites in Milky Way-like Systems

abstract: We develop a comprehensive and flexible model for the connection between satellite galaxies and dark matter subhalos in dark matter-only zoom-in simulations of Milky Way (MW)--mass host halos. We systematically identify the physical and numerical uncertainties in the galaxy--halo connection and simulations underlying our method, including (i) the influence of host halo properties; (ii) the relationship between satellite luminosities and subhalo properties, including the effects of reionization; (iii) the relationship between satellite and subhalo locations; (iv) the relationship between satellite sizes and subhalo properties, including the effects of tidal stripping; (v) satellite and subhalo disruption due to baryonic effects; and (vi) artificial subhalo disruption and orphan satellites. To illustrate our approach, we fit this model to the luminosity distribution of both classical MW satellites and those discovered in the Sloan Digital Sky Survey by performing realistic mock observations that depend on the luminosity, size, and distance of our predicted satellites, and we infer the total satellite population that will be probed by upcoming surveys. We argue that galaxy size and surface brightness modeling will play a key role in interpreting current and future observations, as the expected number of observable satellites depends sensitively on their surface brightness distribution.


title: Arithmetic of (independent) sigma-fields on probability spaces

abstract: This note gathers what is known about, and provides some new results concerning the operations of intersection, of ``generated $\sigma$-field'', and of ``complementation'' for (independent) complete $\sigma$-fields on probability spaces.


title: Quasar Sightline and Galaxy Evolution (QSAGE) Survey - I. The Galaxy Environment of OVI Absorbers up to z=1.4 around PKS 0232-04

abstract: We present the first results from a study of OVI absorption around galaxies at $z<1.44$ using data from a near-infrared grism spectroscopic Hubble Space Telescope Large Program, the Quasar Sightline and Galaxy Evolution (QSAGE) survey. QSAGE is the first grism galaxy survey to focus on the circumgalactic medium at $z\sim1$, providing a blind survey of the galaxy population. Using the first of 12 fields, we provide details of the reduction methods, in particular the handling of the deep grism data which uses multiple position angles to minimise the effects of contamination from overlapping traces. The resulting galaxy sample is H$\alpha$ flux limited ($f({\rm H\alpha}) > 2\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$) at 0.68<z<1.44, corresponding to $\gtrsim0.2-0.8$ M$_\odot$ yr$^{-1}$. We combine the galaxy data with high-resolution STIS and COS spectroscopy of the background quasar to study OVI in the circumgalactic medium. At z>0.68, we find 5 OVI absorption systems along the line of sight with identified galaxies lying at impact parameters of $b\approx100-350$ kpc (proper), whilst we find a further 13 galaxies with no significant associated OVI absorption (i.e. $N({\rm OVI})<10^{13.5-14}$ cm$^{-2}$) in the same impact parameter and redshift range. We find a large scatter in the stellar mass and star-formation rates of the closest galaxies with associated OVI. Whilst one of the OVI absorber systems is found to be associated with a low mass galaxy group at $z\approx1.08$, we infer that the detected OVI absorbers typically lie in the proximity of dark matter halos of masses $10^{11.5} {\rm M_\odot}\lesssim M_{\rm halo}\lesssim10^{12} {\rm M_\odot}$.


title: Super-Strong Coupling NJL model in arbitrary space-time dimensions

abstract: We study chiral symmetry breaking in the Nambu--Jona-Lasinio model regularized in proper-time in arbitrary space-time dimensions through an iterative procedure by writing the gap equation in the form of a discrete dynamical system with the coupling constant as the control parameter. Expectedly, we obtain the critical coupling for chiral symmetry breaking when a nontrivial solution bifurcates away from the trivial one and becomes an attractor. By increasing further the value of the coupling constant, we observe a second bifurcation where the dynamical solution is no longer an attractor, and observation that holds true in all space-time dimensions. In the super-strong coupling regime, the system becomes chaotic.


title: On the Feasibility of Using a Laser Guide Star Adaptive Optics System in the Daytime

abstract: We investigate the use of ultra-narrow band interference filters to enable daytime use of sodium laser guide star adaptive optics systems. Filter performance is explored using theoretical and vendor supplied filter transmission profiles, a modeled daylight sky background, broadband measurements of the daytime sky brightness on Maunakea, and an assumed photon return from the sodium laser guide star and read noise for the wavefront sensor detector. The critical parameters are the bandpass of the filter, the out-of-band rejection, and the peak throughput at the wavelength of the laser guide star light. Importantly, a systematic trade between these parameters leads to potentially simple solutions enabling daytime observations. Finally, we simulated the Mid-Infrared Adaptive Optics (MIRAO) system planned for the Thirty Meter Telescope with an end-to-end simulation, folding in daytime sky counts. We find that MIRAO with five sodium laser guide stars, commercial off-the-shelf filters to suppress the sky background in the laser guide star wavefront sensors, and a near-infrared natural guide star ($K \le 13$) tip/tilt/focus wavefront sensor can attain daytime Strehl ratio values comparable to those at night.


title: Harvesting correlations from thermal and squeezed coherent states

abstract: We study the harvesting of entanglement and mutual information by Unruh-DeWitt particle detectors from thermal and squeezed coherent field states. We prove (for arbitrary spatial dimensions, switching profiles and detector smearings) that while the entanglement harvesting ability of detectors decreases monotonically with the field temperature $T$, harvested mutual information grows linearly with $T$. We also show that entanglement harvesting from a general squeezed coherent state is independent of the coherent amplitude, but depends strongly on the squeezing amplitude. Moreover, we find that highly squeezed states i) allow for detectors to harvest much more entanglement than from the vacuum, and ii) ensure that the entanglement harvested does not decay with their spatial separation. Finally we analyze the spatial inhomogeneity of squeezed states and its influence on harvesting, and investigate how much entanglement one can actually extract from squeezed states when the squeezing is bandlimited.


title: The Simplified Likelihood Framework

abstract: We discuss the simplified likelihood framework as a systematic approximation scheme for experimental likelihoods such as those originating from LHC experiments. We develop the simplified likelihood from the Central Limit Theorem keeping the next-to-leading term in the large $N$ expansion to correctly account for asymmetries. Moreover, we present an efficient method to compute the parameters of the simplified likelihood from Monte Carlo simulations. The approach is validated using a realistic LHC-like analysis, and the limits of the approximation are explored. Finally, we discuss how the simplified likelihood data can be conveniently released in the HepData error source format and automatically built from it, making this framework a convenient tool to transmit realistic experimental likelihoods to the community.


title: Isothermal current-driven insulator-to-metal transition in VO$_\text{2}$ through strong correlation effect

abstract: Electric current has been experimentally demonstrated to be able to drive the insulator-to-metal transition (IMT) in VO$_2$. The main mechanisms involved are believed to be the Joule heating effect and the strong electron-correlation effect. These effects are often entangled with each other in experiments, which complicates the understanding of the essential nature of the observations. We formulate a phase-field model to investigate theoretically in mesoscale the pure correlation effect brought by the current on the IMT in VO$_2$, i.e., the isothermal process under the current. We find that a current with a large density ($\sim 10^1$ nA/nm$^2$) induces a few-nanosecond ultrafast switch in VO$_2$, in agreement with the experiment. The temperature-current phase diagram is further calculated, which reveals that the current may induce the M2 phase at low temperatures. The current is also shown capable of driving domain walls to move. Our work may assist related experiments and provide guidance to the engineering of VO$_2$-based electric switching devices.


title: Efficient Structured Surrogate Loss and Regularization in Structured Prediction

abstract: In this dissertation, we focus on several important problems in structured prediction. In structured prediction, the label has a rich intrinsic substructure, and the loss varies with respect to the predicted label and the true label pair. Structured SVM is an extension of binary SVM to adapt to such structured tasks. In the first part of the dissertation, we study the surrogate losses and its efficient methods. To minimize the empirical risk, a surrogate loss which upper bounds the loss, is used as a proxy to minimize the actual loss. Since the objective function is written in terms of the surrogate loss, the choice of the surrogate loss is important, and the performance depends on it. Another issue regarding the surrogate loss is the efficiency of the argmax label inference for the surrogate loss. Efficient inference is necessary for the optimization since it is often the most time-consuming step. We present a new class of surrogate losses named bi-criteria surrogate loss, which is a generalization of the popular surrogate losses. We first investigate an efficient method for a slack rescaling formulation as a starting point utilizing decomposability of the model. Then, we extend the algorithm to the bi-criteria surrogate loss, which is very efficient and also shows performance improvements. In the second part of the dissertation, another important issue of regularization is studied. Specifically, we investigate a problem of regularization in hierarchical classification when a structural imbalance exists in the label structure. We present a method to normalize the structure, as well as a new norm, namely shared Frobenius norm. It is suitable for hierarchical classification that adapts to the data in addition to the label structure.


title: Non-Matrix Tactile Sensors: How Can Be Exploited Their Local Connectivity For Predicting Grasp Stability?

abstract: Tactile sensors supply useful information during the interaction with an object that can be used for assessing the stability of a grasp. Most of the previous works on this topic processed tactile readings as signals by calculating hand-picked features. Some of them have processed these readings as images calculating characteristics on matrix-like sensors. In this work, we explore how non-matrix sensors (sensors with taxels not arranged exactly in a matrix) can be processed as tactile images as well. In addition, we prove that they can be used for predicting grasp stability by training a Convolutional Neural Network (CNN) with them. We captured over 2500 real three-fingered grasps on 41 everyday objects to train a CNN that exploited the local connectivity inherent on the non-matrix tactile sensors, achieving 94.2% F1-score on predicting stability.


title: Packing colorings of subcubic outerplanar graphs

abstract: Given a graph $G$ and a nondecreasing sequence $S=(s_1,\ldots,s_k)$ of positive integers, the mapping $c:V(G)\longrightarrow \{1,\ldots,k\}$ is called an $S$-packing coloring of $G$ if for any two distinct vertices $x$ and $y$ in $c^{-1}(i)$, the distance between $x$ and $y$ is greater than $s_i$. The smallest integer $k$ such that there exists a $(1,2,\ldots,k)$-packing coloring of a graph $G$ is called the packing chromatic number of $G$, denoted $\chi_{\rho}(G)$. The question of boundedness of the packing chromatic number in the class of subcubic (planar) graphs was investigated in several earlier papers; recently it was established that the invariant is unbounded in the class of all subcubic graphs. In this paper, we prove that the packing chromatic number of any 2-connected bipartite subcubic outerplanar graph is bounded by $7$. Furthermore, we prove that every subcubic triangle-free outerplanar graph has a $(1,2,2,2)$-packing coloring, and that there exists a subcubic outerplanar graph with a triangle that does not admit a $(1,2,2,2)$-packing coloring. In addition, there exists a subcubic triangle-free outerplanar graph that does not admit a $(1,2,2,3)$-packing coloring. A similar dichotomy is shown for bipartite outerplanar graphs: every such graph admits an $S$-packing coloring for $S=(1,3,\ldots,3)$, where $3$ appears $\Delta$ times ($\Delta$ being the maximum degree of vertices), and this property does not hold if one of the integers $3$ is replaced by $4$ in the sequence $S$.


title: On Picard Type Theorems and Entire Solutions of Differential Equations

abstract: We give a connection between the Picard type theorem of Polya-Saxer-Milliox and characterization of entire solutions of a differential equation and then their higher dimensional extensions, which leads further results on both (ordinary and partial) differential equations and Picard type theorems.


title: Quantifying and controlling prethermal nonergodicity in interacting Floquet matter

abstract: The use of periodic driving for synthesizing many-body quantum states depends crucially on the existence of a prethermal regime, which exhibits drive-tunable properties while forestalling the effects of heating. This motivates the search for direct experimental probes of the underlying localized nonergodic nature of the wave function in this metastable regime. We report experiments on a many-body Floquet system consisting of atoms in an optical lattice subjected to ultrastrong sign-changing amplitude modulation. Using a double-quench protocol we measure an inverse participation ratio quantifying the degree of prethermal localization as a function of tunable drive parameters and interactions. We obtain a complete prethermal map of the drive-dependent properties of Floquet matter spanning four square decades of parameter space. Following the full time evolution, we observe sequential formation of two prethermal plateaux, interaction-driven ergodicity, and strongly frequency-dependent dynamics of long-time thermalization. The quantitative characterization of the prethermal Floquet matter realized in these experiments, along with the demonstration of control of its properties by variation of drive parameters and interactions, opens a new frontier for probing far-from-equilibrium quantum statistical mechanics and new possibilities for dynamical quantum engineering.

categories: quant-ph cond-mat.quant-gas cond-mat.stat-mech physics.atom-ph


title: Rigid Body Dynamic Simulation with Multiple Convex Contact Patches

abstract: We present a principled method for dynamic simulation of rigid bodies in intermittent contact with each other where the contact is assumed to be a non-convex contact patch that can be modeled as a union of convex patches. The prevalent assumption in simulating rigid bodies undergoing intermittent contact with each other is that the contact is a point contact. In recent work, we introduced an approach to simulate contacting rigid bodies with convex contact patches (line and surface contact). In this paper, for non-convex contact patches modeled as a union of convex patches, we formulate a discrete-time mixed complementarity problem where we solve the contact detection and integration of the equations of motion simultaneously. Thus, our method is a geometrically-implicit method and we prove that in our formulation, there is no artificial penetration between the contacting rigid bodies. We solve for the equivalent contact point (ECP) and contact impulse of each contact patch simultaneously along with the state, i.e., configuration and velocity of the objects. We provide empirical evidence to show that if the number of contact patches between two objects is less than or equal to three, the state evolution of the bodies is unique, although the contact impulses and ECP may not be unique. We also present simulation results showing that our method can seamlessly capture transition between different contact modes like non-convex patch to point (or line contact) and vice-versa during simulation.


title: Playing With Danger: A Taxonomy and Evaluation of Threats to Smart Toys

abstract: Smart toys have captured an increasing share of the toy market, and are growing ubiquitous in households with children. Smart toys are a subset of Internet of Things (IoT) devices, containing sensors, actuators, and/or artificial intelligence capabilities. They frequently have internet connectivity, directly or indirectly through companion apps, and collect information about their users and environments. Recent studies have found security flaws in many smart toys that have led to serious privacy leaks, or allowed tracking a child's physical location. Some well-publicized discoveries of this nature have prompted actions from governments around the world to ban some of these toys. Compared to other IoT devices, smart toys pose unique risks because of their easily-vulnerable user base, and our work is intended to define these risks and assess a subset of toys against them. We provide a classification of threats specific to smart toys in order to unite and complement existing adhoc analyses, and help comprehensive evaluation of other smart toys. Our threat classification framework addresses the potential security and privacy flaws that can lead to leakage of private information or allow an adversary to control the toy to lure, harm, or distress a child. Using this framework, we perform a thorough experimental analysis of eleven smart toys and their companion apps. Our systematic analysis has uncovered that several current toys still expose children to multiple threats for attackers with physical, nearby, or remote access to the toy.


title: Identification of multi-scale hierarchical brain functional networks using deep matrix factorization

abstract: We present a deep semi-nonnegative matrix factorization method for identifying subject-specific functional networks (FNs) at multiple spatial scales with a hierarchical organization from resting state fMRI data. Our method is built upon a deep semi-nonnegative matrix factorization framework to jointly detect the FNs at multiple scales with a hierarchical organization, enhanced by group sparsity regularization that helps identify subject-specific FNs without loss of inter-subject comparability. The proposed method has been validated for predicting subject-specific functional activations based on functional connectivity measures of the hierarchical multi-scale FNs of the same subjects. Experimental results have demonstrated that our method could obtain subject-specific multi-scale hierarchical FNs and their functional connectivity measures across different scales could better predict subject-specific functional activations than those obtained by alternative techniques.


title: Search for leptoquarks coupled to third-generation quarks in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: Three of the most significant measured deviations from standard model predictions, the enhanced decay rate for B $\to$ D$^{(*)}\tau\nu$, hints of lepton universality violation in B $\to$ K$^{(*)}\ell\ell$ decays, and the anomalous magnetic moment of the muon, can be explained by the existence of leptoquarks (LQs) with large couplings to third-generation quarks and masses at the TeV scale. The existence of these states can be probed at the LHC in high energy proton-proton collisions. A novel search is presented for pair production of LQs coupled to a top quark and a muon using data at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$, recorded by the CMS experiment. No deviation from the standard model prediction has been observed and scalar LQs decaying exclusively into t$\mu$ are excluded up to masses of 1420 GeV. The results of this search are combined with those from previous searches for LQ decays into t$\tau$ and b$\nu$, which excluded scalar LQs below masses of 900 and 1080 GeV. Vector LQs are excluded up to masses of 1190 GeV for all possible combinations of branching fractions to t$\mu$, t$\tau$ and b$\nu$. With this analysis, all relevant couplings of LQs with an electric charge of -1/3 to third-generation quarks are probed for the first time.


title: Retrieval method for the bianisotropic polarizability tensor of Willis acoustic scatterers

abstract: Acoustic materials displaying coupling between pressure and momentum are known as Willis materials. The simplest Willis materials are comprised of sub-wavelength scatterers that couple monopoles to dipoles and {\it vice versa}, with the interaction defined by a polarizability tensor. We propose a method for retrieving the polarizability tensor for sub-wavelength Willis acoustic scatterers using a finite set of scattering amplitudes. We relate the polarizability tensor to standard T-matrix and S-matrix scattering formalisms. This leads to an explicit method for retrieving the components of the polarizability tensor in terms of a small set of scattered pressure data in the near- or far-field. Numerical examples demonstrate the retrieval method for one and two dimensional configurations.


title: Identification of temporal transition of functional states using recurrent neural networks from functional MRI

abstract: Dynamic functional connectivity analysis provides valuable information for understanding brain functional activity underlying different cognitive processes. Besides sliding window based approaches, a variety of methods have been developed to automatically split the entire functional MRI scan into segments by detecting change points of functional signals to facilitate better characterization of temporally dynamic functional connectivity patterns. However, these methods are based on certain assumptions for the functional signals, such as Gaussian distribution, which are not necessarily suitable for the fMRI data. In this study, we develop a deep learning based framework for adaptively detecting temporally dynamic functional state transitions in a data-driven way without any explicit modeling assumptions, by leveraging recent advances in recurrent neural networks (RNNs) for sequence modeling. Particularly, we solve this problem in an anomaly detection framework with an assumption that the functional profile of one single time point could be reliably predicted based on its preceding profiles within stable functional state, while large prediction errors would occur around change points of functional states. We evaluate the proposed method using both task and resting-state fMRI data obtained from the human connectome project and experimental results have demonstrated that the proposed change point detection method could effectively identify change points between different task events and split the resting-state fMRI into segments with distinct functional connectivity patterns.


title: Brain decoding from functional MRI using long short-term memory recurrent neural networks

abstract: Decoding brain functional states underlying different cognitive processes using multivariate pattern recognition techniques has attracted increasing interests in brain imaging studies. Promising performance has been achieved using brain functional connectivity or brain activation signatures for a variety of brain decoding tasks. However, most of existing studies have built decoding models upon features extracted from imaging data at individual time points or temporal windows with a fixed interval, which might not be optimal across different cognitive processes due to varying temporal durations and dependency of different cognitive processes. In this study, we develop a deep learning based framework for brain decoding by leveraging recent advances in sequence modeling using long short-term memory (LSTM) recurrent neural networks (RNNs). Particularly, functional profiles extracted from task functional imaging data based on their corresponding subject-specific intrinsic functional networks are used as features to build brain decoding models, and LSTM RNNs are adopted to learn decoding mappings between functional profiles and brain states. We evaluate the proposed method using task fMRI data from the HCP dataset, and experimental results have demonstrated that the proposed method could effectively distinguish brain states under different task events and obtain higher accuracy than conventional decoding models.


title: On the quantization of the Hall conductivity in the Harper-Hofstadter model

abstract: We study the robustness of the quantization of the Hall conductivity in the Harper-Hofstadter model towards the details of the protocol with which a longitudinal uniform driving force $F_x(t)$ is turned on. In the vector potential gauge, through Peierls substitution, this involves the switching-on of complex time-dependent hopping amplitudes $\mathrm{e}^{-\frac{i}{\hbar}\mathcal{A}_x(t)}$ in the $\hat{\mathbf{x}}$-direction such that $\partial_t \mathcal{A}_x(t)=F_x(t)$. The switching-on can be sudden, $F_x(t)=\theta(t) F$, where $F$ is the steady driving force, or more generally smooth $F_x(t)=f(t/t_{0}) F$, where $f(t/t_{0})$ is such that $f(0)=0$ and $f(1)=1$. We investigate how the time-averaged (steady-state) particle current density $j_y$ in the $\hat{\mathbf{y}}$-direction deviates from the quantized value $j_y \, h/F = n$ due to the finite value of $F$ and the details of the switching-on protocol. Exploiting the time-periodicity of the Hamiltonian $\hat{H}(t)$, we use Floquet techniques to study this problem. In this picture the (Kubo) linear response $F\to 0$ regime corresponds to the adiabatic limit for $\hat{H}(t)$. In the case of a sudden quench $j_y \, h/F$ shows $F^2$ corrections to the perfectly quantized limit. When the switching-on is smooth, the result depends on the switch-on time $t_{0}$: for a fixed $t_{0}$ we observe a crossover force $F^*$ between a quadratic regime for $F<F^*$ and a {\em non-analytic} exponential $\mathrm{e}^{-\gamma/|F|}$ for $F>F^*$. The crossover $F^*$ decreases as $t_{0}$ increases, eventually recovering the topological robustness. These effects are in principle amenable to experimental tests in optical lattice cold atomic systems with synthetic gauge fields.


title: Some exit times estimates for Super-Brownian motion and Fleming-Viot Process

abstract: Estimates for exit time from an interval of length 2r before a prescribed time T are derived for solutions of a class of stochastic partial differential equations used to characterize two population models: super-Brownian motion and Fleming-Viot Process. These types of estimates are then derived for the two population models. The corresponding large deviation results are also applied for the acquired bounds.


title: Duality of maximal-torsion and close-packing in double helices

abstract: It is suggested that there is an interesting transition in double helices. The transition is between a helix, which has maximal torsion for the constituting helical strands, and a double helix which is nearly optimally packed illustrated by corresponding dual points on the rope curve. The possible physical relevance for denaturation and renaturation of DNA is discussed.


title: Hyperuniformity of generalized random organization models

abstract: Studies of random organization models of monodisperse spherical particles have shown that a hyperuniform state is achievable when the system goes through an absorbing phase transition to a critical state. Here we investigate to what extent hyperuniformity is preserved when the model is generalized to particles with a size distribution and/or nonspherical shapes. We begin by examining binary disks in two dimensions and demonstrate that their critical states are hyperuniform as two-phase media, but not hyperuniform nor multihyperuniform as point patterns formed by the particle centroids. We further confirm the generality of our findings by studying particles with a continuous size distribution. Finally, to study the effect of rotational degrees of freedom, we extend our model to noncircular particles, namely, hard rectangles with various aspect ratios, including the hard-needle limit. Although these systems exhibit only short-range orientational order, hyperuniformity is still preserved. Our analysis reveals that the redistribution of the "mass" of the particles rather than the particle centroids is central to this dynamical process. The consideration of the "active volume fraction" of generalized random organization models may help to resolve which universality class they belong to and hence may lead to a deeper theoretical understanding of absorbing-state models. Our results suggest that general particle systems subject to random organization can be a robust way to fabricate a wide class of hyperuniform states of matter by tuning the structures via different particle-size and -shape distributions. This in turn potentially enables the creation of multifunctional hyperuniform materials with desirable optical, transport, and mechanical properties.

categories: cond-mat.soft cond-mat.dis-nn cond-mat.mtrl-sci cond-mat.stat-mech


title: Metric Graph Approximations of Geodesic Spaces

abstract: We study the question of approximating a compact geodesic metric space by metric graphs satisfying a uniform upper bound on their first Betti number. We prove that, up to a suitable multiplicative constant, Reeb graphs of distance functions to a point provide optimal approximation in the Gromov-Hausdsorff sense.


title: Multifidelity Dimension Reduction via Active Subspaces

abstract: We propose a multifidelity dimension reduction method to identify a low-dimensional structure present in many engineering models. The structure of interest arises when functions vary primarily on a low-dimensional subspace of the high-dimensional input space, while varying little along the complementary directions. Our approach builds on the gradient-based methodology of active subspaces, and exploits models of different fidelities to reduce the cost of performing dimension reduction through the computation of the active subspace matrix. We provide a non-asymptotic analysis of the number of gradient evaluations sufficient to achieve a prescribed error in the active subspace matrix, both in expectation and with high probability. We show that the sample complexity depends on a notion of intrinsic dimension of the problem, which can be much smaller than the dimension of the input space. We illustrate the benefits of such a multifidelity dimension reduction approach using numerical experiments with input spaces of up to three thousand dimensions.


title: The imaginary Toda field theory

abstract: We consider the two-dimensional $\mathfrak{sl}_n$ quantum Toda field theory with an imaginary background charge. This conformal field theory has a higher spin symmetry ($W_n$ algebra), a central charge $c \leq n-1$ and a continuous spectrum. Using the conformal bootstrap, we compute structure constants involving two arbitrary scalar fields and a semi-degenerate field of Wyllard type. The solution obtained is not the analytic continuation of the usual Toda three-point function. Non-scalar primary fields and their three-point functions are also discussed. Non-scalar primary fields are classified by conjugacy classes of the permutation group $\mathfrak{S}_n$, and their structure constants are computed explicitly, up to an overall factor.


title: On prime order automorphisms of generalized quadrangles

abstract: In this paper, we study prime order automorphisms of generalized quadrangles. We show that, if $\mathcal{Q}$ is a thick generalized quadrangle of order $(s,t)$, where $s > t$ and $s+1$ is prime, and $\mathcal{Q}$ has an automorphism of order $s+1$, then \[ s \left\lceil \left\lceil \frac{t^2}{s+1}\right\rceil\left(\frac{s+1}{t} \right) \right\rceil \le t(s+t),\] with a similar inequality holding in the dual case when $t > s$, $t+1$ is prime, and $\mathcal{Q}$ is a thick generalized quadrangle of order $(s,t)$ with an automorphism of order $t+1$. In particular, if $s+1$ is prime and if there exists a natural number $n$ such that \[ \frac{t^2}{n+1} + t \le s + 1 < \frac{t^2}{n},\] then a thick generalized quadrangle $\mathcal{Q}$ cannot have an automorphism of order $s+1$, and hence the automorphism group of $\mathcal{Q}$ cannot be transitive on points. These results apply to numerous potential orders for which it is still unknown whether or not generalized quadrangles exist, showing that any examples would necessarily be somewhat asymmetric. Finally, we are able to use the theory we have built up about prime order automorphisms of generalized quadrangles to show that the automorphism group of a potential generalized quadrangle of order $(4,12)$ must necessarily be intransitive on both points and lines.


title: Limit distributions of expanding translates of shrinking submanifolds and non-improvability of Dirichlet's approximation theorem

abstract: On the space $\mathcal{L}_{n+1}$ of unimodular lattices in $\mathbb{R}^{n+1}$, we consider the standard action of $a(t)=\mathrm{diag}(t^n,t^{-1},\ldots,t^{-1})\in \mathrm{SL}(n+1,\mathbb{R})$ for $t>1$. Let $M$ be a nondegenerate submanifold of an expanding horospherical leaf in $\mathcal{L}_{n+1}$. We prove that for all $x\in M\setminus E$ and $t>1$, if $\mu_{x,t}$ denotes the normalized Lebesgue measure on the ball of radius $t^{-1}$ around $x$ in $M$, then the translated measure $a(t)\mu_{x,t}$ get equidistributed $\mathcal{L}_{n+1}$ as $t\to\infty$, where $E$ is a union of countably many lower dimensional submanifolds of $M$. In particular, if $\mu$ is an absolutely continuous probability measure on $M$, then $a(t)\mu$ gets equidistributed in $\mathcal{L}_{n+1}$ as $t\to\infty$. This result implies the non-improvability of Dirichlet's Diophantine approximation theorem for almost every point on a $C^{n+1}$-submanifold of $\mathbb{R}^n$ satisfying a non-degeneracy condition, answering a question arising from the work of Davenport and Schmidt (1969).


title: Models Matter, So Does Training: An Empirical Study of CNNs for Optical Flow Estimation

abstract: We investigate two crucial and closely related aspects of CNNs for optical flow estimation: models and training. First, we design a compact but effective CNN model, called PWC-Net, according to simple and well-established principles: pyramidal processing, warping, and cost volume processing. PWC-Net is 17 times smaller in size, 2 times faster in inference, and 11\% more accurate on Sintel final than the recent FlowNet2 model. It is the winning entry in the optical flow competition of the robust vision challenge. Next, we experimentally analyze the sources of our performance gains. In particular, we use the same training procedure of PWC-Net to retrain FlowNetC, a sub-network of FlowNet2. The retrained FlowNetC is 56\% more accurate on Sintel final than the previously trained one and even 5\% more accurate than the FlowNet2 model. We further improve the training procedure and increase the accuracy of PWC-Net on Sintel by 10\% and on KITTI 2012 and 2015 by 20\%. Our newly trained model parameters and training protocols will be available on https://github.com/NVlabs/PWC-Net


title: Entropic optimal transport is maximum-likelihood deconvolution

abstract: We give a statistical interpretation of entropic optimal transport by showing that performing maximum-likelihood estimation for Gaussian deconvolution corresponds to calculating a projection with respect to the entropic optimal transport distance. This structural result gives theoretical support for the wide adoption of these tools in the machine learning community.


title: Rigidity theorems for circle domains

abstract: A circle domain $\Omega$ in the Riemann sphere is conformally rigid if every conformal map from $\Omega$ onto another circle domain is the restriction of a M\"{o}bius transformation. We show that circle domains satisfying a certain quasihyperbolic condition, which was considered by Jones and Smirnov, are conformally rigid. In particular, H\"{o}lder circle domains and John circle domains are all conformally rigid. This provides new evidence for a conjecture of He and Schramm relating rigidity and conformal removability.


title: Multiple Workflows Scheduling in Multi-tenant Distributed Systems: A Taxonomy and Future Directions

abstract: The workflow is a general notion representing the automated processes along with the flow of data. The automation ensures the processes being executed in the order. Therefore, this feature attracts users from various background to build the workflow. However, the computational requirements are enormous and investing for a dedicated infrastructure for these workflows is not always feasible. To cater to the broader needs, multi-tenant platforms for executing workflows were began to be built. In this paper, we identify the problems and challenges in the multiple workflows scheduling that adhere to the platforms. We present a detailed taxonomy from the existing solutions on scheduling and resource provisioning aspects followed by the survey of relevant works in this area. We open up the problems and challenges to shove up the research on multiple workflows scheduling in multi-tenant distributed systems.


title: Bifurcations in synergistic epidemics on random regular graphs

abstract: The role of cooperative effects (i.e. synergy) in transmission of infection is investigated analytically and numerically for epidemics following the rules of Susceptible-Infected-Susceptible (SIS) model defined on random regular graphs. Non-linear dynamics are shown to lead to bifurcation diagrams for such spreading phenomena exhibiting three distinct regimes: non-active, active and bi-stable. The dependence of bifurcation loci on node degree is studied and interesting effects are found that contrast with the behaviour expected for non-synergistic epidemics.


title: Events Beyond ACE: Curated Training for Events

abstract: We explore a human-driven approach to annotation, curated training (CT), in which annotation is framed as teaching the system by using interactive search to identify informative snippets of text to annotate, unlike traditional approaches which either annotate preselected text or use active learning. A trained annotator performed 80 hours of CT for the thirty event types of the NIST TAC KBP Event Argument Extraction evaluation. Combining this annotation with ACE results in a 6% reduction in error and the learning curve of CT plateaus more slowly than for full-document annotation. 3 NLP researchers performed CT for one event type and showed much sharper learning curves with all three exceeding ACE performance in less than ninety minutes, suggesting that CT can provide further benefits when the annotator deeply understands the system.


title: The frictional force on sliding drops

abstract: The dynamic frictional force between solid surfaces in relative motion differs from the static force needed to initiate motion, but this distinction is not usually thought to occur for liquid drops moving on a solid. Recent experiments [Gao, et al., Nature Phys. {\bf 114}, 191 (2018)] have challenged this view, and claim to observe an analog of solid-on-solid friction for sliding drops. We use molecular dynamic simulations to investigate the forces that moving liquids exert on solids in several situations. In contrast to the indirect techniques required in laboratory experiments, the forces involved in friction are directly accessible in these calculations. We find that, aside from possible inertial effects due to the abrupt initiation of motion and aging effects for unconfined drops, the frictional forces are constant in time.


title: Ensemble Clustering for Graphs

abstract: We propose an ensemble clustering algorithm for graphs (ECG), which is based on the Louvain algorithm and the concept of consensus clustering. We validate our approach by replicating a recently published study comparing graph clustering algorithms over artificial networks, showing that ECG outperforms the leading algorithms from that study. We also illustrate how the ensemble obtained with ECG can be used to quantify the presence of community structure in the graph.


title: Phase space exchange-based bunch compression with reduced CSR effects

abstract: Through an interaction of theory and optimization-guided design, we have developed a novel bunch compression scheme based on an asymmetric double emittance exchange (DEEX) with greatly reduced coherent synchrotron radiation (CSR) effects. First, we discuss the benefits of the phase space exchange-based bunch compressor (BC) in comparison to the chicane-based counterpart in the approximation of linear single particle dynamics and with respect to microbunching instabilities driven by longitudinal space charge forces. We present simulation results for nonlinear dynamics in a symmetric DEEX BC which accounts for the CSR effects of the scheme with direct and mirrored longitudinal phase space at the exit of the beamline. The later scheme demonstrates self-compensation of CSR effects on the longitudinal dynamics. Next, we optimize the parameters of the beam elements with an adaptive feedback, resulting in an asymmetric configuration of the scheme with compensated CSR and nonlinear effects on the transverse dynamics, which we further improve by incorporating several sextupoles relying on an eigen-emittance analysis. Finally, we discuss how this method can be generalized for optimizing nonlinear beam dynamics with collective effects.


title: Assessing Bayes factor surfaces using interactive visualization and computer surrogate modeling

abstract: Bayesian model selection provides a natural alternative to classical hypothesis testing based on p-values. While many papers mention that Bayesian model selection is frequently sensitive to prior specification on the parameters, there are few practical strategies to assess and report this sensitivity. This article has two goals. First, we aim educate the broader statistical community about the extent of potential sensitivity through visualization of the Bayes factor surface. The Bayes factor surface shows the value a Bayes factor takes (usually on the log scale) as a function of user-specified hyperparameters. We provide interactive visualization through an R shiny application that allows the user to explore sensitivity in Bayes factor over a range of hyperparameter settings in a familiar regression setting. We compare the surface with three automatic procedures. Second, we suggest surrogate modeling via Gaussian processes (GPs) to visualize the Bayes factor surface in situations where computation of Bayes factors is expensive. That is, we treat Bayes factor calculation as a computer simulation experiment. In this context, we provide a fully reproducible example using accessible GP libraries to augment an important study of the influence of outliers in empirical finance. We suggest Bayes factor surfaces are valuable for scientific reporting since they (i) increase transparency, making potential instability in Bayes factors easy to visualize, (ii) generalize to simple and more complicated examples, and (iii) provide a path for researchers to assess the impact of prior choice on modeling decisions in a wide variety research areas.


title: HD 38858: a solar type star with a $\sim$10.8 yr activity cycle. Searching for Balmer and metallic line variations

abstract: Context. The detection of chromospheric activity cycles in solar-analogue and twin stars, can be used to put in context the solar cycle. However, there is a relatively low percentage of these stars with activity cycles detected. It is well known that the cores of Ca II H&K lines are modulated by stellar activity. The behaviour of Balmer and other optical lines with stellar activity is not yet completely understood. Aims. To search for Ca II H&K, Balmer and Fe II line variations modulated by the stellar activity. In particular, we apply a novel strategy to detect possible shape variations in the H{\alpha} line. Methods. We analyse activity signatures in HD 38858 by using HARPS and CASLEO spectra obtained between 2003 and 2017. We calculate the Mount Wilson index (S$_{MW}$), log($\mathrm{R}'_\mathrm{HK}$) and the statistical moments of Ca II H&K, Balmer and other optical lines. We search for periodicities by using the generalized Lomb-Scargle periodogram. Results. We detect a long-term activity cycle of 10.8 yr in Ca II H&K and H{\alpha}, being HD 38858 one of the most similar stars to the Sun for which a cycle has been detected. In contrast, this cycle is marginally detected in the Fe II lines. Also, we detect a noticeable radial velocity variation, which seems to be produced by stellar activity. Conclusions. HD 38858 is the second solar-analogue star where we find a clear activity cycle which is replicated by Balmer lines. Spectral indexes based on the shape of H{\alpha} seems to be more realiable than the fluxes in the same line to detect activity variations. On the other hand, the cyclic modulation detected gives place to a radial velocity variation, previously associated to a super-Earth planet. Finally, due to the similarity of HD 38858 with the Sun, we recommend to continue monitoring this star.


title: Loschmidt Amplitude and Work Distribution in Quenches of the Sine-Gordon Model

abstract: The Sine-Gordon - equivalently, the massive Thirring - Hamiltonian is ubiquitous in low-dimensional physics, with applications that range from cold atom and strongly correlated systems to quantum impurities. We study here its non-equilibrium dynamics using the quantum quench protocol - following the system as it evolves under the Sine-Gordon Hamiltonian from initial Mott type states with large potential barriers. By means of the Bethe Ansatz we calculate exactly the Loschmidt amplitude, the fidelity and work distribution characterizing these quenches for different values of the interaction strength. Some universal features are noted as well as an interesting duality relating quenches in different parameter regimes of the model.


title: Image registration and super resolution from first principles

abstract: Image registration is the inference of transformations relating noisy and distorted images. It is fundamental in computer vision, experimental physics, and medical imaging. Many algorithms and analyses exist for inferring shift, rotation, and nonlinear transformations between image coordinates. Even in the simplest case of translation, however, all known algorithms are biased and none have achieved the precision limit of the Cramer Rao bound (CRB). Following Bayesian inference, we prove that the standard method of shifting one image to match another cannot reach the CRB. We show that the bias can be cured and the CRB reached if, instead, we use Super Registration: learning an optimal model for the underlying image and shifting that to match the data. Our theory shows that coarse-graining oversampled images can improve registration precision of the standard method. For oversampled data, our method does not yield striking improvements as measured by eye. In these cases, however, we show our new registration method can lead to dramatic improvements in extractable information, for example, inferring $10\times$ more precise particle positions.


title: Gemini GNIRS near-infrared spectroscopy of 50 quasars at z>~5.7

abstract: We report initial results from a large Gemini program to observe z>~5.7 quasars with GNIRS near-IR spectroscopy. Our sample includes 50 quasars with simultaneous ~0.85-2.5 micron spectra covering the rest-frame ultraviolet and major broad emission lines from Ly-alpha to MgII. We present spectral measurements for these quasars and compare to their lower-redshift counterparts at z=1.5-2.3. We find that when quasar luminosity is matched, there are no significant differences between the rest-UV spectra of z>~5.7 quasars and the low-z comparison sample. High-z quasars have similar continuum and emission line properties and occupy the same region in the black hole mass and luminosity space as the comparison sample, accreting at an average Eddington ratio of ~0.3. There is no evidence for super-Eddington accretion or hypermassive (>10^10 Msun) black holes within our sample. We find a mild excess of quasars with weak CIV lines relative to the control sample. Our results, corroborating earlier studies but with better statistics, demonstrate that these high-z quasars are already mature systems of accreting supermassive black holes operating with the same physical mechanisms as those at lower redshifts.


title: Media Accessibility Policy in Theory and Reality: Empirical Outreach to Audio Description Users in the United States

abstract: Audio description, a form of trans-modal media translation, allows people who are blind or visually impaired access to visually-oriented, socio-cultural, or historical public discourse alike. Although audio description has gained more prominence in media policy and research lately, it rarely has been studied empirically. Yet this paper presents quantitative and qualitative survey data on its challenges and opportunities, through the analysis of responses from 483 participants in a national sample, with 334 of these respondents being blind. Our results give insight into audio description use in broadcast TV, streaming services, for physical media, such as DVDs, and in movie theaters. We further discover a multiplicity of barriers and hindrances which prevent a better adoption and larger proliferation of audio description. In our discussion, we present a possible answer to these problems - the UniDescription Project - a media ecosystem for the creation, curation, and dissemination of audio description for multiple media platforms.


title: Pressure and escape rates for random subshifts of finite type

abstract: In this work we consider several aspects of the thermodynamic formalism in a randomized setting. Let $X$ be a non-trivial mixing shift of finite type, and let $f : X \to \mathbb{R}$ be a H\"older continuous potential with associated Gibbs measure $\mu$. Further, fix a parameter $\alpha \in (0,1)$. For each $n \geq 1$, let $\mathcal{F}_n$ be a random subset of words of length $n$, where each word of length $n$ that appears in $X$ is included in $\mathcal{F}_n$ with probability $1-\alpha$ (and excluded with probability $\alpha$), independently of all other words. Then let $Y_n = Y(\mathcal{F}_n)$ be the random subshift of finite type obtained by forbidding the words in $\mathcal{F}_n$ from $X$. In our first main result, for $\alpha$ sufficiently close to $1$ and $n$ tending to infinity, we show that the pressure of $f$ on $Y_n$ converges in probability to the value $P_X(f) + \log(\alpha)$, where $P_X(f)$ is the pressure of $f$ on $X$. Additionally, let $H_n = H(\mathcal{F}_n)$ be the random hole in $X$ consisting of the union of the cylinder sets of the words in $\mathcal{F}_n$. For our second main result, for $\alpha$ sufficiently close to one and $n$ tending to infinity, we show that the escape rate of $\mu$-mass through $H_n$ converges in probability to the value $-\log(\alpha)$ as $n$ tends to infinity.


title: Parameter Optimization and Uncertainty Analysis of FREYA for Spontaneous Fission

abstract: Background: The complete event fission simulation code FREYA is used to study correlations in fission. To make the best possible simulations, FREYA one must find the optimized values of the five physics-based parameters which characterize the events produced by FREYA. So far this has been done only empirically or by brute force computational techniques. Purpose: We seek to develop a method of parameter optimization for FREYA that would include uncertainty quantification and an analysis of the correlations between the parameters. Methods: We focus on spontaneous fission as a simpler test case of the method. We first check the results of previous optimizations for $^{252}$Cf(sf) and then go on to develop a simulated annealing approach to optimize the parameters. Although $^{252}$Cf(sf) has the most measured observables, we are able to apply the technique to all spontaneously fissioning isotopes in the current version of FREYA. Results: We find optimal values of the five FREYA parameters, with uncertainties, for all isotopes studied. Our results are compared with available data. The parameter values are physically reasonable and in accord with intuition, even in cases where the data for optimization are sparse to non-existent. Conclusions: The simulated annealing method provides a way to determine the best parameters for the phenomenologically-successful fission model FREYA. We have developed an algorithm that can be updated to include additional isotopes or more data sets for the current isotopes as they become available.