title: Phase-field modeling of biomineralization in mollusks and corals: Microstructure vs. formation mechanism

abstract: While biological crystallization processes have been studied on the microscale extensively, models addressing the mesoscale aspects of such phenomena are rare. In this work, we investigate whether the phase-field theory developed in materials science for describing complex polycrystalline structures on the mesoscale can be meaningfully adapted to model crystallization in biological systems. We demonstrate the abilities of the phase-field technique by modeling a range of microstructures observed in mollusk shells and coral skeletons, including granular, prismatic, sheet/columnar nacre, and sprinkled spherulitic structures. We also compare two possible micromechanisms of calcification: the classical route via ion-by-ion addition from a fluid state and a non-classical route, crystallization of an amorphous precursor deposited at the solidification front. We show that with appropriate choice of the model parameters microstructures similar to those found in biomineralized systems can be obtained along both routes, though the timescale of the non-classical route appears to be more realistic. The resemblance of the simulated and natural biominerals suggests that, underneath the immense biological complexity observed in living organisms, the underlying design principles for biological structures may be understood with simple math, and simulated by phase-field theory.


title: Single trial ERP amplitudes reveal the time course of acquiring representations of novel faces in individual participants

abstract: The neural correlates of face individuation - the acquisition of memory representations for novel faces - have been studied only in coarse detail and disregarding individual differences between learners. In their seminal study, (Tanaka, Curran, Porterfield, & Collins, 2006) required the identification of a particular novel face across 70 trials and found that the N250 component in the ERP became more negative from the first to the second half of the experiment, where it reached a similar amplitude as a well-known face. We were unable to directly replicate this finding in our study when we used the original split of trials. However, when we applied a different split of trials we observed very similar changes in N250 amplitude. Then, we developed and applied a new two-step explorative-confirmative non-parametric method based on permutation testing to determine the time course of face individuation in individual participants based on single-trial N250 amplitudes. We show that the assumption of a steep initial increase of N250 amplitude across multiple presentations of the target face, followed by a plateau, yields plausible results in fitting linear trends for most participants. The transition point from initial acquisition to the plateau phase differed strongly between participants and tended to be earlier when performance in target face recognition was better. Hence, face individuation may be accounted for by a biphasic process of early, fast acquisition, followed by a slower, asymptotic consolidation or maintenance phase. The current approach might be fruitfully applied to further investigations into face individuation and their neural correlates


title: Circular polarization in radio pulsar PSR B1451-68: coherent mode transitions and intrabeam interference

abstract: The radio emission of pulsar B1451-68 contains two polarization modes of similar strength, which produce two clear orthogonal polarization angle tracks. When viewed on a Poincare sphere, the emission is composed of two flux patches that rotate meridionally as function of pulse longitude and pass through the Stokes V poles, which results in transitions between orthogonal polarization modes (OPMs). Moreover, the ratio of power in the patches is inversed once within the profile window. It is shown that the meridional circularization is caused by a coherent OPM transition (COMT) produced by a varying mode ratio at a fixed quarter-wave phase lag. The COMTs may be ubiquitous and difficult to detect in radio pulsar data, because they can leave no trace in polarized fractions and they are described by equation similar to the rotating vector model. The circularization, which coincides with flux minima at lower frequency, requires that profile components are formed by radiation with an oscillation phase which increases with longitude in steps of 90 degrees per component. The properties can be understood as an interference pattern involving two pairs of linear orthogonal modes (or two nonorthogonal elliptic waves). The frequency-dependent coherent superposition of coplanar oscillations can produce the minima in the pulse profile, and thereby the illusion of components as separate entities. The orthogonally polarized signal which is left after such negative interference explains the enhancement of polarization degree which is commonly observed in the minima between profile components.


title: Flexo-Elastic Control Factors of Domain Morphology in Core-Shell Ferroelectric Nanoparticles: Soft and Rigid Shells

abstract: Within the framework of the Landau-Ginzburg-Devonshire approach we explore the impact of elastic anisotropy, electrostriction, flexoelectric couplings, and mismatch strain on the domain structure morphology in ferroelectric core-shell nanoparticles of spherical shape. We perform finite element modelling (FEM) for multiaxial ferroelectric nanoparticle cores covered with an elastically-isotropic soft or elastically-anisotropic rigid paraelectric shell, with and without mismatch strains. In the case of a core covered with a soft shell, the FEM results show that at room temperature a single polarization vortex with a dipolar kernel can be stable if the electrostriction coupling is weak. With increasing anisotropic electrostriction coupling, the vortex disappears and is replaced by complex flux-closure structures, which are formed in the equatorial plane and transform into an elongated vortex with a central 180-degree domain wall near the core poles. This complex domain morphology develops in the core due to the anisotropic electrostriction, and the flexoelectric coupling leads to an additional curvature and twist of the polarization isosurfaces. In contrast to this, FEM performed for a core covered with a rigid shell shows that, at room temperature, the anisotropic elastic properties of the shell can stabilize vortex-like structures with three flux-closure domains, which gradually cross in the equatorial plane of the core and transform into 120-degree domains near the core poles. The flexoelectric coupling leads to a noticeable curling of the flux-closure domain walls. A mismatch strain compensates the curling of the flux-closure domains in the core confined by the elastically-anisotropic rigid shell. Our analysis reveals different types of topological defects, namely Bloch point structures (BPS) and Ising lines, that form in a ferroelectric core covered with a soft or rigid shell.


title: Fluid-particle suspension by gas release from a granular bed

abstract: We have studied experimentally particle suspension when injecting a gas at the bottom of an immersed granular layer confined in a Hele-Shaw cell. This work focuses on the dynamics of particles slightly denser than the surrounding fluid. The gas, injected at constant flow-rate, rises through the granular bed then forms bubbles which entrain particles in the above liquid layer. The particles settle down on the edges of the cell, avalanche on the crater formed at the granular bed free surface, and are further entrained by the continuous bubbling at the center. We report the existence of a stationary state, resulting from the competition between particle entrainment and sedimentation. The average solid fraction in the suspension is derived from a simple measurement of the granular bed apparent area. A phenomenological model based on the balance between particle lift by bubbles at the center of the cell and their settling on its sides demonstrates that most of the particles entrained by bubbles come from a global recirculation of the suspension.


title: ML-based Flood Forecasting: Advances in Scale, Accuracy and Reach

abstract: Floods are among the most common and deadly natural disasters in the world, and flood warning systems have been shown to be effective in reducing harm. Yet the majority of the world's vulnerable population does not have access to reliable and actionable warning systems, due to core challenges in scalability, computational costs, and data availability. In this paper we present two components of flood forecasting systems which were developed over the past year, providing access to these critical systems to 75 million people who didn't have this access before.


title: Dynamics-based peptide-MHC binding optimization by a convolutional variational autoencoder: a use-case model for CASTELO

abstract: An unsolved challenge in the development of antigen specific immunotherapies is determining the optimal antigens to target. Comprehension of antigen-MHC binding is paramount towards achieving this goal. Here, we present CASTELO, a combined machine learning-molecular dynamics (ML-MD) approach to design novel antigens of increased MHC binding affinity for a Type 1 diabetes (T1D)-implicated system. We build upon a small molecule lead optimization algorithm by training a convolutional variational autoencoder (CVAE) on MD trajectories of 48 different systems across 4 antigens and 4 HLA serotypes. We develop several new machine learning metrics including a structure-based anchor residue classification model as well as cluster comparison scores. ML-MD predictions agree well with experimental binding results and free energy perturbation-predicted binding affinities. Moreover, ML-MD metrics are independent of traditional MD stability metrics such as contact area and RMSF, which do not reflect binding affinity data. Our work supports the role of structure-based deep learning techniques in antigen specific immunotherapy design.


title: Modified Dorfman procedure for pool tests with dilution -- COVID-19 case study

abstract: The outbreak of the global COVID-19 pandemic results in unprecedented demand for fast and efficient testing of large numbers of patients for the presence of SARS-CoV-2 coronavirus. Beside technical improvements of the cost and speed of individual tests, pool testing may be used to improve efficiency and throughput of a population test. Dorfman pool testing procedure is one of the best known and studied methods of this kind. This procedure is, however, based on unrealistic assumptions that the pool test has perfect sensitivity and the only objective is to minimize the number of tests, and is not well adapted to the case of imperfect pool tests. We propose and analyze a simple modification of this procedure in which test of a pool with negative result is independently repeated up to several times. The proposed procedure is evaluated in a computational study using recent data about dilution effect for SARS-CoV-2 PCR tests, showing that the proposed approach significantly reduces the number of false negatives with a relatively small increase of the number of tests, especially for small prevalence rates. For example, for prevalence rate 0.001 the number of tests could be reduced to 22.1% of individual tests, increasing the expected number of false negatives by no more than 1%, and to 16.8% of individual tests increasing the expected number of false negatives by no more than 10%. At the same time, a similar reduction of the expected number of tests in the standard Dorfman procedure would yield 675% and 821% increase of the expected number of false negatives, respectively. This makes the proposed procedure an interesting choice for screening tests in the case of diseases like COVID-19.


title: Single-photon distributed free-space spectroscopy

abstract: Spectroscopy is a well-established nonintrusive tool that has played an important role in identifying substances and quantifying their compositions, from quantum descriptions to chemical and biomedical diagnostics. Challenges exist in accurate measurements in dynamic environments, especially for understanding chemical reactions in arbitrary free-space. We develop a distributed free-space spectroscopy realized by a comb-referenced frequency-scanning single-photon lidar, providing multidimensional (time-range-spectrum) remote sensing. A continuous field experiment over 72 hours is deployed to obtain the spectra of multiple molecules (CO2 and HDO) in free-space over 6 km, with a range resolution of 60 m and a time resolution of 10 min over a spectrum span of 30 GHz. The CO2 and HDO concentrations are retrieved from the spectra acquired. This distributed free-space spectroscopy holds much promise for increasing knowledge of atmospheric environments and chemistry research, especially for complex molecular spectra evolution in any location over large areas.


