title: Achieving Heisenberg scaling with maximally entangled states: an analytic upper bound for the attainable root mean square error

abstract: In this paper we explore the possibility of performing Heisenberg limited quantum metrology of a phase, without any prior, by employing only maximally entangled states. Starting from the estimator introduced by Higgins et al. in New J. Phys. 11, 073023 (2009), the main result of this paper is to produce an analytical upper bound on the associated Mean Squared Error which is monotonically decreasing as a function of the square of the number of quantum probes used in the process. The analysed protocol is non-adaptive and requires in principle (for distinguishable probes) only separable measurements. We explore also metrology in presence of a limitation on the entanglement size and in presence of loss.


title: On the cone of effective surfaces on $\overline{\mathcal A}_3$

abstract: We determine five extremal effective rays of the four-dimensional cone of effective surfaces on the toroidal compactification $\overline{\mathcal A}_3$ of the moduli space ${\mathcal A}_3$ of complex principally polarized abelian threefolds, and we conjecture that the cone of effective surfaces is generated by these surfaces. As the surfaces we define can be defined in any genus $g\ge 3$, we further conjecture that they generate the cone of effective surfaces on the perfect cone toroidal compactification of ${\mathcal A}_g$ for any $g\ge 3$.


title: A study of radiation tolerance in optical cements

abstract: We study the effect ionizing radiation has on light transmission in the wavelength range 190--1100~nm for a number of optically clear epoxies. We find that the transmittance of traditional, commercially available, optical epoxies show significant degradation for exposures of $1\times10^{12}$~MIPs/cm$^2$. Degradation of light transmission progresses from the shortest wavelengths at low doses to longer wavelengths as the dose increases. In epoxy joints that are 0.1~mm thick, we observe that more than 5\% of the light is lost for wavelengths less than 400~nm for traditional optical epoxies. Our studies have identified an optically clear epoxy that shows little degradation for radiation exposures up to $5.9\times 10^{14}$~MIPs/cm$^2$ ($\approx 220$~kGy).


title: Flow-driven branching in a frangible porous medium

abstract: Channel formation and branching is widely seen in physical systems where movement of fluid through a porous structure causes the spatiotemporal evolution of the medium in response to the flow, in turn causing flow pathways to evolve. We provide a simple theoretical framework that embodies this feedback mechanism in a multi-phase model for flow through a fragile porous medium with a dynamic permeability. Numerical simulations of the model show the emergence of branched networks whose topology is determined by the geometry of external flow forcing. This allows us to delineate the conditions under which splitting and/or coalescing branched network formation is favored, with potential implications for both understanding and controlling branching in soft frangible media.


title: Cometary glycolaldehyde as a source of pre-RNA molecules

abstract: Over 200 molecules have been detected in multiple extraterrestrial environments, including glycolaldehyde (C2(H2O)2, GLA), a two-carbon sugar precursor that has been detected in regions of the interstellar medium. Its recent in situ detection on the nucleus of comet 67P/Churyumov-Gerasimenko and via remote observations in the comae of others, provides tantalizing evidence that it is common on most (if not all) comets. Impact experiments conducted at the Experimental Impact Laboratory at NASA's Johnson Space Center have shown that samples of GLA and GLA mixed with montmorillonite clays can survive impact delivery in the pressure range of 4.5 GPa to 25 GPa. Extrapolated to amounts of GLA observed on individual comets and assuming a monotonic impact rate in the first billion years of solar system history, these experimental results show that up to 10^23 kg of cometary GLA could have survived impact delivery, with substantial amounts of threose, erythrose, glycolic acid, and ethylene glycol also produced or delivered. Importantly, independent of the profile of the impact flux in the early solar system, comet delivery of GLA would have provided (and may continue to provide) a reservoir of starting material for the formose reaction (to form ribose) and the Strecker reaction (to form amino acids). Thus, comets may have been important delivery vehicles for starting molecules necessary for life as we know it.


title: Cerman: Software for simulating streamer propagation in dielectric liquids based on the Townsend-Meek criterion

abstract: We present a software to simulate the propagation of positive streamers in dielectric liquids. Such liquids are commonly used for electric insulation of high-power equipment. We simulate electrical breakdown in a needle-plane geometry, where the needle and the extremities of the streamer are modeled by hyperboloids, which are used to calculate the electric field in the liquid. If the field is sufficiently high, electrons released from anions in the liquid can turn into electron avalanches, and the streamer propagates if an avalanche meets the Townsend-Meek criterion. The software is written entirely in Python and released under an MIT license. We also present a set of model simulations demonstrating the capability and versatility of the software.


title: An Iterative Method for Contour-Based Nonlinear Eigensolvers

abstract: Contour integration techniques have become a popular choice for solving the linear and non-linear eigenvalue problems. They principally include the Sakurai-Sugiura methods, the Beyn's algorithm, the FEAST/NLFEAST algorithms and other rational filtering techniques. While these methods can result in effective 'black-box' approach for solving linear eigenvalue problems, they still present several shortcomings for addressing nonlinear eigenvalue problems which are both mathematically and practically far more challenging. In this paper, we introduce a new hybrid algorithm that advantageously combines the iterative nature of NLFEAST with the effectiveness of Beyn's approach to deal with general non-linearity. In doing so, this NLFEAST-Beyn hybrid algorithm can overcome current limitations of both algorithms taken separately. After presenting its derivation from both a Beyn's and NLFEAST's perspective, several numerical examples are discussed to demonstrate the efficiency of the new technique.


title: Massively Multilingual ASR: 50 Languages, 1 Model, 1 Billion Parameters

abstract: We study training a single acoustic model for multiple languages with the aim of improving automatic speech recognition (ASR) performance on low-resource languages, and over-all simplifying deployment of ASR systems that support diverse languages. We perform an extensive benchmark on 51 languages, with varying amount of training data by language(from 100 hours to 1100 hours). We compare three variants of multilingual training from a single joint model without knowing the input language, to using this information, to multiple heads (one per language cluster). We show that multilingual training of ASR models on several languages can improve recognition performance, in particular, on low resource languages. We see 20.9%, 23% and 28.8% average WER relative reduction compared to monolingual baselines on joint model, joint model with language input and multi head model respectively. To our knowledge, this is the first work studying multilingual ASR at massive scale, with more than 50 languages and more than 16,000 hours of audio across them.


title: CHANG-ES XXI. Transport processes and the X-shaped magnetic field of NGC 4217: off-center superbubble structure

abstract: In order to gain a better understanding of the influence of cosmic rays (CRs) and magnetic fields in the disk-halo interface of edge-on spiral galaxies, we investigate the radio continuum halo, the magnetic field, and the transport processes of the CRs of the edge-on spiral galaxy NGC 4217 using CHANG-ES radio data at two frequencies, 6 GHz (C-band) and 1.5 GHz (L-band), and supplemental LOFAR data at 150 MHz and X-ray Chandra data. NGC 4217 shows a large-scale X-shaped magnetic field structure, covering a major part of the galaxy with a mean total magnetic field strength in the disk of 9 micro Gauss (via equipartition). Using rotation measure synthesis (RM-synthesis) at C-band, we found that the direction of the disk magnetic field is pointing inward. A helical outflow structure is furthermore present in the northwestern part of the galaxy, which is extended nearly 7 kpc into the halo. More polarized emission is observed on the approaching side of the galaxy. With a simplified galaxy disk model, we are able to explain that finding and predict that roughly 75% of edge-on spiral galaxies will show higher polarized intensity on the approaching side. Many loop and shell structures are found throughout the galaxy in total intensity at C-band. A superbubble-like structure is prominent in total and polarized intensity, as well as in Halpha and optical dust filaments, being a possible result of concentrated star formation in the disk. The flux density contribution of the disk in comparison to the halo decreases toward lower frequencies. Total intensity profiles at the three radio frequencies were fit with two-component exponential functions. The frequency dependence of the resulting scale heights between C-band and L-band suggests advection to be the main CR transport process. The 1D CR transport modeling (SPINNAKER) shows that advection appears to be more important than diffusion.


title: From orthocomplementations to locality

abstract: After some background on lattices, the locality framework introduced in earlier work by the authors is extended to cover posets and lattices. We then extend the correspondence between Euclidean structures on vector spaces and orthogonal complementations to a one-one correspondence between a class of locality structures and orthocomplementations on bounded lattices. This recasts in the context of renormalisation classical results in lattice theory.


title: Homotopy theory of curved operads and curved algebras

abstract: Curved algebras are algebras endowed with a predifferential, which is an endomorphism of degree -1 whose square is not necessarily 0. This makes the usual definition of quasi-isomorphism meaningless and therefore the homotopical study of curved algebras cannot follow the same path as differential graded algebras. In this article, we propose to study curved algebras by means of curved operads. We develop the theory of bar and cobar constructions adapted to this new notion as well as Koszul duality theory. To be able to provide meaningful definitions, we work in the context of objects which are filtered and complete and become differential graded after applying the associated graded functor. This setting brings its own difficulties but it nevertheless permits us to define a combinatorial model category structure that we can transfer to the category of curved operads and to the category of algebras over a curved operad using free-forgetful adjunctions. We address the case of curved associative algebras. We recover the notion of curved Aoo-algebras, and we show that the homotopy categories of curved associative algebras and of curved Aoo-algebras are Quillen equivalent.


title: Benchmarking Quantum Annealing Controls with Portfolio Optimization

abstract: Quantum annealing offers a novel approach to finding the optimal solutions for a variety of computational problems, where the quantum annealing controls influence the observed performance and error mechanisms by tuning the underlying quantum dynamics. However, the influence of the available controls is often poorly understood, and methods for evaluating the effects of these controls are necessary to tune quantum computational performance. Here we use portfolio optimization as a case study by which to benchmark quantum annealing controls and their relative effects on computational accuracy. We compare empirical results from the D-Wave 2000Q quantum annealer to the computational ground truth for a variety of portfolio optimization instances. We evaluate both forward and reverse annealing methods and we identify control variations that yield optimal performance in terms of probability of success and probability of chain breaks.


