title: Measuring Adolescents' Well-being: Correspondence of Naive Digital Traces to Survey Data

abstract: Digital traces are often used as a substitute for survey data. However, it is unclear whether and how digital traces actually correspond to the survey-based traits they purport to measure. This paper examines correlations between self-reports and digital trace proxies of depression, anxiety, mood, social integration and sleep among high school students. The study is based on a small but rich multilayer data set (N = 144). The data set contains mood and sleep measures, assessed daily over a 4-month period, along with survey measures at two points in time and information about online activity from VK, the most popular social networking site in Russia. Our analysis indicates that 1) the sentiments expressed in social media posts are correlated with depression; namely, adolescents with more severe symptoms of depression write more negative posts, 2) late-night posting indicates less sleep and poorer sleep quality, and 3) students who were nominated less often as somebody's friend in the survey have fewer friends on VK and their posts receive fewer "likes." However, these correlations are generally weak. These results demonstrate that digital traces can serve as useful supplements to, rather than substitutes for, survey data in studies on adolescents' well-being. These estimates of correlations between survey and digital trace data could provide useful guidelines for future research on the topic.


title: Energy dependent morphology of the pulsar wind nebula HESS J1825-137 with Fermi-LAT

abstract: Taking advantage of more than 11 years of Fermi-LAT data, we perform a new and deep analysis of the pulsar wind nebula (PWN) HESS J1825-137. Combining this analysis with recent H.E.S.S. results we investigate and constrain the particle transport mechanisms at work inside the source as well as the system evolution. The PWN is studied using 11.6 years of Fermi-LAT data between 1 GeV and 1 TeV. In particular, we present the results of the spectral analysis and the first energy-resolved morphological study of the PWN HESS J1825-137 at GeV energies, which provide new insights into the gamma-ray characteristics of the nebula. An optimised analysis of the source returns an extended emission region larger than 2$^{\circ}$, corresponding to an intrinsic size of about 150 pc, making HESS J1825-137 the most extended gamma-ray PWN currently known. The nebula presents a strong energy dependent morphology within the GeV range, moving from a radius of $\sim1.4^\circ$ below 10 GeV to a radius of $\sim$0.8$^\circ$ above 100 GeV, with a shift in the centroid location. Thanks to the large extension and peculiar energy-dependent morphology, it is possible to constrain the particle transport mechanisms inside the PWN HESS J1825-137. Using the variation of the source extension and position, as well as the constraints on the particle transport mechanisms, we present a scheme for the possible evolution of the system. Finally, we provide an estimate of the electron energy density and we discuss its nature in the PWN and TeV halo-like scenario.


title: Stability of solutions for a parabolic problem involving fractional p-Laplacian with logarithmic nonlinearity