title: Topological Learning for Brain Networks

abstract: This paper proposes a novel topological learning framework that integrates networks of different sizes and topology through persistent homology. Such challenging task is made possible through the introduction of a computationally efficient topological loss. The use of the proposed loss bypasses the intrinsic computational bottleneck associated with matching networks. We validate the method in extensive statistical simulations to assess its effectiveness when discriminating networks with different topology. The method is further demonstrated in a twin brain imaging study where we determine if brain networks are genetically heritable. The challenge here is due to the difficulty of overlaying the topologically different functional brain networks obtained from resting-state functional MRI onto the template structural brain network obtained through diffusion MRI.


title: Pulmonary drug delivery and retention: a computational study to identify plausible parameters based on a coupled airway-mucus flow model

abstract: Pulmonary drug delivery systems rely on inhalation of drug-laden aerosols produced from aerosol generators such as inhalers, nebulizers etc. On deposition, the drug molecules diffuse in the mucus layer and are also subjected to mucociliary advection which transports the drugs away from the initial deposition site. The availability of the drug at a particular region of the lung is, thus, determined by a balance between these two phenomena. A mathematical analysis of drug deposition and retention in the lungs is developed through a coupled mathematical model of aerosol transport in air as well as drug molecule transport in the mucus layer. The mathematical model is solved computationally to identify suitable conditions for the transport of drug-laden aerosols to the deep lungs. This study identifies the conditions conducive for delivering drugs to the deep lungs which is crucial for achieving systemic drug delivery. The effect of different parameters on drug retention is also characterized for various regions of the lungs, which is important in determining the availability of the inhaled drugs at a target location. Our analysis confirms that drug delivery efficacy remains highest for aerosols in the size range of 1-5 $\mu$m. Moreover, it is observed that amount of drugs deposited in the deep lung increases by a factor of 2 when the breathing time period is doubled, with respect to normal breathing, suggesting breath control as a means to increase the efficacy of drug delivery to the deep lung. A higher efficacy also reduces the drug load required to be inhaled to produce the same health effects and hence, can help in minimizing the side effects of a drug.


title: Phase-resolved ocean wave forecast with ensemble-based data assimilation

abstract: Through ensemble-based data assimilation (DA), we address one of the most notorious difficulties in phase-resolved ocean wave forecast, regarding the deviation of numerical solution from the true surface elevation due to the chaotic nature of and underrepresented physics in the nonlinear wave models. In particular, we develop a coupled approach of the high-order spectral (HOS) method with the ensemble Kalman filter (EnKF), through which the measurement data can be incorporated into the simulation to improve the forecast performance. A unique feature in this coupling is the mismatch between the predictable zone and measurement region, which is accounted for through a special algorithm to modify the analysis equation in EnKF. We test the performance of the new EnKF-HOS method using both synthetic data and real radar measurements. For both cases (though differing in details), it is shown that the new method achieves much higher accuracy than the HOS-only method, and can retain the phase information of an irregular wave field for arbitrarily long forecast time with sequentially assimilated data.


title: Meta-modeling strategy for data-driven forecasting

abstract: Accurately forecasting the weather is a key requirement for climate change mitigation. Data-driven methods offer the ability to make more accurate forecasts, but lack interpretability and can be expensive to train and deploy if models are not carefully developed. Here, we make use of two historical climate data sets and tools from machine learning, to accurately predict temperature fields. Furthermore, we are able to use low fidelity models that are cheap to train and evaluate, to selectively avoid expensive high fidelity function evaluations, as well as uncover seasonal variations in predictive power. This allows for an adaptive training strategy for computationally efficient geophysical emulation.


title: Using Machine Learning to Calibrate Storm-Scale Probabilistic Guidance of Severe Weather Hazards in the Warn-on-Forecast System

abstract: A primary goal of the National Oceanic and Atmospheric Administration (NOAA) Warn-on-Forecast (WoF) project is to provide rapidly updating probabilistic guidance to human forecasters for short-term (e.g., 0-3 h) severe weather forecasts. Maximizing the usefulness of probabilistic severe weather guidance from an ensemble of convection-allowing model forecasts requires calibration. In this study, we compare the skill of a simple method using updraft helicity against a series of machine learning (ML) algorithms for calibrating WoFS severe weather guidance. ML models are often used to calibrate severe weather guidance since they leverage multiple variables and discover useful patterns in complex datasets. \indent Our dataset includes WoF System (WoFS) ensemble forecasts available every 5 minutes out to 150 min of lead time from the 2017-2019 NOAA Hazardous Weather Testbed Spring Forecasting Experiments (81 dates). Using a novel ensemble storm track identification method, we extracted three sets of predictors from the WoFS forecasts: intra-storm state variables, near-storm environment variables, and morphological attributes of the ensemble storm tracks. We then trained random forests, gradient-boosted trees, and logistic regression algorithms to predict which WoFS 30-min ensemble storm tracks will correspond to a tornado, severe hail, and/or severe wind report. For the simple method, we extracted the ensemble probability of 2-5 km updraft helicity (UH) exceeding a threshold (tuned per severe weather hazard) from each ensemble storm track. The three ML algorithms discriminated well for all three hazards and produced more reliable probabilities than the UH-based predictions. Overall, the results suggest that ML-based calibrations of dynamical ensemble output can improve short term, storm-scale severe weather probabilistic guidance


title: Estimate of the location of the neutron drip line for calcium isotopes from an exact Hamiltonian with continuum pair correlations

abstract: The eastern region of the calcium isotope chain of the nuclei chart is, nowadays, of great activity. The experimental assessment of the limit of stability is of interest to confirm or improve microscopic theoretical models. The goal of this work is to provide the drip line of the calcium isotopes from the exact solution of the pairing Hamiltonian which incorporates explicitly the correlations with the continuum spectrum of energy. The modified Richardson equations, which include correlations with the continuum spectrum of energy modeled by the continuum single particle level density, is used to solve the many-body system. Three models are used, two isospin independent models with core 40Ca and 48Ca, and one isospin dependent model. One and two-neutron separation energies and occupation probabilities for bound and continuum states are calculated from the solution of the Richardson equations. The one particle drip line is found at the nucleus 57Ca, while the two neutron drip line is found at the nucleus 60Ca from the isospin independent model and at 66Ca from the isospin dependent one.


title: On the geometry of the kinematic space in special relativity

abstract: The classifying space of inertial reference frames in special relativity is naturally hyperbolic. There is a remarkable interplay between central elements of hyperbolic geometry and those of special relativity -- which, to a certain extent, have already been observed in the past -- that we present and further discuss in the paper. We aim at a geometrization of special relativity at the level of kinematic space by giving to physical concepts/phenomena purely geometric definitions/descriptions. In this way, the differences between special relativity and classical mechanics can be seen as a manifestation of the distinct geometric natures of their kinematic spaces.


title: Learning Disentangled Latent Factors from Paired Data in Cross-Modal Retrieval: An Implicit Identifiable VAE Approach

abstract: We deal with the problem of learning the underlying disentangled latent factors that are shared between the paired bi-modal data in cross-modal retrieval. Our assumption is that the data in both modalities are complex, structured, and high dimensional (e.g., image and text), for which the conventional deep auto-encoding latent variable models such as the Variational Autoencoder (VAE) often suffer from difficulty of accurate decoder training or realistic synthesis. A suboptimally trained decoder can potentially harm the model's capability of identifying the true factors. In this paper we propose a novel idea of the implicit decoder, which completely removes the ambient data decoding module from a latent variable model, via implicit encoder inversion that is achieved by Jacobian regularization of the low-dimensional embedding function. Motivated from the recent Identifiable VAE (IVAE) model, we modify it to incorporate the query modality data as conditioning auxiliary input, which allows us to prove that the true parameters of the model can be identified under some regularity conditions. Tested on various datasets where the true factors are fully/partially available, our model is shown to identify the factors accurately, significantly outperforming conventional encoder-decoder latent variable models. We also test our model on the Recipe1M, the large-scale food image/recipe dataset, where the learned factors by our approach highly coincide with the most pronounced food factors that are widely agreed on, including savoriness, wateriness, and greenness.


title: Prospect for application of mathematical models in combination cancer treatments

abstract: The long-term efficacy of targeted therapeutics for cancer treatment can be significantly limited by the type of therapy and development of drug resistance, inter alia. Experimental studies indicate that the factors enhancing acquisition of drug resistance in cancer cells include cell heterogeneity, drug target alteration, drug inactivation, DNA damage repair, drug efflux, cell death inhibition, as well as microenvironmental adaptations to targeted therapy, among others. Combination cancer therapies (CCTs) are employed to overcome these molecular and pathophysiological bottlenecks and improve the overall survival of cancer patients. CCTs often utilize multiple combinatorial modes of action and thus potentially constitute a promising approach to overcome drug resistance. Considering the colossal cost, human effort, time and ethical issues involved in clinical drug trials and basic medical research, mathematical modeling and analysis can potentially contribute immensely to the discovery of better cancer treatment regimens. In this article, we review mathematical models on CCTs developed thus far for cancer management. Open questions are highlighted and plausible combinations are discussed based on the level of toxicity, drug resistance, survival benefits, preclinical trials and other side effects.


title: Change of Cropping Patterns of Southeastern Anatolia, Turkey in 2019 and 2020

abstract: A short Technical Report of Multi-year Crop Mapping of Southeastern Anatolia Region of Turkey and its potential impacts on water use estimation.


title: A data-driven approach to the forecasting of ground-level ozone concentration

abstract: The ability to forecast the concentration of air pollutants in an urban region is crucial for decision-makers wishing to reduce the impact of pollution on public health through active measures (e.g. temporary traffic closures). In this study, we present a machine learning approach applied to the forecast of the day-ahead maximum value of the ozone concentration for several geographical locations in southern Switzerland. Due to the low density of measurement stations and to the complex orography of the use case terrain, we adopted feature selection methods instead of explicitly restricting relevant features to a neighbourhood of the prediction sites, as common in spatio-temporal forecasting methods. We then used Shapley values to assess the explainability of the learned models in terms of feature importance and feature interactions in relation to ozone predictions; our analysis suggests that the trained models effectively learned explanatory cross-dependencies among atmospheric variables. Finally, we show how weighting observations helps in increasing the accuracy of the forecasts for specific ranges of ozone's daily peak values.