title: Announcing CzEng 2.0 Parallel Corpus with over 2 Gigawords

abstract: We present a new release of the Czech-English parallel corpus CzEng 2.0 consisting of over 2 billion words (2 "gigawords") in each language. The corpus contains document-level information and is filtered with several techniques to lower the amount of noise. In addition to the data in the previous version of CzEng, it contains new authentic and also high-quality synthetic parallel data. CzEng is freely available for research and educational purposes.


title: Optimal Dynamic Mechanism Design with Stochastic Supply and Flexible Consumers

abstract: We consider the problem of designing an expected-revenue maximizing mechanism for allocating multiple non-perishable goods of $k$ varieties to flexible consumers over $T$ time steps. In our model, a random number of goods of each variety may become available to the seller at each time and a random number of consumers may enter the market at each time. Each consumer is present in the market for one time step and wants to consume one good of one of its desired varieties. Each consumer is associated with a flexibility level that indicates the varieties of the goods it is equally interested in. A consumer's flexibility level and the utility it gets from consuming a good of its desired varieties are its private information. We characterize the allocation rule for a Bayesian incentive compatible, individually rational and expected revenue maximizing mechanism in terms of the solution to a dynamic program. The corresponding payment function is also specified in terms of the optimal allocation function. We leverage the structure of the consumers' flexibility model to simplify the dynamic program and provide an alternative description of the optimal mechanism in terms of thresholds computed by the dynamic program.


title: Beam-Shaping PEC Mirror Phase Corrector Design

abstract: The Perfect Electric Conductor (PEC) mirror phase corrector plays an important role in the beam-shaping mirror system design for Quasi-Optical (QO) mode converter (launcher) in the sub-THz high-power gyrotron. In this article, both the Geometry Optical (GO) method and the phase gradient method have been presented for the PEC mirror phase corrector design. The advantages and disadvantages are discussed for both methods. An efficient algorithm has been proposed for the phase gradient method.


title: On the classification of solutions to a weighted elliptic system involving the Grushin operator

abstract: We investigate here the following weighted degenerate elliptic system \begin{align*} -\Delta_{s} u = \Big(1+\|\mathbf{x}\|^{2(s+1)}\Big)^{\frac{\alpha}{2(s+1)}} v^p, \quad -\Delta_{s} v = \Big(1+\|\mathbf{x}\|^{2(s+1)}\Big)^{\frac{\alpha}{2(s+1)}}u^\theta, \quad u,v>0\quad\mbox{in }\; \mathbb{R}^N:=\mathbb{R}^{N_1}\times \mathbb{R}^{N_2}. \end{align*} where $\Delta_{s}=\Delta_{x}+|x|^{2s}\Delta_{y},$ is the Grushin operator, $s \geq 0,$ $\alpha \geq 0$ and $1<p\leq\theta.$ Here $$\|\mathbf{x}\|=\Big(|x|^{2(s+1)}+|y|^2\Big)^{\frac{1}{2(s+1)}}, \;\mbox{and}\;\; \mathbf{x}:=(x, y)\in \mathbb{R}^N:=\mathbb{R}^{N_1}\times \mathbb{R}^{N_2}.$$ In particular, we establish some new Liouville-type theorems for stable solutions of the system, which recover and considerably improve upon the known results \cite{cow, Hfh, HU, Fa, DP}. As a consequence, we obtain a nonexistence result for the weighted Grushin equation \begin{align*} -\Delta_{s} u =\Big(1+\|\mathbf{x}\|^{2(s+1)}\Big)^{\frac{\alpha}{2(s+1)}} u^p,\;\; \quad u>0 \quad \mbox{in }\;\; \mathbb{R}^N. \end{align*}


title: The origin of the black hole offset in M31

abstract: Using state-of-the-art high-resolution fully GPU N-body simulations, we demonstrate for the first time that the infall of a dark matter rich satellite naturally explains a present black hole offset by sub-parsecs in M31. Observational data of the tidal features provide stringent constraints on the initial conditions of our simulations. The heating of the central region of M31 by the satellite via dynamical friction entails a significant black hole offset after the first pericentric passage. After having reached its maximum offset, the massive black hole sinks towards the M31 centre due to dynamical friction and it is determined to be offset by sub-parsecs as derived by observations.


title: Continuous Maps from Spheres Converging to Boundaries of Convex Hulls

abstract: Given $n$ distinct points $\mathbf{x}_1, \ldots, \mathbf{x}_n$ in $\mathbb{R}^d$, let $K$ denote their convex hull, which we assume to be $d$-dimensional, and $B = \partial K $ its $(d-1)$-dimensional boundary. We construct an explicit one-parameter family of continuous maps $\mathbf{f}_{\varepsilon} \colon \mathbb{S}^{d-1} \to K$ which, for $\varepsilon > 0$, are defined on the $(d-1)$-dimensional sphere and have the property that the images $\mathbf{f}_{\varepsilon}(\mathbb{S}^{d-1})$ are codimension $1$ submanifolds contained in the interior of $K$. Moreover, as the parameter $\varepsilon$ goes to $0^+$, the images $\mathbf{f}_{\varepsilon}(\mathbb{S}^{d-1})$ converge, as sets, to the boundary $B$ of the convex hull. We prove this theorem using techniques from convex geometry of (spherical) polytopes and set-valued homology. We further establish an interesting relationship with the Gauss map of the polytope $B$, appropriately defined. Several computer plots illustrating our results will be presented.


title: Neutron Stars in f(R)-Gravity and Its Extension with a Scalar Axion Field

abstract: We present a brief review of general results about non-rotating neutron stars in simple R-square gravity and its extension with a scalar axion field. Modified Einstein equations are presented for metrics in isotropical coordinates. The mass-radius relation, mass profile and dependence of mass from central density on various equations of state are given in comparison to general relativity.


title: On the Validity of Physical Optics for Narrow-band Beam Scattering and Diffraction from the Open Cylindrical Surface

abstract: The exact formulas for the induced electric surface current (in the scattering phenomenon) and the equivalent electric surface current (in the diffraction phenomenon) on the open cylindrical surface due to an arbitrary narrow-band beam have been shown in their closed-form expressions within the context of the cylindrical harmonics, which gives information about the validity of the Physical Optics (PO) approximation. Both the Electric Field Integral Equation (EFIE) and the Magnetic Field Integral Equation (MFIE) are used to find the induced (equivalent) electric surface currents in the context of the cylindrical harmonics. The numerical example of the scattering and diffraction of the Hermite Gaussian beam from the open cylindrical surface is shown. The result is useful for the evaluation of the validity of the PO approximation in the cylinder-like surface.


title: Topic-based Community Search over Spatial-Social Networks (Technical Report)

abstract: Recently, the community search problem has attracted significant attention, due to its wide spectrum of real-world applications such as event organization, friend recommendation, advertisement in e-commence, and so on. Given a query vertex, the community search problem finds dense subgraph that contains the query vertex. In social networks, users have multiple check-in locations, influence score, and profile information (keywords). Most previous studies that solve the CS problem over social networks usually neglect such information in a community. In this paper, we propose a novel problem, named community search over spatial-social networks (TCS-SSN), which retrieves community with high social influence, small traveling time, and covering certain keywords. In order to tackle the TCS-SSN problem over the spatial-social networks, we design effective pruning techniques to reduce the problem search space. We also propose an effective indexing mechanism, namely social-spatial index, to facilitate the community query, and develop an efficient query answering algorithm via index traversal. We verify the efficiency and effectiveness of our pruning techniques, indexing mechanism, and query processing algorithm through extensive experiments on real-world and synthetic data sets under various parameter settings.


title: Big Data links from Climate to Commodity Production Forecasts and Risk Management

abstract: Frozen concentrated orange juice (FCOJ) is a commodity traded in the International Commodity Exchange. The FCOJ future price volatility is high because the world's orange production is concentrated in a few places, which results in extreme sensitivity to weather and disease. Most of the oranges produced in the United States are from Florida. The United States Department of Agriculture (USDA) issues orange production forecasts on the second week of each month from October to July. The October forecast in particular seems to affect FCOJ price volatility. We assess how a prediction of the directionality and magnitude of the error of the USDA October forecast could affect the decision making process of multiple FCOJ market participants, and if the "production uncertainty" of the forecast could be reduced by incorporating other climate variables. The models developed open up the opportunity to assess the application of the resulting probabilistic forecasts of the USDA production forecast error on the trading decisions of the different FCOJ stakeholders, and to perhaps consider the inclusion of climate predictors in the USDA forecast.


title: Multiple Imputation with Massive Data: An Application to the Panel Study of Income Dynamics

abstract: \Multiple imputation (MI) is a popular and well-established method for handling missing data in multivariate data sets, but its practicality for use in massive and complex data sets has been questioned. One such data set is the Panel Study of Income Dynamics (PSID), a longstanding and extensive survey of household income and wealth in the U.S. Missing data for this survey are currently handled using traditional hot deck methods because of the simple implementation; however, the univariate hot deck results in large random wealth fluctuations. MI is effective but faced with operational challenges. We use a sequential regression/ chained-equation approach, using the software IVEware, to multiply impute cross-sectional wealth data in the 2013 PSID, and compare analyses of the resulting imputed data with those from the current hot deck approach. Practical difficulties, such as non-normally distributed variables, skip patterns, categorical variables with many levels, and multicollinearity, are described together with our approaches to overcoming them. We evaluate the imputation quality and validity with internal diagnostics and external benchmarking data. MI produces improvements over the existing hot deck approach by helping preserve correlation structures, such as the associations between PSID wealth components and the relationships between the household net worth and socio-demographic factors, and facilitates completed data analyses with general purposes. MI incorporates highly predictive covariates into imputation models and increases efficiency. We recommend the practical implementation of MI and expect greater gains when the fraction of missing information is large.


