title: Complexity of quantum state verification in the quantum linear systems problem

abstract: We analyze the complexity of quantum state verification in the context of solving systems of linear equations of the form $A \vec x = \vec b$. We show that any quantum operation that verifies whether a given quantum state is within a constant distance from the solution of the quantum linear systems problem requires $q=\Omega(\kappa)$ uses of a unitary that prepares a quantum state $\left| b \right>$, proportional to $\vec b$, and its inverse in the worst case. Here, $\kappa$ is the condition number of the matrix $A$. For typical instances, we show that $q=\Omega(\sqrt \kappa)$ with high probability. These lower bounds are almost achieved if quantum state verification is performed using known quantum algorithms for the quantum linear systems problem. We also analyze the number of copies of $\left| b \right>$ required by verification procedures of the prepare and measure type. In this case, the lower bounds are quadratically worse, being $\Omega(\kappa^2)$ in the worst case and $\Omega(\kappa)$ in typical instances with high probability. We discuss the implications of our results to known variational and related approaches to this problem, where state preparation, gate, and measurement errors will need to decrease rapidly with $\kappa$ for worst-case and typical instances if error correction is not used, and present some open problems.


title: Coexistence of topological Weyl and nodal-ring states in ferromagnetic and ferrimagnetic double perovskites

abstract: Magnetic topological quantum materials have attracted great attention due to their exotic topological quantum physics induced by the interplay among crystalology, magnetism, and topology, which is of profound importance to fundamental research and technology applications. However, limited materials are experimentally available, most of whom are realized by magnetic impurity doping or heterostructural constructions. In this work, based on the first-principles calculations, we predict that double perovskite Ba2CdReO6 is an intrinsic ferromagnetic topological semi-half-metal, while the ferrimagnetic double perovskite with space group symmetry Fm-3m, such as Ba2FeMoO6, belongs to a topological half-metal. One pair of Weyl points and fully spin-polarized nodal-ring states are found in the vicinity of the Fermi level in Ba2CdReO6. Its two-dimensional nearly flat drumhead surface states are fully spin-polarized. In Ba2FeMoO6, however, there exist four pairs of Weyl points and two fully spin-polarized nodal-rings near the Fermi level. These topological properties are stable in the presence of spin-orbit coupling. This makes these materials be an appropriate platform for studying the emerging intriguing properties, especially for the applications in spintronics, information technology, and topological superconductivity.


title: The Unreasonable Effectiveness of Machine Learning in Moldavian versus Romanian Dialect Identification

abstract: Motivated by the seemingly high accuracy levels of machine learning models in Moldavian versus Romanian dialect identification and the increasing research interest on this topic, we provide a follow-up on the Moldavian versus Romanian Cross-Dialect Topic Identification (MRC) shared task of the VarDial 2019 Evaluation Campaign. The shared task included two sub-task types: one that consisted in discriminating between the Moldavian and Romanian dialects and one that consisted in classifying documents by topic across the two dialects of Romanian. Participants achieved impressive scores, e.g. the top model for Moldavian versus Romanian dialect identification obtained a macro F1 score of 0.895. We conduct a subjective evaluation by human annotators, showing that humans attain much lower accuracy rates compared to machine learning (ML) models. Hence, it remains unclear why the methods proposed by participants attain such high accuracy rates. Our goal is to understand (i) why the proposed methods work so well (by visualizing the discriminative features) and (ii) to what extent these methods can keep their high accuracy levels, e.g. when we shorten the text samples to single sentences or when we use tweets at inference time. A secondary goal of our work is to propose an improved ML model using ensemble learning. Our experiments show that ML models can accurately identify the dialects, even at the sentence level and across different domains (news articles versus tweets). We also analyze the most discriminative features of the best performing models, providing some explanations behind the decisions taken by these models. Interestingly, we learn new dialectal patterns previously unknown to us or to our human annotators. Furthermore, we conduct experiments showing that the machine learning performance on the MRC shared task can be improved through an ensemble based on stacking.


title: MHD accretion-ejection: jets launched by a non-isotropic accretion disk dynamo. I. Validation and application of selected dynamo tensorial components

abstract: Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. The origin of the jet launching magnetic field is one of the open questions for modeling the accretion-ejection process. Here we address the question how to generate the accretion disk magnetization and field structure required for jet launching. Applying the PLUTO code, we present the first resistive MHD simulations of jet launching including a non-scalar accretion disk mean-field $\alpha^2\Omega$-dynamo in the context of large scale disk-jet simulations. Essentially, we find the $\alpha_\phi$-dynamo component determining the amplification of the poloidal magnetic field, which is strictly related to the disk magnetization (and, as a consequence, to the jet speed, mass and collimation), while the $\alpha_R$ and $\alpha_\theta$-dynamo components trigger the formation of multiple, anti-aligned magnetic loops in the disk, with strong consequences on the stability and dynamics of the disk-jet system. In particular, such loops trigger the formation of dynamo inefficient zones, which are characterized by a weak magnetic field, and therefore a lower value of the magnetic diffusivity. The jet mass, speed and collimation are strongly affected by the formation of the dynamo inefficient zones. Moreover, the $\theta$-component of the $\alpha$-dynamo plays a key role when interacting with a non-radial component of the seed magnetic field. We also present correlations between the strength of the disk toy dynamo coefficients and the dynamical parameters of the jet that is launched.


title: X-ray evolution of the nova V959 Mon suggests a delayed ejection and a non-radiative shock

abstract: X-ray observations of shocked gas in novae can provide a useful probe of the dynamics of the ejecta. Here we report on X-ray observations of the nova V959 Mon, which was also detected in GeV gamma-rays with the Fermi satellite. We find that the X-ray spectra are consistent with a two-temperature plasma model with non-solar abundances. We interpret the X-rays as due to shock interaction between the slow equatorial torus and the fast polar outflow that were inferred from radio observations of V959 Mon. We further propose that the hotter component, responsible for most of the flux, is from the reverse shock driven into the fast outflow. We find a systematic drop in the column density of the absorber between Days 60 and 140, consistent with the expectations for such a picture. We present intriguing evidence for a delay of around 40 days in the expulsion of the ejecta from the central binary. Moreover, we infer a relatively small (a few times 10$^{-6}$ Msun) ejecta mass ahead of the shock, considerably lower than the mass of 10$^4$ K gas inferred from radio observations. Finally, we infer that the dominant X-ray shock was likely not radiative at the time of our observations, and that the shock power was considerably higher than the observed X-ray luminosity. It is unclear why high X-ray luminosity, closer to the inferred shock power, is never seen in novae at early times, when the shock is expected to have high enough density to be radiative.


title: Improving Multi-Agent Cooperation using Theory of Mind

abstract: Recent advances in Artificial Intelligence have produced agents that can beat human world champions at games like Go, Starcraft, and Dota2. However, most of these models do not seem to play in a human-like manner: People infer others' intentions from their behaviour, and use these inferences in scheming and strategizing. Here, using a Bayesian Theory of Mind (ToM) approach, we investigated how much an explicit representation of others' intentions improves performance in a cooperative game. We compared the performance of humans playing with optimal-planning agents with and without ToM, in a cooperative game where players have to flexibly cooperate to achieve joint goals. We find that teams with ToM agents significantly outperform non-ToM agents when collaborating with all types of partners: non-ToM, ToM, as well as human players, and that the benefit of ToM increases the more ToM agents there are. These findings have implications for designing better cooperative agents.


title: Job market effects of COVID-19 on urban Ukrainian households

abstract: The employment status of billions of people has been affected by the COVID epidemic around the Globe. New evidence is needed on how to mitigate the job market crisis, but there exists only a handful of studies mostly focusing on developed countries. We fill in this gap in the literature by using novel data from Ukraine, a transition country in Eastern Europe, which enacted strict quarantine policies early on. We model four binary outcomes to identify respondents (i) who are not working during quarantine, (ii) those who are more likely to work from home, (iii) respondents who are afraid of losing a job, and, finally, (iv) survey participants who have savings for 1 month or less if quarantine is further extended. Our findings suggest that respondents employed in public administration, programming and IT, as well as highly qualified specialists, were more likely to secure their jobs during the quarantine. Females, better educated respondents, and those who lived in Kyiv were more likely to work remotely. Working in the public sector also made people more confident about their future employment perspectives. Although our findings are limited to urban households only, they provide important early evidence on the correlates of job market outcomes, expectations, and financial security, indicating potential deterioration of socio-economic inequalities.


title: A note on the class of languages generated by F-systems over regular languages

abstract: An F-system is a computational model that performs a folding operation on words of a given language, following directions coded on words of another given language. This paper considers the case in which both given languages are regular, and it shows that the class of languages generated by such F-systems is a proper subset of the class of linear context-free languages.


title: Machine learning for complete intersection Calabi-Yau manifolds: a methodological study

abstract: We revisit the question of predicting both Hodge numbers $h^{1,1}$ and $h^{2,1}$ of complete intersection Calabi-Yau (CICY) 3-folds using machine learning (ML), considering both the old and new datasets built respectively by Candelas-Dale-Lutken-Schimmrigk / Green-H\"ubsch-Lutken and by Anderson-Gao-Gray-Lee. In real world applications, implementing a ML system rarely reduces to feed the brute data to the algorithm. Instead, the typical workflow starts with an exploratory data analysis (EDA) which aims at understanding better the input data and finding an optimal representation. It is followed by the design of a validation procedure and a baseline model. Finally, several ML models are compared and combined, often involving neural networks with a topology more complicated than the sequential models typically used in physics. By following this procedure, we improve the accuracy of ML computations for Hodge numbers with respect to the existing literature. First, we obtain 97% (resp. 99%) accuracy for $h^{1,1}$ using a neural network inspired by the Inception model for the old dataset, using only 30% (resp. 70%) of the data for training. For the new one, a simple linear regression leads to almost 100% accuracy with 30% of the data for training. The computation of $h^{2,1}$ is less successful as we manage to reach only 50% accuracy for both datasets, but this is still better than the 16% obtained with a simple neural network (SVM with Gaussian kernel and feature engineering and sequential convolutional network reach at best 36%). This serves as a proof of concept that neural networks can be valuable to study the properties of geometries appearing in string theory.