title: Some comments on events associated with falling terrestrial rocks and iron from the sky

abstract: Some examples of the events associated with falls of rocks and iron terrestrial origin from the sky are considered (in scientific publications they are often called as meteor-wrongs or pseudo-meteorites). Their possible connections with other natural phenomena (like a whirlwind and a ball-lightning) are considered. Some compilation of info connected with the event in the US town of Elma in 2003 is presented. Also some possible parallels with the 1908 Tunguska event are given. The author hopes to draw attention to the need for more detailed research of such phenomena.


title: Hey Alexa what did I just type? Decoding smartphone sounds with a voice assistant

abstract: Voice assistants are now ubiquitous and listen in on our everyday lives. Ever since they became commercially available, privacy advocates worried that the data they collect can be abused: might private conversations be extracted by third parties? In this paper we show that privacy threats go beyond spoken conversations and include sensitive data typed on nearby smartphones. Using two different smartphones and a tablet we demonstrate that the attacker can extract PIN codes and text messages from recordings collected by a voice assistant located up to half a meter away. This shows that remote keyboard-inference attacks are not limited to physical keyboards but extend to virtual keyboards too. As our homes become full of always-on microphones, we need to work through the implications.


title: Short term solar energy prediction by machine learning algorithms

abstract: Smooth power generation from solar stations demand accurate, reliable and efficient forecast of solar energy for optimal integration to cater market demand; however, the implicit instability of solar energy production may cause serious problems for the smooth power generation. We report daily prediction of solar energy by exploiting the strength of machine learning techniques to capture and analyze complicated behavior of enormous features effectively. For this purpose, dataset comprising of 98 solar stations has been taken from energy competition of American Meteorological Society (AMS) for predicting daily solar energy. Forecast models of base line regressors including linear, ridge, lasso, decision tree, random forest and artificial neural networks have been implemented on the AMS solar dataset. Grid size is converted into two sections: 16x9 and 10x4 to ascertain attributes contributing more towards the generated power from densely located stations on global ensemble forecast system (GEFS). To evaluate the models, statistical measures of prediction error in terms of RMSE, MAE and R2_score have been analyzed and compared with the existing techniques. It has been observed that improved accuracy is achieved through random forest and ridge regressor for both grid sizes in contrast to all other proposed methods. Stability and reliability of the proposed schemes are evaluated on a single solar station as well as on multiple independent runs.


title: Dynamic Weighted Matching with Heterogeneous Arrival and Departure Rates

abstract: We study a dynamic non-bipartite matching problem. There is a fixed set of agent types, and agents of a given type arrive and depart according to type-specific Poisson processes. Agent departures are not announced in advance. The value of a match is determined by the types of the matched agents. We present an online algorithm that is (1/8)-competitive with respect to the value of the optimal-in-hindsight policy, for arbitrary weighted graphs. Our algorithm treats agents heterogeneously, interpolating between immediate and delayed matching in order to thicken the market while still matching valuable agents opportunistically.


title: Synchronizing subgrid scale models of turbulence to data

abstract: Large Eddy Simulations of turbulent flows are powerful tools used in many engineering and geophysical settings. Choosing the right value of the free parameters for their subgrid scale models is a crucial task for which the current methods present several shortcomings. Using a technique called nudging we show that Large Eddy Simulations can synchronize to data coming from a high-resolution direct numerical simulation of homogeneous and isotropic turbulence. Furthermore, we found that the degree of synchronization is dependent on the value of the parameters of the subgrid scale models utilized, suggesting that nudging can be used as a way to select the best parameters for a model. For example, we show that for the Smagorinsky model synchronization is optimal when its constant takes the usual value of $0.16$. Analyzing synchronization dynamics puts the focus on reconstructing trajectories in phase space, contrary to traditional a posteriori tests of Large Eddy Simulations where the statistics of the flows are compared. These results open up the possibility of utilizing non-statistical analysis in a posteriori tests of Large Eddy Simulations.


title: Tunable Majorana corner modes in noncentrosymmetric superconductors: Tunneling spectroscopy and edge imperfections

abstract: Majorana corner modes appearing in two-dimensional second-order topological superconductors have great potential applications for fault-tolerant topological quantum computations. We demonstrate that in the presence of an in-plane magentic field two-dimensional ($s+p$)-wave superconductors host Majorana corner modes, whose location can be manipulated by the direction of the magnetic field. In addition, we discuss the effects of edge imperfections on the Majorana corner modes. We describe how different edge shapes and edge disorder affect the number and controllability of the Majorana corner modes, which is of relevance for the implementation of topological quantum computations. We also discuss tunneling spectroscopy in the presence of the Majorana corner modes, where a lead-wire is attached to the corner of the noncentrosymmetric superconductor. The zero-bias differential conductance shows a distinct periodicity with respect to the direction of the magnetic field, which demonstrates the excellent controllability of the Majorana corner modes in this setup. Our results lay down the theoretical groundwork for observing and tuning Majoran corner modes in experiments on ($s+p$)-wave superconductors.


title: Phase of Nonlinear Systems

abstract: In this paper, we propose a definition of phase for a class of stable nonlinear systems called semi-sectorial systems, from an input-output perspective. The definition involves the Hilbert transform as a critical instrument to complexify real-valued signals since the notion of phase arises most naturally in the complex domain. The proposed nonlinear system phase, serving as a counterpart of $\mathcal{L}_2$-gain, quantifies the passivity and is highly related to the dissipativity. It also possesses a nice physical interpretation which quantifies the tradeoff between the real energy and reactive energy. A nonlinear small phase theorem is then established for feedback stability analysis of semi-sectorial systems. Additionally, its generalized version is proposed via the use of multipliers. These nonlinear small phase theorems generalize a version of the classical passivity theorem and a recently appeared linear time-invariant small phase theorem.


title: Unilateral and nonreciprocal transmission through bilinear spring systems

abstract: Longitudinal wave propagation is considered in a pair of waveguides connected by bilinear spring systems. The nature of the nonlinearity causes the compressive and tensile force-displacement relations of the bilinear spring to behave in a piece-wise linear manner, and all transmitted and reflected waves scale linearly with the incident wave amplitude. We first concentrate on a single bilinear spring connecting two waveguides. By controlling the bilinear stiffness parameters it is possible to convert a time harmonic incident wave into a transmitted wave of the same period but with particle displacement of a single sign, positive or negative, an effect we call unilateral transmission. Nonreciprocal wave phenomena are obtained by introducing spatial asymmetry. A simple combination of a single bilinear spring with a mass and a linear spring shows significant nonreciprocity with transmission relatively high in one direction and low in the opposite direction.


title: Propagation dynamics of abruptly autofocusing circular Airy-Gaussian vortex beams in the fractional Schr\"odinger equation

abstract: We introduce axisymmetric Airy-Gaussian vortex beams in a model of an optical system based on the (2+1)-dimensional fractional Schr\"odinger equation, characterized by its L\'evy index (LI). By means of numerical methods, we explore propagation dynamics of the beams with vorticities from 0 to 4. The propagation leads to abrupt autofocusing, followed by its reversal (rebound from the center). It is shown that LI, the relative width of the Airy and Gaussian factors, and the vorticity determine properties of the autofocusing dynamics, including the focusing distance, radius of the focal light spot, and peak intensity at the focus. A maximum of the peak intensity is attained at intermediate values of LI, close to LI=1.4 . Dynamics of the abrupt autofocusing of Airy-Gaussian beams carrying vortex pairs (split double vortices) is considered too.


title: The topological non-local braid-group concept of information processing in brain, the different role of the gray and white matter

abstract: The velocity of the action potential transduction along myelinated axons in the peripheral nervous system or in the white matter of brain and spinal cord reaches hundreds of meters per second to assure proper functioning of the body, which exceeds the ability of diffusive ion conduction. We propose the new model of the saltatory conduction based on a plasmon-polariton kinetics in myelinated axons, which excludes, however, the white matter form the information storage and its identification in the brain via e-m response. We propose a nonlocal topological approach to information processing in the cortex of brain in consistence with the ion electricity of the gray matter and a supplementary only communication role of the white matter.


title: Non-trivial Minkowski backgrounds in f(T) gravity

abstract: Boosted and rotated tetrad backgrounds for the Minkowski space are studied in $f(T)$ gravity. We perform Lorentzian perturbations at first order around non-trivial backgrounds and show that some Lorentz modes can exhibit non-trivial dynamics and can propagate in time. This remarkable feature gives evidence of additional mode(s) in the Lorentzian sector which have no precedent in Lorentz-violating modified gravities, however they can cast doubts onto even theoretical viability of these models.


title: Sigma Worksheet: Interactive Construction of OLAP Queries

abstract: The new generation of cloud data warehouses (CDWs) brings large amounts of data and compute power closer to users in enterprises. The ability to directly access the warehouse data, interactively analyze and explore it at scale can empower users to improve their decision making cycles. However, existing tools for analyzing data in CDWs are either limited in ad-hoc transformations or difficult to use for business users, the largest user segment in enterprises. Here we introduce Sigma Worksheet, a new interactive system that enables users to easily perform ad-hoc visual analysis of data in CDWs at scale. For this, Sigma Worksheet provides an accessible spreadsheet-like interface for data analysis through direct manipulation. Sigma Worksheet dynamically constructs matching SQL queries from user interactions on this familiar interface, building on the versatility and expressivity of SQL. Sigma Worksheet executes constructed queries directly on CDWs, leveraging the superior characteristics of the new generation CDWs, including scalability. To evaluate Sigma Worksheet, we first demonstrate its expressivity through two real life use cases, cohort analysis and sessionization. We then measure the performance of the Worksheet generated queries with a set of experiments using the TPC-H benchmark. Results show the performance of our compiled SQL queries is comparable to that of the reference queries of the benchmark. Finally, to assess the usefulness of Sigma Worksheet in deployment, we elicit feedback through a 100-person survey followed by a semi-structured interview study with 70 participants. We find that Sigma Worksheet is easier to use and learn, improving the productivity of users. Our findings also suggest Sigma Worksheet can further improve user experience by providing guidance to users at various steps of data analysis.