title: Confinement controlled bend instability of three-dimensional active fluids

abstract: Spontaneous growth of long-wavelength deformations is a defining feature of active fluids with orientational order. We investigate the effect of biaxial rectangular confinement on the instability of initially shear-aligned 3D isotropic active fluids composed of extensile microtubule bundles and kinesin molecular motors. Under confinement, such fluids exhibit finite-wavelength self-amplifying bend deformations which grow in the plane orthogonal to the direction of the strongest confinement. Both the instability wavelength and the growth rate increase with weakening confinement. These findings are consistent with a minimal hydrodynamic model, which predicts that the fastest growing deformation is set by a balance of active driving and elastic relaxation. Experiments in the highly confined regime confirm that the instability wavelength is set by the balance of active and elastic stresses, which are independently controlled by the concentration of motors and non-motile crosslinkers.


title: 3D EIT Reconstructions from Electrode Data using Direct Inversion D-bar and Calderon Methods

abstract: The first numerical implementation of a D-bar method in 3D using electrode data is presented. Results are compared to Calder\'on's method as well as more common TV and smoothness regularization-based methods. D-bar methods are based on tailor-made non-linear Fourier transforms involving the measured current and voltage data. Low-pass filtering in the non-linear Fourier domain is used to stabilize the reconstruction process. D-bar methods have shown great promise in 2D for providing robust real-time absolute and time-difference conductivity reconstructions but have yet to be used on practical electrode data in 3D, until now. Results are presented for simulated data for conductivity and permittivity with disjoint non-radially symmetric targets on spherical domains and noisy voltage data. The 3D D-bar and Calder\'on methods are demonstrated to provide comparable quality to their 2D CGO counterparts, and hold promise for real-time reconstructions.


title: Incorrect Data in the Widely Used Inside Airbnb Dataset

abstract: Several recently published papers in Decision Support Systems discussed issues related to data quality in Information Systems research. In this short research note, I build on the work introduced in these papers and document two data quality issues discovered in a large open dataset commonly used in research. Inside Airbnb (IA) collects data from places and reviews as posted by users of Airbnb.com. Visitors can effortlessly download data collected by IA for several locations around the globe. While the dataset is widely used in academic research, no thorough investigation of the dataset and its validity has been conducted. This note examines the dataset and explains an issue of incorrect data added to the dataset. Findings suggest that this issue can be attributed to systemic errors in the data collection process. The results suggest that the use of unverified open datasets can be problematic, although the discoveries presented in this work may not be significant enough to challenge all published research that used the IA dataset. Additionally, findings indicate that the incorrect data happens because of a new feature implemented by Airbnb. Thus, unless changes are made, it is likely that the consequences of this issue will only become more severe. Finally, this note explores why reproducibility is a problem when two different releases of the dataset are compared.


title: Deep Contextual Embeddings for Address Classification in E-commerce

abstract: E-commerce customers in developing nations like India tend to follow no fixed format while entering shipping addresses. Parsing such addresses is challenging because of a lack of inherent structure or hierarchy. It is imperative to understand the language of addresses, so that shipments can be routed without delays. In this paper, we propose a novel approach towards understanding customer addresses by deriving motivation from recent advances in Natural Language Processing (NLP). We also formulate different pre-processing steps for addresses using a combination of edit distance and phonetic algorithms. Then we approach the task of creating vector representations for addresses using Word2Vec with TF-IDF, Bi-LSTM and BERT based approaches. We compare these approaches with respect to sub-region classification task for North and South Indian cities. Through experiments, we demonstrate the effectiveness of generalized RoBERTa model, pre-trained over a large address corpus for language modelling task. Our proposed RoBERTa model achieves a classification accuracy of around 90% with minimal text preprocessing for sub-region classification task outperforming all other approaches. Once pre-trained, the RoBERTa model can be fine-tuned for various downstream tasks in supply chain like pincode suggestion and geo-coding. The model generalizes well for such tasks even with limited labelled data. To the best of our knowledge, this is the first of its kind research proposing a novel approach of understanding customer addresses in e-commerce domain by pre-training language models and fine-tuning them for different purposes.


title: Determination of the most representative descriptor among a set of feature vectors for the same object

abstract: On an example of solution of the face recognition problem the approach for estimation of the most representative descriptor among a set of feature vectors for the same face is considered in present study. The estimation is based on robust calculation of the mode-median mixture vector for the set as the descriptor by means of Welsch/Leclerc loss function application in case of very sparse filling of the feature space with feature vectors


title: Renormalisation scale setting for D-mixing

abstract: A naive application of the heavy quark expansion (HQE) yields theory estimates for the decay rate of neutral $D$ mesons that are four orders of magnitude below the experimental determination. It is well known that this huge suppression results from severe GIM cancellations. We find that this mismatch can be solved by individually choosing the renormalisation scale of the different internal quark contributions. For $b$ and $c$ hadron lifetimes, as well as for the decay rate difference of neutral $B$ mesons the effect of our scale setting procedure lies within the previously quoted theory uncertainties, while we get enlarged theory uncertainties for the semi-leptonic CP asymmetries in the $B$ system.


title: Rank inequalities for the Heegaard Floer homology of branched covers

abstract: Given a double cover between 3-manifolds branched along a nullhomologous link, we establish an inequality between the dimensions of their Heegaard Floer homologies. We discuss the relationship with the L-space conjecture and give some other topological applications, as well as an analogous result for sutured Floer homology.


title: On self-gravitating strange dark matter halos around galaxies

abstract: A new family of nonrelativistic, Newtonian, non-quantum equilibrium configurations describing galactic halos is introduced, by considering strange quark matter conglomerates with masses larger than about 8 GeV as new possible components of the dark matter. Originally introduced to explain the state of matter in neutron stars, such conglomerates may also form in the high-density and temperature conditions of the primordial Universe and then decouple from ordinary baryonic matter, providing the fundamental components of dark matter for the formation of pristine gravitational potential wells and the subsequent evolution of cosmic structures. The obtained results for halo mass and radius are consistent with the rotational velocity curve observed in the Galaxy. Additionally, the average density of such dark matter halos is similar to that derived for halos of dwarf spheroidal galaxies, which can therefore be interpreted as downscaled versions of larger dark matter distributions around Milky Way-sized galaxies and hint for a common origin of the two families of cosmic structures.


title: Deep Reinforcement Learning for Cybersecurity Assessment of Wind Integrated Power Systems

abstract: The integration of renewable energy sources (RES) is rapidly increasing in electric power systems (EPS). While the inclusion of intermittent RES coupled with the wide-scale deployment of communication and sensing devices is important towards a fully smart grid, it has also expanded the cyber-threat landscape, effectively making power systems vulnerable to cyberattacks. This paper proposes a cybersecurity assessment approach designed to assess the cyberphysical security of EPS. The work takes into consideration the intermittent generation of RES, vulnerabilities introduced by microprocessor-based electronic information and operational technology (IT/OT) devices, and contingency analysis results. The proposed approach utilizes deep reinforcement learning (DRL) and an adapted Common Vulnerability Scoring System (CVSS) score tailored to assess vulnerabilities in EPS in order to identify the optimal attack transition policy based on N-2 contingency results, i.e., the simultaneous failure of two system elements. The effectiveness of the work is validated via numerical and real-time simulation experiments performed on literature-based power grid test cases. The results demonstrate how the proposed method based on deep Q-network (DQN) performs closely to a graph-search approach in terms of the number of transitions needed to find the optimal attack policy, without the need for full observation of the system. In addition, the experiments present the method's scalability by showcasing the number of transitions needed to find the optimal attack transition policy in a large system such as the Polish 2383 bus test system. The results exhibit how the proposed approach requires one order of magnitude fewer transitions when compared to a random transition policy.


title: Real Constituents of Permutation Characters

abstract: We prove a broad generalization of a theorem of W. Burnside on real characters using permutation characters. Under a necessary hypothesis, We can give some control on multiplicities (a result that needs the Classification of Finite Simple Groups). Along the way, we also give a new characterization of the 2-closed finite groups using odd-order real elements of the group. All this can be seen as a contribution to Brauer's Problem 11. We also obtain similar results for 2-Brauer characters. We also classify finite primitive permutation groups in which every real element has a fixed point.


title: Limits and possibilities of refractive index in atomic systems

abstract: Manipulation of the refractive index has been of growing interest lately. We consider parameters and possibilities of enhancing the absolute-value limit of the linear index in coherent atomic systems. Starting with a review of how two-level transitions - without and with added coherence effects - do not realistically allow for a significantly enhanced index at fixed atomic density, we discuss a possible way around this using wave mixing. This confirms that the only parameters, besides the medium optical depth/density, that can effectively change the value of the attainable index are the frequencies of the involved transitions.


title: Labeling of Multilingual Breast MRI Reports

abstract: Medical reports are an essential medium in recording a patient's condition throughout a clinical trial. They contain valuable information that can be extracted to generate a large labeled dataset needed for the development of clinical tools. However, the majority of medical reports are stored in an unregularized format, and a trained human annotator (typically a doctor) must manually assess and label each case, resulting in an expensive and time consuming procedure. In this work, we present a framework for developing a multilingual breast MRI report classifier using a custom-built language representation called LAMBR. Our proposed method overcomes practical challenges faced in clinical settings, and we demonstrate improved performance in extracting labels from medical reports when compared with conventional approaches.


title: Detailed spectroscopy of doubly magic $^{132}$Sn