abstract: In this paper, we study the following Dirichlet problem for a parabolic equation involving fractional $p$-Laplacian with logarithmic nonlinearity \begin{equation*}\label{eq}\left\{ \begin{array}{llc} u_{t}+(-\Delta)^{s}_{p}u+|u|^{p-2}u=|u|^{p-2}u\log(|u|) & \text{in}\ & \Omega,\;t>0 , u =0 & \text{in} & \mathbb{R}^{N}\backslash \Omega,\;t > 0, u(x,0)=u_{0}(x), & \text{in} &\Omega , \end{array}\right. \end{equation*} where $\Omega \subset \mathbb{R}^N \, ( N\geq 1)$ is a bounded domain with Lipschitz boundary and $2\leq p< \infty$. The local existence will be done by using the Galerkin approximations. By combining the potential well theory with the Nehari manifold we establish the existence of global solutions. Then, by virtue of a differential inequality technique, we prove that the local solutions blow-up in finite time with arbitrary negative initial energy and suitable initial values. Moreover, we give decay estimates of global solutions. The main difficulty here is the lack of logarithmic Sobolev inequality concerning fractional $p$-Laplacian.


title: Rapid late-time X-ray brightening of the tidal disruption event OGLE16aaa

abstract: Stars that pass too close to a super-massive black hole may be disrupted by strong tidal forces. OGLE16aaa is one such tidal disruption event (TDE) which rapidly brightened and peaked in the optical/UV bands in early 2016 and subsequently decayed over the rest of the year. OGLE16aaa was detected in an XMM-Newton X-ray observation on June 9, 2016 with a flux slightly below the Swift/XRT upper limits obtained during the optical light curve peak. Between June 16-21, 2016, Swift/XRT also detected OGLE16aaa and based on the stacked spectrum, we could infer that the X-ray luminosity had jumped up by more than a factor of ten in just one week. No brightening signal was seen in the simultaneous optical/UV data to cause the X-ray luminosity to exceed the optical/UV one. A further XMM-Newton observation on November 30, 2016 showed that almost a year after the optical/UV peak, the X-ray emission was still at an elevated level, while the optical/UV flux decay had already leveled off to values comparable to those of the host galaxy. In all X-ray observations, the spectra were nicely modeled with a 50-70 eV thermal component with no intrinsic absorption, with a weak X-ray tail seen only in the November 30 XMM-Newton observation. The late-time X-ray behavior of OGLE16aaa strongly resembles the tidal disruption events ASASSN-15oi and AT2019azh. We were able to pinpoint the time delay between the initial optical TDE onset and the X-ray brightening to $182 \pm 5$ days, which may possibly represent the timescale between the initial circularization of the disrupted star around the super-massive black hole and the subsequent delayed accretion. Alternatively, the delayed X-ray brightening could be related to a rapid clearing of a thick envelope that covers the central X-ray engine during the first six months.


title: A Mini-Neptune and a Venus-Zone Planet in the Radius Valley Orbiting the Nearby M2-dwarf TOI-1266: Validation with the Habitable-zone Planet Finder

abstract: We report on the validation of two planets orbiting the nearby (36pc) M2 dwarf TOI-1266 observed by the TESS mission. The inner planet is sub-Neptune-sized ($R=2.46 \pm 0.08 R_\oplus$) with an orbital period of 10.9 days. The outer planet has a radius of $1.67_{-0.11}^{+0.09} R_\oplus$ and resides in the exoplanet Radius Valley---the transition region between rocky and gaseous planets. With an orbital period of 18.8 days, the outer planet receives an insolation flux of 2.4 times that of Earth, similar to the insolation of Venus. Using precision near-infrared radial velocities with the Habitable-zone Planet Finder Spectrograph, we place upper mass limits of $15.9 M_\oplus$ and $6.4 M_\oplus$ at 95% confidence for the inner and outer planet, respectively. A more precise mass constraint of planet c, achievable with current RV instruments given the host star brightness (V=12.9, J=9.7), will yield further insights into the dominant processes sculpting the exoplanet Radius Valley.


title: Improving the accuracy of quantum computational chemistry using the transcorrelated method

abstract: Accurately treating electron correlation in the wavefunction is a key challenge for both classical and quantum computational chemistry. Classical methods have been developed which explicitly account for this correlation by incorporating inter-electronic distances into the wavefunction. The transcorrelated method transfers this explicit correlation from the wavefunction to a transformed, non-Hermitian Hamiltonian, whose right-hand eigenvectors become easier to obtain than those of the original Hamiltonian. In this work, we show that the transcorrelated method can reduce the resources required to obtain accurate energies from electronic structure calculations on quantum computers. We overcome the limitations introduced by the non-Hermitian Hamiltonian by using quantum algorithms for imaginary time evolution.


title: $\lambda$-Regularized A-Optimal Design and its Approximation by $\lambda$-Regularized Proportional Volume Sampling

abstract: In this work, we study the $\lambda$-regularized $A$-optimal design problem and introduce the $\lambda$-regularized proportional volume sampling algorithm, generalized from [Nikolov, Singh, and Tantipongpipat, 2019], for this problem with the approximation guarantee that extends upon the previous work. In this problem, we are given vectors $v_1,\ldots,v_n\in\mathbb{R}^d$ in $d$ dimensions, a budget $k\leq n$, and the regularizer parameter $\lambda\geq0$, and the goal is to find a subset $S\subseteq [n]$ of size $k$ that minimizes the trace of $\left(\sum_{i\in S}v_iv_i^\top + \lambda I_d\right)^{-1}$ where $I_d$ is the $d\times d$ identity matrix. The problem is motivated from optimal design in ridge regression, where one tries to minimize the expected squared error of the ridge regression predictor from the true coefficient in the underlying linear model. We introduce $\lambda$-regularized proportional volume sampling and give its polynomial-time implementation to solve this problem. We show its $(1+\frac{\epsilon}{\sqrt{1+\lambda'}})$-approximation for $k=\Omega\left(\frac d\epsilon+\frac{\log 1/\epsilon}{\epsilon^2}\right)$ where $\lambda'$ is proportional to $\lambda$, extending the previous bound in [Nikolov, Singh, and Tantipongpipat, 2019] to the case $\lambda>0$ and obtaining asymptotic optimality as $\lambda\rightarrow \infty$.


title: Evaluation Of Hidden Markov Models Using Deep CNN Features In Isolated Sign Recognition

abstract: Isolated sign recognition from video streams is a challenging problem due to the multi-modal nature of the signs, where both local and global hand features and face gestures needs to be attended simultaneously. This problem has recently been studied widely using deep Convolutional Neural Network (CNN) based features and Long Short-Term Memory (LSTM) based deep sequence models. However, the current literature is lack of providing empirical analysis using Hidden Markov Models (HMMs) with deep features. In this study, we provide a framework that is composed of three modules to solve isolated sign recognition problem using different sequence models. The dimensions of deep features are usually too large to work with HMM models. To solve this problem, we propose two alternative CNN based architectures as the second module in our framework, to reduce deep feature dimensions effectively. After extensive experiments, we show that using pretrained Resnet50 features and one of our CNN based dimension reduction models, HMMs can classify isolated signs with 90.15% accuracy in Montalbano dataset using RGB and Skeletal data. This performance is comparable with the current LSTM based models. HMMs have fewer parameters and can be trained and run on commodity computers fast, without requiring GPUs. Therefore, our analysis with deep features show that HMMs could also be utilized as well as deep sequence models in challenging isolated sign recognition problem.


title: Poisson Learning: Graph Based Semi-Supervised Learning At Very Low Label Rates

abstract: We propose a new framework, called Poisson learning, for graph based semi-supervised learning at very low label rates. Poisson learning is motivated by the need to address the degeneracy of Laplacian semi-supervised learning in this regime. The method replaces the assignment of label values at training points with the placement of sources and sinks, and solves the resulting Poisson equation on the graph. The outcomes are provably more stable and informative than those of Laplacian learning. Poisson learning is efficient and simple to implement, and we present numerical experiments showing the method is superior to other recent approaches to semi-supervised learning at low label rates on MNIST, FashionMNIST, and Cifar-10. We also propose a graph-cut enhancement of Poisson learning, called Poisson MBO, that gives higher accuracy and can incorporate prior knowledge of relative class sizes.

categories: cs.LG cs.CV cs.NA math.AP math.NA stat.ML


title: Desorption energy of soft particles from a fluid interface

abstract: The efficiency of soft particles to stabilize emulsions is examined by measuring their desorption free energy, i.e., the mechanical work required to detach the particle from a fluid interface. Here, we consider rubber-like elastic as well as microgel particles, using coarse-grained molecular dynamics simulations. The energy of desorption is computed for two and three-dimensional configurations by means of the mean thermodynamic integration method. It is shown that the softness affects the particle-interface binding in two opposing directions as compared to rigid particles. On the one hand, a soft particle spreads at the interface and thereby removes a larger unfavorable liquid-liquid contact area compared to rigid particles. On the other hand, softness provides the particle with an additional degree of freedom to get reshaped instead of deforming the interface, resulting in a smaller restoring force during the detachment. It is shown that the first effect prevails so that a soft spherical particle attaches to the fluid interface more strongly than rigid spheres. Finally, we consider microgel particles both in the swollen and in the collapsed state. Surprisingly, we find that the latter has a larger binding energy. All results are rationalised using thermodynamic arguments and thereby offer detailed insights into the desorption energy of soft particles from fluid interfaces.


title: Operator evolution from the similarity renormalization group and the Magnus expansion

abstract: The Magnus expansion is an efficient alternative to solving similarity renormalization group (SRG) flow equations with high-order, memory-intensive ordinary differential equation solvers. The numerical simplifications it offers for operator evolution are particularly valuable for in-medium SRG calculations, though challenges remain for difficult problems involving intruder states. Here we test the Magnus approach in an analogous but more accessible situation, which is the free-space SRG treatment of the spurious bound-states arising from a leading-order chiral effective field theory (EFT) potential with very high cutoffs. We show that the Magnus expansion passes these tests and then use the investigations as a springboard to address various aspects of operator evolution that have renewed relevance in the context of the scale and scheme dependence of nuclear processes. These aspects include SRG operator flow with band- versus block-diagonal generators, universality for chiral EFT Hamiltonians and associated operators with different regularization schemes, and the impact of factorization arising from scale separation. Implications for short-range correlations physics and the possibilities for reconciling high- and low-resolution treatments of nuclear structure and reactions are discussed.


title: The Hermitian Jacobi process: simplified formula for the moments and application to optical fibers MIMO channels

abstract: Using a change of basis in the algebra of symmetric functions, we compute the moments of the Hermitian Jacobi process. After a careful arrangement of the terms and the evaluation of the determinant of an `almost upper-triangular' matrix, we end up with a moment formula which is considerably simpler than the one derived in \cite{Del-Dem}. As an application, we propose the Hermitian Jacobi process as a dynamical model for optical fibers MIMO channels and compute its Shannon capacity for small enough power at the transmitter. Moreover, when the size of the Hermitian Jacobi process is larger than the moment order, our moment formula may be written as a linear combination of balanced terminating ${}_4F_3$-series evaluated at unit argument.


title: Dynamic symmetry-breaking in mutually annihilating fluids with selective interfaces

abstract: The selective entrapment of mutually annihilating species within a phase-changing carrier fluid is explored by both analytical and numerical means. The model takes full account of the dynamic heterogeneity which arises as a result of the coupling between hydrodynamic transport, dynamic phase-transitions and chemical reactions between the participating species, in the presence of a selective droplet interface. Special attention is paid to the dynamic symmetry breaking between the mass of the two species entrapped within the expanding droplet as a function of time. It is found that selective sources are much more effective symmetry breakers than selective diffusion. The present study may be of interest for a broad variety of advection-diffusion-reaction phenomena with selective fluid interfaces, including the problem of electroweak baryogenesis.


title: Cut cotorsion pairs

abstract: We present the concept of cotorsion pairs cut along subcategories of an abelian category. This provides a generalization of complete cotorsion pairs, and represents a general framework to find approximations restricted to certain subcategories. We also exhibit some connections between cut cotorsion pairs and Auslander-Buchweitz approximation theory, by considering relative analogs for Frobenius pairs and Auslander-Buchweitz contexts. Several applications are given in the settings of relative Gorenstein homological algebra, chain complexes and quasi-coherent sheaves, but also to characterize some important results on the Finitistic Dimension Conjecture, the existence of right adjoints of quotient functors by Serre subcategories, and the description of cotorsion pairs in triangulated categories as co-$t$-structures.


title: Solving optimization problems with Rydberg analog quantum computers: Realistic requirements for quantum advantage using noisy simulation and classical benchmarks

abstract: Platforms of Rydberg atoms have been proposed as promising candidates to solve some combinatorial optimization problems. Here, we compute quantitative requirements on the system sizes and noise levels that these platforms must fulfill to reach quantum advantage in approximately solving the Unit-Disk Maximum Independent Set problem. Using noisy simulations of Rydberg platforms of up to 26 atoms interacting through realistic van der Waals interactions, we compute the average approximation ratio that can be attained with a simple quantum annealing-based heuristic within a fixed temporal computational budget. Based on estimates of the correlation lengths measured in the engineered quantum state, we extrapolate the results to large atom numbers and compare them to a simple classical approximation heuristic. We find that approximation ratios of at least $\approx 0.84$ are within reach for near-future noise levels. Not taking into account further classical and quantum algorithmic improvements, we estimate that quantum advantage could be reached by attaining a number of controlled atoms of $\sim8,000$ for a time budget of 2 seconds, and $\sim 1,000-1,200$ for a time budget of 0.2 seconds, provided the coherence levels of the system can be improved by a factor 10 while maintaining a constant repetition rate.


title: Observation of the production of three massive gauge bosons at $\sqrt{s} =$ 13 TeV

abstract: The first observation is reported of the combined production of three massive gauge bosons (VVV with V = W,Z) in proton-proton collisions at a center-of-mass energy of 13 TeV. The analysis is based on a data sample recorded by the CMS experiment at the CERN LHC corresponding to an integrated luminosity of 137 fb$^{-1}$. The searches for individual WWW, WWZ, WZZ, and ZZZ production are performed in final states with three, four, five, and six leptons (electrons or muons), or with two same-sign leptons plus one or two jets. The observed (expected) significance of the combined VVV production signal is 5.7 (5.9) standard deviations and the corresponding measured cross section relative to the standard model prediction is 1.02 $^{+0.26}_{-0.23}$. The significances of the individual WWW and WWZ production are 3.3 and 3.4 standard deviations, respectively. Measured production cross sections for the individual triboson processes are also reported.


title: Modeling and Optimal Control of Hybrid UAVs with Wind Disturbance

abstract: This paper addresses modeling and control of a six-degree-of-freedom unmanned aerial vehicle capable of vertical take-off and landing in the presence of wind disturbances. We design a hybrid vehicle that combines the benefits of both the fixed-wing and the rotary-wing UAVs. A non-linear model for the hybrid vehicle is rapidly built, combining rigid body dynamics, aerodynamics of wing, and dynamics of the motor and propeller. Further, we design an H2 optimal controller to make the UAV robust to wind disturbances. We compare its results against that of PID and LQR-based control. Our proposed controller results in better performance in terms of root mean squared errors and time responses during two scenarios: hover and level-flight.


title: Adaptive feature recombination and recalibration for semantic segmentation with Fully Convolutional Networks

abstract: Fully Convolutional Networks have been achieving remarkable results in image semantic segmentation, while being efficient. Such efficiency results from the capability of segmenting several voxels in a single forward pass. So, there is a direct spatial correspondence between a unit in a feature map and the voxel in the same location. In a convolutional layer, the kernel spans over all channels and extracts information from them. We observe that linear recombination of feature maps by increasing the number of channels followed by compression may enhance their discriminative power. Moreover, not all feature maps have the same relevance for the classes being predicted. In order to learn the inter-channel relationships and recalibrate the channels to suppress the less relevant ones, Squeeze and Excitation blocks were proposed in the context of image classification with Convolutional Neural Networks. However, this is not well adapted for segmentation with Fully Convolutional Networks since they segment several objects simultaneously, hence a feature map may contain relevant information only in some locations. In this paper, we propose recombination of features and a spatially adaptive recalibration block that is adapted for semantic segmentation with Fully Convolutional Networks - the SegSE block. Feature maps are recalibrated by considering the cross-channel information together with spatial relevance. Experimental results indicate that Recombination and Recalibration improve the results of a competitive baseline, and generalize across three different problems: brain tumor segmentation, stroke penumbra estimation, and ischemic stroke lesion outcome prediction. The obtained results are competitive or outperform the state of the art in the three applications.


title: Does Explainable Artificial Intelligence Improve Human Decision-Making?

abstract: Explainable AI provides insight into the "why" for model predictions, offering potential for users to better understand and trust a model, and to recognize and correct AI predictions that are incorrect. Prior research on human and explainable AI interactions has focused on measures such as interpretability, trust, and usability of the explanation. Whether explainable AI can improve actual human decision-making and the ability to identify the problems with the underlying model are open questions. Using real datasets, we compare and evaluate objective human decision accuracy without AI (control), with an AI prediction (no explanation), and AI prediction with explanation. We find providing any kind of AI prediction tends to improve user decision accuracy, but no conclusive evidence that explainable AI has a meaningful impact. Moreover, we observed the strongest predictor for human decision accuracy was AI accuracy and that users were somewhat able to detect when the AI was correct versus incorrect, but this was not significantly affected by including an explanation. Our results indicate that, at least in some situations, the "why" information provided in explainable AI may not enhance user decision-making, and further research may be needed to understand how to integrate explainable AI into real systems.


title: REBD:A Conceptual Framework for Big Data Requirements Engineering

abstract: Requirements engineering (RE), as a part of the project development life cycle, has increasingly been recognized as the key to ensuring on-time, on-budget, and goal-based delivery of software projects;compromising this vital phase is nothing but project failures. RE of big data projects is even more crucial because of the main characteristics of big data, including high volume, velocity, and variety. As the traditional RE methods and tools are user-centric rather than data-centric, employing these methodologies is insufficient to fulfill the RE processes for big data projects. Because of the importance of RE and limitations of traditional RE methodologies in the context of big data software projects, in this paper, a big data requirements engineering framework, named REBD, has been proposed. This conceptual framework describes the systematic plan to carry out big data projects starting from requirements engineering to the development, assuring successful execution, and increased productivity of the big data projects.


title: Interplay Between Game Theory and Control Theory

abstract: An interrelationship between Game Theory and Control Theory is seeked. In this respect two aspects of this relationship are brought up. To establish the direct relationship Control Based Games and to establish the inverse relationship Game Based Control are discussed. In the attempt to establish the direct relationship Control Based Boolean Networks are discussed with the novel technique of application of Semi Tensor Product of matrices in Differential Calculus. For the inverse relationship H Infinity robust optimal control has been discussed with the help of Dynamic Programming and Pontryagin Minimization Principle.


title: Abstract Diagrammatic Reasoning with Multiplex Graph Networks

abstract: Abstract reasoning, particularly in the visual domain, is a complex human ability, but it remains a challenging problem for artificial neural learning systems. In this work we propose MXGNet, a multilayer graph neural network for multi-panel diagrammatic reasoning tasks. MXGNet combines three powerful concepts, namely, object-level representation, graph neural networks and multiplex graphs, for solving visual reasoning tasks. MXGNet first extracts object-level representations for each element in all panels of the diagrams, and then forms a multi-layer multiplex graph capturing multiple relations between objects across different diagram panels. MXGNet summarises the multiple graphs extracted from the diagrams of the task, and uses this summarisation to pick the most probable answer from the given candidates. We have tested MXGNet on two types of diagrammatic reasoning tasks, namely Diagram Syllogisms and Raven Progressive Matrices (RPM). For an Euler Diagram Syllogism task MXGNet achieves state-of-the-art accuracy of 99.8%. For PGM and RAVEN, two comprehensive datasets for RPM reasoning, MXGNet outperforms the state-of-the-art models by a considerable margin.


title: Deterministic chaos and fractal entropy scaling in Floquet CFT

abstract: In this paper, we study 2d Floquet conformal field theory, where the external periodic driving is described by iterated logistic or tent maps. These maps are known to be typical examples of dynamical systems exhibiting the order-chaos transition, and we show that, as a result of such driving, the entanglement entropy scaling develops fractal features when the corresponding dynamical system approaches the chaotic regime. For the driving set by the logistic map, fractal contribution to the scaling dominates, making entanglement entropy highly oscillating function of the subsystem size.

categories: cond-mat.stat-mech cond-mat.mes-hall cond-mat.str-el hep-th math-ph math.MP


title: Modeling Individual and Team Behavior through Spatio-temporal Analysis

abstract: Modeling players' behaviors in games has gained increased momentum in the past few years. This area of research has wide applications, including modeling learners and understanding player strategies, to mention a few. In this paper, we present a new methodology, called Interactive Behavior Analytics (IBA), comprised of two visualization systems, a labeling mechanism, and abstraction algorithms that use Dynamic Time Warping and clustering algorithms. The methodology is packaged in a seamless interface to facilitate knowledge discovery from game data. We demonstrate the use of this methodology with data from two multiplayer team-based games: BoomTown, a game developed by Gallup, and DotA 2. The results of this work show the effectiveness of this method in modeling, and developing human-interpretable models of team and individual behavior.


title: Transitions and Multi-Scaling in Rayliegh-Benard Convection. Small-Scale Universality

abstract: Asymptotically large Reynolds number hydrodynamic turbulence is characterized by multi-scaling of moments of velocity increments and spatial derivatives. With decreasing Reynolds number toward $R_{\lambda}=R^{tr}_{\lambda}\approx 9.0$, the anomalous scaling disappears in favor of the "normal" one and close-to-Gaussian probability densities [Yakhot \& Donzis, {\bf 119}, 044501 (2017)]. The nature of this transition and its universality are subjects of this work. Here we consider Benard convection ( Prandtl number $Pr=1$) between infinite horizontal plates. It is shown that in this system the "competition" between Bolgiano and Kolmogorov processes, results in small-scale velocity fluctuations driven by effective "large-scale" Gaussian random temperature field. Therefore, the intermittent dynamics of velocity derivatives are similar or even identical to that in homogeneous and isotropic turbulence generated by the large-scale random forcing. It is shown that low-Rayleigh number instabilities make the problem much more involved and may lead to transition from Gaussian to exponential PDF of the temperature field. The developed {\it mean-field theory} yielded dimensionless heat flux $Nu\propto Ra^{\beta}$ with $\beta\approx 15/56\approx 0.27$, close to the outcome of Chicago experiment. These results point to an unusual small-scale universality of turbulent flows. It is also shown that at $R_{\lambda}\leq 9.0$, a flow "remembers" its laminar background and, therefore, cannot be universal.


title: Sparse Quantile Regression

abstract: We consider both $\ell _{0}$-penalized and $\ell _{0}$-constrained quantile regression estimators. For the $\ell _{0}$-penalized estimator, we derive an exponential inequality on the tail probability of excess quantile prediction risk and apply it to obtain non-asymptotic upper bounds on the mean-square parameter and regression function estimation errors. We also derive analogous results for the $\ell _{0}$-constrained estimator. The resulting rates of convergence are nearly minimax-optimal and the same as those for $\ell _{1}$-penalized and non-convex penalized estimators. Further, we characterize expected Hamming loss for the $\ell _{0}$-penalized estimator. We implement the proposed procedure via mixed integer linear programming and also a more scalable first-order approximation algorithm. We illustrate the finite-sample performance of our approach in Monte Carlo experiments and its usefulness in a real data application concerning conformal prediction of infant birth weights (with $n\approx 10^{3}$ and up to $p>10^{3}$). In sum, our $\ell _{0}$-based method produces a much sparser estimator than the $\ell _{1}$-penalized and non-convex penalized approaches without compromising precision.


title: Fabrication of nanocrystal superlattice microchannels by soft-lithography for electronic measurements of single crystalline domains

abstract: We report a high-throughput and easy-to-implement approach to fabricate microchannels of nanocrystal superlattices with dimensions of ~4 micrometer^2, thus approaching the size of typical single-crystalline domains. By means of microcontact printing, highly ordered superlattices with microscale dimensions are transferred onto photolithographically prepatterned microelectrodes, obtaining well-defined superlattice microchannels. We present step-by-step guidelines for microfabrication, nanocrystal self-assembly and patterning to archive large quantities of up to 330 microchannels per device for statistically meaningful investigations of charge transport in single-crystalline superlattice domains. As proof-of-concept, we perform conductivity and field-effect transistor measurements on microchannels of PbS nanocrystal superlattices. We find that the electric transport within microchannel superlattices is orders of magnitude more efficient than within conventional large-scale channels, highlighting the advantage of the near single-crystalline microchannels presented in this paper.


title: Modeling unconventional superconductivity at the crossover between strong and weak electronic interactions

abstract: High-temperature superconductivity emerges in a host of different quantum materials, often in a region of the phase diagram where the electronic kinetic energy is comparable in magnitude with the electron-electron Coulomb repulsion. Describing such an intermediate-coupling regime has proven challenging, as standard perturbative approaches are inapplicable. Hence, it is of enormous interest to find models that are amenable to be solved using exact methods. While important advances have been made in elucidating the properties of one such minimal model -- the Hubbard model -- via numerical simulations, the infamous fermionic sign-problem significantly limits the accessible parameter space. Here, we employ Quantum Monte Carlo (QMC) methods to solve a multi-band version of the Hubbard model that does not suffer from the sign-problem and in which only repulsive interband interactions are present. In contrast to previous sign-problem-free QMC studies, this model does not have pre-existing fine-tuned magnetic order, and thus treats superconducting, magnetic, and charge degrees of freedom on an equal footing. We find that, as the electron-electron repulsion increases, a dome of antiferromagnetic order emerges in the intermediate-coupling regime, accompanied by a metal-to-insulator crossover line. Superconductivity is found only near the antiferromagnetic quantum phase transition located on the metallic side of the magnetic dome. Across the antiferromagnetic quantum phase transition we find a change in the dynamical character of the magnetic fluctuations, from slow and overdamped in the metallic side to fast and propagating in the insulating side. Our findings shed new light on the intertwining between superconductivity, magnetism, and charge correlations in quantum materials.


title: Differentially Private Variational Autoencoders with Term-wise Gradient Aggregation

abstract: This paper studies how to learn variational autoencoders with a variety of divergences under differential privacy constraints. We often build a VAE with an appropriate prior distribution to describe the desired properties of the learned representations and introduce a divergence as a regularization term to close the representations to the prior. Using differentially private SGD (DP-SGD), which randomizes a stochastic gradient by injecting a dedicated noise designed according to the gradient's sensitivity, we can easily build a differentially private model. However, we reveal that attaching several divergences increase the sensitivity from O(1) to O(B) in terms of batch size B. That results in injecting a vast amount of noise that makes it hard to learn. To solve the above issue, we propose term-wise DP-SGD that crafts randomized gradients in two different ways tailored to the compositions of the loss terms. The term-wise DP-SGD keeps the sensitivity at O(1) even when attaching the divergence. We can therefore reduce the amount of noise. In our experiments, we demonstrate that our method works well with two pairs of the prior distribution and the divergence.


title: Path Planning in a Riemannian Manifold using Optimal Control

abstract: We consider the motion planning of an object in a Riemannian manifold where the object is steered from an initial point to a final point utilizing optimal control. Considering Pontryagin Minimization Principle we compute the Optimal Controls needed for steering the object from initial to final point. The Optimal Controls were solved with respect to time t and shown to have norm 1 which should be the case when the extremal trajectories, which are the solutions of Pontryagin Principle, are arc length parametrized. The extremal trajectories are supposed to be the geodesics on the Riemannian manifold. So we compute the geodesic curvature and the Gaussian curvature of the Riemannian structure.


title: The uncertainty principle: variations on a theme

abstract: We show how a number of well-known uncertainty principles for the Fourier transform, such as the Heisenberg uncertainty principle, the Donoho--Stark uncertainty principle, and Meshulam's non-abelian uncertainty principle, have little to do with the structure of the Fourier transform itself. Rather, all of these results follow from very weak properties of the Fourier transform (shared by numerous linear operators), namely that it is bounded as an operator $L^1 \to L^\infty$, and that it is unitary. Using a single, simple proof template, and only these (or weaker) properties, we obtain some new proofs and many generalizations of these basic uncertainty principles, to new operators and to new settings, in a completely unified way. Together with our general overview, this paper can also serve as a survey of the many facets of the phenomena known as uncertainty principles.

categories: math.FA cs.IT math-ph math.CO math.GR math.IT math.MP


title: Frustratingly Simple Domain Generalization via Image Stylization

abstract: Convolutional Neural Networks (CNNs) show impressive performance in the standard classification setting where training and testing data are drawn i.i.d. from a given domain. However, CNNs do not readily generalize to new domains with different statistics, a setting that is simple for humans. In this work, we address the Domain Generalization problem, where the classifier must generalize to an unknown target domain. Inspired by recent works that have shown a difference in biases between CNNs and humans, we demonstrate an extremely simple yet effective method, namely correcting this bias by augmenting the dataset with stylized images. In contrast with existing stylization works, which use external data sources such as art, we further introduce a method that is entirely in-domain using no such extra sources of data. We provide a detailed analysis as to the mechanism by which the method works, verifying our claim that it changes the shape/texture bias, and demonstrate results surpassing or comparable to the state of the arts that utilize much more complex methods.


title: Alpha clustering and alpha-capture reaction rate from ab initio symmetry-adapted description of $^{20}$Ne

abstract: We introduce a new framework for studying clustering and for calculating alpha partial widths using ab initio wave functions. We demonstrate the formalism for $^{20}$Ne, by calculating the overlap between the $^{16}$O$+\alpha$ cluster configuration and states in $^{20}$Ne computed in the ab initio symmetry-adapted no-core shell model. We present spectroscopic amplitudes and spectroscopic factors, and compare those to no-core symplectic shell-model results in larger model spaces, to gain insight into the underlying physics that drives alpha-clustering. Specifically, we report on the alpha partial width of the lowest $1^-$ resonance in $^{20}$Ne, which is found to be in good agreement with experiment. We also present first no-core shell-model estimates for asymptotic normalization coefficients for the ground state, as well as for the first excited $4^{+}$ state in $^{20}$Ne that lies in a close proximity to the $^{16}$O$+\alpha$ threshold. This outcome highlights the importance of correlations for developing cluster structures and for describing alpha widths. The widths can then be used to calculate alpha-capture reaction rates for narrow resonances of interest to astrophysics. We explore the reaction rate for the alpha-capture reaction $^{16}$O$(\alpha,\gamma)^{20}$Ne at astrophysically relevant temperatures and determine its impact on simulated X-ray burst abundances.


title: An analytic theory of shallow networks dynamics for hinge loss classification

abstract: Neural networks have been shown to perform incredibly well in classification tasks over structured high-dimensional datasets. However, the learning dynamics of such networks is still poorly understood. In this paper we study in detail the training dynamics of a simple type of neural network: a single hidden layer trained to perform a classification task. We show that in a suitable mean-field limit this case maps to a single-node learning problem with a time-dependent dataset determined self-consistently from the average nodes population. We specialize our theory to the prototypical case of a linearly separable dataset and a linear hinge loss, for which the dynamics can be explicitly solved. This allow us to address in a simple setting several phenomena appearing in modern networks such as slowing down of training dynamics, crossover between rich and lazy learning, and overfitting. Finally, we asses the limitations of mean-field theory by studying the case of large but finite number of nodes and of training samples.

categories: stat.ML cond-mat.dis-nn cond-mat.stat-mech cs.LG


title: Improved Short-Baseline Neutrino Oscillation Search and Energy Spectrum Measurement with the PROSPECT Experiment at HFIR

abstract: We present a detailed report on sterile neutrino oscillation and U-235 antineutrino energy spectrum measurement results from the PROSPECT experiment at the highly enriched High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory. In 96 calendar days of data taken at an average baseline distance of 7.9 m from the center of the 85 MW HFIR core, the PROSPECT detector has observed more than 50,000 interactions of antineutrinos produced in beta decays of U-235 fission products. New limits on the oscillation of antineutrinos to light sterile neutrinos have been set by comparing the detected energy spectra of ten reactor-detector baselines between 6.7 and 9.2 meters. Measured differences in energy spectra between baselines show no statistically significant indication of antineutrinos to sterile neutrino oscillation and disfavor the Reactor Antineutrino Anomaly best-fit point at the 2.5$\sigma$ confidence level. The reported U-235 antineutrino energy spectrum measurement shows excellent agreement with energy spectrum models generated via conversion of the measured U-235 beta spectrum, with a $\chi^2$/DOF of 31/31. PROSPECT is able to disfavor at 2.4$\sigma$ confidence level the hypothesis that U-235 antineutrinos are solely responsible for spectrum discrepancies between model and data obtained at commercial reactor cores. A data-model deviation in PROSPECT similar to that observed by commercial core experiments is preferred with respect to no observed deviation, at a 2.2$\sigma$ confidence level.


title: Rumor source detection with multiple observations under adaptive diffusions

abstract: Recent work, motivated by anonymous messaging platforms, has introduced adaptive diffusion protocols which can obfuscate the source of a rumor: a "snapshot adversary" with access to the subgraph of "infected" nodes can do no better than randomly guessing the entity of the source node. What happens if the adversary has access to multiple independent snapshots? We study this question when the underlying graph is the infinite $d$-regular tree. We show that (1) a weak form of source obfuscation is still possible in the case of two independent snapshots, but (2) already with three observations there is a simple algorithm that finds the rumor source with constant probability, regardless of the adaptive diffusion protocol. We also characterize the tradeoff between local spreading and source obfuscation for adaptive diffusion protocols (under a single snapshot). These results raise questions about the robustness of anonymity guarantees when spreading information in social networks.


title: Some new results on relative entropy production, time reversal, and optimal control of time-inhomogeneous diffusion processes

abstract: This paper studies time-inhomogeneous nonequilibrium diffusion processes, including both Brownian dynamics and Langevin dynamics. We derive upper bounds of the relative entropy production of the time-inhomogeneous process with respect to the transient invariant probability measures. We also study the time reversal of the reverse process in Crooks' fluctuation theorem. We show that the time reversal of the reverse process coincides with the optimally controlled forward process that leads to zero variance importance sampling estimator based on Jarzynski's equality.


title: Shallow-circuit variational quantum eigensolver based on symmetry-inspired Hilbert space partitioning for quantum chemical calculations

abstract: Development of resource-friendly quantum algorithms remains highly desirable for noisy intermediate-scale quantum computing. Based on the variational quantum eigensolver (VQE) with unitary coupled cluster ansatz, we demonstrate that partitioning of the Hilbert space made possible by the point group symmetry of the molecular systems greatly reduces the number of variational operators by confining the variational search within a subspace. In addition, we found that instead of including all subterms for each excitation operator, a single-term representation suffices to reach required accuracy for various molecules tested, resulting in an additional shortening of the quantum circuit. With these strategies, VQE calculations on a noiseless quantum simulator achieve energies within a few meVs of those obtained with the full UCCSD ansatz for $\mathrm{H}_4$ square, $\mathrm{H}_4$ chain and $\mathrm{H}_6$ hexagon molecules; while the number of controlled-NOT (CNOT) gates, a measure of the quantum-circuit depth, is reduced by a factor of as large as 35. Furthermore, we introduced an efficient "score" parameter to rank the excitation operators, so that the operators causing larger energy reduction can be applied first. Using $\mathrm{H}_4$ square and $\mathrm{H}_4$ chain as examples, We demonstrated on noisy quantum simulators that the first few variational operators can bring the energy within the chemical accuracy, while additional operators do not improve the energy since the accumulative noise outweighs the gain from the expansion of the variational ansatz.


title: Spin-glass thermo-remanent magnetization revisited: a numerical and analytical study

abstract: Thermoremanent magnetization data for the 3D Edwards-Anderson spin glass are generated using the Waiting Time Method as simulational tool and interpreted using Record Dynamics. We verify that clusters of contiguous spins are overturned by quakes, non-equilibrium events linked to record sized energy fluctuations and show that quaking is a log-Poisson process, i.e. a Poisson process whose average depends on the logarithm of the system age, counted from the initial quench. Our results compare favorably with experimental thermoremanent magnetization findings and with the spontaneous fluctuation dynamics of the E-A model. The logarithm growth of the size of overturned clusters is related to similar experimental results and to the growing length scale of the spin-spin spatial correlation function. The analysis buttresses the applicability of the Waiting Time Method as a simulational tool and of Record Dynamics as coarse-graining method for aging dynamics.


title: The physical properties and impact of AGN outflows from high to low redshift

abstract: Feedback from active galactic nuclei (AGN) on their host galaxies, in the form of gas outflows capable of quenching star formation, is considered a major player in galaxy evolution. However, clear observational evidence of such major impact is still missing; uncertainties in measuring outflow properties might be partly responsible because of their critical role in comparisons with models and in constraining the impact of outflows on galaxies. Here we briefly review the challenges in measuring outflow physical properties and present an overview of outflow studies from high to low redshift. Finally, we present highlights from our MAGNUM survey of nearby AGN with VLT/MUSE, where the high intrinsic spatial resolution (down to $\sim$10 pc) allows us to accurately measure the physical and kinematic properties of ionised gas outflows.


title: Angular and polarization observables for Majorana-mediated B decays with effective interactions

abstract: We probe the effective field theory extending the Standard Model with a sterile neutrino in B meson decays at B factories and lepton colliders, using angular and polarization observables. We put bounds on different effective operators characterized by their distinct Dirac-Lorentz structure, and probe the $N$-mediated B decays sensitivity to these interactions. We define a Forward-Backward asymmetry $A_{FB}^{\ell \gamma}$ between the muon and photon directions for the $B \to \mu \nu \gamma$ decay, which allows us to separate the SM contribution from the effective lepton number conserving and violating processes, mediated by a near on-shell $N$. Using the most stringent constraints on the effective parameter space from Belle and BaBar we find that a measurement of the final polarization $P_{\tau}$ in the rare $B^-\rightarrow \ell^-_{1}\ell^-_{2} \pi^+$ decays can help us infer the scalar or vector interaction content in the $N$ production or decay vertices. We find that the B meson decays are more sensitive to scalar operators.


title: Causal representation of multi-loop Feynman integrands within the loop-tree duality

abstract: The numerical evaluation of multi-loop scattering amplitudes in the Feynman representation usually requires to deal with both physical (causal) and unphysical (non-causal) singularities. The loop-tree duality (LTD) offers a powerful framework to easily characterise and distinguish these two types of singularities, and then simplify analytically the underling expressions. In this paper, we work explicitly on the dual representation of multi-loop Feynman integrals generated from three parent topologies, which we refer to as Maximal, Next-to-Maximal and Next-to-Next-to-Maximal loop topologies. In particular, we aim at expressing these dual contributions, independently of the number of loops and internal configurations, in terms of causal propagators only. Thus, providing very compact and causal integrand representations to all orders. In order to do so, we reconstruct their analytic expressions from numerical evaluation over finite fields. This procedure implicitly cancels out all unphysical singularities. We also interpret the result in terms of entangled causal thresholds. In view of the simple structure of the dual expressions, we integrate them numerically up to four loops in integer space-time dimensions, taking advantage of their smooth behaviour at integrand level.


title: A Computational Multi-Criteria Optimization Approach to Controller Design for Physical Human-Robot Interaction

abstract: Physical human-robot interaction (pHRI) integrates the benefits of human operator and a collaborative robot in tasks involving physical interaction, with the aim of increasing the task performance. However, the design of interaction controllers that achieve safe and transparent operations is challenging, mainly due to the contradicting nature of these objectives. Knowing that attaining perfect transparency is practically unachievable, controllers that allow better compromise between these objectives are desirable. In this paper, we propose a multi-criteria optimization framework, which jointly optimizes the stability robustness and transparency of a closed-loop pHRI system for a given interaction controller. In particular, we propose a Pareto optimization framework that allows the designer to make informed decisions by thoroughly studying the trade-off between stability robustness and transparency. The proposed framework involves a search over the discretized controller parameter space to compute the Pareto front curve and a selection of controller parameters that yield maximum attainable transparency and stability robustness by studying this trade-off curve. The proposed framework not only leads to the design of an optimal controller, but also enables a fair comparison among different interaction controllers. In order to demonstrate the practical use of the proposed approach, integer and fractional order admittance controllers are studied as a case study and compared both analytically and experimentally. The experimental results validate the proposed design framework and show that the achievable transparency under fractional order admittance controller is higher than that of integer order one, when both controllers are designed to ensure the same level of stability robustness.


title: An Integral Basis for the Universal Enveloping Algebra of the Onsager Algebra

abstract: We construct an integral form for the universal enveloping algebra of the Onsager algebra and an explicit integral basis for this integral form. We also formulate straightening identities among some products of basis elements.


title: Localized Spectral Graph Filter Frames: A Unifying Framework, Survey of Design Considerations, and Numerical Comparison (Extended Cut)

abstract: Representing data residing on a graph as a linear combination of building block signals can enable efficient and insightful visual or statistical analysis of the data, and such representations prove useful as regularizers in signal processing and machine learning tasks. Designing collections of building block signals -- or more formally, dictionaries of atoms -- that specifically account for the underlying graph structure as well as any available representative training signals has been an active area of research over the last decade. In this article, we survey a particular class of dictionaries called localized spectral graph filter frames, whose atoms are created by localizing spectral patterns to different regions of the graph. After showing how this class encompasses a variety of approaches from spectral graph wavelets to graph filter banks, we focus on the two main questions of how to design the spectral filters and how to select the center vertices to which the patterns are localized. Throughout, we emphasize computationally efficient methods that ensure the resulting transforms and their inverses can be applied to data residing on large, sparse graphs. We demonstrate how this class of transform methods can be used in signal processing tasks such as denoising and non-linear approximation, and provide code for readers to experiment with these methods in new application domains.


title: Leveraging Machine Learning to Alleviate Hubbard Model Sign Problems

abstract: Lattice Monte Carlo calculations of interacting systems on non-bipartite lattices exhibit an oscillatory imaginary phase known as the phase or sign problem, even at zero chemical potential. One method to alleviate the sign problem is to analytically continue the integration region of the state variables into the complex plane via holomorphic flow equations. For asymptotically large flow times the state variables approach manifolds of constant imaginary phase known as Lefschetz thimbles. However, flowing such variables and calculating the ensuing Jacobian is a computationally demanding procedure. In this paper we demonstrate that neural networks can be trained to parameterize suitable manifolds for this class of sign problem and drastically reduce the computational cost. We apply our method to the Hubbard model on the triangle and tetrahedron, both of which are non-bipartite. At strong interaction strengths and modest temperatures the tetrahedron suffers from a severe sign problem that cannot be overcome with standard reweighting techniques, while it quickly yields to our method. We benchmark our results with exact calculations and comment on future directions of this work.


title: Valuing the quality option in agricultural commodity futures: a Monte Carlo simulation based approach

abstract: Agricultural commodity futures are often settled by delivery. Quality options that allow the futures short to deliver one of several underlying assets are commonly used in such contracts to prevent manipulation. Inclusion of these options reduces the price of the futures contract and leads to degraded contract hedging performance. Valuation of these options is a first step in assessing the impact of the quality options embedded into a futures contract. This paper demonstrates a Monte Carlo simulation based approach to estimate the value of a quality option. In order to improve simulation efficiency, the technique of antithetic variables is used. This approach can help in the assessment of the impact of embedded quality options.


title: Unified Representation Learning for Efficient Medical Image Analysis

abstract: Medical image analysis typically includes several tasks such as enhancement, segmentation, and classification. Traditionally, these tasks are implemented using separate deep learning models for separate tasks, which is not efficient because it involves unnecessary training repetitions, demands greater computational resources, and requires a relatively large amount of labeled data. In this paper, we propose a multi-task training approach for medical image analysis, where individual tasks are fine-tuned simultaneously through relevant knowledge transfer using a unified modality-specific feature representation (UMS-Rep). We explore different fine-tuning strategies to demonstrate the impact of the strategy on the performance of target medical image tasks. We experiment with different visual tasks (e.g., image denoising, segmentation, and classification) to highlight the advantages offered with our approach for two imaging modalities, chest X-ray and Doppler echocardiography. Our results demonstrate that the proposed approach reduces the overall demand for computational resources and improves target task generalization and performance. Further, our results prove that the performance of target tasks in medical images is highly influenced by the utilized fine-tuning strategy.


title: Shadowing, Finite Order Shifts and Ultrametric Spaces

abstract: Inspired by a recent novel work of Good and Meddaugh, we establish fundamental connections between shadowing, finite order shifts, and ultrametric complete spaces. We develop a theory of shifts of finite type for infinite alphabets. We call them shifts of finite order. We develop the basic theory of the shadowing property in general metric spaces, exhibiting similarities and differences with the theory in compact spaces. We connect these two theories in the setting of zero-dimensional complete spaces, showing that a uniformly continuous map of an ultrametric complete space has the finite shadowing property if, and only if, it is an inverse limit of a system of shifts of finite order satisfying the Mittag-Leffler Condition. Furthermore, in this context, we show that the shadowing property is equivalent to the finite shadowing property and the fulfillment of the Mittag-Leffler Condition in the inverse limit description of the system. As corollaries, we obtain that a variety of maps in ultrametric spaces have the shadowing property, such as similarities and, more generally, maps which themselves, or their inverses, have Lipschitz constant 1. Finally, we apply our results to the dynamics of $p$-adic integers and $p$-adic rationals.


title: Solar neutrino probes of the muon anomalous magnetic moment in the gauged $U(1)_{L_\mu-L_\tau}$

abstract: Models of gauged $U(1)_{L_\mu-L_\tau}$ can provide a solution to the long-standing discrepancy between the theoretical prediction for the muon anomalous magnetic moment and its measured value. The extra contribution is due to a new light vector mediator, which also helps to alleviate an existing tension in the determination of the Hubble parameter. In this article, we explore ways to probe this solution via the scattering of solar neutrinos with electrons and nuclei in a range of experiments and considering high and low solar metallicity scenarios. In particular, we reevaluate Borexino constraints on neutrino-electron scattering, finding them to be more stringent than previously reported, and already excluding a part of the $(g-2)_\mu$ explanation with mediator masses smaller than $2\times10^{-2}$ GeV. We then show that future direct dark matter detectors will be able to probe most of the remaining solution. Due to its large exposure, LUX-ZEPLIN will explore regions with mediator masses up to $5\times10^{-2}$ GeV and DARWIN will be able to extend the search beyond $10^{-1}$ GeV, thereby covering most of the area compatible with $(g-2)_\mu$. For completeness, we have also computed the constraints derived from the recent XENON1T electron recoil search and from the CENNS-10 LAr detector, showing that none of them excludes new areas of the parameter space. Should the excess in the muon anomalous magnetic moment be confirmed, our work suggests that direct detection experiments could provide crucial information with which to test the $U(1)_{L_\mu-L_\tau}$ solution, complementary to efforts in neutrino experiments and accelerators.


title: Gradient descent follows the regularization path for general losses

abstract: Recent work across many machine learning disciplines has highlighted that standard descent methods, even without explicit regularization, do not merely minimize the training error, but also exhibit an implicit bias. This bias is typically towards a certain regularized solution, and relies upon the details of the learning process, for instance the use of the cross-entropy loss. In this work, we show that for empirical risk minimization over linear predictors with arbitrary convex, strictly decreasing losses, if the risk does not attain its infimum, then the gradient-descent path and the algorithm-independent regularization path converge to the same direction (whenever either converges to a direction). Using this result, we provide a justification for the widely-used exponentially-tailed losses (such as the exponential loss or the logistic loss): while this convergence to a direction for exponentially-tailed losses is necessarily to the maximum-margin direction, other losses such as polynomially-tailed losses may induce convergence to a direction with a poor margin.


title: Lookahead Adversarial Learning for Near Real-Time Semantic Segmentation

abstract: Semantic segmentation is one of the most fundamental problems in computer vision with significant impact on a wide variety of applications. Adversarial learning is shown to be an effective approach for improving semantic segmentation quality by enforcing higher-level pixel correlations and structural information. However, state-of-the-art semantic segmentation models cannot be easily plugged into an adversarial setting because they are not designed to accommodate convergence and stability issues in adversarial networks. We bridge this gap by building a conditional adversarial network with a state-of-the-art segmentation model (DeepLabv3+) at its core. To battle the stability issues, we introduce a novel lookahead adversarial learning (LoAd) approach with an embedded label map aggregation module. We focus on semantic segmentation models that run fast at inference for near real-time field applications. Through extensive experimentation, we demonstrate that the proposed solution can alleviate divergence issues in an adversarial semantic segmentation setting and results in considerable performance improvements (+5% in some classes) on the baseline for three standard datasets.


title: Distortion estimates for approximate Bayesian inference

abstract: Current literature on posterior approximation for Bayesian inference offers many alternative methods. Does our chosen approximation scheme work well on the observed data? The best existing generic diagnostic tools treating this kind of question by looking at performance averaged over data space, or otherwise lack diagnostic detail. However, if the approximation is bad for most data, but good at the observed data, then we may discard a useful approximation. We give graphical diagnostics for posterior approximation at the observed data. We estimate a "distortion map" that acts on univariate marginals of the approximate posterior to move them closer to the exact posterior, without recourse to the exact posterior.


title: Aligning Nuclear Cluster Orbits with an Active Galactic Nucleus Accretion Disk

abstract: Active galactic nuclei (AGN) are powered by the accretion of disks of gas onto supermassive black holes (SMBHs). Stars and stellar remnants orbiting the SMBH in the nuclear star cluster (NSC) will interact with the AGN disk. Orbiters plunging through the disk experience a drag force and, through repeated passage, can have their orbits captured by the disk. A population of embedded objects in AGN disks may be a significant source of binary black hole mergers, supernovae, tidal disruption events and embedded gamma-ray bursts. For two representative AGN disk models we use geometric drag and Bondi-Hoyle-Littleton drag to determine the time to capture for stars and stellar remnants. We assume a range of initial inclination angles and semi-major axes for circular Keplerian prograde orbiters. Capture time strongly depends on the density and aspect ratio of the chosen disk model, the relative velocity of the stellar object with respect to the disk, and the AGN lifetime. We expect that for an AGN disk density $\rho \gtrsim 10^{-11}\rm g/cm^3$ and disk lifetime $\geq 1$Myr, there is a significant population of embedded stellar objects, which can fuel mergers detectable in gravitational waves with LIGO-Virgo and LISA.


title: On Power integral bases for certain pure number fields

abstract: Let $K=\mathbb{Q}(\alpha)$ be a number field generated by a complex root $\alpha$ of a monic irreducible polynomial $f(x)=x^{12}-m$, with $m\neq 1$ is a square free rational integer. In this paper, we prove that if $m \equiv 2$ or $3$ (mod 4) and $m\not\equiv \mp 1$ (mod 9), then the number field $K$ is monogenic. If $m \equiv 1$ (mod 8) or $m\equiv \mp 1$ (mod 9), then the number field $K$ is not monogenic.


title: Predicting the phase stability of multi-component high entropy compounds

abstract: A generic method to estimate the relative feasibility of formation of high entropy compounds in a single phase, directly from first principles, is developed. As a first step, the relative formation abilities of 56 multi-component, AO, oxides were evaluated. These were constructed from 5 cation combinations chosen from A={Ca, Co, Cu, Fe, Mg, Mn, Ni, Zn}. Candidates for multi-component oxides are predicted from descriptors related to the enthalpy and configurational entropy obtained from the mixing enthalpies of two component oxides. The utility of this approach is evaluated by comparing the predicted combinations with the experimentally realized entropy stabilized oxide, (MgCoCuNiZn)O. In the second step, Monte Carlo simulations are utilized to investigate the phase composition and local ionic segregation as a function of temperature. This approach allows for the evaluation of potential secondary phases, thereby making realistic predictions of novel multi-component compounds that can be synthesized.


title: Cartesian products of the $g$-topologies are a $g$-topology

abstract: We show that unlike the usual topologies the $g$-topologies are closed with respect to the Cartesian products. Moreover, we bring much detailed explanations some examples of concepts related the statistical metric spaces.


title: Non-repudiable provenance for clinical decision support systems

abstract: Provenance templates are now a recognised methodology for the construction of data provenance records. Each template defines the provenance of a domain-specific action in abstract form, which may then be instantiated as required by a single call to the provenance template service. As data reliability and trustworthiness becomes a critical issue in an increasing number of domains, there is a corresponding need to ensure that the provenance of that data is non-repudiable. In this paper we contribute two new, complementary modules to our template model and implementation to produce non-repudiable data provenance. The first, a module that traces the operation of the provenance template service itself, and records a provenance trace of the construction of an object-level document, at the level of individual service calls. The second, a non-repudiation module that generates evidence for the data recorded about each call, annotates the service trace accordingly, and submits a representation of that evidence to a provider-agnostic notary service. We evaluate the applicability of our approach in the context of a clinical decision support system. We first define a policy to ensure the non-repudiation of evidence with respect to a security threat analysis in order to demonstrate the suitability of our solution. We then select three use cases from within a particular system, Consult, with contrasting data provenance recording requirements and analyse the subsequent performance of our prototype implementation against three different notary providers.


title: Semi-Discriminative Representation Loss for Online Continual Learning

abstract: The use of episodic memory in continual learning has demonstrated effectiveness for alleviating catastrophic forgetting. In recent studies, gradient-based approaches have been developed to make more efficient use of compact episodic memory. Such approaches refine the gradients resulting from new samples by those from memorized samples, aiming to reduce the diversity of gradients from different tasks. In this paper, we clarify the relation between diversity of gradients and discriminativeness of representations, showing shared as well as conflicting interests between Deep Metric Learning and continual learning, thus demonstrating pros and cons of learning discriminative representations in continual learning. Based on these findings, we propose a simple method -- Semi-Discriminative Representation Loss (SDRL) -- for continual learning. In comparison with state-of-the-art methods, SDRL shows better performance with low computational cost on multiple benchmark tasks in the setting of online continual learning.


title: Constraining Early Dark Energy with Large-Scale Structure

abstract: An axion-like field comprising $\sim 10\%$ of the energy density of the universe near matter-radiation equality is a candidate to resolve the Hubble tension; this is the "early dark energy" (EDE) model. However, as shown in Hill et al. (2020), the model fails to simultaneously resolve the Hubble tension and maintain a good fit to both cosmic microwave background (CMB) and large-scale structure (LSS) data. Here, we use redshift-space galaxy clustering data to sharpen constraints on the EDE model. We perform the first EDE analysis using the full-shape power spectrum likelihood from the Baryon Oscillation Spectroscopic Survey (BOSS), based on the effective field theory (EFT) of LSS. The inclusion of this likelihood in the EDE analysis yields a $25\%$ tighter error bar on $H_0$ compared to primary CMB data alone, yielding $H_0 = 68.54^{+0.52}_{-0.95}$ km/s/Mpc ($68\%$ CL). In addition, we constrain the maximum fractional energy density contribution of the EDE to $f_{\rm EDE} < 0.072$ ($95\%$ CL). We explicitly demonstrate that the EFT BOSS likelihood yields much stronger constraints on EDE than the standard BOSS likelihood. Including further information from photometric LSS surveys,the constraints narrow by an additional $20\%$, yielding $H_0 = 68.73^{+0.42}_{-0.69}$ km/s/Mpc ($68\%$ CL) and $f_{\rm EDE}<0.053$ ($95\%$ CL). These bounds are obtained without including local-universe $H_0$ data, which is in strong tension with the CMB and LSS, even in the EDE model. We also refute claims that MCMC analyses of EDE that omit SH0ES from the combined dataset yield misleading posteriors. Finally, we demonstrate that upcoming Euclid/DESI-like spectroscopic galaxy surveys can greatly improve the EDE constraints. We conclude that current data preclude the EDE model as a resolution of the Hubble tension, and that future LSS surveys can close the remaining parameter space of this model.


title: Using Deep Learning to Predict Beam-Tunable Pareto Optimal Dose Distribution for Intensity Modulated Radiation Therapy

abstract: We propose to develop deep learning models that can predict Pareto optimal dose distributions by using any given set of beam angles, along with patient anatomy, as input to train the deep neural networks. We implement and compare two deep learning networks that predict with two different beam configuration modalities. We generated Pareto optimal plans for 70 patients with prostate cancer. We used fluence map optimization to generate 500 IMRT plans that sampled the Pareto surface for each patient, for a total of 35,000 plans. We studied and compared two different models, Model I and Model II. Model I directly uses beam angles as a second input to the network as a binary vector. Model II converts the beam angles into beam doses that are conformal to the PTV. Our deep learning models predicted voxel-level dose distributions that precisely matched the ground truth dose distributions. Quantitatively, Model I prediction error of 0.043 (confirmation), 0.043 (homogeneity), 0.327 (R50), 2.80% (D95), 3.90% (D98), 0.6% (D50), 1.10% (D2) was lower than that of Model II, which obtained 0.076 (confirmation), 0.058 (homogeneity), 0.626 (R50), 7.10% (D95), 6.50% (D98), 8.40% (D50), 6.30% (D2). Treatment planners who use our models will be able to use deep learning to control the tradeoffs between the PTV and OAR weights, as well as the beam number and configurations in real time. Our dose prediction methods provide a stepping stone to building automatic IMRT treatment planning.


title: 2020 Global reassessment of the neutrino oscillation picture

abstract: We present an updated global fit of neutrino oscillation data in the simplest three-neutrino framework. In the present study we include up-to-date analyses from a number of experiments. Concerning the atmospheric and solar sectors, we give updated analyses of DeepCore and SNO data, respectively. We have also included the latest electron antineutrino data collected by the Daya Bay and RENO reactor experiments, and the long-baseline T2K and NO$\nu$A measurements. These new analyses result in more accurate measurements of $\theta_{13}$, $\theta_{12}$, $\Delta m_{21}^2$ and $|\Delta m_{31}^2|$. The best fit value for the atmospheric angle $\theta_{23}$ lies in the second octant, but first octant solutions remain allowed at $\sim2.4\sigma$. Regarding CP violation measurements, the preferred value of $\delta$ we obtain is 1.08$\pi$ (1.58$\pi$) for normal (inverted) neutrino mass ordering. The global analysis prefers normal neutrino mass ordering with 2.5$\sigma$. This preference is milder than the one found in previous global analyses. The new results should be regarded as robust due to the agreement found between our Bayesian and frequentist approaches. Taking into account only oscillation data, there is a weak/moderate preference for the normal neutrino mass ordering of $2.00\sigma$. While adding neutrinoless double beta decay from the latest Gerda, CUORE and KamLAND-Zen results barely modifies this picture, cosmological measurements raise the preference to $2.68\sigma$ within a conservative approach. A more aggressive data set combination of cosmological observations leads to a similar preference, namely $2.70\sigma$. This very same cosmological data set provides $2\sigma$ upper limits on the total neutrino mass corresponding to $\sum\nu<0.12$ ($0.15$)~eV for normal (inverted) neutrino mass ordering.


title: Geometry of polarised varieties

abstract: In this paper we investigate the geometry of projective varieties polarised by ample and more generally nef and big Weil divisors. First we study birational boundedness of linear systems. We show that if $X$ is a projective variety of dimension $d$ with $\epsilon$-lc singularities for $\epsilon>0$, and if $N$ is a nef and big Weil divisor on $X$ such that $N-K_X$ is pseudo-effective, then the linear system $|mN|$ defines a birational map for some natural number $m$ depending only on $d,\epsilon$. This is key to proving various other results. For example, it implies that if $N$ is a big Weil divisor (not necessarily nef) on a klt Calabi-Yau variety of dimension $d$, then the linear system $|mN|$ defines a birational map for some natural number $m$ depending only on $d$. It also gives new proofs of some known results, for example, if $X$ is an $\epsilon$-lc Fano variety of dimension $d$ then taking $N=-K_X$ we recover birationality of $|-mK_X|$ for bounded $m$. We prove similar birational boundedness results for nef and big Weil divisors $N$ on projective klt varieties $X$ when both $K_X$ and $N-K_X$ are pseudo-effective (here $X$ is not assumed $\epsilon$-lc). Using the above, we show boundedness of polarised varieties under some natural conditions. We extend these to boundedness of semi-log canonical Calabi-Yau pairs polarised by effective ample Weil divisors not containing lc centres. We will briefly discuss applications to existence of projective coarse moduli spaces of such polarised Calabi-Yau pairs.


title: Denoising Diffusion Probabilistic Models

abstract: We present high quality image synthesis results using diffusion probabilistic models, a class of latent variable models inspired by considerations from nonequilibrium thermodynamics. Our best results are obtained by training on a weighted variational bound designed according to a novel connection between diffusion probabilistic models and denoising score matching with Langevin dynamics, and our models naturally admit a progressive lossy decompression scheme that can be interpreted as a generalization of autoregressive decoding. On the unconditional CIFAR10 dataset, we obtain an Inception score of 9.46 and a state-of-the-art FID score of 3.17. On 256x256 LSUN, we obtain sample quality similar to ProgressiveGAN. Our implementation is available at https://github.com/hojonathanho/diffusion


title: Optimal control for the nonlinear Fisher-Kolmogorov system with applications to aquatic plant management

abstract: Spatiotemporal dynamics of populations may be described by the reaction-diffusion Fisher-Kolmogorov model. In this work we have proposed a new formulation for a control problem of aquatic plants in a temporal dynamics. The solution of this problem is extended to a spatiotemporal Fisher-Kolmogorov system with multiple species of plants interacting in the same place. The control consists on human intervention as a strategy for management of the aquatic plants. In our applications, one plant and two plants cases have been considered. Simulation results are presented to show the effectiveness of the proposed control strategies.


title: A simple convergence proof for the lace expansion

abstract: We use the lace expansion to give a simple proof that the critical two-point function for weakly self-avoiding walk on $\mathbb{Z}^d$ has decay $|x|^{-(d-2)}$ in dimensions $d>4$. The proof uses elementary Fourier analysis and the Riemann--Lebesgue Lemma.


abstract: Consistency Guided Scene Flow Estimation (CGSF) is a self-supervised framework for the joint reconstruction of 3D scene structure and motion from stereo video. The model takes two temporal stereo pairs as input, and predicts disparity and scene flow. The model self-adapts at test time by iteratively refining its predictions. The refinement process is guided by a consistency loss, which combines stereo and temporal photo-consistency with a geometric term that couples disparity and 3D motion. To handle inherent modeling error in the consistency loss (e.g. Lambertian assumptions) and for better generalization, we further introduce a learned, output refinement network, which takes the initial predictions, the loss, and the gradient as input, and efficiently predicts a correlated output update. In multiple experiments, including ablation studies, we show that the proposed model can reliably predict disparity and scene flow in challenging imagery, achieves better generalization than the state-of-the-art, and adapts quickly and robustly to unseen domains.


title: Hidden Photon Dark Matter in the Light of XENON1T and Stellar Cooling

abstract: The low-energy electronic recoil spectrum in XENON1T provides an intriguing hint for potential new physics. At the same time, observations of horizontal branch stars favor the existence of a small amount of extra cooling compared to the one expected from the Standard Model particle content. In this note, we argue that a hidden photon with a mass of $\sim 2.5$ keV and a kinetic mixing of $\sim 10^{-15}$ allows for a good fit to both of these excesses. In this scenario, the signal detected in XENON1T is due to the absorption of hidden photon dark matter particles, whereas the anomalous cooling of horizontal branch stars arises from resonant production of hidden photons in the stellar interior.


title: Counting Risk Increments to Make Decisions During an Epidemic

abstract: I propose a smartphone app that will allow people to participate in the management of their own safety during an epidemic or pandemic such as COVID-19 by enabling them to view, in advance, the risks they would take if they visit some given venue (a cafe, the gym, the workplace, the park,...) and, furthermore, track the accumulation of such risks during the course of any given day or week. This idea can be presented to users of the app as counting points. One point represents some constant probability, $p_\text{point}$, of infection. Then the app would work in a similar way to a calorie counting app (instead of counting calories we count probability increments of being infected). Government could set a maximum recommended number of daily (or weekly) points available to each user in accord with its objectives (bringing the disease under control, allowing essential workers to work, protecting vulnerable individuals, ...). It is posited that this, along with other proposed "levers" would allow government to manage a gradual transition to normalcy. I discuss a circuit framework with wires running between boxes. In this framework the wires represent possible sources of infection, namely individuals and the venues themselves (through deposits of pathogens left at the venue). The boxes represent interactions of these sources (when individuals visit a venue). This circuit framework allows (i) calculation of points cost for visiting venues and (ii) probabilistic contact tracing. The points systems proposed here could complement existing contact tracing apps by adding functionality to permit users to participate in decision making up front.


title: List decoding of Convolutional Codes over integer residue rings

abstract: A convolutional code $\C$ over $\ZZ[D]$ is a $\ZZ[D]$-submodule of $\ZZN[D]$ where $\ZZ[D]$ stands for the ring of polynomials with coefficients in $\ZZ$. In this paper, we study the list decoding problem of these codes when the transmission is performed over an erasure channel, that is, we study how much information one can recover from a codeword $w\in \C$ when some of its coefficients have been erased. We do that using the $p$-adic expansion of $w$ and particular representations of the parity-check polynomial matrix of the code. From these matrix polynomial representations we recursively select certain equations that $w$ must satisfy and have only coefficients in the field $p^{r-1}\ZZ$. We exploit the natural block Toeplitz structure of the sliding parity-check matrix to derive a step by step methodology to obtain a list of possible codewords for a given corrupted codeword $w$, that is, a list with the closest codewords to $w$.


title: Laser-induced ultrafast demagnetization and perpendicular magnetic anisotropy reduction in a Co$_{88}$Tb$_{12}$ thin film with stripe domains

abstract: We use time-resolved x-ray resonant magnetic scattering (tr-XRMS) at the Co M$_{2,3}$- and Tb O$_1$-edges to study ultrafast demagnetization in an amorphous Co$_{88}$Tb$_{12}$ alloy with stripe domains. Combining the femtosecond temporal with nanometer spatial resolution of our experiment, we demonstrate that the equilibrium spin texture of the thin film remains unaltered by the optical pump-pulse on ultrashort timescales ($<$1 ps). However, after $\simeq$ 4 ps, we observe the onset of a significant domain wall broadening, which we attribute to a reduction of the uniaxial magnetic anisotropy of the system, due to energy transfer to the lattice. Static temperature dependent magnetometry measurements combined with analytical modeling of the magnetic structure of the thin film corroborate this interpretation.


title: Cycle completion times probe interactions with environment

abstract: Recent measurements of durations of non-equilibrium processes provide valuable information on microscopic mechanisms and energetics. Comprehensive theory for corresponding experiments so far is well developed for single-particle systems only. Little is known for interacting systems in non-equilibrium environments. Here, we introduce and study a basic model for cycle processes interacting with an environment that can exhibit a net particle flow. We find a surprising richness of cycle time variations with environmental conditions. This manifests itself in unequal cycle times $\tau^+$ and $\tau^-$ in forward and backward cycle direction with both asymmetries $\tau^-<\tau^+$ and $\tau^->\tau^+$, speeding up of backward cycles by interactions, and dynamical phase transitions, where cycle times become multimodal functions of the bias. The model allows us to relate these effects to specific microscopic mechanisms, which can be helpful for interpreting experiments.


title: An Alon-Boppana theorem for powered graphs and generalized Ramanujan graphs

abstract: The r-th power of a graph modifies a graph by connecting every vertex pair within distance r. This paper gives a generalization of the Alon-Boppana Theorem for the r-th power of graphs, including irregular graphs. This leads to a generalized notion of Ramanujan graphs, those for which the powered graph has a spectral gap matching the derived Alon-Boppana bound. In particular, we show that certain graphs that are not good expanders due to local irregularities, such as Erdos-Renyi random graphs, become almost Ramanujan once powered. A different generalization of Ramanujan graphs can also be obtained from the nonbacktracking operator. We next argue that the powering operator gives a more robust notion than the latter: Sparse Erdos-Renyi random graphs with an adversary modifying a subgraph of log(n)^c$ vertices are still almost Ramanujan in the powered sense, but not in the nonbacktracking sense. As an application, this gives robust community testing for different block models.


title: Dehn surgery and non-separating two-spheres

abstract: When can surgery on a null-homologous knot K in a rational homology sphere produce a non-separating sphere? We use Heegaard Floer homology to give sufficient conditions for K to be unknotted. We also discuss some applications to homology cobordism, concordance, and Mazur manifolds.


title: Light new physics in XENON1T

abstract: We examine the recently-reported low-energy electron recoil spectrum observed at the XENON1T underground dark matter direct detection experiment, in the context of new interactions with solar neutrinos. In particular we show that scalar and vector mediators with masses $\lesssim 50$ keV coupled to leptons could already leave a visible signature in the XENON1T experiment, similar to the observed peak below 7 keV. This signals that dark matter detectors are already competing with neutrino scattering experiments such as GEMMA, CHARM-II and Borexino. If these results from XENON1T are interpreted as a new signal of such physics, the parameters which fit the excess face challenges from astrophysics which seem very difficult to overcome. If they are rather viewed as a constraint on new couplings, they herald the start of an era of novel precise probes of physics beyond the standard model with dark matter detectors.


title: Halving spaces and lower bounds in real enumerative geometry

abstract: We develop the theory of halving spaces to obtain lower bounds in real enumerative geometry. Halving spaces are topological spaces with an action of a Lie group $\Gamma$ with additional cohomological properties. For $\Gamma=\mathbb{Z}_2$ we recover the conjugation spaces of Hausmann, Holm and Puppe. For $\Gamma=\mathrm{U}(1)$ we obtain the circle spaces. We show that real even and quaternionic partial flag manifolds are circle spaces leading to non-trivial lower bounds for even real and quaternionic Schubert problems. To prove that a given space is a halving space, we generalize results of Borel and Haefliger on the cohomology classes of real subvarieties and their complexifications. The novelty is that we are able to obtain results in rational cohomology instead of modulo 2. The equivariant extension of the theory of circle spaces leads to generalizations of the results of Borel and Haefliger on Thom polynomials.


title: Generalized Lewis-Riesenfeld invariance for dynamical effective mass in jammed granullar media under a potential well in non-commutative space

abstract: Consideration of the asteroid belt (Kuiper belt) as a jammed-granular media establishes a bridge between condensed matter physics and astrophysics. It opens up an experimental possibility to determine the deformation parameters for noncommutative space-time. Dynamics of the Kuiper belt can be simplified as dynamics of a dynamical effective mass for a jammed granular media under a gravitational well. Alongside, if one considers the space-time to be noncommutative, then an experimental model for the determination of the deformation parameters for noncommutative space-time can be done. The construction of eigenfunctions and invariance for this model is in general a tricky problem. We have utilized the Lewis-Riesenfeld invariant method to determine the invariance for this time-dependent quantum system. In this article, we have shown that a class of generalized time-dependent Lewis-Riesenfeld invariant operators exist for the system with dynamical effective mass in jammed granular media under a potential well in noncommutative space. To keep the discussion fairly general, we have considered both position-position and momentum-momentum noncommutativity. Since, up to a time-dependent phase-factor, the eigenfunctions of the invariant operator will satisfy the time-dependent Schr\"{o}dinger equation for the time-dependent Hamiltonian of the system, the construction of the invariant operator fairly solve the problem mathematically, the results of which can be utilized to demonstrate an experiment.


title: Operator expansions, layer susceptibility and two-point functions in BCFT

abstract: We show that in boundary CFTs, there exists a one-to-one correspondence between the boundary operator expansion of the two-point correlation function and a power series expansion of the layer susceptibility. This general property allows the direct identification of the boundary spectrum and expansion coefficients from the layer susceptibility and opens a new way for efficient calculations of two-point correlators in BCFTs. To show how it works we derive an explicit expression for the correlation function $\langle\phi_i \phi^i\rangle$ of the O(N) model at the extraordinary transition in 4-$\epsilon$ dimensional semi-infinite space to order $O(\epsilon)$. The bulk operator product expansion of the two-point function gives access to the spectrum of the bulk CFT. In our example, we obtain the averaged anomalous dimensions of scalar composite operators of the O(N) model to order $O(\epsilon^2)$. These agree with the known results both in $\epsilon$ and large-N expansions.


title: Physics of Psychophysics: two coupled square lattices of spiking neurons have huge dynamic range at criticality

abstract: Psychophysics try to relate physical input magnitudes to psychological or neural correlates. Microscopic models to account for macroscopic psychophysical laws, in the sense of statistical physics, are an almost unexplored area. Here we examine a sensory epithelium composed of two connected square lattices of stochastic integrate-and-fire cells. With one square lattice we obtain a Stevens's law $\rho \propto h^m$ with Stevens's exponent $m = 0.254$ and a sigmoidal saturation, where $\rho$ is the neuronal network activity and $h$ is the input intensity (external field). We relate Stevens's power law exponent with the field critical exponent as $m = 1/\delta_h = \beta/\sigma$. We also show that this system pertains to the Directed Percolation (DP) universality class (or perhaps the Compact-DP class). With stacked two layers of square lattices, and a fraction of connectivity between the first and second layer, we obtain at the output layer $\rho_ 2 \propto h^{m_2}$, with $m_2 = 0.08 \approx m^2$, which corresponds to a huge dynamic range. This enhancement of the dynamic range only occurs when the layers are close to their critical point.


title: A Revisit of Chen-Teboulle's Proximal-based Decomposition Method

abstract: In this paper, we show that Chen-Teboulle's proximal-based decomposition method can be interpreted as a proximal augmented Lagrangian method. More precisely, it coincides with a linearized augmented Lagrangian method. We then proposed three generalized methods based on this interpretation. By invoking recent work (He et al., IMA J. Numer. Anal., 32 (2020), pp. 227--245), we show that the step size condition of Chen-Teboulle's method can be relaxed without adding any further assumptions. Our analysis offers a new insight into this proximal-based decomposition method.


title: Large-scale parallel server system with multi-component jobs

abstract: A broad class of parallel server systems is considered, for which we prove the steady-state asymptotic independence of server workloads, as the number of servers goes to infinity, while the system load remains sub-critical. Arriving jobs consist of multiple components. There are multiple job classes, and each class may be of one of two types, which determines the rule according to which the job components add workloads to the servers. The model is broad enough to include as special cases some popular queueing models with redundancy, such as cancel-on-start and cancel-on-completion redundancy. Our analysis uses mean-field process representation and the corresponding mean-field limits. In essence, our approach relies almost exclusively on three fundamental properties of the model: (a) monotonicity, (b) work conservation, (c) the property that, on average, "new arriving workload prefers to go to servers with lower workloads."


title: Gravitational Waves Created During the EWPT

abstract: We study gravitational waves generated by bubble expansion created during the Cosmological Electroweak Phase Transition (EWPT). The energy radiated via gravitational waves is produced by the stress-energy tensor created by the magnetic field produced by bubble collisions during the EWPT, which occured about $10{-11}$ seconds after the Big Bang.


title: Tuning the structural and antiferromagnetic phase transitions in UCr$_{2}$Si$_2$: hydrostatic pressure and chemical substitution

abstract: Structural phase transitions in $f$-electron materials have attracted sustained attention both for practical and basic science reasons, including that they offer an environment to directly investigate relationships between structure and the $f$-state. Here we present results for UCr$_2$Si$_2$, where structural (tetragonal $\rightarrow$ monoclinic) and antiferromagnetic phase transitions are seen at $T_{\rm{S}}$ $=$ 205 K and $T_{\rm{N}}$ $=$ 25 K, respectively. We also provide evidence for an additional second order phase transition at $T_{\rm{X}}$ = 280 K. We show that $T_{\rm{X}}$, $T_{\rm{S}}$, and $T_{\rm{N}}$ respond in distinct ways to the application of hydrostatic pressure and Cr $\rightarrow$ Ru chemical substitution. In particular, hydrostatic compression increases the structural ordering temperature, eventually causes it to merge with $T_{\rm{X}}$ and destroys the antiferromagnetism. In contrast, chemical substitution in the series UCr$_{2-x}$Ru$_x$Si$_2$ suppresses both $T_{\rm{S}}$ and $T_{\rm{N}}$, causing them to approach zero temperature near $x$ $\approx$ 0.16 and 0.08, respectively. The distinct $T-P$ and $T-x$ phase diagrams are related to the evolution of the rigid Cr-Si and Si-Si substructures, where applied pressure semi-uniformly compresses the unit cell and Cr $\rightarrow$ Ru substitution results in uniaxial lattice compression along the tetragonal $c$-axis and an expansion in the $ab$-plane. These results provide insights into an interesting class of strongly correlated quantum materials where degrees of freedom associated with $f$-electron magnetism, strong electronic correlations, and structural instabilities are readily controlled.


title: Learning to Prove from Synthetic Theorems

abstract: A major challenge in applying machine learning to automated theorem proving is the scarcity of training data, which is a key ingredient in training successful deep learning models. To tackle this problem, we propose an approach that relies on training with synthetic theorems, generated from a set of axioms. We show that such theorems can be used to train an automated prover and that the learned prover transfers successfully to human-generated theorems. We demonstrate that a prover trained exclusively on synthetic theorems can solve a substantial fraction of problems in TPTP, a benchmark dataset that is used to compare state-of-the-art heuristic provers. Our approach outperforms a model trained on human-generated problems in most axiom sets, thereby showing the promise of using synthetic data for this task.