title: Data-driven optimal control of a SEIR model for COVID-19

abstract: We present a data-driven optimal control approach which integrates the reported partial data with the epidemic dynamics for COVID-19. We use a basic Susceptible-Exposed-Infectious-Recovered (SEIR) model, the model parameters are time-varying and learned from the data. This approach serves to forecast the evolution of the outbreak over a relatively short time period and provide scheduled controls of the epidemic. We provide efficient numerical algorithms based on a generalized Pontryagin Maximum Principle associated with the optimal control theory. Numerical experiments demonstrate the effective performance of the proposed model and its numerical approximations.


title: Lithium-cesium slow beam from a two-dimensional magneto-optical trap

abstract: We present the creation of a lithium-cesium slow beam using a two-dimensional magneto-optical trap. The two-species atomic beam is directed to load a three-dimensional magneto-optical trap in ultrahigh vacuum. We achieve a loading rate of 1.3$\times$10$^{7}$ lithium atoms/s with the lithium oven temperature at 370~$^{\circ}$C and 2.2$\times$10$^{7}$ cesium atoms/s with the Cs oven temperature at 20~$^{\circ}$C. The maximum numbers of lithium and cesium atoms in the trap are 1.0$\times$10$^{9}$ and 1.4$\times$10$^{8}$, respectively. Our results show that the simple and compact two-dimensional magneto-optical trap is suitable for producing an atomic beam of two species that have a large mass ratio and very different volatilities.


title: Gapless to gapless phase transitions in quantum spin chains

abstract: We investigate spin chains with bilinear-biquadratic (BLBQ) spin interactions as a function of an applied magnetic field $h$. At the Uimin-Lai-Sutherland (ULS) critical point we find a gapless to gapless transition revealed by the dynamical structure factor $S(q,\omega)$ as a function of $h$. At $h=0$, the envelope of the lowest energy excitations goes soft at {\it two points} $q_1=2\pi/3$ and $q_2=4\pi/3$, dubbed the phase A. With increasing field, the spectral peaks at each of the gapless points bifurcate, making in total {\it four} soft modes, and combine to form a new set of excitations that soften at a {\it single} point $q=\pi$ at $h_{c1}\approx 0.94$. Beyond $h_{c1}$ the system enters another gapless B-phase until the transition at $h_{c2}=4$ to the fully polarized phase. We compare the ULS model results with those for the AKLT model as a representative of gapped Haldane phase. We explain the mechanism of the gapless to gapless transition in the ULS model using its conserved charges and a spinon band picture. We also discuss the universality of central charges of the BLBQ family of models subjected to a magnetic field.


title: Spontaneous symmetry breaking in free theories with boundary potentials

abstract: Patterns of symmetry breaking induced by potentials at the boundary of free $O(N)$ models in $d=3- \epsilon$ dimensions are studied. We show that the spontaneous symmetry breaking in these theories leads to a boundary RG flow ending with $N - 1$ Neumann modes in the IR. The possibility of fluctuation-induced symmetry breaking is examined and we derive a general formula for computing one-loop effective potentials at the boundary. Using the $\epsilon-$expansion we test these ideas in an $O(N)\oplus O(N)$ model with boundary interactions. We determine the RG flow diagram of this theory and find that it has an IR-stable critical point satisfying conformal boundary conditions. The leading correction to the effective potential is computed and we argue the existence of a phase boundary separating the region flowing to the symmetric fixed point from the region flowing to a symmetry-broken phase with a combination of Neumann and Dirchlet boundary conditions.


title: Curing the unphysical behaviour of NLO quarkonium production at the LHC and its relevance to constrain the gluon PDF at low scales

abstract: We address the unphysical energy dependence of quarkonium-hadroproduction cross sections at Next-to-Leading Order (NLO) in alpha_S which we attribute to an over-subtraction in the factorisation of the collinear singularities inside the PDFs in the MSbar scheme. Such over- or under-subtractions have a limited phenomenological relevance in most of the scattering processes in particle physics. On the contrary, it is particularly harmful for P_T-integrated charmonium hadroproduction which renders a wide class of NLO results essentially unusable. Indeed, in such processes, alphaS is not so small, the PDFs are not evolved much and can be rather flat for the corresponding momentum fractions and, finally, some process-dependent NLO pieces are either too small or too large. We propose a scale-fixing criterion which avoids such an over-subtraction. We demonstrate its efficiency for eta(c,b) but also for a fictitious light elementary scalar boson. Having provided stable NLO predictions for eta(c,b) P_T-integrated cross sections, sigma^NLO(eta(Q)), and discussed the options to study eta(b) hadroproduction, we argue that their measurement at the LHC can help better determine the gluon PDF at low scales and tell whether the local minimum in conventional NLO gluon PDFs around x=0.001 at scales below 2 GeV is physical or not.


title: Modular Hamiltonians for the massless Dirac field in the presence of a boundary

abstract: We study the modular Hamiltonians of an interval for the massless Dirac fermion on the half-line. The most general boundary conditions ensuring the global energy conservation lead to consider two phases, where either the vector or the axial symmetry is preserved. In these two phases we derive the corresponding modular Hamiltonian in explicit form. Its density involves a bi-local term localised in two points of the interval, one conjugate to the other. The associated modular flows are also established. Depending on the phase, they mix fields with different chirality or charge that follow different modular trajectories. Accordingly, the modular flow preserves either the vector or the axial symmetry. We compute the two-point correlation functions along the modular flow and show that they satisfy the Kubo-Martin-Schwinger condition in both phases. The entanglement entropies are also derived.

categories: hep-th cond-mat.stat-mech math-ph math.MP


title: Lower Bounds for Semialgebraic Range Searching and Stabbing Problems

abstract: In the semialgebraic range searching problem, we are to preprocess $n$ points in $\mathbb{R}^d$ s.t. for any query range from a family of constant complexity semialgebraic sets, all the points intersecting the range can be reported or counted efficiently. When the ranges are composed of simplices, the problem can be solved using $S(n)$ space and with $Q(n)$ query time with $S(n)Q^d(n) = \tilde{O}(n^d)$ and this trade-off is almost tight. Consequently, there exists low space structures that use $\tilde{O}(n)$ space with $O(n^{1-1/d})$ query time and fast query structures that use $O(n^d)$ space with $O(\log^{d} n)$ query time. However, for the general semialgebraic ranges, only low space solutions are known, but the best solutions match the same trade-off curve as the simplex queries. It has been conjectured that the same could be done for the fast query case but this open problem has stayed unresolved. Here, we disprove this conjecture. We give the first nontrivial lower bounds for semilagebraic range searching and related problems. We show that any data structure for reporting the points between two concentric circles with $Q(n)$ query time must use $S(n)=\Omega(n^{3-o(1)}/Q(n)^5)$ space, meaning, for $Q(n)=O(\log^{O(1)}n)$, $\Omega(n^{3-o(1)})$ space must be used. We also study the problem of reporting the points between two polynomials of form $Y=\sum_{i=0}^\Delta a_i X^i$ where $a_0, \cdots, a_\Delta$ are given at the query time. We show $S(n)=\Omega(n^{\Delta+1-o(1)}/Q(n)^{\Delta^2+\Delta})$. So for $Q(n)=O(\log^{O(1)}n)$, we must use $\Omega(n^{\Delta+1-o(1)})$ space. For the dual semialgebraic stabbing problems, we show that in linear space, any data structure that solves 2D ring stabbing must use $\Omega(n^{2/3})$ query time. This almost matches the linearization upper bound. For general semialgebraic slab stabbing problems, again, we show an almost tight lower bounds.


title: Ultrafast transfer and transient entrapment of photoexcited Mg electron in Mg@C60

abstract: Electron relaxation is studied in endofullerene Mg@C60, after an initial localized photoexcitation in Mg, by nonadiabtic molecular dynamics simulations. To ensure reliability, two methods are used: i) an independent particle approach with a DFT description of the ground state and ii) HF ground state with many-body effects for the excited state dynamics. Both methods exhibit similar relaxation times leading to an ultrafast decay and charge transfer from Mg to C60 within tens of femtoseconds. Method (i) further elicits a robust transient-trap of the transferred electron that can delay the electron-hole recombination. Results shall motivate experiments to probe these ultrafast processes by two-photon transient absorption spectroscopy in gas phase, in solution, or as thin films.


title: Mixed precision matrix interpolative decompositions for model reduction

abstract: Renewed interest in mixed-precision algorithms has emerged due to growing data capacity and bandwidth concerns, as well as the advancement of GPUs, which enable significant speedup for low precision arithmetic. In light of this, we propose a mixed-precision algorithm to generate a double-precision accurate matrix interpolative decomposition approximation under a given set of criteria. Though low precision arithmetic suffers from quicker accumulation of round-off error, for many data-rich applications we nevertheless attain viable approximation accuracy, as the error incurred using low precision arithmetic is dominated by the error inherent to low-rank approximation. We then conduct several simulated numerical tests to demonstrate the efficacy of the algorithms and the corresponding error estimates. Finally, we present the application of our algorithms to a problem in model reduction for particle-laden turbulent flow.


title: A note on the limit of Orlicz norms

abstract: We generalize the well-known inequality that the limit of the $L^p$ norm of a function as $p\rightarrow\infty$ is the $L^\infty$ norm to the scale of Orlicz spaces.


title: Mutual Information Constraints for Monte-Carlo Objectives