title: Solar: $L_0$ solution path averaging for fast and accurate variable selection in high-dimensional data

abstract: We propose a new variable selection algorithm, subsample-ordered least-angle regression (solar), and its coordinate descent generalization, solar-cd. Solar re-constructs lasso paths using the $L_0$ norm and averages the resulting solution paths across subsamples. Path averaging retains the ranking information of the informative variables while averaging out sensitivity to high dimensionality, improving variable selection stability, efficiency, and accuracy. We prove that: (i) with a high probability, path averaging perfectly separates informative variables from redundant variables on the average $L_0$ path; (ii) solar variable selection is consistent and accurate; and (iii) the probability that solar omits weak signals is controllable for finite sample size. We also demonstrate that: (i) solar yields, with less than $1/3$ of the lasso computation load, substantial improvements over lasso in terms of the sparsity (64-84\% reduction in redundant variable selection) and accuracy of variable selection; (ii) compared with the lasso safe/strong rule and variable screening, solar largely avoids selection of redundant variables and rejection of informative variables in the presence of complicated dependence structures; (iii) the sparsity and stability of solar conserves residual degrees of freedom for data-splitting hypothesis testing, improving the accuracy of post-selection inference on weak signals with limited $n$; (iv) replacing lasso with solar in bootstrap selection (e.g., bolasso or stability selection) produces a multi-layer variable ranking scheme that improves selection sparsity and ranking accuracy with the computation load of only one lasso realization; and (v) given the computation resources, solar bootstrap selection is substantially faster (98\% lower computation time) than the theoretical maximum speedup for parallelized bootstrap lasso (confirmed by Amdahl's law).


title: The PDR structure and kinematics around the compact HII regions S235A and S235C with [CII], [13CII], [OI] and HCO+ line profiles

abstract: The aim of the present work is to study structure and gas kinematics in the photodissociation regions (PDRs) around the compact HII regions S235A and S235C. We observe the [CII], [13CII] and [OI] line emission, using SOFIA/upGREAT and complement them by data of HCO+ and CO. We use the [13CII] line to measure the optical depth of the [CII] emission, and find that the [CII] line profiles are influenced by self-absorption, while the [13CII] line remains unaffected by these effects. Hence, for dense PDRs, [13CII] emission is a better tracer of gas kinematics. The optical depth of the [CII] line is up to 10 in S235A. We find an expanding motion of the [CII]-emitting layer of the PDRs into the front molecular layer in both regions. Comparison of the gas and dust columns shows that gas components visible neither in the [CII] nor in low-J CO lines may contribute to the total column across S235A. We test whether the observed properties of the PDRs match the predictions of spherical models of expanding HII region + PDR + molecular cloud. Integrated intensities of the [13CII], [CII] and [OI] lines are well-represented by the model, but the models do not reproduce the double-peaked [CII] line profiles due to an insufficient column density of C+. The model predicts that the [OI] line could be a more reliable tracer of gas kinematics, but the foreground self-absorbing material does not allow using it in the considered regions.


title: An Asymptotic Lower Bound for Online Vector Bin Packing

abstract: We consider the online vector bin packing problem where $n$ items specified by $d$-dimensional vectors must be packed in the fewest number of identical $d$-dimensional bins. Azar et al. (STOC'13) showed that for any online algorithm $A$, there exist instances I, such that $A(I)$, the number of bins used by $A$ to pack $I$, is $\Omega(d/\log^2 d)$ times $OPT(I)$, the minimal number of bins to pack $I$. However in those instances, $OPT(I)$ was only $O(\log d)$, which left open the possibility of improved algorithms with better asymptotic competitive ratio when $OPT(I) \gg d$. We rule this out by showing that for any arbitrary function $q(\cdot)$ and any randomized online algorithm $A$, there exist instances $I$ such that $ E[A(I)] \geq c\cdot d/\log^3d \cdot OPT(I) + q(d)$, for some universal constant $c$.


title: Privacy Enhancing Machine Learning via Removal of Unwanted Dependencies

abstract: The rapid rise of IoT and Big Data has facilitated copious data driven applications to enhance our quality of life. However, the omnipresent and all-encompassing nature of the data collection can generate privacy concerns. Hence, there is a strong need to develop techniques that ensure the data serve only the intended purposes, giving users control over the information they share. To this end, this paper studies new variants of supervised and adversarial learning methods, which remove the sensitive information in the data before they are sent out for a particular application. The explored methods optimize privacy preserving feature mappings and predictive models simultaneously in an end-to-end fashion. Additionally, the models are built with an emphasis on placing little computational burden on the user side so that the data can be desensitized on device in a cheap manner. Experimental results on mobile sensing and face datasets demonstrate that our models can successfully maintain the utility performances of predictive models while causing sensitive predictions to perform poorly.


title: Detecting Distrust Towards the Skills of a Virtual Assistant Using Speech

abstract: Research has shown that trust is an essential aspect of human-computer interaction directly determining the degree to which the person is willing to use the system. An automatic prediction of the level of trust that a user has on a certain system could be used to attempt to correct potential distrust by having the system take relevant actions like, for example, explaining its actions more thoroughly. In this work, we explore the feasibility of automatically detecting the level of trust that a user has on a virtual assistant (VA) based on their speech. We use a dataset collected for this purpose, containing human-computer speech interactions where subjects were asked to answer various factual questions with the help of a virtual assistant, which they were led to believe was either very reliable or unreliable. We find that the subject's speech can be used to detect which type of VA they were using, which could be considered a proxy for the user's trust toward the VA's abilities, with an accuracy up to 76\%, compared to a random baseline of 50\%. These results are obtained using features that have been previously found useful for detecting speech directed to infants and non-native speakers.


title: Description of the UPPAAL Models for SRP and CSRP and Verification of their Termination and Consistency Properties

abstract: The IEEE Audio Video Bridging (AVB) Task Group (TG) was created to provide Ethernet with soft real-time guarantees. Later on, the TG was renamed to Time-Sensitive Networking (TSN) and its scope broadened to support hard real-time and critical applications. The Stream Reservation Protocol (SRP) is a key work of the TGs as it allows reserving resources in the network, guaranteeing the required quality of service (QoS). AVB's SRP is based on a distributed architecture, while TSN's is based on centralized ones. The distributed version of SRP is supported and used in TSN. Nevertheless, it was not designed to provide properties that are important for critical applications. Therefore, we propose a new version of SRP with enhanced services called Consistent Stream Reservation Protocol (CSRP). In this document we describe the SRP and CSRP UPPAAL models we developed and the queries we used to verify their termination and consistency properties.


title: Protocol for Optically Pumping AlH$^+$ to a Pure Quantum State

abstract: We propose an optical pumping scheme to prepare trapped $\mathrm{AlH}^+$ molecules in a pure state, the stretched hyperfine state $\lvert F=\frac{7}{2},\, m_F=\frac{7}{2}\rangle$ of the rovibronic ground manifold $\lvert \mathrm{X}^2\Sigma^+,\, v=0,\, N=0\rangle$. Our scheme utilizes linearly-polarized and circularly-polarized fields of a broadband pulsed laser to cool the rotational degree of freedom and drive the population to the hyperfine state, respectively. We simulate the population dynamics by solving a representative system of rate equations that accounts for the laser fields, blackbody radiation, and spontaneous emission. In order to model the hyperfine structure, new hyperfine constants of the $\mathrm{A}^2\Pi$ excited state were computed using a RASSCF wavefunction. We find that adding an infrared laser to drive the $1 \,-\; 0$ vibrational transition within the $ \mathrm{X}^2\Sigma^+$ manifold accelerates the cooling process. The results show that under optimum conditions, the population in the target state of the rovibronic ground manifold can reach 63 $\%$ after 68 $\mathrm{\mu}$s (330 ms) and 95 $\%$ after 25 ms (1.2 s) with (without) the infrared laser.


title: An oscillation-free fully partitioned scheme for the numerical modeling of cardiac active mechanics

abstract: In silico models of cardiac electromechanics couple together mathematical models describing different physics. One instance is represented by the model describing the generation of active force, coupled with the one of tissue mechanics. For the numerical solution of the coupled model, partitioned schemes, that foresee the sequential solution of the two subproblems, are often used. However, this approach may be unstable. For this reason, the coupled model is commonly solved as a unique system using Newton type algorithms, at the price, however, of high computational costs. In light of this motivation, in this paper we propose a new numerical scheme, that is numerically stable and accurate, yet within a fully partitioned (i.e. segregated) framework. Specifically, we introduce, with respect to standard segregated scheme, a numerically consistent stabilization term, capable of removing the nonphysical oscillations otherwise present in the numerical solution of the commonly used segregated scheme. Our new method is derived moving from a physics-based analysis on the microscale energetics of the force generation dynamics. By considering a model problem of active mechanics we prove that the proposed scheme is unconditionally absolutely stable (i.e. it is stable for any time step size), unlike the standard segregated scheme, and we also provide an interpretation of the scheme as a fractional step method. We show, by means of several numerical tests, that the proposed stabilization term successfully removes the nonphysical numerical oscillations characterizing the non stabilized segregated scheme solution. Our numerical tests are carried out for several force generation models available in the literature, namely the Niederer-Hunter-Smith model, the model by Land and coworkers, and the mean-field force generation model that we have recently proposed. Finally, we apply the proposed scheme [...]


title: Doppler beaming factors for white dwarfs, main sequence stars, and giant stars Limb-darkening coefficients for 3D (DA and DB) white dwarf models

abstract: We computed Doppler beaming factors for DA, DB, and DBA white dwarf models, as well as for main sequence and giant stars covering the transmission curves of the Sloan, UBVRI, HiPERCAM, Kepler, TESS, and Gaia photometric systems. The calculations of the limb-darkening coefficients for 3D models were carried out using the least-squares method for these photometric systems. The beaming factor calculations, which use realistic models of stellar atmospheres, show that the black body approximation is not accurate, particularly for the filters $u$, $u'$, $U$, $g$, $g'$, and $B$. The black body approach is only valid for high effective temperatures and/or long effective wavelengths. Therefore, for more accurate analyses of light curves, we recommend the use of the beaming factors presented in this paper. Concerning limb-darkening, the distribution of specific intensities for 3D models indicates that, in general, these models are less bright toward the limb than their 1D counterparts, which implies steeper profiles. To describe these intensities better, we recommend the use of the four-term law (also for 1D models) given the level of precision that is being achieved with Earth-based instruments and space missions such as Kepler and TESS (and PLATO in the future).


title: Derivations and automorphisms of locally matrix algebras

abstract: We describe derivations and automorphisms of infinite tensor products of matrix algebras. Using this description we show that for a countable--dimensional locally matrix algebra $A$ over a field $\mathbb{F}$ the dimension of the Lie algebra of outer derivations of $A$ and the order of the group of outer automorphisms of $A$ are both equal to $|\mathbb{F}|^{\aleph_0},$ where $|\mathbb{F}|$ is the cardinality of the field $\mathbb{F}.$


title: Physics in precision-dependent normal neighborhoods

abstract: We introduce a procedure to determine the size and shape of normal neighborhoods in any spacetimes and their dependence on the precision of the measurements performed by arbitrary observers. As an example, we consider the Schwarzschild geometry in Riemann and Fermi normal coordinates and determine the size and shape of normal neighborhoods in the vicinity of the event horizon. Depending on the observers, normal neighborhoods extend to the event horizon and even beyond into the black hole interior. It is shown that the causal structure supported by normal neighborhoods across an event horizon is consistent with general relativity. In particular, normal neighborhoods reaching over an event horizon are void of the Schwarzschild coordinate singularity. In addition, we introduce a new variant of normal coordinates which we call Fermi normal coordinates around a point, unifying features of Riemann and Fermi normal coordinates, and analyze their neighborhoods.


title: The $\eta$-pseudo-hermitic generator in the deformed Woods Saxons potential

abstract: In this paper, we present a general method to solve non-hermetic potentials with PT symmetry using the definition of two $\eta$-pseudo-hermetic and first-order operators. This generator applies to the Dirac equation which consists of two spinor wave functions and non-hermetic potentials, with position that mass is considered a constant value and also Hamiltonian hierarchy method and the shape invariance property are used to perform calculations. Furthermore, we show the correlation between the potential parameters with transmission probabilities that $\eta$-pseudo-hermetic using the change of focal points on Hamiltonian can be formalized based on Schr\"{o}dinger-like equation. By using this method for some solvable potentials such as deformed Woods Saxon's potential, it can be shown that these real potentials can be decomposed into complex potentials consisting of eigenvalues of a class of $\eta$-pseudo-hermitic.


title: The Entropic Dynamics of Spin

abstract: In the Entropic Dynamics (ED) approach the essence of quantum theory lies in its probabilistic nature while the Hilbert space structure plays a secondary and ultimately optional role. The dynamics of probability distributions is driven by the maximization of an entropy subject to constraints that carry the relevant physical information -- directionality, correlations, gauge interactions, etc. The challenge is to identify those constraints and to establish a criterion for how the constraints themselves are updated. In this paper the ED framework is extended to describe a spin-1/2 point particle. In ED spin is neither modelled as a rotating body, nor through the motion of a point particle; it is an epistemic property of the wave function. The constraint that reflects the peculiar rotational properties of spin is most effectively expressed in the language of geometric algebra. The updating of all constraints is carried out in a way that stresses the central importance of symmetry principles. First we identify the appropriate symplectic and metric structures in the phase space of probabilities, their conjugate momenta, and the spin variables. This construction yields a derivation of the Fubini-Study metric for a spin-1/2 particle which highlights its deep connection to information geometry. Then we construct an ED that preserves both the symplectic structure (a Hamiltonian flow) and the metric structure (a Killing flow). We show that generic Hamiltonian-Killing flows are linear in the wave function. Imposing further that the Hamiltonian be the generator of an entropic evolution in time leads to an entropic dynamics described by the Pauli equation. We conclude with a discussion of the new interpretation of the formalism which yields a physical picture that is significantly different from that provided by other interpretations.


title: Algebraic 3D Graphic Statics: reciprocal constructions

abstract: The recently developed 3D graphic statics (3DGS) lacks a rigorous mathematical definition relating the geometrical and topological properties of the reciprocal polyhedral diagrams as well as a precise method for the geometric construction of these diagrams. This paper provides a fundamental algebraic formulation for 3DGS by developing equilibrium equations around the edges of the primal diagram and satisfying the equations by the closeness of the polygons constructed by the edges of the corresponding faces in the dual/reciprocal diagram. The research provides multiple numerical methods for solving the equilibrium equations and explains the advantage of using each technique. The approach of this paper can be used for compression-and-tension combined form-finding and analysis as it allows constructing both the form and force diagram based on the interpretation of the input diagram. Besides, the paper expands on the geometric/static degrees of (in)determinacies of the diagrams using the algebraic formulation and shows how these properties can be used for the constrained manipulation of the polyhedrons in an interactive environment without breaking the reciprocity between the two.


title: Dimension Groups and Dynamical Systems

abstract: We give a description of the link between topological dynamical systems and their dimension groups. The focus is on minimal systems and, in particular, on substitution shifts. We describe in detail the various classes of systems including Sturmian shifts and interval exchange shifts. This is a preliminary version of a book which will be published by Cambridge University Press. Any comments are of course welcome.


title: Dynamic Transitions of the Swift-Hohenberg Equation with Third-Order Dispersion

abstract: The Swift-Hohenberg equation is ubiquitous in the study of bistable dynamics. In this paper, we study the dynamic transitions of the Swift-Hohenberg equation with a third-order dispersion term in one spacial dimension with a periodic boundary condition. As a control parameter crosses a critical value, the trivial stable equilibrium solution will lose its stability, and undergoes a dynamic transition to a new physical state, described by a local attractor. The main result of this paper is to fully characterize the type and detailed structure of the transition using dynamic transition theory. In particular, employing techniques from center manifold theory, we reduce this infinite dimensional problem to a finite one since the space on which the exchange of stability occurs is finite dimensional. The problem then reduces to analysis of single or double Hopf bifurcations, and we completely classify the possible phase changes depending on the dispersion for every spacial period.


title: Orbital angular momentum beam generation using a free-electron laser oscillator

abstract: With wavelength tunability, free-electron lasers (FELs) are well-suited for generating orbital angular momentum (OAM) beams in a wide photon energy range. We report the first experimental demonstration of OAM beam generation using an oscillator FEL. Lasing around 458 nm, we have produced the four lowest orders of coherently mixed OAM beams with good beam quality, excellent stability, and substantial intracavity power. We have also developed a pulsed mode operation of the OAM beam with a highly reproducible temporal structure for a range of modulation frequencies from 1 to 30 Hz. This development can be extended to short wavelengths, for example to x-rays using a future x-ray FEL oscillator. The operation of such an OAM FEL also paves the way for the generation of OAM gamma-ray beams via Compton scattering.


title: MAPPER: Multi-Agent Path Planning with Evolutionary Reinforcement Learning in Mixed Dynamic Environments

abstract: Multi-agent navigation in dynamic environments is of great industrial value when deploying a large scale fleet of robot to real-world applications. This paper proposes a decentralized partially observable multi-agent path planning with evolutionary reinforcement learning (MAPPER) method to learn an effective local planning policy in mixed dynamic environments. Reinforcement learning-based methods usually suffer performance degradation on long-horizon tasks with goal-conditioned sparse rewards, so we decompose the long-range navigation task into many easier sub-tasks under the guidance of a global planner, which increases agents' performance in large environments. Moreover, most existing multi-agent planning approaches assume either perfect information of the surrounding environment or homogeneity of nearby dynamic agents, which may not hold in practice. Our approach models dynamic obstacles' behavior with an image-based representation and trains a policy in mixed dynamic environments without homogeneity assumption. To ensure multi-agent training stability and performance, we propose an evolutionary training approach that can be easily scaled to large and complex environments. Experiments show that MAPPER is able to achieve higher success rates and more stable performance when exposed to a large number of non-cooperative dynamic obstacles compared with traditional reaction-based planner LRA* and the state-of-the-art learning-based method.


title: Convexifying Multilinear Sets with Cardinality Constraints: Structural Properties, Nested Case and Extensions

abstract: The problem of minimizing a multilinear function of binary variables is a well-studied NP-hard problem. The set of solutions of the standard linearization of this problem is called the multilinear set. We study a cardinality constrained version of it with upper and lower bounds on the number of nonzero variables. We call the set of solutions of the standard linearization of this problem a multilinear set with cardinality constraints. We characterize a set of conditions on these multilinear terms (called properness) and observe that under these conditions the convex hull description of the set is tractable via an extended formulation. We then give an explicit polyhedral description of the convex hull when the multilinear terms have a nested structure. Our description has an exponential number of inequalities which can be separated in polynomial time. Finally, we generalize these inequalities to obtain valid inequalities for the general case.


title: Dynamic Hypergames for Synthesis of Deceptive Strategies with Temporal Logic Objectives

abstract: In this paper, we study the use of deception for strategic planning in adversarial environments. We model the interaction between the agent (player 1) and the adversary (player 2) as a two-player concurrent game in which the adversary has incomplete information about the agent's task specification in temporal logic. During the online interaction, the adversary can infer the agent's intention from observations and adapt its strategy so as to prevent the agent from satisfying the task. To plan against such an adaptive opponent, the agent must leverage its knowledge about the adversary's incomplete information to influence the behavior of the opponent, and thereby being deceptive. To synthesize a deceptive strategy, we introduce a class of hypergame models that capture the interaction between the agent and its adversary given asymmetric, incomplete information. A hypergame is a hierarchy of games, perceived differently by the agent and its adversary. We develop the solution concept of this class of hypergames and show that the subjectively rationalizable strategy for the agent is deceptive and maximizes the probability of satisfying the task in temporal logic. This deceptive strategy is obtained by modeling the opponent evolving perception of the interaction and integrating the opponent model into proactive planning. Following the deceptive strategy, the agent chooses actions to influence the game history as well as to manipulate the adversary's perception so that it takes actions that benefit the goal of the agent. We demonstrate the correctness of our deceptive planning algorithm using robot motion planning examples with temporal logic objectives and design a detection mechanism to notify the agent of potential errors in modeling of the adversary's behavior.


title: Cumulated burden of Covid-19 in Spain from a Bayesian perspective

abstract: The main goal of this work is to estimate the actual number of cases of Covid-19 in Spain in the period 01-31-2020 / 06-01-2020 by Autonomous Communities. Based on these estimates, this work allows us to accurately re-estimate the lethality of the disease in Spain, taking into account unreported cases. A hierarchical Bayesian model recently proposed in the literature has been adapted to model the actual number of Covid-19 cases in Spain. The results of this work show that the real load of Covid-19 in Spain in the period considered is well above the data registered by the public health system. Specifically, the model estimates show that, cumulatively until June 1st, 2020, there were 2,425,930 cases of Covid-19 in Spain with characteristics similar to those reported (95\% credibility interval: 2,148,261 - 2,813,864), from which were actually registered only 518,664. Considering the results obtained from the second wave of the Spanish seroprevalence study, which estimates 2,350,324 cases of Covid-19 produced in Spain, in the period of time considered, it can be seen that the estimates provided by the model are quite good. This work clearly shows the key importance of having good quality data to optimize decision-making in the critical context of dealing with a pandemic.


title: Revisiting non-linear functional brain co-activations: directed, dynamic and delayed

abstract: The center stage of neuro-imaging is currently occupied by studies of functional correlations between brain regions. These correlations define the brain functional networks, which are the most frequently used framework to represent and interpret a variety of experimental findings. In previous work we first demonstrated that the relatively stronger BOLD activations contain most of the information relevant to understand functional connectivity and subsequent work confirmed that a large compression of the original signals can be obtained without significant loss of information. In this work we revisit the correlation properties of these epochs to define a measure of nonlinear dynamic directed functional connectivity (nldFC) across regions of interest. We show that the proposed metric provides at once, without extensive numerical complications, directed information of the functional correlations, as well as a measure of temporal lags across regions, overall offering a different perspective in the analysis of brain co-activation patterns. In this paper we provide for a proof of concept, based on replicating and completing existing results on an Autism database, to discuss the main features and advantages of the proposed strategy for the study of brain functional correlations. These results show new interpretations of the correlations found on this sample.


title: Phenomenology and scaling of optimal flapping wing kinematics

abstract: Biological flapping wing fliers operate efficiently and robustly in a wide range of flight conditions and are a great source of inspiration to engineers. The unsteady aerodynamics of flapping-wings are dominated by large-scale vortical structures that augment the aerodynamic performance but are sensitive to minor changes in the wing actuation. We experimentally optimise the pitch angle kinematics of a flapping wing system in hover to maximise the stroke average lift and hovering efficiency using a evolutionary algorithm and in-situ force and torque measurements at the wing root. Additional flow field measurements are conducted to link the vortical flow structures to the aerodynamic performance for the Pareto-optimal kinematics. The optimised pitch angle profiles yielding maximum stroke-average lift coefficients have trapezoidal shapes and high average angles of attack. These kinematics create strong leading-edge vortices early in the cycle which enhance the force production on the wing. The most efficient pitch angle kinematics resemble sinusoidal evolutions and have lower average angles of attack. The leading-edge vortex grows slower and stays close-bound to the wing throughout the majority of the stroke-cycle. This requires less aerodynamic power and increases the hovering efficiency by 93% but sacrifices 43% of the maximum lift. In all cases, a leading-edge vortex is fed by vorticity through the leading edge shear-layer which makes the shear-layer velocity a good indicator for the growth of the vortex and its impact on the aerodynamic forces. We estimate the shear-layer velocity at the leading edge solely from the input kinematics and use it to scale the average and the time-resolved evolution of the circulation and the aerodynamic forces.


title: Decoherence-free radiofrequency dressed subspaces

abstract: We study the spectral signatures and coherence properties of radiofrequency dressed hyperfine Zeeman sub-levels of 87Rb. Experimentally, we engineer combinations of static and RF magnetic fields to modify the response of the atomic spin states to environmental magnetic field noise. We demonstrate analytically and experimentally the existence of 'magic' dressing conditions where decoherence due to electromagnetic field noise is strongly suppressed. Building upon this result, we propose a bi-chromatic dressing configuration that reduces the global sensitivity of the atomic ground states to low-frequency noise, and enables the simultaneous protection of multiple transitions between the two ground hyperfine manifolds of atomic alkali species. Our methods produce protected transitions between any pair of hyperfine sub-levels at arbitrary (low) DC-magnetic fields.


title: Is there something I'm missing? Topic Modeling in eDiscovery

abstract: In legal eDiscovery, the parties are required to search through their electronically stored information to find documents that are relevant to a specific case. Negotiations over the scope of these searches are often based on a fear that something will be missed. This paper continues an argument that discovery should be based on identifying the facts of a case. If a search process is less than complete (if it has Recall less than 100%), it may still be complete in presenting all of the relevant available topics. In this study, Latent Dirichlet Allocation was used to identify 100 topics from all of the known relevant documents. The documents were then categorized to about 80% Recall (i.e., 80% of the relevant documents were found by the categorizer, designated the hit set and 20% were missed, designated the missed set). Despite the fact that less than all of the relevant documents were identified by the categorizer, the documents that were identified contained all of the topics derived from the full set of documents. This same pattern held whether the categorizer was a na\"ive Bayes categorizer trained on a random selection of documents or a Support Vector Machine trained with Continuous Active Learning (which focuses evaluation on the most-likely-to-be-relevant documents). No topics were identified in either categorizer's missed set that were not already seen in the hit set. Not only is a computer-assisted search process reasonable (as required by the Federal Rules of Civil Procedure), it is also complete when measured by topics.


title: $\mathcal{PT}$-symmetry in compact phase space for a linear Hamiltonian

abstract: We study the time evolution of a PT-symmetric, non-Hermitian quantum system for which the associated phase space is compact. We focus on the simplest non-trivial example of such a Hamiltonian, which is linear in the angular momentum operators. In order to describe the evolution of the system, we use a particular disentangling decomposition of the evolution operator, which remains numerically accurate even in the vicinity of the Exceptional Point. We then analyze how the non-Hermitian part of the Hamiltonian affects the time evolution of two archetypical quantum states, coherent and Dicke states. For that purpose we calculate the Husimi distribution or Q function and study its evolution in phase space. For coherent states, the characteristics of the evolution equation of the Husimi function agree with the trajectories of the corresponding angular momentum expectation values. This allows to consider these curves as the trajectories of a classical system. For other types of quantum states, e.g. Dicke states, the equivalence of characteristics and trajectories of expectation values is lost.


title: Mean-square convergence rates of implicit Milstein type methods for SDEs with non-Lipschitz coefficients

abstract: A class of implicit Milstein type methods is introduced and analyzed in the present article for stochastic differential equations (SDEs) with non-globally Lipschitz drift and diffusion coefficients. By incorporating a pair of method parameters $\theta, \eta \in [0, 1]$ into both the drift and diffusion parts, the new schemes are indeed a kind of drift-diffusion double implicit methods. Within a general framework, we offer upper mean-square error bounds for the proposed schemes, based on certain error terms only getting involved with the exact solution processes. Such error bounds help us to easily analyze mean-square convergence rates of the schemes, without relying on a priori high-order moment estimates of numerical approximations. Putting further globally polynomial growth condition, we successfully recover the expected mean-square convergence rate of order one for the considered schemes with $\theta \in [\tfrac12, 1], \eta \in [0, 1]$. Also, some of the proposed schemes are applied to solve three SDE models evolving in the positive domain $(0, \infty)$. More specifically, the particular drift-diffusion implicit Milstein method ($ \theta = \eta = 1 $) is utilized to approximate the Heston $\tfrac32$-volatility model and the stochastic Lotka-Volterra competition model. The semi-implicit Milstein method ($\theta =1, \eta = 0$) is used to solve the Ait-Sahalia interest rate model. Thanks to the previously obtained error bounds, we reveal the optimal mean-square convergence rate of the positivity preserving schemes under more relaxed conditions, compared with existing relevant results in the literature. Numerical examples are also reported to confirm the previous findings.


title: On the inverse scattering problem for radially-symmetric domains in two dimensions

abstract: In the present paper we describe a method for solving inverse problems for the Helmholtz equation in radially-symmetric domains given multi-frequency data. Our approach is based on the construction of suitable trace formulas which relate the impedance of the total field at multiple frequencies to derivatives of the potential. Using this trace formula we obtain a system of coupled differential equations which can be solved to obtain the potential in a stable manner. Finally, the performance of the reconstruction algorithm is illustrated with several numerical examples.


title: Discovery of a Nearby Young Brown Dwarf Disk

abstract: We report the discovery of the youngest brown dwarf with a disk at 102 pc from the Sun, WISEA~J120037.79-784508.3 (W1200-7845), via the Disk Detective citizen science project. We establish that W1200-7845 is located in the 3.7$\substack{+4.6 \\ -1.4}$ Myr-old $\varepsilon$~Cha association. Its spectral energy distribution (SED) exhibits clear evidence of an infrared (IR) excess, indicative of the presence of a warm circumstellar disk. Modeling this warm disk, we find the data are best fit using a power-law description with a slope $\alpha = -0.94$, which suggests it is a young, Class II type disk. Using a single blackbody disk fit, we find $T_{eff, disk} = 521 K$ and $L_{IR}/L_{*} = 0.14$. The near-infrared spectrum of W1200-7845 matches a spectral type of M6.0$\gamma \pm 0.5$, which corresponds to a low surface gravity object, and lacks distinctive signatures of strong Pa$\beta$ or Br$\gamma$ accretion. Both our SED fitting and spectral analysis indicate the source is cool ($T_{eff} = $2784-2850 K), with a mass of 42-58 $M_{Jup}$, well within the brown dwarf regime. The proximity of this young brown dwarf disk makes the system an ideal benchmark for investigating the formation and early evolution of brown dwarfs.


title: Curating Long-term Vector Maps

abstract: Autonomous service mobile robots need to consistently, accurately, and robustly localize in human environments despite changes to such environments over time. Episodic non-Markov Localization addresses the challenge of localization in such changing environments by classifying observations as arising from Long-Term, Short-Term, or Dynamic Features. However, in order to do so, EnML relies on an estimate of the Long-Term Vector Map (LTVM) that does not change over time. In this paper, we introduce a recursive algorithm to build and update the LTVM over time by reasoning about visibility constraints of objects observed over multiple robot deployments. We use a signed distance function (SDF) to filter out observations of short-term and dynamic features from multiple deployments of the robot. The remaining long-term observations are used to build a vector map by robust local linear regression. The uncertainty in the resulting LTVM is computed via Monte Carlo resampling the observations arising from long-term features. By combining occupancy-grid based SDF filtering of observations with continuous space regression of the filtered observations, our proposed approach builds, updates, and amends LTVMs over time, reasoning about all observations from all robot deployments in an environment. We present experimental results demonstrating the accuracy, robustness, and compact nature of the extracted LTVMs from several long-term robot datasets.


title: Heaviside Set Constrained Optimization: Optimality and Newton Method

abstract: Data in the real world frequently involve binary status: truth or falsehood, positiveness or negativeness, similarity or dissimilarity, spam or non-spam, and to name a few, with applications into the regression, classification problems and so on. To characterize the binary status, one of the ideal functions is the Heaviside step function that returns one for one status and zero for the other. Hence, it is of dis-continuity. Because of this, the conventional approaches to deal with the binary status tremendously benefit from its continuous surrogates. In this paper, we target the Heaviside step function directly and study the Heaviside set constrained optimization: calculating the tangent and normal cones of the feasible set, establishing several first-order sufficient and necessary optimality conditions, as well as developing a Newton type method that enjoys locally quadratic convergence and excellent numerical performance.


title: Joint DOD and DOA Estimation in Slow-Time MIMO Radar via PARAFAC Decomposition

abstract: We develop a new tensor model for slow-time multiple-input multiple output (MIMO) radar and apply it for joint direction-of-departure (DOD) and direction-of-arrival (DOA) estimation. This tensor model aims to exploit the independence of phase modulation matrix and receive array in the received signal for slow-time MIMO radar. Such tensor can be decomposed into two tensors of different ranks, one of which has identical structure to that of the conventional tensor model for MIMO radar, and the other contains all phase modulation values used in the transmit array. We then develop a modification of the alternating least squares algorithm to enable parallel factor decomposition of tensors with extra constants. The Vandermonde structure of the transmit and receive steering matrices (if both arrays are uniform and linear) is then utilized to obtain angle estimates from factor matrices. The multi-linear structure of the received signal is maintained to take advantage of tensor-based angle estimation algorithms, while the shortage of samples in Doppler domain for slow-time MIMO radar is mitigated. As a result, the joint DOD and DOA estimation performance is improved as compared to existing angle estimation techniques for slow-time MIMO radar. Simulation results verify the effectiveness of the proposed method.


title: Hearing What You Cannot See: Acoustic Vehicle Detection Around Corners

abstract: This work proposes to use passive acoustic perception as an additional sensing modality for intelligent vehicles. We demonstrate that approaching vehicles behind blind corners can be detected by sound before such vehicles enter in line-of-sight. We have equipped a research vehicle with a roof-mounted microphone array, and show on data collected with this sensor setup that wall reflections provide information on the presence and direction of occluded approaching vehicles. A novel method is presented to classify if and from what direction a vehicle is approaching before it is visible, using as input Direction-of-Arrival features that can be efficiently computed from the streaming microphone array data. Since the local geometry around the ego-vehicle affects the perceived patterns, we systematically study several environment types, and investigate generalization across these environments. With a static ego-vehicle, an accuracy of 0.92 is achieved on the hidden vehicle classification task. Compared to a state-of-the-art visual detector, Faster R-CNN, our pipeline achieves the same accuracy more than one second ahead, providing crucial reaction time for the situations we study. While the ego-vehicle is driving, we demonstrate positive results on acoustic detection, still achieving an accuracy of 0.84 within one environment type. We further study failure cases across environments to identify future research directions.


title: A scaling improved inner-outer decomposition of near-wall turbulent motions

abstract: Near-wall turbulent velocities in turbulent channel flows are decomposed into small-scale and large-scale components at $y^+<100$ by improving the predictive inner-outer model of Baars et al. [Phys. Rev. Fluids 1, 054406 (2016)], where $y^+$ is the viscous-normalized wall-normal height. The small-scale one is obtained by reducing the outer reference height (a parameter in the model) from the center of the logarithmic layer to $y^+=100$, which can fully remove outer influences. On the other hand, the large-scale one represents the near-wall footprints of outer energy-containing motions. We present plenty of evidences that demonstrate that the small-scale motions are Reynolds-number invariant with the viscous scaling, at friction Reynolds numbers between 1000 and 5200. At lower Reynolds numbers from 180 to 600, the small scales can not be scaled by the viscous units, and the vortical structures are progressively strengthened as Reynolds number increases, which is proposed as a possible mechanism responsible for the anomalous scaling behavior. Finally, it is found that a small-scale part of the outer large-scale footprint can be well scaled by the viscous units.


title: Using GSM SMS controller alarm Configurator to develop cost effective, intelligent fire safety system in a developing country

abstract: Electricity supply to facilities is essential, but can cause fires when care is not taken, and can destroy the properties within some few minutes. The need for fire protection is therefore essential. This paper address-es this pertinent issue facing Ghana as a developing country. At the fire outbreak, the designed system re-sponds to the smoke and cuts off the electricity supply. When a fire is detected, an alarm is turned on and Short Messaging Service (SMS) alert is sent to the owner. After ten to fifteen seconds, the system, resends SMS to the owner and the National Fire Service (NFS) safety officer giving the exact location. The call alert is sent to the owner of a facility and fire safety personnel respectively, when the system is not reset. A microcontroller serves as the command center of the cost effective, intelligent fire safety system. The circuit is designed and simulated using Proteus software, and programing was done using Global System for Mo-bile Communications (GSM) SMS Controller Alarm Configurator software. The prototype was constructed and tested in real-time. The proposed system is cost effective and will help policy makers, and save lives and property when implemented.


title: High-Energy Neutrinos and Gamma-Rays from Non-Relativistic Shock-Powered Transients

abstract: Shock interaction has been argued to play a role in powering a range of optical transients, including supernovae (particularly the superluminous class), classical novae, stellar mergers, tidal disruption events, and fast blue optical transients. These same shocks can accelerate relativistic ions, generating high-energy neutrino and gamma-ray emission via hadronic pion production. The recent discovery of time-correlated optical and gamma-ray emission in classical novae has revealed the important role of radiative shocks in powering these events, enabling an unprecedented view of the properties of ion acceleration, including its efficiency and energy spectrum, under similar physical conditions to shocks in extragalactic transients. Here we introduce a model for connecting the radiated optical fluence of non-relativistic transients to their maximal neutrino and gamma-ray fluence. We apply this technique to a wide range of extragalactic transient classes in order to place limits on their contributions to the cosmological high-energy gamma-ray and neutrino backgrounds. Based on a simple model for diffusive shock acceleration at radiative shocks, calibrated to novae, we demonstrate that several of the most luminous transients can accelerate protons up to energies $E_{\rm max} \gtrsim 10^{16}$ eV, sufficient to contribute to the IceCube astrophysical background. Furthermore, several of the considered sources$-$particularly hydrogen-poor supernovae$-$may serve as "hidden" gamma-ray sources due to the high gamma-ray opacity of their ejecta, evading constraints imposed by the non-blazar Fermi-LAT background. However, adopting an ion acceleration efficiency $\sim$ 0.3-1$\%$ motivated by nova observations, we find that currently known classes of non-relativistic, potentially shock-powered transients contribute at most a few percent of the total IceCube background.


title: Distribution-Free Models of Social Networks

abstract: The structure of large-scale social networks has predominantly been articulated using generative models, a form of average-case analysis. This chapter surveys recent proposals of more robust models of such networks. These models posit deterministic and empirically supported combinatorial structure rather than a specific probability distribution. We discuss the formal definitions of these models and how they relate to empirical observations in social networks, as well as the known structural and algorithmic results for the corresponding graph classes.


title: Nonlinear Optics of optomagnetics: Quantum and Classical Treatments

abstract: Optomagnetics emerges as a growing field of research cross-linking optics, magnetism and material science. Here, we provide a microscopic quantum mechanical and a macroscopic classical models to describe optomagnetic effects from nonlinear optics point of view. Our self-consistent quantum mechanical formulation considers all orders of perturbing field and results not only in finding generalized Pitaevskii's relationship, where photoinduced magnetization can be expanded in terms of light power, but also provide compact and analytical expressions for optical gyration vector coefficients. classical treatment is then developed based on the anharmonic Drude-Lorentz model showing that the photo-induced DC magnetization is proportional to odd harmonics of the light power. The difference in quantum and classical results are revealed and discussed. Having a pomp-probe setup in mind, we describe how a probe light signal can propagate down an optomagnetic medium, i.e. a medium that is magnetized by intense circularly-polarized pump light, via its permittivity tensor and find light propagation characteristics. Inverse Faraday and Cotton-Mouton Effects are discussed as a result of circular and linear birefringences and their Verdet constants have been analytically found.


title: On Hyperparameter Optimization of Machine Learning Algorithms: Theory and Practice

abstract: Machine learning algorithms have been used widely in various applications and areas. To fit a machine learning model into different problems, its hyper-parameters must be tuned. Selecting the best hyper-parameter configuration for machine learning models has a direct impact on the model's performance. It often requires deep knowledge of machine learning algorithms and appropriate hyper-parameter optimization techniques. Although several automatic optimization techniques exist, they have different strengths and drawbacks when applied to different types of problems. In this paper, optimizing the hyper-parameters of common machine learning models is studied. We introduce several state-of-the-art optimization techniques and discuss how to apply them to machine learning algorithms. Many available libraries and frameworks developed for hyper-parameter optimization problems are provided, and some open challenges of hyper-parameter optimization research are also discussed in this paper. Moreover, experiments are conducted on benchmark datasets to compare the performance of different optimization methods and provide practical examples of hyper-parameter optimization. This survey paper will help industrial users, data analysts, and researchers to better develop machine learning models by identifying the proper hyper-parameter configurations effectively.


title: Laser2Vec: Similarity-based Retrieval for Robotic Perception Data

abstract: As mobile robot capabilities improve and deployment times increase, tools to analyze the growing volume of data are becoming necessary. Current state-of-the-art logging, playback, and exploration systems are insufficient for practitioners seeking to discover systemic points of failure in robotic systems. This paper presents a suite of algorithms for similarity-based queries of robotic perception data and implements a system for storing 2D LiDAR data from many deployments cheaply and evaluating top-k queries for complete or partial scans efficiently. We generate compressed representations of laser scans via a convolutional variational autoencoder and store them in a database, where a light-weight dense network for distance function approximation is run at query time. Our query evaluator leverages the local continuity of the embedding space to generate evaluation orders that, in expectation, dominate full linear scans of the database. The accuracy, robustness, scalability, and efficiency of our system is tested on real-world data gathered from dozens of deployments and synthetic data generated by corrupting real data. We find our system accurately and efficiently identifies similar scans across a number of episodes where the robot encountered the same location, or similar indoor structures or objects.


title: Probing Reionization and Early Cosmic Enrichment with the MgII Forest

abstract: Because the same massive stars that reionized the intergalactic medium (IGM) inevitably exploded as supernovae that polluted the Universe with metals, the history of cosmic reionization and enrichment are intimately intertwined. While the overly sensitive Ly-alpha transition completely saturates in a neutral IGM, strong low-ionization metal lines like the MgII 2796,2804 doublet will give rise to a detectable `metal-line forest' if the metals produced during reionization (Z ~ 10^{-3}Z_sol) permeate the neutral IGM. We simulate the MgII forest for the first time by combining a large hydrodynamical simulation with a semi-numerical reionization topology, assuming a simple enrichment model where the IGM is uniformly suffused with metals. In contrast to the traditional approach of identifying discrete absorbers, we treat the absorption as a continuous random field and measure its two-point correlation function, leveraging techniques from precision cosmology. We show that a realistic mock dataset of 10 JWST spectra can simultaneously determine the Mg abundance, [Mg/H], with a 1sigma precision of 0.02 dex and measure the global neutral fraction <x_HI> to 5% for a Universe with <x_HI> = 0.74 and [Mg/H] = -3.7. Alternatively, if the IGM is pristine, a null-detection of the MgII forest would set a stringent upper limit on the IGM metallicity of [Mg/H] < -4.4 at 95% credibility, assuming <x_HI> > 0.5 from another probe. Concentrations of metals in the circumgalactic environs of galaxies can significantly contaminate the IGM signal, but we demonstrate how these discrete absorbers can be easily identified and masked such that their impact on the correlation function is negligible. The MgII forest thus has tremendous potential to precisely constrain the reionization and enrichment history of the Universe.


title: Descending the Stable Matching Lattice: How many Strategic Agents are required to turn Pessimality to Optimality?

abstract: The set of stable matchings induces a distributive lattice. The supremum of the stable matching lattice is the boy-optimal (girl-pessimal) stable matching and the infimum is the girl-optimal (boy-pessimal) stable matching. The classical boy-proposal deferred-acceptance algorithm returns the supremum of the lattice, that is, the boy-optimal stable matching. In this paper, we study the smallest group of girls, called the {\em minimum winning coalition of girls}, that can act strategically, but independently, to force the boy-proposal deferred-acceptance algorithm to output the girl-optimal stable matching. We characterize the minimum winning coalition in terms of stable matching rotations and show that its cardinality can take on any value between $0$ and $\left\lfloor \frac{n}{2}\right\rfloor$, for instances with $n$ boys and $n$ girls. Our main result is that, for the random matching model, the expected cardinality of the minimum winning coalition is $(\frac{1}{2}+o(1))\log{n}$. This resolves a conjecture of Kupfer \cite{Kup18}.


title: On multi-symmetric functions and transportation polytopes

abstract: We present a study of the transportation polytopes appearing in the product rule of elementary multi-symmetric functions introduced by F. Vaccarino.


title: Optical clocks based on molecular vibrations as probes of variation of the proton-to-electron mass ratio

abstract: Some new physics models of quantum gravity or dark matter predict drifts or oscillations of the fundamental constants. A relatively simple model relates molecular vibrations to the proton-to-electron mass ratio $\mu$. Many vibrational transitions are at optical frequencies with prospects for use as highly accurate optical clocks. We give a brief summary of new physics models that lead to changes in $\mu$ and the current limits on drifts and oscillation amplitudes. After an overview of laboratory procedures, we give examples of molecules with experiments currently in development or underway. These experiments' projected systematic and statistical uncertainties make them leading candidates in next-generation searches for time-variation of $\mu$.


title: From calibration to parameter learning: Harnessing the scaling effects of big data in geoscientific modeling

abstract: The behaviors and skills of models in many geosciences (e.g., hydrology and ecosystem sciences) strongly depend on spatially-varying parameters that need calibration. A well-calibrated model can reasonably propagate information from observations to unobserved variables via model physics, but traditional calibration is highly inefficient and results in non-unique solutions. Here we propose a novel differentiable parameter learning (dPL) framework that efficiently learns a global mapping between inputs (and optionally responses) and parameters. Crucially, dPL exhibits beneficial scaling curves not previously demonstrated to geoscientists: as training data increases, dPL achieves better performance, more physical coherence, and better generalizability (across space and uncalibrated variables), all with orders-of-magnitude lower computational cost. We demonstrate examples that learned from soil moisture and streamflow, where dPL drastically outperformed existing evolutionary and regionalization methods, or required only ~12.5% of the training data to achieve similar performance. The generic scheme promotes the integration of deep learning and process-based models, without mandating reimplementation.


title: Materials Graph Transformer predicts the outcomes of inorganic reactions with reliable uncertainties

abstract: A common bottleneck for materials discovery is synthesis. While recent methodological advances have resulted in major improvements in the ability to predicatively design novel materials, researchers often still rely on trial-and-error approaches for determining synthesis procedures. In this work, we develop a model that predicts the major product of solid-state reactions. The cardinal feature of this approach is the construction of fixed-length, learned representations of reactions. Precursors are represented as nodes on a `reaction graph', and message-passing operations between nodes are used to embody the interactions between precursors in the reaction mixture. Through an ablation study, it is shown that this framework not only outperforms less physically-motivated baseline methods but also more reliably assesses the uncertainty in its predictions.


title: Bose-Einstein condensate fluctuations versus an interparticle interaction

abstract: We calculate the Bose-Einstein condensate (BEC) occupation statistics vs. the interparticle interaction in a dilute gas with a nonuniform condensate in a box trap within the Bogoliubov approach. The results are compared against the previously found BEC-occupation statistics in (i) an ideal gas and (ii) a weakly interacting gas with a uniform condensate. In particular, we reveal and explicitly describe an appearance of a nontrivial transition from the ideal gas to the Thomas-Fermi regime. The results include finding the main regimes of the BEC statistics - the anomalous non-Gaussian thermally-dominated fluctuations and the Gaussian quantum-dominated fluctuations - as well as a crossover between them and their manifestations in a mesoscopic system. Remarkably, we show that the effect of the boundary conditions, imposed at the box trap, on the BEC fluctuations does not vanish in the thermodynamic limit of a macroscopic system even in the presence of the interparticle interactions. Finally, we discuss a challenging problem of an experimental verification of the theory of the BEC fluctuations addressing a much deeper level of the many-body statistical physics than usually studied quantities related to the mean condensate occupation.


title: Bounding non-rationality of divisors on 3-fold Fano fibrations

abstract: In this paper we investigate non-rationality of divisors on 3-fold log Fano fibrations $(X,B)\to Z$ under mild conditions. We show that if $D$ is a component of $B$ with coefficient $\ge t>0$ which is contracted to a point on $Z$, then $D$ is birational to $\mathbb{P}^1\times C$ where $C$ is a smooth projective curve with gonality bounded depending only on $t$. Moreover, if $t>\frac{1}{2}$, then genus of $C$ is bounded depending only on $t$.


title: Blending Controllers via Multi-Objective Bandits

abstract: Safety and performance are often two competing objectives in sequential decision-making problems. Existing performant controllers, such as controllers derived from reinforcement learning algorithms, often fall short of safety guarantees. On the contrary, controllers that guarantee safety, such as those derived from classical control theory, require restrictive assumptions and are often conservative in performance. Our goal is to blend a performant and a safe controller to generate a single controller that is safer than the performant and accumulates higher rewards than the safe controller. To this end, we propose a blending algorithm using the framework of contextual multi-armed multi-objective bandits. At each stage, the algorithm observes the environment's current context alongside an immediate reward and cost, which is the underlying safety measure. The algorithm then decides which controller to employ based on its observations. We demonstrate that the algorithm achieves sublinear Pareto regret, a performance measure that models coherence with an expert that always avoids picking the controller with both inferior safety and performance. We derive an upper bound on the loss in individual objectives, which imposes no additional computational complexity. We empirically demonstrate the algorithm's success in blending a safe and a performant controller in a safety-focused testbed, the Safety Gym environment. A statistical analysis of the blended controller's total reward and cost reflects two key takeaways: The blended controller shows a strict improvement in performance compared to the safe controller, and it is safer than the performant controller.


title: Marcus-Hush-Chidsey Kinetics at Electrode-Electrolyte Interfaces

abstract: Electrochemical kinetics at electrode-electrolyte interfaces limit performance of devices including fuel cells and batteries. While the importance of moving beyond Butler-Volmer kinetics and incorporating the effect of electronic density of states of the electrode have been recognized, a unified framework that incorporates these aspects directly into electrochemical performance models is still lacking. In this work, we explicitly account for the DFT-calculated density of states numerically in calculating electrochemical reaction rates for a variety of electrode-electrolyte interfaces. We first show the utility of this for two cases related to Li metal electrodeposition and stripping on a Li surface and a Cu surface (anode-free configuration). The deviation in reaction rates is minor for cases with flat densities of states such as Li, but is significant for Cu due to nondispersive d-bands creating large variation. Finally, we consider a semiconducting case of a solid-electrolyte interphase (SEI) consisting of LiF and Li$_2$CO$_3$ and note the importance of the Fermi level at the interface, pinned by the redox reaction occuring there. We identify the asymmetry in reaction rates as a function of discharge/charge naturally within this approach. The analysis code used in this work is available open-source on Github.


title: Unidentified Floating Object detection in maritime environment using dictionary learning

abstract: Maritime domain is one of the most challenging scenarios for object detection due to the complexity of the observed scene. In this article, we present a new approach to detect unidentified floating objects in the maritime environment. The proposed approach is capable of detecting floating objects without any prior knowledge of their visual appearance, shape or location. The input image from the video stream is denoised using a visual dictionary learned from a K-SVD algorithm. The denoised image is made of self-similar content. Later, we extract the residual image, which is the difference between the original image and the denoised (self-similar) image. Thus, the residual image contains noise and salient structures (objects). These salient structures can be extracted using an a contrario model. We demonstrate the capabilities of our algorithm by testing it on videos exhibiting varying maritime scenarios.


title: Eulerian dynamics in multi-dimensions with radial symmetry

abstract: We study the global wellposedness of pressure-less Eulerian dynamics in multi-dimensions, with radially symmetric data. Compared with the 1D system, a major difference in multi-dimensional Eulerian dynamics is the presence of the spectral gap, which is difficult to control in general. We propose a new pair of scalar quantities that provides a significant better control of the spectral gap. Two applications are presented. (i) the Euler-Poisson equations: we show a sharp threshold condition on initial data that distinguish global regularity and finite time blowup; (ii) the Euler-alignment equations: we show a large subcritical region of initial data that leads to global smooth solutions.


title: The Program with a Personality: Analysis of Elk Cloner, the First Personal Computer Virus

abstract: Although self-replicating programs and viruses have existed since the 1960s and 70s, Elk Cloner was the first virus to circulate among personal computers in the wild. Despite its historical significance, it received comparatively little attention when it first appeared in 1982. In this paper, we: present the first detailed examination of the operation and structure of Elk Cloner; discuss the effect of environmental characteristics on its virulence; and provide supporting evidence for several hypotheses about why its release was largely ignored in the early 1980s.


title: On precision holography for the circular Wilson loop in $AdS_{5}\times S^{5}$

abstract: The string theory calculation of the 1/2-BPS circular Wilson loop of N=4 SYM in the planar limit at next to leading order at strong coupling is revisited in the ratio of its semiclassical string partition function and the one dual to a latitude Wilson loop with trivial expectation value. After applying a conformal transformation from the disk to the cylinder, this problem can be approached by means of the Gel'fand-Yaglom formalism. Using results from the literature and the exclusion of zero modes from a modified Gel'fand-Yaglom formula, we obtain matching with the known field theory result. As seen in the phaseshift method computation, non-zero mode contributions cancel and the end result comes from the zero mode degeneracies of the latitude Wilson loop.


title: DART Mission Determination of Momentum Transfer: Model of Ejecta Plume Observations

abstract: The NASA Double Asteroid Redirection Test (DART) spacecraft will impact the secondary member of the [65803] Didymos binary in order to perform the first demonstration of asteroid deflection by kinetic impact. Determination of the momentum transfer to the target body from the kinetic impact is a primary planetary defense objective, using ground-based telescopic observations of the orbital period change of Didymos and imaging of the DART impact ejecta plume by the LICIACube cubesat, along with modeling and simulation of the DART impact. LICIACube, contributed by the Italian Space Agency, will perform a flyby of Didymos a few minutes after the DART impact, to resolve the ejecta plume spatial structure and to study the temporal evolution. LICIACube ejecta plume images will help determine the vector momentum transfer from the DART impact, by determining or constraining the direction and the magnitude of the momentum carried by ejecta. A model is developed for the impact ejecta plume optical depth, using a point source scaling model of the DART impact. The model is applied to expected LICIACube plume images and shows how plume images enable characterization of the ejecta mass versus velocity distribution. The ejecta plume structure, as it evolves over time, is determined by the amount of ejecta that has reached a given altitude at a given time. The evolution of the plume optical depth profiles determined from LICIACube images can distinguish between strength-controlled and gravity-controlled impacts, by distinguishing the respective mass versus velocity distributions. LICIACube plume images discriminate the differences in plume structure and evolution that result from different target physical properties, mainly strength and porosity, thereby allowing inference of these properties to improve the determination of momentum transfer.


title: Deformations of Higher-Page Analogues of $\partial\bar\partial$-Manifolds

abstract: We extend the notion of essential deformations from the case of the Iwasawa manifold, for which they were introduced recently by the first-named author, to the general case of page-$1$-$\partial\bar\partial$-manifolds that were jointly introduced very recently by all three authors. We go on to obtain an analogue of the unobstructedness theorem of Bogomolov, Tian and Todorov for page-$1$-$\partial\bar\partial$-manifolds. As applications of this discussion, we study the small deformations of certain Nakamura solvmanifolds and reinterpret the cases of the Iwasawa manifold and its $5$-dimensional analogue from this standpoint.


title: A fast adaptive algorithm for scattering from a two dimensional radially-symmetric potential

abstract: In the present paper we describe a simple black box algorithm for efficiently and accurately solving scattering problems related to the scattering of time-harmonic waves from radially-symmetric potentials in two dimensions. The method uses FFTs to convert the problem into a set of decoupled second-kind Fredholm integral equations for the Fourier coefficients of the scattered field. Each of these integral equations are solved using scattering matrices, which exploit certain low-rank properties of the integral operators associated with the integral equations. The performance of the algorithm is illustrated with several numerical examples including scattering from singular and discontinuous potentials. Finally, the above approach can be easily extended to time-dependent problems. After outlining the necessary modifications we show numerical experiments illustrating the performance of the algorithm in this setting.


title: The Abundance Pattern of $\alpha$ elements in the Triangulum-Andromeda Overdensity

abstract: The close relationship between the nature of the Triangulum-Andromeda (TriAnd) overdensity and the Galactic disk has become increasingly evident in recent years. However, the chemical pattern of this overdensity (R$_{GC}$ = 20 - 30 kpc) is unique and differs from what we know of the local disk. In this study, we analyze the chemical abundances of five $\alpha$ elements (Mg, O, Si, Ca, and Ti) in a sample of stars belonging to the TriAnd overdensity, including stars with [Fe/H] $<$ $-$1.2, to investigate the evolution of the $\alpha$ elements with metallicity. High-resolution spectra from Gemini North with GRACES were analyzed. Overall, the TriAnd population presents an $\alpha$-element pattern that differs from that of the local disk; the TriAnd stars fall in between the local disk and the dwarf galaxies in the [X/Fe] vs. [Fe/H] plane. The high [Mg/Fe] ratios obtained for the lower metallicity TriAnd stars may indicate a roughly parallel sequence to the Milky Way local disk at lower values of [Fe/H], revealing a 'knee' shifted towards lower metallicities for the TriAnd population. Similar behavior is also exhibited in the [Ca/Fe] and [Si/Fe] ratios. However, for O and Ti the behavior of the [X/Fe] ratios shows a slight decay with decreasing metallicity. Our results reinforce the TriAnd overdensity as a unique stellar population of the Milky Way, with an abundance pattern that is different from all stellar populations studied to date. The complete understanding of the complex TriAnd population will require high-resolution spectroscopic observations of a larger sample of TriAnd stars.


title: On asymptotic expansions for the fractional infinity Laplacian

abstract: We propose two asymptotic expansions of two interrelated integral-type averages, in the context of the fractional $\infty$-Laplacian $\Delta_\infty^s$ for $s\in (\frac{1}{2},1)$. This operator has been introduced and first studied in [Bjorland, C., Caffarelli, L. and Figalli, A., \textsl{Nonlocal Tug-of-War and the inifnity fractional Laplacian}, Comm. Pure Appl. Math., \textbf{65}, pp. 337--380, (2012)]. Our expansions are parametrised by the radius of the removed singularity $\epsilon$, and allow for the identification of $\Delta_\infty^s\phi(x)$ as the $\epsilon^{2s}$-order coefficient of the deviation of the $\epsilon$-average from the value $\phi(x)$, in the limit $\epsilon\to 0+$. The averages are well posed for functions $\phi$ that are only Borel regular and bounded.


title: Additive interaction modelling using I-priors

abstract: Additive regression models with interactions are widely studied in the literature, using methods such as splines or Gaussian process regression. However, these methods can pose challenges for estimation and model selection, due to the presence of many smoothing parameters and the lack of suitable criteria. We propose to address these challenges by extending the I-prior methodology (Bergsma, 2020) to multiple covariates, which may be multidimensional. The I-prior methodology has some advantages over other methods, such as Gaussian process regression and Tikhonov regularization, both theoretically and practically. In particular, the I-prior is a proper prior, is based on minimal assumptions, yields an admissible posterior mean, and estimation of the scale (or smoothing) parameters can be done using an EM algorithm with simple E and M steps. Moreover, we introduce a parsimonious specification of models with interactions, which has two benefits: (i) it reduces the number of scale parameters and thus facilitates the estimation of models with interactions, and (ii) it enables straightforward model selection (among models with different interactions) based on the marginal likelihood.


title: The Saturn Ring Skimmer Mission Concept: The next step to explore Saturn's rings, atmosphere, interior, and inner magnetosphere

abstract: The innovative Saturn Ring Skimmer mission concept enables a wide range of investigations that address fundamental questions about Saturn and its rings, as well as giant planets and astrophysical disk systems in general. This mission would provide new insights into the dynamical processes that operate in astrophysical disk systems by observing individual particles in Saturn's rings for the first time. The Ring Skimmer would also constrain the origin, history, and fate of Saturn's rings by determining their compositional evolution and material transport rates. In addition, the Ring Skimmer would reveal how the rings, magnetosphere, and planet operate as an inter-connected system by making direct measurements of the ring's atmosphere, Saturn's inner magnetosphere and the material owing from the rings into the planet. At the same time, this mission would clarify the dynamical processes operating in the planet's visible atmosphere and deep interior by making extensive high-resolution observations of cloud features and repeated measurements of the planet's extremely dynamic gravitational field. Given the scientific potential of this basic mission concept, we advocate that it be studied in depth as a potential option for the New Frontiers program.


title: Decomposition of the Uniform Projection of the Weil Character

abstract: In this paper we obtain a decomposition formula of the uniform projection of the Weil character of a finite reductive dual pair consisting of a symplectic group and an even orthogonal group. This is the first and major step to give an explicit description of the Howe correspondence on unipotent characters confirming a conjecture by Aubert-Michel-Rouquier.


title: Instrument variable detection with graph learning : an application to high dimensional GIS-census data for house pricing

abstract: Endogeneity bias and instrument variable validation have always been important topics in statistics and econometrics. In the era of big data, such issues typically combine with dimensionality issues and, hence, require even more attention. In this paper, we merge two well-known tools from machine learning and biostatistics---variable selection algorithms and probablistic graphs---to estimate house prices and the corresponding causal structure using 2010 data on Sydney. The estimation uses a 200-gigabyte ultrahigh dimensional database consisting of local school data, GIS information, census data, house characteristics and other socio-economic records. Using "big data", we show that it is possible to perform a data-driven instrument selection efficiently and purge out the invalid instruments. Our approach improves the sparsity of variable selection, stability and robustness in the presence of high dimensionality, complicated causal structures and the consequent multicollinearity, and recovers a sparse and intuitive causal structure. The approach also reveals an efficiency and effectiveness in endogeneity detection, instrument validation, weak instrument pruning and the selection of valid instruments. From the perspective of machine learning, the estimation results both align with and confirms the facts of Sydney house market, the classical economic theories and the previous findings of simultaneous equations modeling. Moreover, the estimation results are consistent with and supported by classical econometric tools such as two-stage least square regression and different instrument tests. All the code may be found at \url{https://github.com/isaac2math/solar_graph_learning}.


title: Effect of Orientation of Cation CH3NH3PbI3 on Ambipolarity, Open-Circuit Voltage, and Excitons Lifetime

abstract: The effect of the orientation of the cation MA (i.e., CH3NH3) on the physical properties of MAPbI3 are investigated using the density functional theory. We report that the Fermi energy level exhibits a large variation with MA orientation, and it is this Fermi level variation that makes this material unique and so different from the non-organic perovskites. The Fermi level variation with orientation of MA is proposed to be responsible for the experimentally observed intrinsic open-circuit voltage V_OC, ambipolarity, and long lifetime of the excitons. Based on our results, ferroelectric domains in the low-temperature orthorhombic phase or clusters of MA molecules not rotating in unison (as a consequence of e.g. impurities and/or defects) in the higher-temperature tetragonal and cubic phases are proposed to be responsible for ambipolarity. This is because any given two adjacent domains or clusters with different MA s average orientations may have very different Fermi levels, causing them to behave as an effective diode. Also, the measured significant lifetime of excitons in MAPbI3 finds a natural explanation using this effective diode. Thus, our work provides a unified manner for explaining at least these three important properties in terms of the rotation-induced Fermi level variations. We used the density functional theory to get our results.


title: Deformed Quantum Mechanics and the Landau Problem

abstract: A deformation of the Landau problem based on a modification of Fock algebra is considered. Systems with Hamiltonians f(H) where H is the Landau Hamiltonian in the lowest level are discussed. The case $f(H) = {\alpha} H + b H^2$ is studied and it is shown that in this particular example, parameters of the problem can be fixed by using the quadratic Zeeman effect data and the Breit- Rabi formula. The proposed approach allows to solve exactly Landau-like families of problems not previously discussed in the literature.


title: On the Cohen-Macaulay property of the Rees algebra of the module of differentials

abstract: Let $R$ be an algebra essentially of finite type over a field $k$ and let $\Omega_k(R)$ be its module of K\"ahler differentials over $k$. If $R$ is a homogeneous complete intersection and $\mathrm{char}(k)=0$, we prove that $\Omega_k(R)$ is of linear type whenever its Rees algebra is Cohen-Macaulay and locally at every homogeneous prime $\mathfrak{p}$ the embedding dimension of $R_{\mathfrak{p}}$ is at most twice its dimension.


title: Mass-concentration of low-regularity blow-up solutions to the focusing 2D modified Zakharov-Kuznetsov equation

abstract: We consider the focusing modified Zakharov-Kuznetsov (mZK) equation in two space dimensions. We prove that solutions which blow up in finite time in the $H^1(\R^{2})$ norm have the property that they concentrate a non-trivial portion of their mass (more precisely, at least the amount equal to the mass of the ground state) at blow-up time. For finite-time blow-up solutions in the $H^s(\R^2)$ norm for $\frac{17}{18} < s < 1$, we prove a slightly weaker result. Moreover, we prove that the stronger concentration result can be extended to the range $ \frac{17}{18} < s \le 1$ under an additional assumption on the upper bound of the blow-up rate of the solution. The main tools used here are the $I$-method and a profile decomposition theorem for a bounded family of $H^1(\R^{2})$ functions.


title: Reconstruction of Protein Structures from Single-Molecule Time Series

abstract: Single-molecule experimental techniques track the real-time dynamics of molecules by recording a small number of experimental observables. Following these observables provides a coarse-grained, low-dimensional representation of the conformational dynamics but does not furnish an atomistic representation of the instantaneous molecular structure. Takens' Delay Embedding Theorem asserts that, under quite general conditions, these low-dimensional time series can contain sufficient information to reconstruct the full molecular configuration of the system up to an a priori unknown transformation. By combining Takens' Theorem with tools from statistical thermodynamics, manifold learning, artificial neural networks, and rigid graph theory, we establish an approach Single-molecule TAkens Reconstruction (STAR) to learn this transformation and reconstruct molecular configurations from time series in experimentally-measurable observables such as intramolecular distances accessible to single molecule F\"orster resonance energy transfer. We demonstrate the approach in applications to molecular dynamics simulations of a C24H50 polymer chain and the artificial mini-protein Chignolin. The trained models reconstruct molecular configurations from synthetic time series data in the head-to-tail molecular distances with atomistic root mean squared deviation accuracies better than 0.2 nm. This work demonstrates that it is possible to accurately reconstruct protein structures from time series in experimentally-measurable observables and establishes the theoretical and algorithmic foundations to do so in applications to real experimental data.


title: Alerts in High-resolution TEM characterization of perovskite material

abstract: High-resolution TEM (HRTEM) is a powerful tool for structure characterization. However, methylammonium lead iodide (MAPbI3) perovskite is highly sensitive to electron beams and easily decompose into lead iodide (PbI2). Universal misidentifications that PbI2 is incorrectly labeled as perovskite are widely exist in HRTEM characterization, which would negatively affect the development of perovskite research field. Here misidentifications in MAPbI3 perovskite calibration are summarized, classified and corrected based on corresponding electron diffraction (ED) simulations. Corresponding crystallographic parameters of intrinsic tetragonal MAPbI3 and the confusable hexagonal PbI2 are also presented clearly. Finally, the method of proper phase identification and some ways to control the radiation damage in HRTEM are provided. This work paves the way to avoid misleadings in HRTEM characterization of perovskite and other electron beam-sensitive materials in the future.