abstract: The structure of the doubly magic $^{132}_{50}$Sn$_{82}$ has been investigated at the ISOLDE facility at CERN, populated both by the $\beta^-$decay of $^{132}$In and $\beta^-$-delayed neutron emission of $^{133}$In. The level scheme of $^{132}$Sn is greatly expanded with the addition of 68 $\gamma$-transitions and 17 levels observed for the first time in the $\beta$ decay. The information on the excited structure is completed by new $\gamma$-transitions and states populated in the $\beta$-n decay of $^{133}$In. Improved delayed neutron emission probabilities are obtained both for $^{132}$In and $^{133}$In. Level lifetimes are measured via the Advanced Time-Delayed $\beta\gamma\gamma$(t) fast-timing method. An interpretation of the level structure is given based on the experimental findings and the particle-hole configurations arising from core excitations both from the \textit{N} = 82 and \textit{Z} = 50 shells, leading to positive and negative parity particle-hole multiplets. The experimental information provides new data to challenge the theoretical description of $^{132}$Sn.


title: Time-energy uncertainty does not create particles

abstract: In this contribution in honour of Paul Busch, we criticise the claims of many expositions that the time-energy uncertainty principle allows both a violation of energy conservation, and particle creation, provided that this happens for a sufficiently short time. But we agree that there are grains of truth in these claims: which we make precise and justify using perturbation theory.


title: Theoretical and Practical Limits of Kolmogorov-Zurbenko Periodograms with DiRienzo-Zurbenko Algorithm Smoothing in the Spectral Analysis of Time Series Data

abstract: The Kolomogorov-Zurbenko periodogram with DiRienzo-Zurbenko algorithm smoothing is the state-of-the-art method for spectral analysis of time series data. Because this approach assumes that a sinusoidal model underlies time-series data and because its algorithms are adaptive in nature, it is superior to traditional use of autoregressive integral moving average (ARIMA) algorithms. This article begins with a presentation of its statistical derivation and development followed by instructions for accessing and utilizing this approach within the R statistical program platform. The discussion then turns to a presentation of its theoretical and practical limits with regard to sensitivity (i.e., ability to detect weak signals), accuracy (i.e., ability to correctly identify signal frequencies), resolution (i.e., ability to resolve signals with close frequencies), and robustness with respect to missing data (i.e., sensitivity and accuracy despite high levels of missingness). Next using a simulated time series in which two signals close in frequency are embedded in significant amounts of random noise, the predictive power of this approach is compared to the traditional ARIMA approach, with support also garnered for its being robust even in the face of significant levels of missing data. The article concludes with brief descriptions of studies across a range of scientific disciplines that have capitalized on the power of the Kolmogorov-Zurbenko periodogram with DiRienzo-Zurbenko algorithm smoothing.


title: Continual Learning in Human Activity Recognition: an Empirical Analysis of Regularization

abstract: Given the growing trend of continual learning techniques for deep neural networks focusing on the domain of computer vision, there is a need to identify which of these generalizes well to other tasks such as human activity recognition (HAR). As recent methods have mostly been composed of loss regularization terms and memory replay, we provide a constituent-wise analysis of some prominent task-incremental learning techniques employing these on HAR datasets. We find that most regularization approaches lack substantial effect and provide an intuition of when they fail. Thus, we make the case that the development of continual learning algorithms should be motivated by rather diverse task domains.


title: Renormalization of Galilean Electrodynamics

abstract: We study the quantum properties of a Galilean-invariant abelian gauge theory coupled to a Schr\"odinger scalar in 2+1 dimensions. At the classical level, the theory with minimal coupling is obtained from a null-reduction of relativistic Maxwell theory coupled to a complex scalar field in 3+1 dimensions and is closely related to the Galilean electromagnetism of Le-Bellac and L\'evy-Leblond. Due to the presence of a dimensionless, gauge-invariant scalar field in the Galilean multiplet of the gauge-field, we find that at the quantum level an infinite number of couplings is generated. We explain how to handle the quantum corrections systematically using the background field method. Due to a non-renormalization theorem, the beta function of the gauge coupling is found to vanish to all orders in perturbation theory, leading to a continuous family of fixed points where the non-relativistic conformal symmetry is preserved.


title: Nonlinear Transform Coding

abstract: We review a class of methods that can be collected under the name nonlinear transform coding (NTC), which over the past few years have become competitive with the best linear transform codecs for images, and have superseded them in terms of rate--distortion performance under established perceptual quality metrics such as MS-SSIM. We assess the empirical rate--distortion performance of NTC with the help of simple example sources, for which the optimal performance of a vector quantizer is easier to estimate than with natural data sources. To this end, we introduce a novel variant of entropy-constrained vector quantization. We provide an analysis of various forms of stochastic optimization techniques for NTC models; review architectures of transforms based on artificial neural networks, as well as learned entropy models; and provide a direct comparison of a number of methods to parameterize the rate--distortion trade-off of nonlinear transforms, introducing a simplified one.


title: Robust Technique for Representative Volume Element Identification in Noisy Microtomography Images of Porous Materials Based on Pores Morphology and Their Spatial Distribution

abstract: Microtomography is a powerful method of materials investigation. It enables to obtain physical properties of porous media non-destructively that is useful in studies. One of the application ways is a calculation of porosity, pore sizes, surface area, and other parameters of metal-ceramic (cermet) membranes which are widely spread in the filtration industry. The microtomography approach is efficient because all of those parameters are calculated simultaneously in contrast to the conventional techniques. Nevertheless, the calculations on Micro-CT reconstructed images appear to be time-consuming, consequently representative volume element should be chosen to speed them up. This research sheds light on representative elementary volume identification without consideration of any physical parameters such as porosity, etc. Thus, the volume element could be found even in noised and grayscale images. The proposed method is flexible and does not overestimate the volume size in the case of anisotropic samples. The obtained volume element could be used for computations of the domain's physical characteristics if the image is filtered and binarized, or for selections of optimal filtering parameters for denoising procedure.


title: A Discrete Ion Stochastic Continuum Overdamped Solvent Algorithm for Modeling Electrolytes

abstract: In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the Fluctuating Immersed Boundary (FIB) approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both cases the Immersed Boundary (IB) method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long range electrostatic interactions are computed by solving the Poisson equation, with short range corrections included using a novel immersed-boundary variant of the classical Particle-Particle Particle-Mesh (P3M) technique. Also included is a short range repulsive force based on the Weeks-Chandler-Andersen (WCA) potential. The new methodology is validated by comparison to Debye-H{\"u}ckel theory for ion-ion pair correlation functions, and Debye-H{\"u}ckel-Onsager theory for conductivity, including the Wein effect for strong electric fields. In each case good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.


title: Curve counting and S-duality

abstract: We work on a projective threefold $X$ which satisfies the Bogomolov-Gieseker conjecture of Bayer-Macr\`i-Toda, such as $\mathbb P^3$ or the quintic threefold. We prove certain moduli spaces of 2-dimensional torsion sheaves on $X$ are smooth bundles over Hilbert schemes of ideal sheaves of curves and points in $X$. When $X$ is Calabi-Yau this gives a simple wall crossing formula expressing curve counts (and so ultimately Gromov-Witten invariants) in terms of counts of D4-D2-D0 branes. These latter invariants are predicted to have modular properties which we discuss from the point of view of S-duality and Noether-Lefschetz theory.


title: Stillman's question for twisted commutative algebras

abstract: Let $\mathbf{A}_{n, m}$ be the polynomial ring $\text{Sym}(\mathbf{C}^n \otimes \mathbf{C}^m)$ with the natural action of $\mathbf{GL}_m(\mathbf{C})$. We construct a family of $\mathbf{GL}_m(\mathbf{C})$-stable ideals $J_{n, m}$ in $\mathbf{A}_{n, m}$, each equivariantly generated by one homogeneous polynomial of degree $2$. Using the Ananyan-Hochster principle, we show that the regularity of this family is unbounded. This negatively answers a question raised by Erman-Sam-Snowden on a generalization of Stillman's conjecture.


title: Streaming Verification for Graph Problems: Optimal Tradeoffs and Nonlinear Sketches

abstract: We study graph computations in an enhanced data streaming setting, where a space-bounded client reading the edge stream of a massive graph may delegate some of its work to a cloud service. We seek algorithms that allow the client to verify a purported proof sent by the cloud service that the work done in the cloud is correct. A line of work starting with Chakrabarti et al. (ICALP 2009) has provided such algorithms, which we call schemes, for several statistical and graph-theoretic problems, many of which exhibit a tradeoff between the length of the proof and the space used by the streaming verifier. This work designs new schemes for a number of basic graph problems---including triangle counting, maximum matching, topological sorting, and single-source shortest paths---where past work had either failed to obtain smooth tradeoffs between these two key complexity measures or only obtained suboptimal tradeoffs. Our key innovation is having the verifier compute certain nonlinear sketches of the input stream, leading to either new or improved tradeoffs. In many cases, our schemes in fact provide optimal tradeoffs up to logarithmic factors. Specifically, for most graph problems that we study, it is known that the product of the verifier's space cost $v$ and the proof length $h$ must be at least $\Omega(n^2)$ for $n$-vertex graphs. However, matching upper bounds are only known for a handful of settings of $h$ and $v$ on the curve $h \cdot v=\tilde{\Theta}(n^2)$. For example, for counting triangles and maximum matching, schemes with costs lying on this curve are only known for $(h=\tilde{O}(n^2), v=\tilde{O}(1))$, $(h=\tilde{O}(n), v=\tilde{O}(n))$, and the trivial $(h=\tilde{O}(1), v=\tilde{O}(n^2))$. A major message of this work is that by exploiting nonlinear sketches, a significant ``portion'' of costs on the tradeoff curve $h \cdot v = n^2$ can be achieved.


title: Near-Linear Time Edit Distance for Indel Channels

abstract: We consider the following model for sampling pairs of strings: $s_1$ is a uniformly random bitstring of length $n$, and $s_2$ is the bitstring arrived at by applying substitutions, insertions, and deletions to each bit of $s_1$ with some probability. We show that the edit distance between $s_1$ and $s_2$ can be computed in $O(n \ln n)$ time with high probability, as long as each bit of $s_1$ has a mutation applied to it with probability at most a small constant. The algorithm is simple and only uses the textbook dynamic programming algorithm as a primitive, first computing an approximate alignment between the two strings, and then running the dynamic programming algorithm restricted to entries close to the approximate alignment. The analysis of our algorithm provides theoretical justification for alignment heuristics used in practice such as BLAST, FASTA, and MAFFT, which also start by computing approximate alignments quickly and then find the best alignment near the approximate alignment. Our main technical contribution is a partitioning of alignments such that the number of the subsets in the partition is not too large and every alignment in one subset is worse than an alignment considered by our algorithm with high probability. Similar techniques may be of interest in the average-case analysis of other problems commonly solved via dynamic programming.


title: Parallel Algorithms for Successive Convolution

abstract: In this work, we consider alternative discretizations for PDEs which use expansions involving integral operators to approximate spatial derivatives. These constructions use explicit information within the integral terms, but treat boundary data implicitly, which contributes to the overall speed of the method. This approach is provably unconditionally stable for linear problems and stability has been demonstrated experimentally for nonlinear problems. Additionally, it is matrix-free in the sense that it is not necessary to invert linear systems and iteration is not required for nonlinear terms. Moreover, the scheme employs a fast summation algorithm that yields a method with a computational complexity of $\mathcal{O}(N)$, where $N$ is the number of mesh points along a direction. While much work has been done to explore the theory behind these methods, their practicality in large scale computing environments is a largely unexplored topic. In this work, we explore the performance of these methods by developing a domain decomposition algorithm suitable for distributed memory systems along with shared memory algorithms. As a first pass, we derive an artificial CFL condition that enforces a nearest-neighbor communication pattern and briefly discuss possible generalizations. We also analyze several approaches for implementing the parallel algorithms by optimizing predominant loop structures and maximizing data reuse. Using a hybrid design that employs MPI and Kokkos for the distributed and shared memory components of the algorithms, respectively, we show that our methods are efficient and can sustain an update rate $> 1\times10^8$ DOF/node/s. We provide results that demonstrate the scalability and versatility of our algorithms using several different PDE test problems, including a nonlinear example, which employs an adaptive time-stepping rule.

categories: physics.comp-ph cs.CE cs.DC cs.NA math.NA


title: Proca-stinated Cosmology I: A N-body code for the vector Galileon

abstract: We investigate the nonlinear growth of large-scale structure in the generalised Proca theory, in which a self-interacting massive vector field plays the role of driving the acceleration of the cosmic expansion. Focusing to the Proca Lagrangian at cubic order -- the cubic vector Galileon model -- we derive the simplified equations for gravity as well as the longitudinal and transverse modes of the vector field under the weak-field and quasi-static approximations, and implement them in a modified version of the ECOSMOG N-body code. Our simulations incorporate the Vainshtein screening effect, which reconciles the fifth force propagated by the longitudinal mode of the cubic vector Galileon model with local tests of gravity. The results confirm that for all scales probed by the simulation, the transverse mode has a negligible impact on structure formation in a realistic cosmological setup. It is well known that in this model the strength of the fifth force is controlled by a free model parameter, which we denote as beta_3. By running a suite of cosmological simulations for different values of beta_3, we show that this parameter also determines the effectiveness of the Vainshtein screening. The model behaves identically to the cubic scalar Galileon for beta_3 going to zero, in which the fifth force is strong in unscreened regions but is efficiently screened in high-density regions. In the opposite limit, beta_3 going to infinity, the model approaches its `quintessence' counterpart, which has a vanishing fifth force but a modified expansion history compared to LambdaCDM. This endows the model with rich phenomenology, which will be investigated in future works.


title: Criterion for the functional dissipativity of second order differential operators with complex coefficients

abstract: In the present paper we consider the Dirichlet problem for the second order differential operator $E=\nabla(A \nabla)$,where $A$ is a matrix with complex valued $L^\infty$ entries. We introduce the concept of dissipativity of $E$ with respect to a given function $\varphi:R^+ \to R^+$. Under the assumption that the $Im\, A$ is symmetric, we prove that the condition $|s\, \varphi'(s)| \, | \langle Im\, A (x)\, \xi,\xi\rangle |\leq 2\, \sqrt{\varphi(s)\, [s\, \varphi(s)]'}\, \langle Re\, A(x) \, \xi,\xi\rangle $ (for almost every $x\in\Omega\subset R^N$ and for any $s>0$, $\xi\in R^N$) is necessary and sufficient for the functional dissipativity of $E$.