abstract: A common failure mode of density models trained as variational autoencoders is to model the data without relying on their latent variables, rendering these variables useless. Two contributing factors, the underspecification of the model and the looseness of the variational lower bound, have been studied separately in the literature. We weave these two strands of research together, specifically the tighter bounds of Monte-Carlo objectives and constraints on the mutual information between the observable and the latent variables. Estimating the mutual information as the average Kullback-Leibler divergence between the easily available variational posterior $q(z|x)$ and the prior does not work with Monte-Carlo objectives because $q(z|x)$ is no longer a direct approximation to the model's true posterior $p(z|x)$. Hence, we construct estimators of the Kullback-Leibler divergence of the true posterior from the prior by recycling samples used in the objective, with which we train models of continuous and discrete latents at much improved rate-distortion and no posterior collapse. While alleviated, the tradeoff between modelling the data and using the latents still remains, and we urge for evaluating inference methods across a range of mutual information values.


title: Study of Vortex Dynamics and Phase Transitions in Superconducting Thin Films

abstract: The work reported in my doctoral thesis is an experimental study of vortex dynamics and phase transitions in thin films of type II superconductors using scanning tunneling spectroscopy, low frequency ac susceptibility measurements and complimentary transport measurements. Chapter I and II cover basics of superconductivity and methodologies used in the thesis. Chapter III discusses the effect of periodic pinning centers on the geometry and dynamics of vortex lattice in NbN thin films. Consequent study of dynamic transition of vortex Mott-like to vortex metal-like state is described. In Chapter IV, effect of strong disorder on vortex lattice in NbN thin films is studied. Here we show magnetic field induced granularity gives rise to pseudogap phase which is utilized to explain superconductor to insulator-like transition in stronger disorder. Chapter V contains the study of 2-dimensional vortex lattice melting in a-MoGe thin films. We discuss the observation of hexatic vortex fluid phase and the BKTHNY two-step melting in this scenario. Effect of sample thickness on this phases is also described. In Chapter VI, we study the possibility of quantum fluctuation of vortices in weakly pinned a-MoGe thin films. Effect of pinning on this fluctuation and a possibility of a quantum to thermal crossover is also explored.


title: Angle-Independent Plasmonic Substrates for Multi-Mode Vibrational Strong Coupling with Molecular Thin Films

abstract: Vibrational strong coupling of molecules to optical cavities based on plasmonic resonances has been explored recently, because plasmonic near-fields can provide strong coupling in sub-diffraction limited volumes. Such field localization maximizes coupling strength, which is crucial for modifying the vibrational response of molecules and, thereby, manipulating chemical reactions. Here, we demonstrate an angle-independent plasmonic nanodisk substrate that overcomes limitations of traditional Fabry-Perot optical cavities, because the design can strongly couple with all molecules on the surface of the substrate regardless of molecular orientation. We demonstrate that the plasmonic substrate provides strong coupling with the C=O vibrational stretch of deposited films of PMMA. We also show that the large linewidths of the plasmon resonance allows for simultaneous strong coupling to two, orthogonal water symmetric and asymmetric vibrational modes in a thin film of copper sulfate monohydrate deposited on the substrate surface. A three-coupled-oscillator model is developed to analyze the coupling strength of the plasmon resonance with these two water modes. With precise control over the nanodisk diameter, the plasmon resonance is tuned systematically through the modes, with the coupling strength to both modes varying as a function of the plasmon frequency, and with strong coupling to both modes achieved simultaneously for a range of diameters. This work may aid further studies into manipulation of the ground-state chemical landscape of molecules by perturbing multiple vibrational modes simultaneously and increasing the coupling strength in sub-diffraction limited volumes.


title: Learning with Knowledge of Structure: A Neural Network-Based Approach for MIMO-OFDM Detection

abstract: In this paper, we explore neural network-based strategies for performing symbol detection in a MIMO-OFDM system. Building on a reservoir computing (RC)-based approach towards symbol detection, we introduce a symmetric and decomposed binary decision neural network to take advantage of the structure knowledge inherent in the MIMO-OFDM system. To be specific, the binary decision neural network is added in the frequency domain utilizing the knowledge of the constellation. We show that the introduced symmetric neural network can decompose the original $M$-ary detection problem into a series of binary classification tasks, thus significantly reducing the neural network detector complexity while offering good generalization performance with limited training overhead. Numerical evaluations demonstrate that the introduced hybrid RC-binary decision detection framework performs close to maximum likelihood model-based symbol detection methods in terms of symbol error rate in the low SNR regime with imperfect channel state information (CSI).


title: Complex powers of the wave operator and the spectral action on Lorentzian scattering spaces

abstract: We consider perturbations of Minkowski space as well as more general spacetimes on which the wave operator $\square_g$ is known to be essentially self-adjoint. We define complex powers $(\square_g-i\varepsilon)^{-\alpha}$ by functional calculus, and show that the trace density exists as a meromorphic function of $\alpha$. We relate its poles to geometric quantities, in particular to the scalar curvature. The results allow us to formulate a spectral action principle which serves as a simple Lorentzian model for the bosonic part of the Chamseddine-Connes action. Our proof combines microlocal resolvent estimates, including radial propagation estimates, with uniform estimates for the Hadamard parametrix. The arguments operate in Lorentzian signature directly and do not rely on a transition from the Euclidean setting. The results hold also true in the case of ultrastatic spacetimes.


title: On the Lagrangian formulation of the double copy to cubic order

abstract: We investigate the Lagrangian formulation of the double-copy correspondence between gauge theories and gravity, up to the cubic order. Building on the definition of the double-copy field as a convolution of two vectors, we obtain free gravitational Lagrangians as products of two Yang-Mills Lagrangians, in a form amenable to be easily extended to the massive case. We discuss the origin of these results from tensionless strings and show the existence of gauge fixings that mix the two spin-one sectors and lead to an alternative, especially simple, version of the free Lagrangian. We then construct cubic vertices for the full double-copy multiplet, comprising a graviton, a two-form and a scalar particle, by means of the Noether procedure. Both at the free and at the cubic level the result gets uniquely fixed only upon imposing, on top of gauge invariance, a left-right Lorentz symmetry ruling contraction of indices among double-copy fields. Whereas the outcome nicely matches the cubic interactions of $\mathcal{N}=0$ supergravity, including the gauge-invariant coupling between the scalar particle and the two-form, such a twofold Lorentz symmetry seems to conflict with the perturbative reconstruction of space-time geometry.


title: Debiasing Evaluations That are Biased by Evaluations

abstract: It is common to evaluate a set of items by soliciting people to rate them. For example, universities ask students to rate the teaching quality of their instructors, and conference organizers ask authors of submissions to evaluate the quality of the reviews. However, in these applications, students often give a higher rating to a course if they receive higher grades in a course, and authors often give a higher rating to the reviews if their papers are accepted to the conference. In this work, we call these external factors the "outcome" experienced by people, and consider the problem of mitigating these outcome-induced biases in the given ratings when some information about the outcome is available. We formulate the information about the outcome as a known partial ordering on the bias. We propose a debiasing method by solving a regularized optimization problem under this ordering constraint, and also provide a carefully designed cross-validation method that adaptively chooses the appropriate amount of regularization. We provide theoretical guarantees on the performance of our algorithm, as well as experimental evaluations.


title: On symbolic factors of $\mathcal{S}$-adic subshifts of finite topological rank

abstract: This paper studies several aspects of symbolic factors of $\mathcal{S}$-adic subshifts of finite alphabet rank. First, we address a problem raised in [DDPM20] about the topological rank of symbolic factors of $\mathcal{S}$-adic subshifts and prove that this rank is at most the one of the extension system, improving results from [E20] and [GH2020]. As a consequence of our methods, we prove that finite topological rank systems are coalescent. Second, we investigate the structure of fibers $\pi^{-1}(y)$ of factor maps $\pi\colon(X,T)\to(Y,T)$ between minimal $\mathcal{S}$-adic subshifts of finite alphabet rank and show that they have the same finite cardinality for all $y$ in a residual subset of $Y$. Finally, we prove that the number of symbolic factors (up to conjugacy) of a fixed subshift of finite topological rank is finite, thus extending Durand's similar theorem on linearly recurrent subshifts.


title: On decay-surge population models

abstract: We consider continuous space-time decay-surge population models which are semi- stochastic processes for which deterministically declining populations, bound to fade away, are rein- vigorated at random times by bursts or surges of random sizes. In a particular separable framework (in a sense made precise below) we provide explicit formulae for the scale (or harmonic) function and the speed measure of the process. The behavior of the scale function at infinity allows to formulate conditions under which such processes either explode or are transient at infinity, or Harris recurrent. A description of the structures of both the discrete-time embedded chain and extreme record chain of such continuous-time processes is supplied.


title: Programmability of covariant quantum channels