title: High-responsivity graphene photodetectors integrated on silicon microring resonators

abstract: Graphene integrated photonics provides several advantages over conventional Si photonics. Single layer graphene (SLG) enables fast, broadband, and energy-efficient electro-optic modulators, optical switches and photodetectors (GPDs), and is compatible with any optical waveguide. The last major barrier to SLG-based optical receivers lies in the low responsivity - electrical output per optical input - of GPDs compared to conventional PDs. Here we overcome this shortfall by integrating a photo-thermoelectric GPD with a Si microring resonator. Under critical coupling, we achieve $>$90% light absorption in a $\sim$6 $\mu$m SLG channel along the Si waveguide. Exploiting the cavity-enhanced light-matter interaction, causing carriers in SLG to reach $\sim$400 K for an input power of $\sim$0.6 mW, we get a voltage responsivity $\sim$90 V/W, demonstrating the feasibility of our approach. Our device is capable of detecting data rates up to 20 Gbit/s, with a receiver sensitivity enabling it to operate at a 10$^{-9}$ bit-error rate, on par with mature semiconductor technology. The natural generation of a voltage rather than a current, removes the need for transimpedance amplification, with a reduction of the energy-per-bit cost and foot-print, when compared to a traditional semiconductor-based receiver.


title: RAPTOR II: Polarized radiative transfer in curved spacetime

abstract: Accreting supermassive black holes are sources of polarized radiation that propagates through highly curved spacetime before reaching the observer. In order to help interpret observations of such polarized emission, accurate and efficient numerical schemes for polarized radiative transfer in curved spacetime are needed. In this manuscript we extend our publicly available radiative transfer code RAPTOR to include polarization. We provide a brief review of different codes and methods for covariant polarized radiative transfer available in the literature and existing codes, and present an efficient new scheme. For the spacetime-propagation aspect of the computation, we develop a compact, Lorentz-invariant representation of a polarized ray. For the plasma-propagation aspect of the computation, we perform a formal analysis of the stiffness of the polarized radiative-transfer equation with respect to our explicit integrator, and develop a hybrid integration scheme that switches to an implicit integrator in case of stiffness, in order to solve the equation with optimal speed and accuracy for all possible values of the local optical/Faraday thickness of the plasma. We perform a comprehensive code verification by solving a number of well-known test problems using RAPTOR and comparing its output to exact solutions. We also demonstrate convergence with existing polarized radiative-transfer codes in the context of complex astrophysical problems. RAPTOR is capable of performing polarized radiative transfer in arbitrary, highly curved spacetimes. This capability is crucial for interpreting polarized observations of accreting black holes, which can yield information about the magnetic-field configuration in such accretion flows. The efficient formalism implemented in RAPTOR is computationally light and conceptually simple. The code is publicly available.


title: Machine Learning with the Sugeno Integral: The Case of Binary Classification

abstract: In this paper, we elaborate on the use of the Sugeno integral in the context of machine learning. More specifically, we propose a method for binary classification, in which the Sugeno integral is used as an aggregation function that combines several local evaluations of an instance, pertaining to different features or measurements, into a single global evaluation. Due to the specific nature of the Sugeno integral, this approach is especially suitable for learning from ordinal data, that is, when measurements are taken from ordinal scales. This is a topic that has not received much attention in machine learning so far. The core of the learning problem itself consists of identifying the capacity underlying the Sugeno integral. To tackle this problem, we develop an algorithm based on linear programming. The algorithm also includes a suitable technique for transforming the original feature values into local evaluations (local utility scores), as well as a method for tuning a threshold on the global evaluation. To control the flexibility of the classifier and mitigate the problem of overfitting the training data, we generalize our approach toward $k$-maxitive capacities, where $k$ plays the role of a hyper-parameter of the learner. We present experimental studies, in which we compare our method with competing approaches on several benchmark data sets.


title: Leveraging Class Hierarchies with Metric-Guided Prototype Learning

abstract: In many classification tasks, the set of target classes can be organized into a hierarchy. This structure induces a semantic distance between classes, and can be summarised under the form of a cost matrix, which defines a finite metric on the class set. In this paper, we propose to model the hierarchical class structure by integrating this metric in the supervision of a prototypical network. Our method relies on jointly learning a feature-extracting network and a set of class prototypes whose relative arrangement in the embedding space follows an hierarchical metric. We show that this approach allows for a consistent improvement of the error rate weighted by the cost matrix when compared to traditional methods and other prototype-based strategies. Furthermore, when the induced metric contains insight on the data structure, our method improves the overall precision as well. Experiments on four different public datasets - from agricultural time series classification to depth image semantic segmentation - validate our approach.


title: PHELP: Pixel Heating Experiment Learning Platform for Education and Research on IAI-based Smart Control Engineering

abstract: Thermal processes are one of the most common systems in the industry, making its understanding a mandatory skill for control engineers. So, multiple efforts are focused on developing low-cost and portable experimental training rigs recreating the thermal process dynamics and controls, usually limited to SISO or low order 2x2 MIMO systems. This paper presents PHELP, a low-cost, portable, and high order MIMO educational platform for uniformity temperature control training. The platform is composed of an array of 16 Peltier modules as heating elements, with a lower heating and cooling times, resulting in a 16x16 high order MIMO system. A low-cost real-time infrared thermal camera is employed as a temperature feedback sensor instead of a standard thermal sensor, ideal for high order MIMO system sensing and temperature distribution tracking. The control algorithm is developed in Matlab/Simulink and employs an Arduino board in hardware in the loop configuration to apply the control action to each Peltier module in the array. A temperature control experiment is performed, showing that the platform is suitable for teaching and training experiences not only in the classroom but also for engineers in the industry. Furthermore, various abnormal conditions can be introduced so that smart control engineering features can be tested.