abstract: A programmable quantum processor uses the states of a program register to specify one element of a set of quantum channels which is applied to an input register. It is well-known that such a device is impossible with a finite-dimensional program register for any set that contains infinitely many unitary quantum channels (Nielsen and Chuang's No-Programming Theorem), meaning that a universal programmable quantum processor does not exist. The situation changes if the system has symmetries. Indeed, here we consider group-covariant channels. If the group acts irreducibly on the channel input, these channels can be implemented exactly by a programmable quantum processor with finite program dimension (via teleportation simulation, which uses the Choi-Jamiolkowski state of the channel as a program). Moreover, by leveraging the representation theory of the symmetry group action, we show how to remove redundancy in the program and prove that the resulting program register has minimum Hilbert space dimension. Furthermore, we provide upper and lower bounds on the program register dimension of a processor implementing all group-covariant channels approximately.


title: Simulating Surface Wave Dynamics with Convolutional Networks

abstract: We investigate the performance of fully convolutional networks to simulate the motion and interaction of surface waves in open and closed complex geometries. We focus on a U-Net architecture and analyse how well it generalises to geometric configurations not seen during training. We demonstrate that a modified U-Net architecture is capable of accurately predicting the height distribution of waves on a liquid surface within curved and multi-faceted open and closed geometries, when only simple box and right-angled corner geometries were seen during training. We also consider a separate and independent 3D CNN for performing time-interpolation on the predictions produced by our U-Net. This allows generating simulations with a smaller time-step size than the one the U-Net has been trained for.


title: Gull's theorem revisited

abstract: Steve Gull, in unpublished work available on his Cambridge University homepage, has outlined a proof of Bell's theorem using Fourier theory. Gull's philosophy is that Bell's theorem (or perhaps a key lemma in its proof) can be seen as a no-go theorem for a project in distributed computing with classical, not quantum, computers. We present his argument, correcting misprints and filling gaps. In his argument, there were two completely separated computers in the network. We need three in order to fill all the gaps in his proof: a third computer supplies a stream of random numbers to the two computers representing the two measurement stations in Bell's work. One could also imagine that computer replaced by a cloned, virtual computer, generating the same pseudo-random numbers within each of Alice and Bob's computers. Either way, we need an assumption of the presence of shared i.i.d. randomness in the form of a synchronised sequence of realisations of i.i.d. hidden variables underlying the otherwise deterministic physics of the sequence of trials. Gull's proof then just needs a third step: rewriting an expectation as the expectation of a conditional expectation given the hidden variables.


title: Fully Convolutional Networks for Panoptic Segmentation

abstract: In this paper, we present a conceptually simple, strong, and efficient framework for panoptic segmentation, called Panoptic FCN. Our approach aims to represent and predict foreground things and background stuff in a unified fully convolutional pipeline. In particular, Panoptic FCN encodes each object instance or stuff category into a specific kernel weight with the proposed kernel generator and produces the prediction by convolving the high-resolution feature directly. With this approach, instance-aware and semantically consistent properties for things and stuff can be respectively satisfied in a simple generate-kernel-then-segment workflow. Without extra boxes for localization or instance separation, the proposed approach outperforms previous box-based and -free models with high efficiency on COCO, Cityscapes, and Mapillary Vistas datasets with single scale input. Our code is made publicly available at https://github.com/Jia-Research-Lab/PanopticFCN.


title: Targeted particle delivery via vortex ring reconnection

abstract: A conceptual model for targeted particle delivery is proposed using controlled vortex ring reconnection. Entrained particles can be efficiently transported within the core of the vortex ring which is propelled via self-induction. A pair of these particle-transporting vortices travelling in the streamwise direction along parallel trajectories will mutually interact resulting in vortex reconnection. The reconnection causes a topological change to the vortex rings which is accompanied by a rapid repulsion in a perpendicular plane to the direction of travel; effectively, transporting the particles to the desired delivery site on the sidewalls. This work explains the physics of the process and the considerations for targeted delivery.


title: Measuring and assessing economic uncertainty

abstract: This paper evaluates the dynamic response of economic activity to shocks in uncertainty as percieved by agents.The study focuses on the comparison between the perception of economic uncertainty by manufacturers and consumers.Since uncertainty is not directly observable, we approximate it using the geometric discrepancy indicator of Claveria et al.(2019).This approach allows us quantifying the proportion of disagreement in business and consumer expectations of eleven European countries and the Euro Area.First, we compute three independent indices of discrepancy corresponding to three dimensions of uncertainty (economic, inflation and employment) and we average them to obtain aggregate disagreement measures for businesses and for consumers.Next, we use a bivariate Bayesian vector autoregressive framework to estimate the impulse response functions to innovations in disagreement in every country.We find that the effect on economic activity of shocks to the perception of uncertainty differ markedly between manufacturers and consumers.On the one hand, shocks to consumer discrepancy tend to be of greater magnitude and duration than those to manufacturer discrepancy.On the other hand, innovations in disagreement between the two collectives have an opposite effect on economic activity:shocks to manufacturer discrepancy lead to a decrease in economic activity, as opposed to shocks to consumer discrepancy.This finding is of particular relevance to researchers when using cross-sectional dispersion of survey-based expectations, since the effect on economic growth of shocks to disagreement depend on the type of agent.


title: Magnetoelastic coupling and Gr\"uneisen scaling in NdB$_4$

abstract: We report high-resolution capacitance dilatometry studies on the uniaxial length changes in a NdB$_4$ single crystal. The evolution of magnetically ordered phases below $T_{\rm N}$= 17.2~K (commensurate antiferromagnetic phase, cAFM), $T_{\rm IT}$= 6.8~K (intermediate incommensurate phase, IT), and $T_{\rm LT}$= 4.8~K (low-temperature phase, LT) is associated with pronounced anomalies in the thermal expansion coefficients. The data imply significant magneto-elastic coupling and evidence of a structural phase transition at $T_{\rm LT}$ . While both cAFM and LT favor structural anisotropy $\delta$ between in-plane and out-of-plane length changes, it competes with the IT-type of order, i.e., $\delta$ is suppressed in that phase. Notably, finite anisotropy well above $T_{\rm N}$ indicates short-range correlations which are, however, of neither cAFM, IT, nor LT-type. Gr\"uneisen analysis of the ratio of thermal expansion coefficient and specific heat enables the derivation of uniaxial as well as hydrostatic pressure dependencies. While $\alpha$/$c_{\rm p}$ evidences a single dominant energy scale in LT, our data imply precursory fluctuations of a competing phase in IT and cAFM, respectively. Our results suggest the presence of orbital degrees of freedom competing with cAFM and successive evolution of a magnetically and orbitally ordered ground state.


title: Assessing and Accelerating Coverage in Deep Reinforcement Learning

abstract: Current deep reinforcement learning (DRL) algorithms utilize randomness in simulation environments to assume complete coverage in the state space. However, particularly in high dimensions, relying on randomness may lead to gaps in coverage of the trained DRL neural network model, which in turn may lead to drastic and often fatal real-world situations. To the best of the author's knowledge, the assessment of coverage for DRL is lacking in current research literature. Therefore, in this paper, a novel measure, Approximate Pseudo-Coverage (APC), is proposed for assessing the coverage in DRL applications. We propose to calculate APC by projecting the high dimensional state space on to a lower dimensional manifold and quantifying the occupied space. Furthermore, we utilize an exploration-exploitation strategy for coverage maximization using Rapidly-Exploring Random Tree (RRT). The efficacy of the assessment and the acceleration of coverage is demonstrated on standard tasks such as Cartpole, highway-env.


title: Regular multidimensional stationary time series

abstract: The aim of this paper is to give a simpler, more usable sufficient condition to the regularity of generic weakly stationary time series. Also, this condition is used to show how regular processes satisfying these sufficient conditions can be approximated by a lower rank \emph{regular} process. The relevance of these issues is shown by the ever increasing presence of high-dimensional data in many fields lately, and because of this, low rank processes and low rank approximations are becoming more important. Moreover, regular processes are the ones which are completely influenced by random innovations, so they are primary targets both in the theory and applications.


title: RAFT-3D: Scene Flow using Rigid-Motion Embeddings

abstract: We address the problem of scene flow: given a pair of stereo or RGB-D video frames, estimate pixelwise 3D motion. We introduce RAFT-3D, a new deep architecture for scene flow. RAFT-3D is based on the RAFT model developed for optical flow but iteratively updates a dense field of pixelwise SE3 motion instead of 2D motion. A key innovation of RAFT-3D is rigid-motion embeddings, which represent a soft grouping of pixels into rigid objects. Integral to rigid-motion embeddings is Dense-SE3, a differentiable layer that enforces geometric consistency of the embeddings. Experiments show that RAFT-3D achieves state-of-the-art performance. On FlyingThings3D, under the two-view evaluation, we improved the best published accuracy (d < 0.05) from 34.3% to 83.7%. On KITTI, we achieve an error of 5.77, outperforming the best published method (6.31), despite using no object instance supervision. Code is available at https://github.com/princeton-vl/RAFT-3D.


title: Spatial Distribution of OVI Covering Fractions in the Simulated Circumgalactic Medium

abstract: We use adaptive mesh refinement cosmological simulations to study the spatial distribution and covering fraction of OVI absorption in the circumgalactic medium (CGM) as a function of projected virial radius and azimuthal angle. We compare these simulations to an observed sample of 53 galaxies from the Multiphase Galaxy Halos Survey. Using Mockspec, an absorption line analysis pipeline, we generate synthetic quasar absorption line observations of the simulated CGM. To best emulate observations, we studied the averaged properties of 15,000 "mock samples" each of 53 sightlines having a distribution of $D/R_{vir}$ and sightline orientation statistically consistent with the observations. We find that the OVI covering fraction obtained for the simulated galaxies agrees well with the observed value for the inner halo ($D/R_{vir} \leq 0.375$) and is within $1.1\sigma$ in the outer halo ($D/R_{vir} > 0.75$), but is underproduced within $0.375 < D/R_{vir} \leq 0.75$. The observed bimodal distribution of OVI covering fraction with azimuthal angle, showing higher frequency of absorption along the projected major and minor axes of galaxies, is not reproduced in the simulations. Further analysis reveals the spatial-kinematic distribution of OVI absorbing gas is dominated by outflows in the inner halo mixed with a inflowing gas that originates from further out in the halo. Though the CGM of the individual simulated galaxies exhibit spatial structure, the flat azimuthal distribution occurs because the individual simulated galaxies do not develop a CGM structure that is universal from galaxy to galaxy.


title: Intrinsic analysis for dual word embedding space models

abstract: Recent word embeddings techniques represent words in a continuous vector space, moving away from the atomic and sparse representations of the past. Each such technique can further create multiple varieties of embeddings based on different settings of hyper-parameters like embedding dimension size, context window size and training method. One additional variety appears when we especially consider the Dual embedding space techniques which generate not one but two-word embeddings as output. This gives rise to an interesting question - "is there one or a combination of the two word embeddings variety, which works better for a specific task?". This paper tries to answer this question by considering all of these variations. Herein, we compare two classical embedding methods belonging to two different methodologies - Word2Vec from window-based and Glove from count-based. For an extensive evaluation after considering all variations, a total of 84 different models were compared against semantic, association and analogy evaluations tasks which are made up of 9 open-source linguistics datasets. The final Word2vec reports showcase the preference of non-default model for 2 out of 3 tasks. In case of Glove, non-default models outperform in all 3 evaluation tasks.


title: mlOSP: Towards a Unified Implementation of Regression Monte Carlo Algorithms

abstract: We introduce mlOSP, a computational template for Machine Learning for Optimal Stopping Problems. The template is implemented in the R statistical environment and publicly available via a GitHub repository. mlOSP presents a unified numerical implementation of Regression Monte Carlo (RMC) approaches to optimal stopping, providing a state-of-the-art, open-source, reproducible and transparent platform. Highlighting its modular nature, we present multiple novel variants of RMC algorithms, especially in terms of constructing simulation designs for training the regressors, as well as in terms of machine learning regression modules. Furthermore, mlOSP nests most of the existing RMC schemes, allowing for a consistent and verifiable benchmarking of extant algorithms. The article contains extensive R code snippets and figures, and serves as a vignette to the underlying software package.


title: Homological Dehn functions of groups of type $FP_2$

abstract: We prove foundational results for homological Dehn functions of groups of type $FP_2$ such as superadditivity and the invariance under quasi-isometry. We then study the homological Dehn functions of Leary's groups $G_L(S)$ providing methods to obtain uncountably many groups with a given homological Dehn function. This allows us to show that there exist groups of type $FP_2$ with quartic homological Dehn function and unsolvable word problem.


title: Large eddy simulation/probability density function simulations of the Cambridge turbulent stratified flame series under swirling conditions

abstract: The LES/PDF methodology is applied to the Cambridge/Sandia turbulent stratified flame series under swirling conditions. The methane/air chemistry is represented by a 16-species reduced ARM1 mechanism, and the in situ adaptive tabulation method is adopted to accelerate the chemistry calculation. Effects of differential diffusion and heat loss through the bluffbody surface are taken into account. The simulations are conducted for premixed (SwB3), moderately and highly stratified (SwB7 and SwB11, respectively) cases under swirling conditions. The results from LES/PDF simulations are compared with experimental measurements. The computed mean and r.m.s. profiles of velocity, temperature, equivalence ratio and mass fractions of species are in very good agreement with the measurements for all three conditions. Scatter plots and conditional means of species mole fractions and temperature are compared with the experimental data, where overall good consistency with the measurements is achieved. The recirculation zones are five times longer than those obtained under non-swirling conditions. A low-equivalence-ratio high-temperature region near the bluffbody is not captured in the computation, which is attributed to the insuffcient entrainment of downstream mixtures into the recirculation zone. The parametric studies show that the differential diffusion has a negligible effect on the mean and r.m.s. results, whereas the heat loss has a considerable effect on the temperature and CO profiles close to the bluffbody.


title: Convergence and Sample Complexity of SGD in GANs

abstract: We provide theoretical convergence guarantees on training Generative Adversarial Networks (GANs) via SGD. We consider learning a target distribution modeled by a 1-layer Generator network with a non-linear activation function $\phi(\cdot)$ parametrized by a $d \times d$ weight matrix $\mathbf W_*$, i.e., $f_*(\mathbf x) = \phi(\mathbf W_* \mathbf x)$. Our main result is that by training the Generator together with a Discriminator according to the Stochastic Gradient Descent-Ascent iteration proposed by Goodfellow et al. yields a Generator distribution that approaches the target distribution of $f_*$. Specifically, we can learn the target distribution within total-variation distance $\epsilon$ using $\tilde O(d^2/\epsilon^2)$ samples which is (near-)information theoretically optimal. Our results apply to a broad class of non-linear activation functions $\phi$, including ReLUs and is enabled by a connection with truncated statistics and an appropriate design of the Discriminator network. Our approach relies on a bilevel optimization framework to show that vanilla SGDA works.


title: Multistage games and Bell scenarios with communication

abstract: Bell nonlocality is a cornerstone of quantum theory with applications in information processing ranging from cryptography to distributed computing and game theory. Indeed, it is known that Bell's theorem can be formally linked to Bayesian games, allowing the use of nonlocal correlations to advise players and thereby achieve new points of equilibrium that are unavailable classically. Here we generalize this link, proving the connection between multistage games of incomplete information with Bell scenarios involving the communication of measurement outcomes between the parties. We apply the general framework for a few cases of interest and analyze the equilibrium reached by quantum nonlocal correlations.


title: Spectral decomposition and decay to grossly determined solutions for a simplified BGK model

abstract: Extending work of Carty, we show that $H^1$ solutions of a simplified 1D BGK model decay exponentially in $L^2$ to a subclass of the class of grossly determined solutions as defined by Truesdell and Muncaster. In the process, we determine the spectrum and generalized eigenfunctions of the associated non-selfadjoint linearized operator and derive the associated generalized Fourier transform and Parseval's identity. Notably, our analysis makes use of rigged space techniques originating from quantum mechanics, as adapted by Ljance and others to the nonselfadjoint case.


title: Sol-gel synthesis of Li3VO4/C composites as anode materials for lithium-ion batteries

abstract: Li3VO4/C composites have been synthesized by a sol-gel method and post-annealing at 650 {\deg}C for 1 h in N2 flow using either tartaric acid, malic acid, or glucose as both chelating agents and carbon source. The presence of these organic additives crucially affects morphology and crystallite size of the final product. It is found that the electrochemical properties of Li3VO4/C as anode material for Li-ion batteries (LIBs) are influenced by the morphology, texture and carbon content of the material. When using carboxylic acids as carbon source composites with mesoporous structure and a high surface area are obtained that display an enhanced electrochemical activity. Initially, reversible capacity of about 400 mAh g-1 is obtained. In contrast, Li3VO4/C synthesized with glucose outperforms in terms of cycling stability. It exhibits a discharge capacity of 299 mAh g-1 after 100 cycles corresponding to an excellent capacity retention of 96 %. The favorable effect of carbon composites on the electrochemical performance of Li3VO4 is shown.


title: Infinite-Dimensional Programmable Quantum Processors

abstract: A universal programmable quantum processor uses program quantum states to apply an arbitrary quantum channel to an input state. We generalize the concept of a finite-dimensional programmable quantum processor to infinite dimension assuming an energy constraint on the input and output of the target quantum channels. By proving reductions to and from finite-dimensional processors, we obtain upper and lower bounds on the program dimension required to approximately implement energy-limited quantum channels. In particular, we consider the implementation of Gaussian channels. Due to their practical relevance, we investigate the resource requirements for gauge-covariant Gaussian channels. Additionally, we give upper and lower bounds on the program dimension of a processor implementing all Gaussian unitary channels. These lower bounds rely on a direct information-theoretic argument, based on the generalization from finite to infinite dimension of a certain replication lemma for unitaries.


title: Conservation of radial actions in time-dependent spherical potentials

abstract: In slowly evolving spherical potentials, $\Phi(r,t)$, radial actions are typically assumed to remain constant. Here, we construct dynamical invariants that allow us to derive the evolution of radial actions in spherical central potentials with an arbitrary time dependence. We show that to linear order, radial actions oscillate around a constant value with an amplitude $\Delta J_r \propto \dot{\Phi}/\Phi\,P(E,L)$. Using this result, we develop a diffusion theory that describes the evolution of the radial action distribution of ensembles of tracer particles orbiting in generic time-dependent spherical potentials. Tests against restricted $N$-body simulations in a varying Kepler potential indicate that our linear theory is accurate in regions of phase-space in which the diffusion coefficient $\tilde{D}(J_r) < 0.01\,J_r^2$. For illustration, we apply our theory to two astrophysical processes. We show that the median mass accretion rate of a Milky Way (MW) dark matter (DM) halo leads to slow global time-variation of the gravitational potential, in which the evolution of radial actions is linear (i.e. either adiabatic or diffusive) for $\sim 84$ per cent of the DM halo at redshift $z=0$. This fraction grows considerably with lookback time, suggesting that diffusion may be relevant to the modelling of several Gyr-old tidal streams in action-angle space. As a second application, we show that dynamical tracers in a self-interacting DM (SIDM) dwarf halo (with $\sigma/m_\chi = 1\,{\rm cm^2g^{-1}}$) have invariant radial actions during the formation of a cored density profile.


title: Searching, Sorting, and Cake Cutting in Rounds

abstract: We study searching and sorting in rounds motivated by a fair division question: given a cake cutting problem with $n$ players, compute a fair allocation in at most $k$ rounds of interaction with the players. Rounds interpolate between the simultaneous and the fully adaptive settings, also capturing parallel complexity. We find that proportional cake cutting in rounds is equivalent to sorting with rank queries in rounds. We design a protocol for proportional cake cutting in rounds, while lower bounds for sorting in rounds with rank queries were given by Alon and Azar. Inspired by the rank query model, we then consider two basic search problems: ordered and unordered search. In unordered search, we get an array $\vec{x}=(x_1, \ldots, x_n)$ and an element $z$ promised to be in $\vec{x}$. We have access to an oracle that receives queries of the form "Is $z$ at location $i$?" and answers "Yes" or "No". The goal is to find the location of $z$ with success probability at least $p$ in at most $k$ rounds of interaction with the oracle. We show the expected query complexity of randomized algorithms on a worst case input is $np\bigl(\frac{k+1}{2k}\bigr) \pm O(1)$, while that of deterministic algorithms on a worst case input distribution is $np \bigl(1 - \frac{k-1}{2k}p \bigr) \pm O(1)$. These bounds apply even to fully adaptive unordered search, where the ratio between the two complexities converges to $2-p$ as the size of the array grows. In ordered search, we get sorted array $\vec{x}=(x_1, \ldots, x_n)$ and element $z$ promised to be in $\vec{x}$. We have access to an oracle that gets comparison queries. Here we find that the expected query complexity of randomized algorithms on a worst case input and deterministic algorithms on a worst case input distribution is essentially the same: $p k \cdot n^{\frac{1}{k}} \pm O(1+pk)$.


title: GLEAN: Generative Latent Bank for Large-Factor Image Super-Resolution

abstract: We show that pre-trained Generative Adversarial Networks (GANs), e.g., StyleGAN, can be used as a latent bank to improve the restoration quality of large-factor image super-resolution (SR). While most existing SR approaches attempt to generate realistic textures through learning with adversarial loss, our method, Generative LatEnt bANk (GLEAN), goes beyond existing practices by directly leveraging rich and diverse priors encapsulated in a pre-trained GAN. But unlike prevalent GAN inversion methods that require expensive image-specific optimization at runtime, our approach only needs a single forward pass to generate the upscaled image. GLEAN can be easily incorporated in a simple encoder-bank-decoder architecture with multi-resolution skip connections. Switching the bank allows the method to deal with images from diverse categories, e.g., cat, building, human face, and car. Images upscaled by GLEAN show clear improvements in terms of fidelity and texture faithfulness in comparison to existing methods.


title: MYSTIKO : : Cloud-Mediated, Private, Federated Gradient Descent

abstract: Federated learning enables multiple, distributed participants (potentially on different clouds) to collaborate and train machine/deep learning models by sharing parameters/gradients. However, sharing gradients, instead of centralizing data, may not be as private as one would expect. Reverse engineering attacks on plaintext gradients have been demonstrated to be practically feasible. Existing solutions for differentially private federated learning, while promising, lead to less accurate models and require nontrivial hyperparameter tuning. In this paper, we examine the use of additive homomorphic encryption (specifically the Paillier cipher) to design secure federated gradient descent techniques that (i) do not require addition of statistical noise or hyperparameter tuning, (ii) does not alter the final accuracy or utility of the final model, (iii) ensure that the plaintext model parameters/gradients of a participant are never revealed to any other participant or third party coordinator involved in the federated learning job, (iv) minimize the trust placed in any third party coordinator and (v) are efficient, with minimal overhead, and cost effective.


title: A converse to Lieb-Robinson bounds in one dimension using index theory

abstract: Unitary dynamics with a strict causal cone (or "light cone") have been studied extensively, under the name of quantum cellular automata (QCAs). In particular, QCAs in one dimension have been completely classified by an index theory. Physical systems often exhibit only approximate causal cones; Hamiltonian evolutions on the lattice satisfy Lieb-Robinson bounds rather than strict locality. This motivates us to study approximately locality preserving unitaries (ALPUs). We show that the index theory is robust and completely extends to one-dimensional ALPUs. As a consequence, we achieve a converse to the Lieb-Robinson bounds: any ALPU of index zero can be exactly generated by some time-dependent, quasi-local Hamiltonian in constant time. For the special case of finite chains with open boundaries, any unitary satisfying the Lieb-Robinson bound may be generated by such a Hamiltonian. We also discuss some results on the stability of operator algebras which may be of independent interest.

categories: quant-ph cond-mat.str-el math-ph math.MP


title: Spectral Analysis of Word Statistics

abstract: Given a random text over a finite alphabet, we study the frequencies at which fixed-length words occur as subsequences. As the data size grows, the joint distribution of word counts exhibits a rich asymptotic structure. We investigate all linear combinations of subword statistics, and fully characterize their different orders of magnitude using diverse algebraic tools. Moreover, we establish the spectral decomposition of the space of word statistics of each order. We provide explicit formulas for the eigenvectors and eigenvalues of the covariance matrix of the multivariate distribution of these statistics. Our techniques include and elaborate on a set of algebraic word operators, recently studied and employed by Dieker and Saliola (Adv Math, 2018). Subword counts find applications in Combinatorics, Statistics, and Computer Science. We revisit special cases from the combinatorial literature, such as intransitive dice, random core partitions, and questions on random walk. Our structural approach describes in a unified framework several classical statistical tests. We propose further potential applications to data analysis and machine learning.


title: Adaptive Neural Architectures for Recommender Systems

abstract: Deep learning has proved an effective means to capture the non-linear associations of user preferences. However, the main drawback of existing deep learning architectures is that they follow a fixed recommendation strategy, ignoring users' real time-feedback. Recent advances of deep reinforcement strategies showed that recommendation policies can be continuously updated while users interact with the system. In doing so, we can learn the optimal policy that fits to users' preferences over the recommendation sessions. The main drawback of deep reinforcement strategies is that are based on predefined and fixed neural architectures. To shed light on how to handle this issue, in this study we first present deep reinforcement learning strategies for recommendation and discuss the main limitations due to the fixed neural architectures. Then, we detail how recent advances on progressive neural architectures are used for consecutive tasks in other research domains. Finally, we present the key challenges to fill the gap between deep reinforcement learning and adaptive neural architectures. We provide guidelines for searching for the best neural architecture based on each user feedback via reinforcement learning, while considering the prediction performance on real-time recommendations and the model complexity.


title: A Framework and Dataset for Abstract Art Generation via CalligraphyGAN

abstract: With the advancement of deep learning, artificial intelligence (AI) has made many breakthroughs in recent years and achieved superhuman performance in various tasks such as object detection, reading comprehension, and video games. Generative Modeling, such as various Generative Adversarial Networks (GAN) models, has been applied to generate paintings and music. Research in Natural Language Processing (NLP) also had a leap forward in 2018 since the release of the pre-trained contextual neural language models such as BERT and recently released GPT3. Despite the exciting AI applications aforementioned, AI is still significantly lagging behind humans in creativity, which is often considered the ultimate moonshot for AI. Our work is inspired by Chinese calligraphy, which is a unique form of visual art where the character itself is an aesthetic painting. We also draw inspirations from paintings of the Abstract Expressionist movement in the 1940s and 1950s, such as the work by American painter Franz Kline. In this paper, we present a creative framework based on Conditional Generative Adversarial Networks and Contextual Neural Language Model to generate abstract artworks that have intrinsic meaning and aesthetic value, which is different from the existing work, such as image captioning and text-to-image generation, where the texts are the descriptions of the images. In addition, we have publicly released a Chinese calligraphy image dataset and demonstrate our framework using a prototype system and a user study.


title: Double machine learning for sample selection models

abstract: This paper considers the evaluation of discretely distributed treatments when outcomes are only observed for a subpopulation due to sample selection or outcome attrition. For identification, we combine a selection-on-observables assumption for treatment assignment with either selection-on-observables or instrumental variable assumptions concerning the outcome attrition/sample selection process. We also consider dynamic confounding, meaning that covariates that jointly affect sample selection and the outcome may (at least partly) be influenced by the treatment. To control in a data-driven way for a potentially high dimensional set of pre- and/or post-treatment covariates, we adapt the double machine learning framework for treatment evaluation to sample selection problems. We make use of (a) Neyman-orthogonal, doubly robust, and efficient score functions, which imply the robustness of treatment effect estimation to moderate regularization biases in the machine learning-based estimation of the outcome, treatment, or sample selection models and (b) sample splitting (or cross-fitting) to prevent overfitting bias. We demonstrate that the proposed estimators are asymptotically normal and root-n consistent under specific regularity conditions concerning the machine learners and investigate their finite sample properties in a simulation study. We also apply our proposed methodology to the Job Corps data for evaluating the effect of training on hourly wages which are only observed conditional on employment. The estimator is available in the causalweight package for the statistical software R.


title: Projectors on invariant subspaces of representations $\operatorname{ad}^{\otimes 2}$ of Lie algebras $so(N)$ and $sp(2r)$ and Vogel parametrization

abstract: Explicit formulae for the projectors onto invariant subspaces of the $\operatorname{ad}^{\otimes 2}$ representation of the Lie algebras $so(N)$ and $sp(2r)$ have been found by means of the split Casimir operator. These projectors have also been considered from the viewpoint of the universal complex simple Lie algebra description by using the Vogel parametrisation.


title: The Rapid ASKAP Continuum Survey I: Design and First Results

abstract: The Rapid ASKAP Continuum Survey (RACS) is the first large-area survey to be conducted with the full 36-antenna Australian Square Kilometre Array Pathfinder (ASKAP) telescope. RACS will provide a shallow model of the ASKAP sky that will aid the calibration of future deep ASKAP surveys. RACS will cover the whole sky visible from the ASKAP site in Western Australia, and will cover the full ASKAP band of $700-1800$ MHz. The RACS images are generally deeper than the existing NRAO VLA Sky Survey (NVSS) and Sydney University Molonglo Sky Survey (SUMSS) radio surveys and have better spatial resolution. All RACS survey products will be public, including radio images (with $\sim 15$ arcsecond resolution) and catalogues of about three million source components with spectral index and polarisation information. In this paper, we present a description of the RACS survey and the first data release of 903 images covering the sky south of declination $+41^\circ$ made over a 288 MHz band centred at 887.5 MHz.


title: On The Gaussian Approximation To Bayesian Posterior Distributions

abstract: The present article derives the minimal number $N$ of observations needed to consider a Bayesian posterior distribution as Gaussian. Two examples are presented. Within one of them, a chi-squared distribution, the observable $x$ as well as the parameter $\xi$ are defined all over the real axis, in the other one, the binomial distribution, the observable $x$ is an entire number while the parameter $\xi$ is defined on a finite interval of the real axis. The required minimal $N$ is high in the first case and low for the binomial model. In both cases the precise definition of the measure $\mu$ on the scale of $\xi$ is crucial.


title: Keplers's Equation and Angular Momentum: Historical Perspective, Critical Analysis and Implications for Development of the Orbital Mechanics/Dynamics, Mathematics and Physics

abstract: After some more than four centuries from the formulation and publication (in Astronomia Nova) of the Kepler's Equation, which relates the eccentric (and, intermediately, the true) anomaly of the planetary trajectories to the uniformly flowing time, in accordance with his Second ("Area") law, the subsequently -- in course of development of Orbital Mechanics -- to the 2nd law related and formally derived non-existent (zero-valued) transverse acceleration is questioned. Certain implications to Elliptic Integration, Symplectic Integration, Symplectic Geometry/Topology, as well as the connection between physical and mathematical continua in the context of the multi-level, scale-invariant mechanics/dynamics (with the augmented central and torquing forces) are also briefly hinted to.