abstract: The authors and co-signers of the Terrestrial Planets Comparative Climatology (TPCC) mission concept white paper advocate that planetary science in the next decade would greatly benefit from comparatively studying the fundamental behavior of the atmospheres of Venus and Mars, contemporaneously and with the same instrumentation, to capture atmospheric response to the same solar forcing, and with a minimum of instrument-related variability. Thus, this white paper was created for the 2023-2032 Planetary Science Decadal Survey process. It describes the science rationale for such a mission, and a mission concept that could achieve such a mission.


title: Tidal driving of inertial modes of Maclaurin spheroids

abstract: We examine the inertial modes of vibration of a Maclaurin spheroid and determine how they are driven by an external tidal field, either Newtonian (gravitoelectric) or post-Newtonian (gravitomagnetic). The context and motivation for this work come from the realization that inertial modes of rotating neutron stars can be resonantly excited in binary inspirals, with a measurable effect on the phasing of the emitted gravitational waves. We aim to provide additional insights into this phenomenon. We calculate how the frequencies of the relevant inertial modes, and their overlap integrals with the tidal forces, vary as the star's rotation rate increases, spanning almost the full range between zero rotation and the mass-shedding limit. We prove that a single inertial mode couples to a Newtonian tidal field; overlap integrals with all other modes vanish. We prove also that four inertial modes couple to a post-Newtonian, gravitomagnetic tidal field; overlap integrals with all other modes vanish. Finally, we determine the rather extreme conditions under which the gravitoelectric driving of inertial modes dominates over the gravitomagnetic driving.


title: Carbontracker: Tracking and Predicting the Carbon Footprint of Training Deep Learning Models

abstract: Deep learning (DL) can achieve impressive results across a wide variety of tasks, but this often comes at the cost of training models for extensive periods on specialized hardware accelerators. This energy-intensive workload has seen immense growth in recent years. Machine learning (ML) may become a significant contributor to climate change if this exponential trend continues. If practitioners are aware of their energy and carbon footprint, then they may actively take steps to reduce it whenever possible. In this work, we present Carbontracker, a tool for tracking and predicting the energy and carbon footprint of training DL models. We propose that energy and carbon footprint of model development and training is reported alongside performance metrics using tools like Carbontracker. We hope this will promote responsible computing in ML and encourage research into energy-efficient deep neural networks.


title: Learning to Segment Anatomical Structures Accurately from One Exemplar

abstract: Accurate segmentation of critical anatomical structures is at the core of medical image analysis. The main bottleneck lies in gathering the requisite expert-labeled image annotations in a scalable manner. Methods that permit to produce accurate anatomical structure segmentation without using a large amount of fully annotated training images are highly desirable. In this work, we propose a novel contribution of Contour Transformer Network (CTN), a one-shot anatomy segmentor including a naturally built-in human-in-the-loop mechanism. Segmentation is formulated by learning a contour evolution behavior process based on graph convolutional networks (GCNs). Training of our CTN model requires only one labeled image exemplar and leverages additional unlabeled data through newly introduced loss functions that measure the global shape and appearance consistency of contours. We demonstrate that our one-shot learning method significantly outperforms non-learning-based methods and performs competitively to the state-of-the-art fully supervised deep learning approaches. With minimal human-in-the-loop editing feedback, the segmentation performance can be further improved and tailored towards the observer desired outcomes. This can facilitate the clinician designed imaging-based biomarker assessments (to support personalized quantitative clinical diagnosis) and outperforms fully supervised baselines.


title: Benefiting from Bicubically Down-Sampled Images for Learning Real-World Image Super-Resolution

abstract: Super-resolution (SR) has traditionally been based on pairs of high-resolution images (HR) and their low-resolution (LR) counterparts obtained artificially with bicubic downsampling. However, in real-world SR, there is a large variety of realistic image degradations and analytically modeling these realistic degradations can prove quite difficult. In this work, we propose to handle real-world SR by splitting this ill-posed problem into two comparatively more well-posed steps. First, we train a network to transform real LR images to the space of bicubically downsampled images in a supervised manner, by using both real LR/HR pairs and synthetic pairs. Second, we take a generic SR network trained on bicubically downsampled images to super-resolve the transformed LR image. The first step of the pipeline addresses the problem by registering the large variety of degraded images to a common, well understood space of images. The second step then leverages the already impressive performance of SR on bicubically downsampled images, sidestepping the issues of end-to-end training on datasets with many different image degradations. We demonstrate the effectiveness of our proposed method by comparing it to recent methods in real-world SR and show that our proposed approach outperforms the state-of-the-art works in terms of both qualitative and quantitative results, as well as results of an extensive user study conducted on several real image datasets.


title: An encryption-decryption framework for validating single-particle imaging

abstract: We propose an encryption-decryption framework for validating diffraction intensity volumes reconstructed using single-particle imaging (SPI) with x-ray free-electron lasers (XFELs) when the ground truth volume is absent. This framework exploits each reconstructed volumes' ability to decipher latent variables (e.g. orientations) of unseen sentinel diffraction patterns. Using this framework, we quantify novel measures of orientation disconcurrence, inconsistency, and disagreement between the decryptions by two independently reconstructed volumes. We also study how these measures can be used to define data sufficiency and its relation to spatial resolution, and the practical consequences of focusing XFEL pulses to smaller foci. This framework overcomes critical ambiguities in using Fourier Shell Correlation (FSC) as a validation measure for SPI. Finally, we show how this encryption-decryption framework naturally leads to an information-theoretic reformulation of the resolving power of XFEL-SPI, which we hope will lead to principled frameworks for experiment and instrument design.


title: A Unified Approach to Enhanced Sampling

abstract: The sampling problem lies at the heart of atomistic simulations and over the years many different enhanced sampling methods have been suggested towards its solution. These methods are often grouped into two broad families. On the one hand methods such as umbrella sampling and metadynamics that build a bias potential based on few order parameters or collective variables. On the other hand, tempering methods such as replica exchange that combine different thermodynamic ensembles in one single expanded ensemble. We instead adopt a unifying perspective, focusing on the target probability distribution sampled by the different methods. This allows us to introduce a new class of collective-variables-based bias potentials that can be used to sample any of the expanded ensembles normally sampled via replica exchange. We also provide a practical implementation, by properly adapting the iterative scheme of the recently developed on-the-fly probability enhanced sampling method [Invernizzi and Parrinello, J. Phys. Chem. Lett. 11.7 (2020)], which was originally introduced for metadynamics-like sampling. The resulting method is very general and can be used to achieve different types of enhanced sampling. It is also reliable and simple to use, since it presents only few and robust external parameters and has a straightforward reweighting scheme. Furthermore, it can be used with any number of parallel replicas. We show the versatility of our approach with applications to multicanonical and multithermal-multibaric simulations, thermodynamic integration, umbrella sampling, and combinations thereof.


title: VPN: Learning Video-Pose Embedding for Activities of Daily Living

abstract: In this paper, we focus on the spatio-temporal aspect of recognizing Activities of Daily Living (ADL). ADL have two specific properties (i) subtle spatio-temporal patterns and (ii) similar visual patterns varying with time. Therefore, ADL may look very similar and often necessitate to look at their fine-grained details to distinguish them. Because the recent spatio-temporal 3D ConvNets are too rigid to capture the subtle visual patterns across an action, we propose a novel Video-Pose Network: VPN. The 2 key components of this VPN are a spatial embedding and an attention network. The spatial embedding projects the 3D poses and RGB cues in a common semantic space. This enables the action recognition framework to learn better spatio-temporal features exploiting both modalities. In order to discriminate similar actions, the attention network provides two functionalities - (i) an end-to-end learnable pose backbone exploiting the topology of human body, and (ii) a coupler to provide joint spatio-temporal attention weights across a video. Experiments show that VPN outperforms the state-of-the-art results for action classification on a large scale human activity dataset: NTU-RGB+D 120, its subset NTU-RGB+D 60, a real-world challenging human activity dataset: Toyota Smarthome and a small scale human-object interaction dataset Northwestern UCLA.


title: Approximation algorithms for car-sharing problems

abstract: We consider several variants of a car-sharing problem. Given are a number of requests each consisting of a pick-up location and a drop-off location, a number of cars, and nonnegative, symmetric travel times that satisfy the triangle inequality. Each request needs to be served by a car, which means that a car must first visit the pick-up location of the request, and then visit the drop-off location of the request. Each car can serve two requests. One problem is to serve all requests with the minimum total travel time (called $\CS_{sum}$), and the other problem is to serve all requests with the minimum total latency (called $\CS_{lat}$). We also study the special case where the pick-up and drop-off location of a request coincide. We propose two basic algorithms, called the match and assign algorithm and the transportation algorithm. We show that the best of the resulting two solutions is a $ 2$-approximation for $\CS_{sum}$ (and a $7/5$-approximation for its special case), and a $5/3 $-approximation for $\CS_{lat}$ (and a $3/2$-approximation for its special case); these ratios are better than the ratios of the individual algorithms. Finally, we indicate how our algorithms can be applied to more general settings where each car can serve more than two requests, or where cars have distinct speeds.


title: Proposal on Application of the Multi-Wire Proportional Chambers of the LHCb MUON Detector at Very High Rates for the Future Upgrades

abstract: The MUON Detector (MD) of LHCb is one of the largest instruments of this kind worldwide, and one of the most irradiated. It has performed exceptionally well during the RUN1 and RUN2 of the LHC at an instantaneous luminosity of 4$\times$10$^{32}$ cm$^{-2}$s$^{-1}$, with tracking inefficiencies at the level of 1$\%$ and 2.6$\%$, respectively. Looking forward for the future LHCb Upgrade 2 (U2) planned in 2031 and aiming in running the detector at increased luminosity by a factor $\sim$50, and at the same time keeping a very high ($\sim$99$\%$) detection efficiency, an option with reuse significant part of the present Multi-Wire Proportional Chambers (MWPC) in a new Muon System is presented. In addition, the first idea of new Front End Electronics (FEE) and an existing test setup applicable for designing both: new MWPCs with a higher granularity of the cathode readout pads and new FEE are described.


title: Guided Fine-Tuning for Large-Scale Material Transfer

abstract: We present a method to transfer the appearance of one or a few exemplar SVBRDFs to a target image representing similar materials. Our solution is extremely simple: we fine-tune a deep appearance-capture network on the provided exemplars, such that it learns to extract similar SVBRDF values from the target image. We introduce two novel material capture and design workflows that demonstrate the strength of this simple approach. Our first workflow allows to produce plausible SVBRDFs of large-scale objects from only a few pictures. Specifically, users only need take a single picture of a large surface and a few close-up flash pictures of some of its details. We use existing methods to extract SVBRDF parameters from the close-ups, and our method to transfer these parameters to the entire surface, enabling the lightweight capture of surfaces several meters wide such as murals, floors and furniture. In our second workflow, we provide a powerful way for users to create large SVBRDFs from internet pictures by transferring the appearance of existing, pre-designed SVBRDFs. By selecting different exemplars, users can control the materials assigned to the target image, greatly enhancing the creative possibilities offered by deep appearance capture.


title: Perverse sheaves and finite-dimensional algebras

abstract: Let $X$ be a topologically stratified space, $p$ be any perversity on $X$, and $k$ be a field. We show that the category of $p$-perverse sheaves on $X$, constructible with respect to the stratification and with coefficients in $k$, is equivalent to the category of finite-dimensional modules over a finite-dimensional algebra if and only if $X$ has finitely many strata and the same holds for the category of local systems on each of these. The main component in the proof is a construction of projective covers for simple perverse sheaves.


title: Four-fold anisotropy of the parallel upper critical magnetic field in a pure layered d-wave superconductor at $T = 0$

abstract: It is well known that a four-fold symmetry of the parallel upper critical magnetic field disappears in the Ginzburg-Landau (GL) region in quasi-two-dimensional (Q2D) $d$-wave superconductors. Therefore, it has been accurately calculated so far as a correction to the GL results, which is valid close to superconducting transition temperature and is expected to be stronger at low temperatures. As to the case $T=0$, some approximated methods have been used, which are good only for closed electron orbits and unappropriate for the open orbits which exist in a parallel magnetic field in Q2D superconductors. For the first time, we accurately calculate the four-fold anisotropy of the parallel upper critical magnetic field in a pure Q2D $d$-wave superconductor at $T=0$, where it has the highest possible value. Our results are applicable to Q2D $d$-wave high-Tc and organic superconductors.


title: Fast convex optimization via a third-order in time evolution equation: TOGES-V an improved version of TOGES

abstract: In a Hilbert space setting H, for convex optimization, we analyze the fast convergence properties as t tends to infinity of the trajectories generated by a third-order in time evolution system. The function f to minimize is supposed to be convex, continuously differentiable, with a nonempty set of minimizers. It enters into the dynamic through its gradient. Based on this new dynamical system, we improve the results obtained by [Attouch, Chbani, Riahi: Fast convex optimization via a third-order in time evolution equation, Optimization 2020]. As a main result, when the damping parameter $\alpha$ satisfies $\alpha > 3$, we show that the convergence of the values at the order 1/t3 as t goes to infinity, as well as the convergence of the trajectories. We complement these results by introducing into the dynamic an Hessian driven damping term, which reduces the oscillations. In the case of a strongly convex function f, we show an autonomous evolution system of the third order in time with an exponential rate of convergence. All these results have natural extensions to the case of a convex lower semicontinuous function with extended real values. Just replace f with its Moreau envelope.


title: ARC-Net: Activity Recognition Through Capsules

abstract: Human Activity Recognition (HAR) is a challenging problem that needs advanced solutions than using handcrafted features to achieve a desirable performance. Deep learning has been proposed as a solution to obtain more accurate HAR systems being robust against noise. In this paper, we introduce ARC-Net and propose the utilization of capsules to fuse the information from multiple inertial measurement units (IMUs) to predict the activity performed by the subject. We hypothesize that this network will be able to tune out the unnecessary information and will be able to make more accurate decisions through the iterative mechanism embedded in capsule networks. We provide heatmaps of the priors, learned by the network, to visualize the utilization of each of the data sources by the trained network. By using the proposed network, we were able to increase the accuracy of the state-of-the-art approaches by 2%. Furthermore, we investigate the directionality of the confusion matrices of our results and discuss the specificity of the activities based on the provided data.


title: Maximizing five-cycles in $K_r$-free graphs

abstract: The Erd\H{o}s Pentagon problem asks to find an $n$-vertex triangle-free graph that is maximizing the number of $5$-cycles. The problem was solved using flag algebras by Grzesik and independently by Hatami, Hladk\'{y}, Kr\'{a}l', Norin, and Razborov. Recently, Palmer suggested the general problem of maximizing the number of $5$-cycles in $K_{k+1}$-free graphs. Using flag algebras, we show that every $K_{k+1}$-free graph of order $n$ contains at most \[\frac{1}{10k^4}(k^4 - 5k^3 + 10k^2 - 10k + 4)n^5 + o(n^5)\] copies of $C_5$ for any $k \geq 3$, with the Tur\'an graph begin the extremal graph for large enough $n$.


title: A Boolean Control Network Approach to the Formal Verification of Feedback Context-Aware Pervasive Systems

abstract: The emergence of Context-aware systems in the domains of autonomic, monitoring, and safety-critical applications asks for the definition of methods to formally assess their correctness and dependability properties. Many of these properties are common to Automatic Control systems, a field that developed well established analysis and design techniques to formalize and investigate them. In this paper, we use Boolean Control Networks, to discuss some properties of a feedback Context-aware system in a case study based on a healthcare management example.


title: Ultrafast viscosity measurement with ballistic optical tweezers

abstract: Viscosity is an important property of out-of-equilibrium systems such as active biological materials and driven non-Newtonian fluids, and for fields ranging from biomaterials to geology, energy technologies and medicine. However, noninvasive viscosity measurements typically require integration times of seconds. Here we demonstrate a four orders-of-magnitude improvement in speed, down to twenty microseconds, with uncertainty dominated by fundamental thermal noise for the first time. We achieve this using the instantaneous velocity of a trapped particle in an optical tweezer. To resolve the instantaneous velocity we develop a structured-light detection system that allows particle tracking with megahertz bandwidths. Our results translate viscosity from a static averaged property, to one that may be dynamically tracked on the timescales of active dynamics. This opens a pathway to new discoveries in out-of-equilibrium systems, from the fast dynamics of phase transitions, to energy dissipation in motor molecule stepping, to violations of fluctuation laws of equilibrium thermodynamics.

categories: cond-mat.soft physics.bio-ph physics.ins-det physics.optics quant-ph


title: A Physics-informed Assembly of Feed-Forward Neural Network Engines to Predict Inelasticity in Cross-Linked Polymers

abstract: In solid mechanics, Data-driven approaches are widely considered as the new paradigm that can overcome the classic problems of constitutive models such as limiting hypothesis, complexity, and high dependence on training data. However, implementation of machine-learned approaches in material modeling has been modest due to the high-dimensionality of the data space, significant size of missing data, and limited convergence. This work proposes a framework to hire concepts from polymer science, statistical physics, and continuum mechanics to provide super-constrained machine-learning techniques of reduced-order to overcome many of the existing difficulties. Using a sequential order-reduction, we have simplified the 3D stress-strain tensor mapping problem into a limited number of super-constrained 1D mapping problems. Next, we introduce an assembly of multiple replicated Neural Network agents to systematically classify those mapping problems into a few categories, all of which are replications of a few distinct agent types. By capturing all loading modes through a simplified set of disperse experimental data, the proposed hybrid assembly of agents provides a new generation of machine learned approaches that simply outperforms most constitutive laws in training data volume, training speed, and accuracy even in complicated loading scenarios. Also, it avoids low interpretability of conventional AI-based models.


title: Dielectrophoretic Equilibrium of Complex Particles

abstract: In contrast to the commonly used spherical Janus particles, here we used engineered Janus particles that are fabricated using photolithography technique for precise control over their geometry and coated regions. Specifically, we studied a lollipop-shaped complex particle where its head is coated with gold while its tail is left bare. Due to their distinct electrical properties (i.e. electrical polarizability) the particle exhibits force equilibrium where opposite dielectrophoretic forces acting on its head and tail exactly cancel each other to yield a stable equilibrium position. This was realized in a quadrupolar electrode array where the equilibrium position of the engineered particle could be tuned by the frequency. This stands in contrast to the standard dielectrophoretic behavior where the particle shifts positions from either the center of the quad to the very edge of the electrodes when shifting from a negative to positive dielectrophoretic response, respectively. This opens new opportunities for positioning control of such complex particles for self-assembly, biosensing, biomimetic spermatozoa and more.


title: Outcome-Driven Dynamic Refugee Assignment with Allocation Balancing

abstract: This study proposes two new dynamic assignment algorithms to match refugees and asylum seekers to geographic localities within a host country. The first, currently implemented in a multi-year randomized control trial in Switzerland, seeks to maximize the average predicted employment level (or any measured outcome of interest) of refugees through a minimum-discord online assignment algorithm. The performance of this algorithm is tested on real refugee resettlement data from both the US and Switzerland, where we find that it is able to achieve near-optimal expected employment compared to the hindsight-optimal solution, and is able to improve upon the status quo procedure by 40-50%. However, pure outcome maximization can result in a periodically imbalanced allocation to the localities over time, leading to implementation difficulties and an undesirable workflow for resettlement resources and agents. To address these problems, the second algorithm balances the goal of improving refugee outcomes with the desire for an even allocation over time. We find that this algorithm can achieve near-perfect balance over time with only a small loss in expected employment compared to the employment-maximizing algorithm. In addition, the allocation balancing algorithm offers a number of ancillary benefits compared to pure outcome maximization, including robustness to unknown arrival flows and greater exploration.


title: Novel current actuated piezoelectric composite model with fully dynamic electromagnetic field

abstract: In this paper we propose a novel current actuated fully dynamic piezoelectric composite model. This new model complements the existing modeling framework of deriving piezoelectric models. We show that the novel composite model is well-posed. Furthermore, we consider two approximation methods, FEM and MFEM and investigate the stabilizability properties. Finally, we include the closed-loop behaviour in the numerical results.


title: Deep Partial Updating: Towards Communication Efficient Updating for On-device Inference

abstract: Emerging edge intelligence applications require the server to retrain and update deep neural networks deployed on remote edge nodes to leverage newly collected data samples. Unfortunately, it may be impossible in practice to continuously send fully updated weights to these edge nodes due to the highly constrained communication resource. In this paper, we propose the weight-wise deep partial updating paradigm, which smartly selects a small subset of weights to update in each server-to-edge communication round, while achieving a similar performance compared to full updating. Our method is established through analytically upper-bounding the loss difference between partial updating and full updating, and only updates the weights which make the largest contributions to the upper bound. Extensive experimental results demonstrate the efficacy of our partial updating methodology which achieves a high inference accuracy while updating a rather small number of weights.


title: On real solutions of the Dirac equation for a one-dimensional Majorana particle

abstract: We construct general solutions of the time-dependent Dirac equation in (1+1) dimensions with a Lorentz scalar potential, subject to the so-called Majorana condition, in the Majorana representation. In this situation, these solutions are real-valued and describe a one-dimensional Majorana single particle. We specifically obtain solutions for the following cases: a Majorana particle at rest inside a box, a free (i.e., in a penetrable box with the periodic boundary condition), in an impenetrable box with no potential (here we only have four boundary conditions), and in a linear potential. All these problems are treated in a very detailed and systematic way. In addition, we obtain and discuss various results related to real wave functions. Finally, we also wish to point out that, in choosing the Majorana representation, the solutions of the Dirac equation with a Lorentz scalar potential can be chosen to be real but do not need to be real. In fact, complex solutions for this equation can also be obtained. Thus, a Majorana particle cannot be described only with the Dirac equation in the Majorana representation without explicitly imposing the Majorana condition.


title: Generative Model-Based Loss to the Rescue: A Method to Overcome Annotation Errors for Depth-Based Hand Pose Estimation

abstract: We propose to use a model-based generative loss for training hand pose estimators on depth images based on a volumetric hand model. This additional loss allows training of a hand pose estimator that accurately infers the entire set of 21 hand keypoints while only using supervision for 6 easy-to-annotate keypoints (fingertips and wrist). We show that our partially-supervised method achieves results that are comparable to those of fully-supervised methods which enforce articulation consistency. Moreover, for the first time we demonstrate that such an approach can be used to train on datasets that have erroneous annotations, i.e. "ground truth" with notable measurement errors, while obtaining predictions that explain the depth images better than the given "ground truth".


title: Kernel Stein Generative Modeling

abstract: We are interested in gradient-based Explicit Generative Modeling where samples can be derived from iterative gradient updates based on an estimate of the score function of the data distribution. Recent advances in Stochastic Gradient Langevin Dynamics (SGLD) demonstrates impressive results with energy-based models on high-dimensional and complex data distributions. Stein Variational Gradient Descent (SVGD) is a deterministic sampling algorithm that iteratively transports a set of particles to approximate a given distribution, based on functional gradient descent that decreases the KL divergence. SVGD has promising results on several Bayesian inference applications. However, applying SVGD on high dimensional problems is still under-explored. The goal of this work is to study high dimensional inference with SVGD. We first identify key challenges in practical kernel SVGD inference in high-dimension. We propose noise conditional kernel SVGD (NCK-SVGD), that works in tandem with the recently introduced Noise Conditional Score Network estimator. NCK is crucial for successful inference with SVGD in high dimension, as it adapts the kernel to the noise level of the score estimate. As we anneal the noise, NCK-SVGD targets the real data distribution. We then extend the annealed SVGD with an entropic regularization. We show that this offers a flexible control between sample quality and diversity, and verify it empirically by precision and recall evaluations. The NCK-SVGD produces samples comparable to GANs and annealed SGLD on computer vision benchmarks, including MNIST and CIFAR-10.


title: A Term-Rewriting Semantics for Imperative Style Programming

abstract: Term rewriting systems have a simple syntax and semantics and facilitate proofs of correctness. However, they are not as popular in industry or academia as imperative languages. We define a term rewriting based abstract programming language with an imperative style and a precise semantics allowing programs to be translatable into efficient imperative languages, to obtain proofs of correctness together with efficient execution. This language is designed to facilitate translations into correct programs in imperative languages with assignment statements, iteration, recursion, arrays, pointers, and side effects. It can also be used in place of a pseudo-programming language to specify algorithms.


title: Modeling R Coronae Borealis Stars: Effects of He-Burning Shell Temperature and Metallicity

abstract: The R Coronae Borealis (RCB) stars are extremely hydrogen-deficient carbon stars which produce large amounts of dust, causing sudden deep declines in brightness. They are believed to be formed primarily through white dwarf mergers. In this paper, we use MESA to investigate how post-merger objects with a range of initial He-burning shell temperatures from 2.1 - 5.4 *10^8 K with solar and subsolar metallicities evolve into RCB stars. The most successful model of these has subsolar metallicity and an initial temperature near 3 *10^8 K. We find a strong dependence on initial He-burning shell temperature for surface abundances of elements involved in the CNO cycle, as well as differences in effective temperature and radius of RCBs. Elements involved in nucleosynthesis present around 1 dex diminished surface abundances in the 10% solar metallicity models, with the exception of carbon and lithium which are discussed in detail. Models with subsolar metallicities also exhibit longer lifetimes than their solar counterparts. Additionally, we find that convective mixing of the burned material occurs only in the first few years of post-merger evolution, after which the surface abundances are constant during and after the RCB phase, providing evidence for why these stars show a strong enhancement of partial He-burning products.


title: Transient compressible flow in a compliant viscoelastic tube

abstract: Motivated by problems arising in the pneumatic actuation of controllers for micro-electromechanical systems (MEMS), labs-on-a-chip or biomimetic soft robots, and the study of microrheology of both gases and soft solids, we analyze the transient fluid--structure interaction (FSIs) between a viscoelastic tube conveying compressible flow at low Reynolds number. We express the density of the fluid as a linear function of the pressure, and we use the lubrication approximation to further simplify the fluid dynamics problem. On the other hand, the structural mechanics is governed by a modified Donnell shell theory accounting for Kelvin--Voigt-type linearly viscoelastic mechanical response. The fluid and structural mechanics problems are coupled through the tube's radial deformation and the hydrodynamic pressure. For small compressibility numbers and weak coupling, the equations are solved analytically via a perturbation expansion. Three illustrative problems are analyzed. First, we obtain exact (but implicit) solutions for the pressure for steady flow conditions. Second, we solve the transient problem of impulsive pressurization of the tube's inlet. Third, we analyze the transient response to an oscillatory inlet pressure. We show that an oscillatory inlet pressure leads to acoustic streaming in the tube, attributed to the nonlinear pressure gradient induced by the interplay of FSI and compressibility. Furthermore, we demonstrate an enhancement in the volumetric flow rate due to FSI coupling. The hydrodynamic pressure oscillations are shown to exhibit a low-pass frequency response (when averaging over the period of oscillations), while the frequency response of the tube deformation is similar to that of a band-pass filter.


title: A direct approach to nonuniqueness and failure of compactness for the SQG equation

abstract: We give an alternative proof of the nonuniqueness of weak solutions to the surface quasigeostrophic equation (SQG) first shown in [Buckmaster-Shkoller-Vicol, '16]. Our approach proceeds directly at the level of the scalar field. Furthermore, we prove that every smooth scalar field with compact support that conserves the integral can be realized as a weak limit of solutions to SQG.


title: A search for rotation periods in 1000 TESS objects of interest

abstract: The high quality light curves from the Transiting Exoplanet Survey Satellite (TESS) represent a unique laboratory for the study of stellar rotation, a fundamental observable driving stellar and planetary evolution, including planetary atmospheres and impacting on habitability conditions and the genesis of life around stars. As of April 14th 2020, this mission delivered public light curves for 1000 TESS Objects of Interest (TOIs), observed with 2 minute cadence during the first 20 months of the mission. Here, we present a search for rotation signatures in these TOIs, using Fast Fourier Transform, Lomb-Scargle, and wavelet techniques, accompanied by a rigorous visual inspection. This effort revealed 163 targets with rotation signatures, 131 of which present unambiguous rotation periods ranging from 0.321 and 13.219 days, whereas 32 of them present dubious rotation periodicities. One hundred and nine of these stars show flux fluctuations whose root-cause is not clearly identified. For 714 TOIs, the light curves show a noisy behavior, corresponding to typically low-amplitude signals. Our analysis has also revealed 10 TOI stars with pulsation periodicities ranging from 0.049 to 2.995 days and four eclipsing binaries. With upcoming TESS data releases, our periodicity analysis will be expanded to almost all TOI stars, thereby contributing in defining criteria for follow-up strategy itself, and the study of star-planet interactions, surface dynamic of host stars and habitability conditions in planets, among other aspects. In this context, a living catalog is maintained on the Filtergraph visualization portal at the URL https://filtergraph.com/tess_rotation_tois.


title: Comultiplication in the Serre Spectral Sequence

abstract: We show the homological Serre spectral sequence with coefficients in a field is a spectral sequence of coalgebras. We also identify the comultiplication on the $E^2$ page of the spectral sequence as being induced by the usual comultiplication in homology. At the end, we provide some example computations highlighting the use the co-Leibniz rule.