title: Real-time simulation of viscoelastic tissue behavior with physics-guided deep learning

abstract: Finite element methods (FEM) are popular approaches for simulation of soft tissues with elastic or viscoelastic behavior. However, their usage in real-time applications, such as in virtual reality surgical training, is limited by computational cost. In this application scenario, which typically involves transportable simulators, the computing hardware severely constrains the size or the level of details of the simulated scene. To address this limitation, data-driven approaches have been suggested to simulate mechanical deformations by learning the mapping rules from FEM generated datasets. Herein, we propose a deep learning method for predicting displacement fields of soft tissues with viscoelastic properties. The main contribution of this work is the use of a physics-guided loss function for the optimization of the deep learning model parameters. The proposed deep learning model is based on convolutional (CNN) and recurrent layers (LSTM) to predict spatiotemporal variations. It is augmented with a mass conservation law in the lost function to prevent the generation of physically inconsistent results. The deep learning model is trained on a set of FEM datasets that are generated from a commercially available state-of-the-art numerical neurosurgery simulator. The use of the physics-guided loss function in a deep learning model has led to a better generalization in the prediction of deformations in unseen simulation cases. Moreover, the proposed method achieves a better accuracy over the conventional CNN models, where improvements were observed in unseen tissue from 8% to 30% depending on the magnitude of external forces. It is hoped that the present investigation will help in filling the gap in applying deep learning in virtual reality simulators, hence improving their computational performance (compared to FEM simulations) and ultimately their usefulness.

categories: cs.LG cond-mat.soft physics.comp-ph q-bio.TO


title: Estimating entanglement in 2D Heisenberg model in the strong rung-coupling limit

abstract: In this paper, we calculate entanglement in the isotropic Heisenberg model in a magnetic field on a two-dimensional rectangular zig-zag lattice in the strong rung-coupling limit, using the one-dimensional XXZ model as a proxy. Focusing on the leading order in perturbation, for arbitrary size of the lattice, we show how the one-dimensional effective description emerges. We point out specific states in the low-energy sector of the two-dimensional model that are well-approximated by the one-dimensional spin-1/2 XXZ model. We propose a systematic approach for mapping matrix-elements of operators defined on the two-dimensional model to their low-energy counterparts on the one-dimensional XXZ model. We also show that partial trace-based description of entanglement in the two-dimensional model can be satisfactorily approximated using the one-dimensional XXZ model as a substitute. We further show numerically that the one-dimensional XXZ model performs well in estimating entanglement quantified using a measurement-based approach in the two-dimensional model for specific choices of measured Hermitian operators.


title: Neutrino Oscillation and Lattice QCD

abstract: Next generation high-precision neutrino scattering experiments have the goal of measuring the as-of-yet unknown parameters governing neutrino oscillation. This effort is hampered by the use of large nuclear targets: secondary interactions within a nucleus can confuse the interpretation of experimental data, leading to ambiguities about the initial neutrino interaction in scattering events. The distribution of energies for neutrino events must instead be inferred from the responses of a sum of dissimilar event topologies. For this reason, precise neutrino cross sections on nucleon targets are of vital importance to the neutrino oscillation experimental program. On the other hand, the necessary experimental data for neutrino scattering with elementary targets are scarce because of the weak interaction cross section, which leads to poorly-constrained nucleon and nuclear cross sections. Lattice QCD is uniquely positioned to provide the requisite nucleon amplitudes needed to enable high-precision oscillation experiments. In particular, LQCD has the ability to probe axial matrix elements that are challenging to isolate or completely inaccessible to experiments. In these proceedings, I will discuss some of my work to quantify neutrino cross sections with realistic uncertainty estimates, primarily focusing on neutrino quasielastic scattering and the nucleon axial form factor. I will also outline how the needs of next-generation neutrino oscillation experimental programs can be met with modern dedicated LQCD computations.


title: Observed polarized scattered light phase functions of planet-forming disks

abstract: Dust particles are the building blocks from which planetary bodies are made. A major goal of the studies of planet-forming disks is to constrain the properties of dust particles and aggregates in order to trace their origin, structure, and the associated growth and mixing processes in the disk. Observations of scattering and/or emission of dust in a location of the disk often lead to degenerate information about the kind of particles, such as size, porosity, or fractal dimension of aggregates. Progress can be made by deriving the full (polarizing) scattering phase function of such particles at multiple wavelengths. This has now become possible by careful extraction from scattered light images. Such an extraction requires knowledge about the shape of the scattering surface in the disk and we discuss how to obtain such knowledge as well as the associated uncertainties. We use a sample of disk images from observations with VLT/SPHERE to, for the first time, extract the phase functions of a whole sample of disks with broad phase angle coverage. We find that polarized phase functions come in two categories. Comparing the extracted functions with theoretical predictions from rigorous T-Matrix computations of aggregates, we show that one category can be linked back to fractal, porous aggregates, while the other is consistent with more compact, less porous aggregates. We speculate that the more compact particles become visible in disks where embedded planets trigger enhanced vertical mixing.


title: Pair production of heavy quarkonia in the color evaporation model

abstract: In the article, we study pair production of heavy quarkonia $(J/\psi J/\psi, \Upsilon\Upsilon, \Upsilon J/\psi)$ in the improved color evaporation model via the parton Reggeization approach. The last one is based on high-energy factorization of hard processes in multi-Regge kinematics, the Kimber-Martin-Ryskin-Watt model for unintegrated parton distribution functions, and the Lipatov effective field theory of Reggezied gluons and quarks. We compare contributions from the single and double parton scattering mechanisms in the pair production of heavy quarkonia. The numerical calculations are performed with the Monte-Carlo event generator KaTie.


title: TinyHD: Efficient Video Saliency Prediction with Heterogeneous Decoders using Hierarchical Maps Distillation

abstract: Video saliency prediction has recently attracted attention of the research community, as it is an upstream task for several practical applications. However, current solutions are particularly computationally demanding, especially due to the wide usage of spatio-temporal 3D convolutions. We observe that, while different model architectures achieve similar performance on benchmarks, visual variations between predicted saliency maps are still significant. Inspired by this intuition, we propose a lightweight model that employs multiple simple heterogeneous decoders and adopts several practical approaches to improve accuracy while keeping computational costs low, such as hierarchical multi-map knowledge distillation, multi-output saliency prediction, unlabeled auxiliary datasets and channel reduction with teacher assistant supervision. Our approach achieves saliency prediction accuracy on par or better than state-of-the-art methods on DFH1K, UCF-Sports and Hollywood2 benchmarks, while enhancing significantly the efficiency of the model. Code is on https://github.com/feiyanhu/tinyHD


title: AdaptSLAM: Edge-Assisted Adaptive SLAM with Resource Constraints via Uncertainty Minimization

abstract: Edge computing is increasingly proposed as a solution for reducing resource consumption of mobile devices running simultaneous localization and mapping (SLAM) algorithms, with most edge-assisted SLAM systems assuming the communication resources between the mobile device and the edge server to be unlimited, or relying on heuristics to choose the information to be transmitted to the edge. This paper presents AdaptSLAM, an edge-assisted visual (V) and visual-inertial (VI) SLAM system that adapts to the available communication and computation resources, based on a theoretically grounded method we developed to select the subset of keyframes (the representative frames) for constructing the best local and global maps in the mobile device and the edge server under resource constraints. We implemented AdaptSLAM to work with the state-of-the-art open-source V- and VI-SLAM ORB-SLAM3 framework, and demonstrated that, under constrained network bandwidth, AdaptSLAM reduces the tracking error by 62% compared to the best baseline method.


title: Entropy of different phases formed by soft rods

abstract: Computation of entropy in liquids and liquid crystal phases is a big challenge in statistical physics. In this work, we extend the two-phase thermodynamic model (2PT) to shape anisotropic soft repulsive spherocylinders (SRSs) and report the absolute values of entropy for different liquid crystal (LC) phases for a range of aspect ratios L/D = 2-5. We calculate the density of states (DoS) for different LC phases and decompose it into contributions arising from translational and rotational degrees of freedom. The translational and rotational modes are further partitioned into diffusive, gas-like, and non-diffusive, solid-like components using a fluidicity factor. In the dilute limit, the entropy values obtained from the 2PT method match exactly those of an ideal rigid rotor. We find that, for a given packing fraction, the magnitude of the total entropy is roughly equal regardless of the different LC phases associated with different aspect ratios. We also compute the excess entropy (for L/D = 5) and compare those with the values obtained using the standard integration approach of molecular dynamics (MD) or Monte Carlo (MC) equation of state (EOS) of SRS. The values obtained using both approaches match very well. The rotational and translational fluidicity factors are further used to determine the phase boundaries of different liquid crystal phases for the respective aspect ratios.


title: Ion filling of a one-dimensional nanofluidic channel in the interaction confinement regime

abstract: Ion transport measurements are widely used as an indirect probe for various properties of confined electrolytes. It is generally assumed that the ion concentration in a nanoscale channel is equal to the ion concentration in the macroscopic reservoirs it connects to, with deviations arising only in the presence of surface charges on the channel walls. Here, we show that this assumption may break down even in a neutral channel, due to electrostatic correlations between the ions arising in the regime of interaction confinement, where Coulomb interactions are reinforced due to the presence of the channel walls. We focus on a one-dimensional channel geometry, where an exact evaluation of the electrolyte's partition function is possible with a transfer operator approach. Our exact solution reveals that in nanometre-scale channels, the ion concentration is generally lower than in the reservoirs, and depends continuously on the bulk salt concentration, in contrast to conventional mean-field theory that predicts an abrupt filling transition. We develop a modified mean-field theory taking into account the presence of ion pairs that agrees quantitatively with the exact solution and provides predictions for experimentally-relevant observables such as the ionic conductivity. Our results will guide the interpretation of nanoscale ion transport measurements.

categories: cond-mat.soft cond-mat.mes-hall cond-mat.stat-mech physics.chem-ph


title: Enhancing ResNet Image Classification Performance by using Parameterized Hypercomplex Multiplication

abstract: Recently, many deep networks have introduced hypercomplex and related calculations into their architectures. In regard to convolutional networks for classification, these enhancements have been applied to the convolution operations in the frontend to enhance accuracy and/or reduce the parameter requirements while maintaining accuracy. Although these enhancements have been applied to the convolutional frontend, it has not been studied whether adding hypercomplex calculations improves performance when applied to the densely connected backend. This paper studies ResNet architectures and incorporates parameterized hypercomplex multiplication (PHM) into the backend of residual, quaternion, and vectormap convolutional neural networks to assess the effect. We show that PHM does improve classification accuracy performance on several image datasets, including small, low-resolution CIFAR 10/100 and large high-resolution ImageNet and ASL, and can achieve state-of-the-art accuracy for hypercomplex networks.


title: Typical Correlation Length of Sequentially Generated Tensor Network States

abstract: The complexity of quantum many-body systems is manifested in the vast diversity of their correlations, making it challenging to distinguish the generic from the atypical features. This can be addressed by analyzing correlations through ensembles of random states, chosen to faithfully embody the relevant physical properties. Here, we focus on spins with local interactions, whose correlations are extremely well captured by tensor network states. Adopting an operational perspective, we define ensembles of random tensor network states in one and two spatial dimensions that admit a sequential generation. As such, they directly correspond to outputs of quantum circuits with a sequential architecture and random gates. In one spatial dimension, the ensemble explores the entire family of matrix product states, while in two spatial dimensions, it corresponds to random isometric tensor network states. We extract the scaling behavior of the average correlations between two subsystems as a function of their distance. Using elementary concentration results, we then deduce the typical case for measures of correlation such as the von Neumann mutual information and a measure arising from the Hilbert-Schmidt norm. We find for all considered cases that the typical behavior is an exponential decay (for both one and two spatial dimensions). We observe the consistent emergence of a correlation length that depends only on the underlying spatial dimension and not the considered measure. Remarkably, increasing the bond dimension leads to a higher correlation length in one spatial dimension but has the opposite effect in two spatial dimensions.


title: Multivariate Regression via Enhanced Response Envelope: Envelope Regularization and Double Descent

abstract: The envelope model provides substantial efficiency gains over the standard multivariate linear regression by identifying the material part of the response to the model and by excluding the immaterial part. In this paper, we propose the enhanced response envelope by incorporating a novel envelope regularization term in its formulation. It is shown that the enhanced response envelope can yield better prediction risk than the original envelope estimator. The enhanced response envelope naturally handles high-dimensional data for which the original response envelope is not serviceable without necessary remedies. In an asymptotic high-dimensional regime where the ratio of the number of predictors over the number of samples converges to a non-zero constant, we characterize the risk function and reveal an interesting double descent phenomenon for the first time for the envelope model. A simulation study confirms our main theoretical findings. Simulations and real data applications demonstrate that the enhanced response envelope does have significantly improved prediction performance over the original envelope method.


title: Deep Axial Hypercomplex Networks

abstract: Over the past decade, deep hypercomplex-inspired networks have enhanced feature extraction for image classification by enabling weight sharing across input channels. Recent works make it possible to improve representational capabilities by using hypercomplex-inspired networks which consume high computational costs. This paper reduces this cost by factorizing a quaternion 2D convolutional module into two consecutive vectormap 1D convolutional modules. Also, we use 5D parameterized hypercomplex multiplication based fully connected layers. Incorporating both yields our proposed hypercomplex network, a novel architecture that can be assembled to construct deep axial-hypercomplex networks (DANs) for image classifications. We conduct experiments on CIFAR benchmarks, SVHN, and Tiny ImageNet datasets and achieve better performance with fewer trainable parameters and FLOPS. Our proposed model achieves almost 2% higher performance for CIFAR and SVHN datasets, and more than 3% for the ImageNet-Tiny dataset and takes six times fewer parameters than the real-valued ResNets. Also, it shows state-of-the-art performance on CIFAR benchmarks in hypercomplex space.


title: Improved Approximations for Extremal Eigenvalues of Sparse Hamiltonians

abstract: We give a classical $1/(qk+1)$-approximation for the maximum eigenvalue of a $k$-sparse fermionic Hamiltonian with strictly $q$-local terms, as well as a $1/(4k+1)$-approximation when the Hamiltonian has both $2$-local and $4$-local terms. More generally we obtain a $1/O(qk^2)$-approximation for $k$-sparse fermionic Hamiltonians with terms of locality at most $q$. Our techniques also yield analogous approximations for $k$-sparse, $q$-local qubit Hamiltonians with small hidden constants and improved dependence on $q$.


title: Face Attribute Editing with Disentangled Latent Vectors

abstract: We propose an image-to-image translation framework for facial attribute editing with disentangled interpretable latent directions. Facial attribute editing task faces the challenges of targeted attribute editing with controllable strength and disentanglement in the representations of attributes to preserve the other attributes during edits. For this goal, inspired by the latent space factorization works of fixed pretrained GANs, we design the attribute editing by latent space factorization, and for each attribute, we learn a linear direction that is orthogonal to the others. We train these directions with orthogonality constraints and disentanglement losses. To project images to semantically organized latent spaces, we set an encoder-decoder architecture with attention-based skip connections. We extensively compare with previous image translation algorithms and editing with pretrained GAN works. Our extensive experiments show that our method significantly improves over the state-of-the-arts. Project page: https://yusufdalva.github.io/vecgan


title: On some identities for confluent hypergeometric functions and Bessel functions

abstract: Mathematical functions, which often appear in mathematical analysis, are referred to as special functions and have been studied over hundreds of years. Many books and dictionaries are available that describe their properties and serve as a foundation of current science. In this paper, we find a new integral representation of the Whittaker function of the first kind and show a relevant summation formula for Kummer's confluent hypergeometric functions. We also perform the specifications of our identities to link to known and new results.


title: ShadowNav: Crater-Based Localization for Nighttime and Permanently Shadowed Region Lunar Navigation

abstract: There has been an increase in interest in missions that drive significantly longer distances per day than what has currently been performed. Further, some of these proposed missions require autonomous driving and absolute localization in darkness. For example, the Endurance A mission proposes to drive 1200km of its total traverse at night. The lack of natural light available during such missions limits what can be used as visual landmarks and the range at which landmarks can be observed. In order for planetary rovers to traverse long ranges, onboard absolute localization is critical to the ability of the rover to maintain its planned trajectory and avoid known hazardous regions. Currently, to accomplish absolute localization, a ground in the loop (GITL) operation is performed wherein a human operator matches local maps or images from onboard with orbital images and maps. This GITL operation limits the distance that can be driven in a day to a few hundred meters, which is the distance that the rover can maintain acceptable localization error via relative methods. Previous work has shown that using craters as landmarks is a promising approach for performing absolute localization on the moon during the day. In this work we present a method of absolute localization that utilizes craters as landmarks and matches detected crater edges on the surface with known craters in orbital maps. We focus on a localization method based on a perception system which has an external illuminator and a stereo camera. We evaluate (1) both monocular and stereo based surface crater edge detection techniques, (2) methods of scoring the crater edge matches for optimal localization, and (3) localization performance on simulated Lunar surface imagery at night. We demonstrate that this technique shows promise for maintaining absolute localization error of less than 10m required for most planetary rover missions.


title: Deep Residual Axial Networks

abstract: While convolutional neural networks (CNNs) demonstrate outstanding performance on computer vision tasks, their computational costs remain high. Several techniques are used to reduce these costs, like reducing channel count, and using separable and depthwise separable convolutions. This paper reduces computational costs by introducing a novel architecture, axial CNNs, which replaces spatial 2D convolution operations with two consecutive depthwise separable 1D operations. The axial CNNs are predicated on the assumption that the dataset supports approximately separable convolution operations with little or no loss of training accuracy. Deep axial separable CNNs still suffer from gradient problems when training deep networks. We modify the construction of axial separable CNNs with residual connections to improve the performance of deep axial architectures and introduce our final novel architecture namely residual axial networks (RANs). Extensive benchmark evaluation shows that RANs achieve at least 1% higher performance with about 77%, 86%, 75%, and 34% fewer parameters and about 75%, 80%, 67%, and 26% fewer flops than ResNets, wide ResNets, MobileNets, and SqueezeNexts on CIFAR benchmarks, SVHN, and Tiny ImageNet image classification datasets. Moreover, our proposed RANs improve deep recursive residual networks performance with 94% fewer parameters on the image super-resolution dataset.


title: Federated Learning under Heterogeneous and Correlated Client Availability

abstract: The enormous amount of data produced by mobile and IoT devices has motivated the development of federated learning (FL), a framework allowing such devices (or clients) to collaboratively train machine learning models without sharing their local data. FL algorithms (like FedAvg) iteratively aggregate model updates computed by clients on their own datasets. Clients may exhibit different levels of participation, often correlated over time and with other clients. This paper presents the first convergence analysis for a FedAvg-like FL algorithm under heterogeneous and correlated client availability. Our analysis highlights how correlation adversely affects the algorithm's convergence rate and how the aggregation strategy can alleviate this effect at the cost of steering training toward a biased model. Guided by the theoretical analysis, we propose CA-Fed, a new FL algorithm that tries to balance the conflicting goals of maximizing convergence speed and minimizing model bias. To this purpose, CA-Fed dynamically adapts the weight given to each client and may ignore clients with low availability and large correlation. Our experimental results show that CA-Fed achieves higher time-average accuracy and a lower standard deviation than state-of-the-art AdaFed and F3AST, both on synthetic and real datasets.


title: Accelerating Machine Learning Inference with GPUs in ProtoDUNE Data Processing

abstract: We study the performance of a cloud-based GPU-accelerated inference server to speed up event reconstruction in neutrino data batch jobs. Using detector data from the ProtoDUNE experiment and employing the standard DUNE grid job submission tools, we attempt to reprocess the data by running several thousand concurrent grid jobs, a rate we expect to be typical of current and future neutrino physics experiments. We process most of the dataset with the GPU version of our processing algorithm and the remainder with the CPU version for timing comparisons. We find that a 100-GPU cloud-based server is able to easily meet the processing demand, and that using the GPU version of the event processing algorithm is two times faster than processing these data with the CPU version when comparing to the newest CPUs in our sample. The amount of data transferred to the inference server during the GPU runs can overwhelm even the highest-bandwidth network switches, however, unless care is taken to observe network facility limits or otherwise distribute the jobs to multiple sites. We discuss the lessons learned from this processing campaign and several avenues for future improvements.


title: Street-View Image Generation from a Bird's-Eye View Layout

abstract: Bird's-Eye View (BEV) Perception has received increasing attention in recent years as it provides a concise and unified spatial representation across views and benefits a diverse set of downstream driving applications. At the same time, data-driven simulation for autonomous driving has been a focal point of recent research but with few approaches that are both fully data-driven and controllable. Instead of using perception data from real-life scenarios, an ideal model for simulation would generate realistic street-view images that align with a given HD map and traffic layout, a task that is critical for visualizing complex traffic scenarios and developing robust perception models for autonomous driving. In this paper, we propose BEVGen, a conditional generative model that synthesizes a set of realistic and spatially consistent surrounding images that match the BEV layout of a traffic scenario. BEVGen incorporates a novel cross-view transformation with spatial attention design which learns the relationship between cameras and map views to ensure their consistency. We evaluate the proposed model on the challenging NuScenes and Argoverse 2 datasets. After training, BEVGen can accurately render road and lane lines, as well as generate traffic scenes with diverse different weather conditions and times of day.


title: On The Determination of Sets By Their Subset Sums

abstract: Let $A$ be a multiset with elements in an abelian group. Let $FS(A)$ be the multiset containing the $2^{|A|}$ sums of all subsets of $A$. We study the reconstruction problem ``Given $FS(A)$, is it possible to identify $A$?'', and we give a satisfactory answer for all abelian groups. We prove that, up to identifying multisets through a natural equivalence relation, the function $A \mapsto FS(A)$ is injective (and thus the reconstruction problem is solvable) if and only if every order $n$ of a torsion element of the abelian group satisfies a certain number-theoretical property linked to the multiplicative group $(\mathbb{Z} / n\mathbb{Z})^*$. The core of the proof relies on a delicate study of the structure of cyclotomic units. Moreover, as a tool, we develop an inversion formula for a novel discrete Radon transform on finite abelian groups that might be of independent interest.


title: Unavoidable structures in infinite tournaments

abstract: We prove a strong dichotomy result for countably-infinite oriented graphs; that is, we prove that for all countably-infinite oriented graphs $G$, either (i) there is a countably-infinite tournament $K$ such that $G\not\subseteq K$, or (ii) every countably-infinite tournament contains a \emph{spanning} copy of $G$. Furthermore, we are able to give a concise characterization of such oriented graphs. Our characterization becomes even simpler in the case of transitive acyclic oriented graphs (i.e. strict partial orders). For uncountable oriented graphs, we are able to extend the dichotomy result mentioned above to all regular cardinals $\kappa$; however, we are only able to provide a concise characterization in the case when $\kappa=\aleph_1$.


title: A Systematic Study of Ia-CSM Supernovae from the ZTF Bright Transient Survey

abstract: Among the supernovae (SNe) that show strong interaction with the circumstellar medium, there is a rare subclass of Type Ia supernovae, SNe Ia-CSM, that show strong narrow hydrogen emission lines much like SNe IIn but on top of a diluted over-luminous Type Ia spectrum. In the only previous systematic study of this class (Silverman et al. 2013), 16 objects were identified, 8 historic and 8 from the Palomar Transient Factory (PTF). Now using the successor survey to PTF, the Zwicky Transient Facility (ZTF), we have classified 12 additional objects of this type through the systematic Bright Transient Survey (BTS). In this study, we present and analyze the optical and mid-IR light curves, optical spectra, and host galaxy properties of this sample. Consistent with previous studies, we find the objects to have slowly evolving light curves compared to normal SNe Ia with peak absolute magnitudes between -19.1 and -21, spectra having weak H$\beta$, large Balmer decrements of ~7 and strong Ca NIR emission. Out of 10 SNe from our sample observed by NEOWISE, 9 have $3\sigma$ detections, along with some showing a clear reduction in red-wing of H$\alpha$, indicative of newly formed dust. We do not find our SN Ia-CSM sample to have a significantly different distribution of equivalent width of He I $\lambda5876$ than SNe IIn as observed in Silverman et al. 2013. The hosts tend to be late-type galaxies with recent star formation. We also derive a rate estimate of 29$^{+27}_{-21}$ Gpc$^{-3}$ yr$^{-1}$ for SNe Ia-CSM which is ~0.02--0.2 % of the SN Ia rate. This work nearly doubles the sample of well-studied Ia-CSM objects in Silverman et al. 2013, increasing the total number to 28.


title: Dynamics of a data-driven low-dimensional model of turbulent minimal Couette flow

abstract: Because the Navier-Stokes equations are dissipative, the long-time dynamics of a flow in state space are expected to collapse onto a manifold whose dimension may be much lower than the dimension required for a resolved simulation. On this manifold, the state of the system can be exactly described in a coordinate system parameterizing the manifold. Describing the system in this low-dimensional coordinate system allows for much faster simulations and analysis. We show, for turbulent Couette flow, that this description of the dynamics is possible using a data-driven manifold dynamics modeling method. This approach consists of an autoencoder to find a low-dimensional manifold coordinate system and a set of ordinary differential equations defined by a neural network. Specifically, we apply this method to minimal flow unit turbulent plane Couette flow at $\textit{Re}=400$, where a fully resolved solutions requires $\mathcal{O}(10^5)$ degrees of freedom. Using only data from this simulation we build models with fewer than $20$ degrees of freedom that quantitatively capture key characteristics of the flow, including the streak breakdown and regeneration cycle. At short-times, the models track the true trajectory for multiple Lyapunov times, and, at long-times, the models capture the Reynolds stress and the energy balance. For comparison, we show that the models outperform POD-Galerkin models with $\sim$2000 degrees of freedom. Finally, we compute unstable periodic orbits from the models. Many of these closely resemble previously computed orbits for the full system; additionally, we find nine orbits that correspond to previously unknown solutions in the full system.


title: Experimental verification of the temperature coefficient of resistivity

abstract: We have created an experimental procedure for determining the temperature coefficient of resistivity, $\alpha_R$, for introductory physics laboratories. This method examines the relationship between temperature and resistivity to establish $\alpha_R$ within 10% of the accepted value.


title: Properties of the multi-index special function $\mathcal{W}^{\left(\bar{\alpha},\bar{\nu}\right)}(z)$

abstract: In this paper, we investigate some properties related to a multi-index special function $\mathcal{W}^{\left(\bar{\alpha},\bar{\nu}\right)}$ that arose from an eigenvalue problem for a multi-order fractional hyper-Bessel operator, involving Caputo fractional derivatives. We show that for particular values of the parameters involved in this special function $\mathcal{W}^{\left(\bar{\alpha},\bar{\nu}\right)}$, this leads to the hyper-Bessel function of Delerue. The Laplace transform of the $\mathcal{W}^{\left(\bar{\alpha},\bar{\nu}\right)}$ is discussed obtaining, in particular cases, the well-known functional relation between hyper-Bessel function and multi-index Mittag-Leffler function, or, quite simply, between classical Wright and Mittag-Leffler functions. Moreover, it is shown that the multi-index special function satisfies the recurrence relation involving fractional derivatives. In a particular case, we derive, to the best of our knowledge, a new differential recurrence relation for the Mittag-Leffler function. We also provide derivatives of the 3-parameters function $\mathcal{W}_{\alpha,\beta,\nu}$ with respect to parameters, leading to infinite power series with coefficients being quotients of digamma and gamma functions.


title: Plug-in Channel Estimation with Dithered Quantized Signals in Spatially Non-Stationary Massive MIMO Systems

abstract: As the array dimension of massive MIMO systems increases to unprecedented levels, two problems occur. First, the spatial stationarity assumption along the antenna elements is no longer valid. Second, the large array size results in an unacceptably high power consumption if high-resolution analog-to-digital converters are used. To address these two challenges, we consider a Bussgang linear minimum mean square error (BLMMSE)-based channel estimator for large scale massive MIMO systems with one-bit quantizers and a spatially non-stationary channel. Whereas other works usually assume that the channel covariance is known at the base station, we consider a plug-in BLMMSE estimator that uses an estimate of the channel covariance and rigorously analyze the distortion produced by using an estimated, rather than the true, covariance. To cope with the spatial non-stationarity, we introduce dithering into the quantized signals and provide a theoretical error analysis. In addition, we propose an angular domain fitting procedure which is based on solving an instance of non-negative least squares. For the multi-user data transmission phase, we further propose a BLMMSE-based receiver to handle one-bit quantized data signals. Our numerical results show that the performance of the proposed BLMMSE channel estimator is very close to the oracle-aided scheme with ideal knowledge of the channel covariance matrix. The BLMMSE receiver outperforms the conventional maximum-ratio-combining and zero-forcing receivers in terms of the resulting ergodic sum rate.


title: Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

abstract: Quantum mechanical calculations for material modelling using Kohn-Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for $N$ smallest eigenvector-eigenvalue pairs with $N$ proportional to the number of electrons in the material system. These calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn-Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh-Ritz projection step of the discretized Kohn-Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node CPU and GPU architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behaviour and energy efficiency.


title: tieval: An Evaluation Framework for Temporal Information Extraction Systems

abstract: Temporal information extraction (TIE) has attracted a great deal of interest over the last two decades, leading to the development of a significant number of datasets. Despite its benefits, having access to a large volume of corpora makes it difficult when it comes to benchmark TIE systems. On the one hand, different datasets have different annotation schemes, thus hindering the comparison between competitors across different corpora. On the other hand, the fact that each corpus is commonly disseminated in a different format requires a considerable engineering effort for a researcher/practitioner to develop parsers for all of them. This constraint forces researchers to select a limited amount of datasets to evaluate their systems which consequently limits the comparability of the systems. Yet another obstacle that hinders the comparability of the TIE systems is the evaluation metric employed. While most research works adopt traditional metrics such as precision, recall, and $F_1$, a few others prefer temporal awareness -- a metric tailored to be more comprehensive on the evaluation of temporal systems. Although the reason for the absence of temporal awareness in the evaluation of most systems is not clear, one of the factors that certainly weights this decision is the necessity to implement the temporal closure algorithm in order to compute temporal awareness, which is not straightforward to implement neither is currently easily available. All in all, these problems have limited the fair comparison between approaches and consequently, the development of temporal extraction systems. To mitigate these problems, we have developed tieval, a Python library that provides a concise interface for importing different corpora and facilitates system evaluation. In this paper, we present the first public release of tieval and highlight its most relevant features.


abstract: Does progress on ImageNet transfer to real-world datasets? We investigate this question by evaluating ImageNet pre-trained models with varying accuracy (57% - 83%) on six practical image classification datasets. In particular, we study datasets collected with the goal of solving real-world tasks (e.g., classifying images from camera traps or satellites), as opposed to web-scraped benchmarks collected for comparing models. On multiple datasets, models with higher ImageNet accuracy do not consistently yield performance improvements. For certain tasks, interventions such as data augmentation improve performance even when architectures do not. We hope that future benchmarks will include more diverse datasets to encourage a more comprehensive approach to improving learning algorithms.


title: Vertical projections in the Heisenberg group for sets of dimension greater than 3

abstract: It is shown that vertical projections in the Heisenberg group of sets of dimension strictly greater than 3 almost surely have positive area. The proof uses the point-plate incidence method introduced by F\"assler and Orponen, and also uses a similar approach to a recent maximal inequality of Zahl for fractal families of tubes. It relies on the endpoint trilinear Kakeya inequality in $\mathbb{R}^3$. Some related results are given on generic intersections with horizontal lines.


title: HAWC Detection of a TeV Halo Candidate Surrounding a Radio-quiet pulsar

abstract: Extended very-high-energy (VHE; 0.1-100 TeV) $\gamma$-ray emission has been observed around several middle-aged pulsars and referred to as ``TeV halos". Their formation mechanism remains under debate. It is also unknown whether they are ubiquitous or related to certain subgroup of pulsars. With 2321 days of observation, the High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory detected VHE $\gamma$-ray emission at the location of the radio-quiet pulsar PSR J0359+5414 with $>6\sigma$ significance. By performing likelihood tests with different spectral and spatial models and comparing the TeV spectrum with multi-wavelength observations of nearby sources, we show that this excess is consistent with a TeV halo associated with PSR J0359+5414, though future observation of HAWC and multi-wavelength follow-ups are needed to confirm this nature. This new halo candidate is located in a non-crowded region in the outer Galaxy. It shares similar properties to the other halos but its pulsar is younger and radio-quiet. Our observation implies that TeV halos could commonly exist around pulsars and their formation does not depend on the configuration of the pulsar magnetosphere.


title: EXIF as Language: Learning Cross-Modal Associations Between Images and Camera Metadata

abstract: We learn a visual representation that captures information about the camera that recorded a given photo. To do this, we train a multimodal embedding between image patches and the EXIF metadata that cameras automatically insert into image files. Our model represents this metadata by simply converting it to text and then processing it with a transformer. The features that we learn significantly outperform other self-supervised and supervised features on downstream image forensics and calibration tasks. In particular, we successfully localize spliced image regions "zero shot" by clustering the visual embeddings for all of the patches within an image.


title: Head-Free Lightweight Semantic Segmentation with Linear Transformer

abstract: Existing semantic segmentation works have been mainly focused on designing effective decoders; however, the computational load introduced by the overall structure has long been ignored, which hinders their applications on resource-constrained hardwares. In this paper, we propose a head-free lightweight architecture specifically for semantic segmentation, named Adaptive Frequency Transformer. It adopts a parallel architecture to leverage prototype representations as specific learnable local descriptions which replaces the decoder and preserves the rich image semantics on high-resolution features. Although removing the decoder compresses most of the computation, the accuracy of the parallel structure is still limited by low computational resources. Therefore, we employ heterogeneous operators (CNN and Vision Transformer) for pixel embedding and prototype representations to further save computational costs. Moreover, it is very difficult to linearize the complexity of the vision Transformer from the perspective of spatial domain. Due to the fact that semantic segmentation is very sensitive to frequency information, we construct a lightweight prototype learning block with adaptive frequency filter of complexity $O(n)$ to replace standard self attention with $O(n^{2})$. Extensive experiments on widely adopted datasets demonstrate that our model achieves superior accuracy while retaining only 3M parameters. On the ADE20K dataset, our model achieves 41.8 mIoU and 4.6 GFLOPs, which is 4.4 mIoU higher than Segformer, with 45% less GFLOPs. On the Cityscapes dataset, our model achieves 78.7 mIoU and 34.4 GFLOPs, which is 2.5 mIoU higher than Segformer with 72.5% less GFLOPs. Code is available at https://github.com/dongbo811/AFFormer.


title: The alignment of galaxies at the Baryon Acoustic Oscillation scale

abstract: Massive elliptical galaxies align pointing their major axis towards each other in the structure of the Universe. Such alignments are well-described at large scales through a linear relation with respect to the tidal field of the large-scale structure. At such scales, galaxy alignments are sensitive to the presence of baryon acoustic oscillations (BAO). The shape of the BAO feature in galaxy alignment correlations differs from the traditional peak in the clustering correlation function. Instead, it appears as a trough feature at the BAO scale. In this work, we show that this feature can be explained by a simple toy model of tidal fields from a spherical shell of matter. This helps give a physical insight for the feature and highlights the need for tailored template-based identification methods for the BAO in alignment statistics. We also discuss the impact of projection baselines and photometric redshift uncertainties for identifying the BAO in intrinsic alignment measurements.


title: Geometry-biased Transformers for Novel View Synthesis

abstract: We tackle the task of synthesizing novel views of an object given a few input images and associated camera viewpoints. Our work is inspired by recent 'geometry-free' approaches where multi-view images are encoded as a (global) set-latent representation, which is then used to predict the color for arbitrary query rays. While this representation yields (coarsely) accurate images corresponding to novel viewpoints, the lack of geometric reasoning limits the quality of these outputs. To overcome this limitation, we propose 'Geometry-biased Transformers' (GBTs) that incorporate geometric inductive biases in the set-latent representation-based inference to encourage multi-view geometric consistency. We induce the geometric bias by augmenting the dot-product attention mechanism to also incorporate 3D distances between rays associated with tokens as a learnable bias. We find that this, along with camera-aware embeddings as input, allows our models to generate significantly more accurate outputs. We validate our approach on the real-world CO3D dataset, where we train our system over 10 categories and evaluate its view-synthesis ability for novel objects as well as unseen categories. We empirically validate the benefits of the proposed geometric biases and show that our approach significantly improves over prior works.


title: 20736-node Weighted Max-Cut Problem Solving by Quadrature Photonic Spatial Ising Machine

abstract: To tackle challenging combinatorial optimization problems, analog computing machines based on the nature-inspired Ising model are attracting increasing attentions in order to disruptively overcome the impending limitations on conventional electronic computers. Photonic spatial Ising machine has become an unique and primitive solution with all-to-all connections to solve large-scale Max-cut problems. However, spin configuration and flipping requires two independent sets of spatial light modulators (SLMs) for amplitude and phase modulation, which will lead to tremendous engineering difficulty of optical alignment and coupling. We report a novel quadrature photonic spatial-Euler Ising machine to realize large-scale and flexible spin-interaction configuration and spin-flip in a single spatial light modulator, and develop a noise enhancement approach by adding digital white noise onto detected optical signals. We experimentally show that such proposal accelerates solving (un)weighted, (non)fully connected, 20736-node Max-cut problems, which offers obvious advantages over simulation and heuristic algorithm results in digital computers.


title: Estimate Deformation Capacity of Non-Ductile RC Shear Walls using Explainable Boosting Machine

abstract: Machine learning is becoming increasingly prevalent for tackling challenges in earthquake engineering and providing fairly reliable and accurate predictions. However, it is mostly unclear how decisions are made because machine learning models are generally highly sophisticated, resulting in opaque black-box models. Machine learning models that are naturally interpretable and provide their own decision explanation, rather than using an explanatory, are more accurate in determining what the model actually computes. With this motivation, this study aims to develop a fully explainable machine learning model to predict the deformation capacity of non-ductile reinforced concrete shear walls based on experimental data collected worldwide. The proposed Explainable Boosting Machines (EBM)-based model is an interpretable, robust, naturally explainable glass-box model, yet provides high accuracy comparable to its black-box counterparts. The model enables the user to observe the relationship between the wall properties and the deformation capacity by quantifying the individual contribution of each wall property as well as the correlations among them. The mean coefficient of determination R2 and the mean ratio of predicted to actual value based on the test dataset are 0.92 and 1.05, respectively. The proposed predictive model stands out with its overall consistency with scientific knowledge, practicality, and interpretability without sacrificing high accuracy.


title: Optirank: classification for RNA-Seq data with optimal ranking reference genes

abstract: Classification algorithms using RNA-Sequencing (RNA-Seq) data as input are used in a variety of biological applications. By nature, RNA-Seq data is subject to uncontrolled fluctuations both within and especially across datasets, which presents a major difficulty for a trained classifier to generalize to an external dataset. Replacing raw gene counts with the rank of gene counts inside an observation has proven effective to mitigate this problem. However, the rank of a feature is by definition relative to all other features, including highly variable features that introduce noise in the ranking. To address this problem and obtain more robust ranks, we propose a logistic regression model, optirank, which learns simultaneously the parameters of the model and the genes to use as a reference set in the ranking. We show the effectiveness of this method on simulated data. We also consider real classification tasks, which present different kinds of distribution shifts between train and test data. Those tasks concern a variety of applications, such as cancer of unknown primary classification, identification of specific gene signatures, and determination of cell type in single-cell RNA-Seq datasets. On those real tasks, optirank performs at least as well as the vanilla logistic regression on classical ranks, while producing sparser solutions. In addition, to increase the robustness against dataset shifts, we propose a multi-source learning scheme and demonstrate its effectiveness when used in combination with rank-based classifiers.


title: Classical and Quantum Elliptical Billiards: Mixed Phase Space and Short Correlations in Singlets and Doublets

abstract: Billiards are flat cavities where a particle is free to move between elastic collisions with the boundary. In chaos theory these systems are simple prototypes, their conservative dynamics of a billiard may vary from regular to chaotic, depending only on the border. The results reported here seek to shed light on the quantization of classically chaotic systems. We present numerical results on classical and quantum properties in two bi-parametric families of Billiards, Elliptical Stadium Billiard (ESB) and Elliptical-$C_3$ Billiards (E-$C_3$B). Both are elliptical perturbations of chaotic billiards with originally circular sectors on their borders. Our numerical calculations show evidence that the elliptical families can present a mixed classical phase space, identified by a parameter $\rho_\text{c} < 1$, which we use to guide our analysis of quantum spectra. We explored the short correlations through nearest neighbor spacing distribution $p(s)$, which showed that in the mixed region of the classical phase space, $p(s)$ is well described by the Berry-Robnik-Brody (BRB) distributions for the ESB. In agreement with the expected from the so-called ergodic parameter $\alpha = t_\text{H}/t_\text{T}$, the ratio between the Heisenberg time and the classical diffusive-like transport time signals the possibility of quantum dynamical localization when $\alpha < 1$. For the E-$C_3$B family, the eigenstates can be split into singlets and doublets. BRB describes $p(s)$ for singlets as the previous family in the mixed region. However, the $p(s)$ for doublets are described by new distributions recently introduced in the literature but only tested in a few cases for $\rho_\text{c} < 1$. We observed that as $\rho_\text{c}$ decreases, the $p(s)$'s tend to move away simultaneously from the GOE (singlets) and GUE (doublets) distributions.


title: ChatGPT is not all you need. A State of the Art Review of large Generative AI models

abstract: During the last two years there has been a plethora of large generative models such as ChatGPT or Stable Diffusion that have been published. Concretely, these models are able to perform tasks such as being a general question and answering system or automatically creating artistic images that are revolutionizing several sectors. Consequently, the implications that these generative models have in the industry and society are enormous, as several job positions may be transformed. For example, Generative AI is capable of transforming effectively and creatively texts to images, like the DALLE-2 model; text to 3D images, like the Dreamfusion model; images to text, like the Flamingo model; texts to video, like the Phenaki model; texts to audio, like the AudioLM model; texts to other texts, like ChatGPT; texts to code, like the Codex model; texts to scientific texts, like the Galactica model or even create algorithms like AlphaTensor. This work consists on an attempt to describe in a concise way the main models are sectors that are affected by generative AI and to provide a taxonomy of the main generative models published recently.


title: Towards a population synthesis of discs and planets. II. Confronting disc models and observations at the population level

abstract: Aims. We want to find the distribution of initial conditions that best reproduces disc observations at the population level. Methods. We first ran a parameter study using a 1D model that includes the viscous evolution of a gas disc, dust, and pebbles, coupled with an emission model to compute the millimetre flux observable with ALMA. This was used to train a machine learning surrogate model that can compute the relevant quantity for comparison with observations in seconds. This surrogate model was used to perform parameter studies and synthetic disc populations. Results. Performing a parameter study, we find that internal photoevaporation leads to a lower dependency of disc lifetime on stellar mass than external photoevaporation. This dependence should be investigated in the future. Performing population synthesis, we find that under the combined losses of internal and external photoevaporation, discs are too short lived. Conclusions. To match observational constraints, future models of disc evolution need to include one or a combination of the following processes: infall of material to replenish the discs, shielding of the disc from internal photoevaporation due to magnetically driven disc winds, and extinction of external high-energy radiation. Nevertheless, disc properties in low-external-photoevaporation regions can be reproduced by having more massive and compact discs. Here, the optimum values of the $\alpha$ viscosity parameter lie between $3\times10^{-4}$ and $10^{-3}$ and with internal photoevaporation being the main mode of disc dispersal.


title: Quantum spin ice in three-dimensional Rydberg atom arrays

abstract: Quantum spin liquids are exotic phases of matter whose low-energy physics is described as the deconfined phase of an emergent gauge theory. With recent theory proposals and an experiment showing preliminary signs of $\mathbb{Z}_2$ topological order [G. Semeghini et al., Science 374, 1242 (2021)], Rydberg atom arrays have emerged as a promising platform to realize a quantum spin liquid. In this work, we propose a way to realize a $U(1)$ quantum spin liquid in three spatial dimensions, described by the deconfined phase of $U(1)$ gauge theory in a pyrochlore lattice Rydberg atom array. We study the ground state phase diagram of the proposed Rydberg system as a function of experimentally relevant parameters. Within our calculation, we find that by tuning the Rabi frequency, one can access both the confinement-deconfinement transition driven by a proliferation of "magnetic" monopoles and the Higgs transition driven by a proliferation of "electric" charges of the emergent gauge theory. We suggest experimental probes for distinguishing the deconfined phase from ordered phases. This work serves as a proposal to access a confinement-deconfinement transition in three spatial dimensions on a Rydberg-based quantum simulator.

categories: cond-mat.quant-gas cond-mat.str-el hep-lat physics.atom-ph quant-ph


title: Many-body resonances in the avalanche instability of many-body localization

abstract: Many-body localized (MBL) systems fail to reach thermal equilibrium under their own dynamics, even though they are interacting, nonintegrable, and in an extensively excited state. One instability towards thermalization of MBL systems is the so-called ``avalanche'', where a locally thermalizing rare region is able to spread thermalization through the full system. The spreading of the avalanche may be modeled and numerically studied in finite one-dimensional MBL systems by weakly coupling an infinite-temperature bath to one end of the system. We find that the avalanche spreads primarily via strong many-body resonances between rare near-resonant eigenstates of the closed system. Thus we find and explore a detailed connection between many-body resonances and avalanches in MBL systems.


title: The JWST Resolved Stellar Populations Early Release Science Program II. Survey Overview

abstract: We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come.


title: Anomalies, Representations, and Self-Supervision

abstract: We develop a self-supervised method for density-based anomaly detection using contrastive learning, and test it using event-level anomaly data from CMS ADC2021. The AnomalyCLR technique is data-driven and uses augmentations of the background data to mimic non-Standard-Model events in a model-agnostic way. It uses a permutation-invariant Transformer Encoder architecture to map the objects measured in a collider event to the representation space, where the data augmentations define a representation space which is sensitive to potential anomalous features. An AutoEncoder trained on background representations then computes anomaly scores for a variety of signals in the representation space. With AnomalyCLR we find significant improvements on performance metrics for all signals when compared to the raw data baseline.


title: Kondo Phase in Twisted Bilayer Graphene -- A Unified Theory for Distinct Experiments

abstract: A number of interesting physical phenomena have been discovered in magic-angle twisted bilayer graphene (MATBG), such as superconductivity, correlated gapped and gapless phases, etc. The gapped phases are believed to be symmetry-breaking states described by mean-field theories, whereas gapless phases exhibit features beyond mean field. This work, combining poor man's scaling, numerical renormalization group, and dynamic mean-field theory, demonstrates that the gapless phases are the heavy Fermi liquid state with some symmetries broken and the others preserved. We adopt the recently proposed topological heavy fermion model for MATBG with effective local orbitals around AA-stacking regions and Dirac fermions surrounding them. At zero temperature and most non-integer fillings, the ground states are found to be heavy Fermi liquids and exhibit Kondo resonance peaks. The Kondo temperature $T_K$ is found at the order of 1meV. A higher temperature than $T_K$ will drive the system into a metallic LM phase where disordered LM's and a Fermi liquid coexist. At integer fillings $\pm1,\pm2$, $T_K$ is suppressed to zero or a value weaker than RKKY interaction, leading to Mott insulators or symmetry-breaking states. This theory offers a unified explanation for several experimental observations, such as zero-energy peaks and quantum-dot-like behaviors in STM, the Pomeranchuk effect, and the saw-tooth feature of inverse compressibility, etc. For future experimental verification, we predict that the Fermi surface in the gapless phase will shrink upon heating - as a characteristic of the heavy Fermi liquid. We also conjecture that the heavy Fermi liquid is the parent state of the observed unconventional superconductivity because the Kondo screening reduces the overwhelming Coulomb interaction (~60meV) to a rather small effective interaction (~1meV) comparable to possible weak attractive interactions.


title: SN 2020bio: A Double-peaked, H-poor Type IIb Supernova with Evidence of Circumstellar Interaction

abstract: We present photometric and spectroscopic observations of SN 2020bio, a double-peaked Type IIb supernova (SN) discovered within a day of explosion, primarily obtained by Las Cumbres Observatory and Swift. SN 2020bio displays a rapid and long-lasting initial decline throughout the first week of its light curve, similarly to other well-studied Type IIb SNe. This early-time emission is thought to originate from the cooling of the extended outer hydrogen-rich (H-rich) envelope of the progenitor star that is shock heated by the SN explosion. We compare SN 2020bio to a sample of other double-peaked Type IIb SNe in order to investigate its progenitor properties. Analytical model fits to the early-time emission give progenitor radius ($\approx$ 100--1500 $R_\odot$) and H-rich envelope mass ($\approx$ 0.01--0.5 $M_\odot$) estimates that are consistent with other Type IIb SNe. However, SN 2020bio displays several peculiarities, including: (1) weak H spectral features indicating a greater amount of mass loss than other Type IIb progenitors; (2) an underluminous secondary light-curve peak that implies a small amount of synthesized $^{56}$Ni ($M_{\text{Ni}}$ $\approx$ 0.02 $M_\odot$); and (3) low-luminosity nebular [O I] and interaction-powered nebular features. These observations are more consistent with a lower-mass progenitor ($M_{\text{ZAMS}} \approx$ 12 $M_\odot$) that was stripped of most of its H-rich envelope before exploding. This study adds to the growing diversity in the observed properties of Type IIb SNe and their progenitors.


title: The OGLE Collection of Variable Stars. Over 2600 Delta Scuti Stars in the Small Magellanic Cloud

abstract: We present the first-ever collection of delta Scuti stars found over the entire area of the Small Magellanic Cloud (SMC). The sample consists of 2810 variables of which over 2600 objects belong to the SMC while the remaining stars are most likely members of the Milky Way's halo. The sample has been divided into 2733 singlemode and 77 multimode pulsators. We provide observational parameters (pulsation periods, mean magnitudes, amplitudes, Fourier coefficients) of all delta Sct stars and the long-term I- and V-band time-series photometric measurements collected during the fourth phase of the Optical Gravitational Lensing Experiment (OGLE-IV).


title: Stellar mass dependence of galaxy size-dark matter halo radius relation probed by Subaru-HSC survey weak lensing measurements

abstract: We investigate the stellar mass-dependence of the galaxy size-dark matter halo radius relation for low redshift galaxies using weak gravitational lensing measurements. Our sample consists of $\sim$38,000 galaxies more massive than $10^{8}{\rm M}_{\odot}h^{-2}$ and within $z<0.3$ drawn from the overlap of GAMA survey DR4 and HSC-SSP PDR2. We divide our sample into a number of stellar mass bins and measure stacked weak lensing signals. We model the signals using a conditional stellar mass function to infer the stellar mass-halo mass relation. We fit a single S\'ersic model to HSC $i$-band images of our galaxies and obtain their three-dimensional half-light radii. We use these measurements to construct a median galaxy size-mass relation. We then combine these relations to infer the galaxy size-halo radius relation. We confirm that this relation appears linear given the statistical errors, i.e. the ratio of galaxy size to halo radius remains constant over two orders of magnitudes in stellar mass above $\sim 10^{9} {\rm M}_{\odot}h^{-2}$. Extrapolating the stellar mass-halo mass relation below this limit, we see an indication of a decreasing galaxy size-halo radius ratio with the decline in stellar mass. At stellar mass $\sim 10^{8} {\rm M}_{\odot}h^{-2}$ the ratio becomes 30% smaller than its value in linear regime. The possible existence of a such trend in dwarf galaxy sectors calls for either modification in models employing a constant fraction of halo angular momentum transferred to explain sizes of dwarfs or else points towards our lack of knowledge about dark matter haloes of low-mass galaxies.


title: A DZ white dwarf with a 30 MG magnetic field

abstract: Magnetic white dwarfs with field strengths below 10 MG are easy to recognise since the Zeeman splitting of spectral lines appears proportional to the magnetic field strength. For fields $\geq 100$ MG, however, transition wavelengths become chaotic, requiring quantum-chemical predictions of wavelengths and oscillator strengths with a non-perturbative treatment of the magnetic field. While highly accurate calculations have previously been performed for hydrogen and helium, the variational techniques employed become computationally intractable for systems with more than three to four electrons. Modern computational techniques, such as finite-field coupled-cluster theory, allow the calculation of many-electron systems in arbitrarily strong magnetic fields. Because around 25 percent of white dwarfs have metal lines in their spectra, and some of those are also magnetic, the possibility arises for some metals to be observed in very strong magnetic fields, resulting in unrecognisable spectra. We have identified SDSSJ114333.48+661531.83 as a magnetic DZ white dwarf, with a spectrum exhibiting many unusually shaped lines at unknown wavelengths. Using atomic data calculated from computational finite-field coupled-cluster methods, we have identified some of these lines arising from Na, Mg, and Ca. Surprisingly, we find a relatively low field strength of 30 MG, where the large number of overlapping lines from different elements make the spectrum challenging to interpret at a much lower field strength than for DAs and DBs. Finally we model the field structure of SDSSJ1143+6615 finding the data are consistent with an offset dipole.


title: First Detection of Radio Emission Associated with a Classical Cepheid

abstract: We report the detection of 15 GHz radio continuum emission associated with the classical Cepheid variable star delta Cephei based on observations with the Karl G. Jansky Very Large Array. Our results constitute the first probable detection of radio continuum emission from a classical Cepheid. We observed the star at pulsation phase phi~0.43 (corresponding to the phase of maximum radius and minimum temperature) during three pulsation cycles in late 2018 and detected statistically significant emission (>5 sigma) during one of the three epochs. The observed radio emission appears to be variable at a >~10% level on timescales of days to weeks. We also present an upper limit on the 10 GHz flux density at pulsation phase phi=0.31 from an observation in 2014. We discuss possible mechanisms that may produce the observed 15 GHz emission, but cannot make a conclusive identification from the present data. The emission does not appear to be consistent with originating from a close-in, late-type dwarf companion, although this scenario cannot yet be strictly excluded. Previous X-ray observations have shown that delta Cephei undergoes periodic increases in X-ray flux during pulsation phase phi~0.43. The lack of radio detection in two out of three observing epochs at phi~0.43 suggests that either the radio emission is not linked with a particular pulsation phase, or else that the strength of the generated radio emission in each pulsation cycle is variable.


title: Unravelling the mass spectrum of destroyed dwarf galaxies with the metallicity distribution function

abstract: Accreted stellar populations are comprised of the remnants of destroyed galaxies, and often dominate the `stellar haloes' of galaxies such as the Milky Way (MW). This ensemble of external contributors is a key indicator of the past assembly history of a galaxy. We introduce a novel statistical method that uses the unbinned metallicity distribution function (MDF) of a stellar population to estimate the mass spectrum of its progenitors. Our model makes use of the well-known mass-metallicity relation of galaxies and assumes Gaussian MDF distributions for individual progenitors: the overall MDF is thus a mixture of MDFs from smaller galaxies. We apply the method to the stellar halo of the MW, as well as the classical MW satellite galaxies. The stellar components of the satellite galaxies have relatively small sample sizes, but we do not find any evidence for accreted populations with L > L_host/100. We find that the MW stellar halo has N~1-3 massive progenitors (L > 10^8 L_Sun) within 10 kpc, and likely several hundred progenitors in total. We also test our method on simulations of MW-mass haloes, and find that our method is able to recover the true accreted population within a factor of two. Future datasets will provide MDFs with orders of magnitude more stars, and this method could be a powerful technique to quantify the accreted populations down to the ultra-faint dwarf mass-scale for both the MW and its satellites.


title: Trapped Ion Quantum Computing using Optical Tweezers and the Magnus Effect

abstract: We consider the implementation of quantum logic gates in trapped ions using tightly focused optical tweezers. Strong polarization gradients near the tweezer focus lead to qubit-state dependent forces on the ion. We show that these may be used to implement quantum logic gates on pairs of ion qubits in a crystal. The qubit-state dependent forces generated by this effect live on the plane perpendicular to the direction of propagation of the laser beams opening new ways of coupling to motional modes of an ion crystal. The proposed gate does not require ground state cooling of the ions and does not rely on the Lamb-Dicke approximation, although the waist of the tightly focused beam needs to be comparable with its wavelength in order to achieve the needed field curvature. Furthermore, the gate can be performed on both ground state and magnetic field insensitive clock state qubits without the need for counter-propagating laser fields. This simplifies the setup and eliminates errors due to phase instabilities between the gate laser beams. Finally, we show that imperfections in the gate execution, in particular pointing errors $<30$ nm in the tweezers reduce the gate fidelity from $\mathcal F\gtrsim 0.99998$ to $\gtrsim 0.999$.


title: Kondo Lattice Model of Magic-Angle Twisted-Bilayer Graphene: Hund's Rule, Local-Moment Fluctuations, and Low-Energy Effective Theory

abstract: We apply a generalized Schrieffer-Wolff transformation to the extended Anderson-like topological heavy fermion (THF) model for the magic-angle $(\theta=1.05^{\circ})$ twisted bilayer graphene (MATBLG) (Phys. Rev. Lett. 129, 047601 (2022)), to obtain its Kondo Lattice limit. In this limit localized $f$-electrons on a triangular lattice interact with topological conduction $c$-electrons. By solving the exact limit of the THF model, we show that the integer fillings $\nu=0, \pm 1, \pm 2$ are controlled by the heavy $f$-electrons, while $\nu = \pm 3$ is at the border of a phase transition between two $f$-electron fillings. For $\nu=0, \pm 1, \pm 2$, we then calculate the RKKY interactions between the $f$-moments in the full model and analytically prove the SU(4) Hund's rule for the ground state which maintains that two $f$-electrons fill the same valley-spin flavor. Our (ferromagnetic interactions in the) spin model dramatically differ from the usual Heisenberg antiferromagnetic interactions expected at strong coupling. We show the ground state in some limits can be found exactly by employing a positive semidefinite "bond-operators" method. We then compute the excitation spectrum of the $f$-moments in the ordered ground state, prove the stability of the ground state favored by RKKY interactions, and discuss the properties of the Goldstone modes, the (reason for the accidental) degeneracy of (some of) the excitation modes, and the physics of their phase stiffness. We develop a low-energy effective theory for the $f$-moments and obtain analytic expressions for the dispersion of the collective modes. We discuss the relevance of our results to the spin-entropy experiments in twisted bilayer graphene.


title: The effect of stellar encounters on the dark matter annihilation signal from prompt cusps

abstract: Prompt cusps are the densest quasi-equilibrium dark matter objects; one forms at the instant of collapse within every isolated peak of the initial cosmological density field. They have power-law density profiles, $\rho \propto r^{-1.5}$ with central phase-space density set by the primordial velocity dispersion of the dark matter. At late times they account for $\sim 1\%$ of the dark matter mass but for $>90\%$ of its annihilation luminosity in all but the densest regions, where they are tidally disrupted. Here we demonstrate that individual stellar encounters, rather than the mean galactic tide, are the dominant disruptors of prompt cusps within galaxies. Their cumulative effect is fully (though stochastically) characterised by an impulsive shock strength $B_* = 2\pi G\int\rho_*({\bf x}(t))\, \mathrm{d}t$ where $\rho_*$, the total mass density in stars, is integrated over a cusp's entire post-formation trajectory. Stellar encounters and mean tides have only a small effect on the halo annihilation luminosity seen by distant observers, but this is not true for the Galactic halo because of the Sun's position. For a 100 GeV WIMP, Earth-mass prompt cusps are predicted, and stellar encounters suppress their mean annihilation luminosity by a factor of two already at 20 kpc, so that their annihilation emission is predicted to appear almost uniform over the sky. The Galactic Center $\gamma$-ray Excess is thus unaffected by cusps. If it is indeed dark matter annihilation radiation, then prompt cusps in the outer Galactic halo and beyond must account for 20-80% of the observed isotropic $\gamma$-ray background in the 1 to 10 GeV range.


title: Circuit Complexity through phase transitions: consequences in quantum state preparation

abstract: In this paper, we analyze the circuit complexity for preparing ground states of quantum many-body systems. In particular, how this complexity grows as the ground state approaches a quantum phase transition. We discuss different definitions of complexity, namely the one following the Fubini-Study metric or the Nielsen complexity. We also explore different models: Ising, ZZXZ or Dicke. In addition, different forms of state preparation are investigated: analytic or exact diagonalization techniques, adiabatic algorithms (with and without shortcuts), and Quantum Variational Eigensolvers. We find that the divergence (or lack thereof) of the complexity near a phase transition depends on the non-local character of the operations used to reach the ground state. For Fubini-Study based complexity, we extract the universal properties and their critical exponents. In practical algorithms, we find that the complexity depends crucially on whether or not the system passes close to a quantum critical point when preparing the state. For both VQE and Adiabatic algorithms, we provide explicit expressions and bound the growth of complexity with respect to the system size and the execution time, respectively.


title: Clues on the presence and segregation of very massive stars in the Sunburst Lyman-continuum cluster at z=2.37

abstract: We report the identification of very massive stars (VMS; mass $> 100$\,\msun) that may be segregated in the center of the young massive star cluster at $z$=2.37 hosted in the lensed galaxy called {\tt Sunburst} galaxy. This result is based on two pieces of evidence: (1) VLT/MUSE spectra of several multiple images of the same star cluster show key spectral signatures of VMS, such as the \heii\ broad emission, \nivblue\ emission, and an \niv\ P-Cygni profile. In particular, \heii\ is broad ($\sim1610\pm300$ \kms), with an equivalent width of 3\AA,\ and asymmetric profile. These features require an extremely young ($\sim2.5$ Myr) stellar population component in which the masses of the stars exceed 100~\msun. When a Salpeter initial mass function and BPASS models for normal massive stars are assumed, the observed spectral features require $\sim$400 VMS. (2) The same star cluster is detected at a signal-to-noise ratio of~$\sim100$ in the Lyman continuum domain ($\lambda < 900$\AA). The Lyman continuum emission emerges from a region with a radius that is at least twice smaller than what is observed at 1700\AA~(independently of magnification) and is located in the center of the cluster. After delensing, the effective radii in absolute scales are R$_{\tt eff}[{\tt LyC}]\sim4.7 \pm 1.5$ pc and R$_{\tt eff}[1700]= 7.8 \pm 1.4$ pc. The Lyman continuum radiation is mainly produced by hot and massive stars, which implies that their spatial distribution (including that of VMS) is preferentially more confined in the central parts of the cluster. Approximately 400 VMS hosted by a cluster of $\sim 10^7$ \msun\ produce $\sim$15\% of the escaping Lyman continuum photons, and the remaining photons are produced by other massive early-type stars.


title: Symmetric Kondo Lattice States in Doped Strained Twisted Bilayer Graphene

abstract: We use the topological heavy fermion (THF) model and its Kondo Lattice (KL) formulation to study the symmetric Kondo state in twisted bilayer graphene. Via a large-N approximation, we find a symmetric Kondo (SK) state in KL mode at fillings $\nu=0,\pm 1,\pm 2$. In the SK state, all symmetries are preserved and the local moments are Kondo screened by the conduction electrons. At the mean-field level of the THF model at $\nu=0,\pm 1, \pm 2, \pm 3$, we also find a similar symmetric state. We study the stability of the symmetric state by comparing its energy with the ordered states and find the ordered states to have lower energy. However, moving away from integer fillings by doping holes to the light bands, we find the energy difference is reduced, which suggests the loss of ordering and a tendency towards Kondo screening. In order to include many-body effects beyond the mean-field approximation, we perform dynamical mean-field theory (DMFT) calculations on the THF model. We find the spin susceptibility follows a Curie behavior at $\nu=0, \pm 1,\pm 2$ down to $\sim 2\text{K}$ where the onset of screening of the local moment becomes visible. This hints to very low Kondo temperatures at these fillings, in agreement with the outcome of our mean-field calculations. At non-integer filling $\nu=\pm 0.5,\pm 0.8,\pm 1.2$ DMFT shows deviations from a $1/T$-susceptibility at much higher temperatures, suggesting a more effective screening of local moments with doping. Finally, we study the effect of a $C_{3z}$-rotational-symmetry-breaking strain via mean-field approaches and find that a symmetric phase (that only breaks $C_{3z}$ symmetry) can be stabilized at sufficiently large strain at $\nu=0,\pm 1, \pm 2$. Our results suggest that a symmetric Kondo phase is strongly suppressed at integer fillings, but could be stabilized either at non-integer fillings or by applying strain.


title: Late vacuum choice and slow roll approximation in gravitational particle production during reheating

abstract: In the transition between inflation and reheating, the curvature scalar typically undergoes oscillations which have significant impact on the density of gravitationally produced particles. The commonly used adiabatic vacuum prescription for the extraction of produced particle spectra becomes a non-reliable definition of vacuum in the regimes for which this oscillatory behavior is important. In this work, we study particle production for a scalar field non-minimally coupled to gravity, taking into account the complete dynamics of spacetime during inflation and reheating. We derive an approximation for the solution to the mode equation during the slow-roll of the inflaton and analyze the importance of Ricci scalar oscillations in the resulting spectra. Additionally, we propose a prescription for the vacuum that allows to safely extrapolate the result to the present, given that the test field interacts only gravitationally. Lastly, we calculate the abundance of dark matter this mechanism yields and compare it to observations.


title: Systematic design of a robust half-W1 photonic crystal waveguide for interfacing slow light and trapped cold atoms

abstract: Novel platforms interfacing trapped cold atoms and guided light in nanoscale waveguides are a promising route to achieve a regime of strong coupling between light and atoms in single pass, with applications to quantum non-linear optics and quantum simulation. A strong challenge for the experimental development of this emerging waveguide-QED field of research is to combine facilitated optical access for atom transport, atom trapping via guided modes and robustness to inherent nanofabrication imperfections. In this endeavor, here we propose to interface Rubidium atoms with a photonic-crystal waveguide based on a large-index GaInP slab. With a specifically tailored half-W1 design, we show that a large chiral coupling to the waveguide can be obtained and guided modes can be used to form two-color dipole traps for atoms at 116~nm from the edge of the structure. This optimized device should greatly improve the level of experimental control and facilitate the atom integration.


title: A viable $A_4$ 3HDM theory of quark mass matrices

abstract: It is known that a three Higgs doublet model (3HDM)symmetric under an exact $A_4$ symmetry is not compatible with nonzero quark masses and/or non-block-diagonal CKM matrix. We show that a 3HDM with softly broken $A_4$ terms in the scalar potential does allow for a fit of quark mass matrices. Moreover, the result is consistent with $m_h=125\textrm{GeV}$ and the $h \rightarrow WW, ZZ$ signal. We also checked numerically that, for each point that passes all the constraints, the minimum is a global minimum of the potential.


title: Path integrals for classical-quantum dynamics

abstract: Consistent dynamics which couples classical and quantum degrees of freedom exists. This dynamics is linear in the hybrid state, completely positive and trace preserving. Starting from completely positive classical-quantum master equations, we derive a general path integral representation for such dynamics in terms of a classical-quantum action, which includes the necessary and sufficient conditions for complete positivity and trace preservation. The path integral we study is a generalization of the Feynman path integral for quantum systems, and the stochastic path integral used to study classical stochastic processes, allowing for interaction between the classical and quantum systems. When the classical-quantum Hamiltonian is at most quadratic in the momenta we are able to derive a configuration space path integral, providing a map between master equations and covariant classical-quantum path integrals.


title: Representation Stability for Disks in a Strip

abstract: We consider the ordered configuration space of $n$ open unit-diameter disks in the infinite strip of width $w$. In the spirit of Arnol'd and Cohen, we provide a finite presentation for the rational homology groups of this ordered configuration space as a twisted algebra. We use this presentation to prove that the ordered configuration space of open unit-diameter disks in the infinite strip of width $w$ exhibits a notion of first-order representation stability similar to Church--Ellenberg--Farb and Miller--Wilson's first-order representation stability for the ordered configuration space of points in a manifold. In addition, we prove that for large $w$ this disk configuration space exhibits notions of second- (and higher) order representation stability.


title: Coupling multi-fluid dynamics equipped with Landau closures to the particle-in-cell method

abstract: The particle-in-cell (PIC) method is successfully used to study magnetized plasmas. However, this requires large computational costs and limits simulations to short physical run-times and often to setups in less than three spatial dimensions. Traditionally, this is circumvented either via hybrid-PIC methods (adopting massless electrons) or via magneto-hydrodynamic-PIC methods (modelling the background plasma as a single charge-neutral magneto-hydrodynamical fluid). Because both methods preclude modelling important plasma-kinetic effects, we introduce a new fluid-PIC code that couples a fully explicit and charge-conservative multi-fluid solver to the PIC code SHARP through a current-coupling scheme and solve the full set of Maxwell's equations. This avoids simplifications typically adopted for Ohm's Law and enables us to fully resolve the electron temporal and spatial scales while retaining the versatility of initializing any number of ion, electron, or neutral species with arbitrary velocity distributions. The fluid solver includes closures emulating Landau damping so that we can account for this important kinetic process in our fluid species. Our fluid-PIC code is second-order accurate in space and time. The code is successfully validated against several test problems, including the stability and accuracy of shocks and the dispersion relation and damping rates of waves in unmagnetized and magnetized plasmas. It also matches growth rates and saturation levels of the gyro-scale and intermediate-scale instabilities driven by drifting charged particles in magnetized thermal background plasmas in comparison to linear theory and PIC simulations. This new fluid-SHARP code is specially designed for studying high-energy cosmic rays interacting with thermal plasmas over macroscopic timescales.

categories: astro-ph.HE astro-ph.GA physics.comp-ph physics.flu-dyn physics.plasm-ph


title: Formation of Rocky Super-Earths From A Narrow Ring of Planetesimals

abstract: The formation of super-Earths, the most abundant planets in the Galaxy, remains elusive. These planets have masses that typically exceed that of the Earth by a factor of a few; appear to be predominantly rocky, although often surrounded by H/He atmospheres; and frequently occur in multiples. Moreover, planets that encircle the same star tend to have similar masses and radii, whereas those belonging to different systems exhibit remarkable overall diversity. Here, we advance a theoretical picture for rocky planet formation that satisfies the aforementioned constraints: building upon recent work - which demonstrates that planetesimals can form rapidly at discrete locations in the disk - we propose that super-Earths originate inside rings of silicate-rich planetesimals at approximately ~1 AU. Within the context of this picture, we show that planets grow primarily through pairwise collisions among rocky planetesimals, until they achieve terminal masses that are regulated by isolation and orbital migration. We quantify our model with numerical simulations and demonstrate that our synthetic planetary systems bear a close resemblance to compact, multi-resonant progenitors of the observed population of short-period extrasolar planets. Our results thus indicate that the absence of short-period super-Earths within the solar system can simply be attributed to the comparatively low mass of the primordial planetesimal ring within the protosolar nebula.


title: Stability of standard Einstein metrics on homogeneous spaces of non-simple Lie groups

abstract: The classification of compact homogeneous spaces of the form $M=G/K$, where $G$ is a non-simple Lie group, such that the standard metric is Einstein is still open. The only known examples are $4$ infinite families and $3$ isolated spaces found by Nikonorov and Rodionov in the 90s. In this paper, we prove that most of these standard Einstein metrics are unstable as critical points of the scalar curvature functional on the manifold of all unit volume $G$-invariant metrics on $M$, providing a lower bound for the coindex in the case of Ledger-Obata spaces. On the other hand, examples of stable (in particular, local maxima) invariant Einstein metrics on certain homogeneous spaces of non-simple Lie groups are also given.


title: Joint k-TE Space Image Reconstruction and Data Fitting for T2 Mapping

abstract: Objectives: To develop a joint k-TE reconstruction algorithm to reconstruct the T2-weighted (T2W) images and T2 map simultaneously. Materials and Methods: The joint k-TE reconstruction model was formulated as an optimization problem subject to a self-consistency condition of the exponential decay relationship between the T2W images and T2 map. The objective function included a data fidelity term enforcing the agreement between the solution and the measured k-space data, together with a spatial regularization term on image properties of the T2W images. The optimization problem was solved using Alternating-Direction Method of Multipliers (ADMM). We tested the joint k-TE method in phantom data and healthy volunteer scans with fully-sampled and under-sampled k-space lines. Image quality of the reconstructed T2W images and T2 map, and the accuracy of T2 measurements derived by the joint k- TE and the conventional signal fitting method were compared. Results: The proposed method improved image quality with reduced noise and less artifacts on both T2W images and T2 map, and increased measurement consistency in T2 relaxation time measurements compared with the conventional method in all data sets. Conclusions: The proposed reconstruction method outperformed the conventional magnitude image-based signal fitting method in image quality and stability of quantitative T2 measurements


title: On Helly numbers of exponential lattices

abstract: Given a set $S \subseteq \mathbb{R}^2$, define the \emph{Helly number of $S$}, denoted by $H(S)$, as the smallest positive integer $N$, if it exists, for which the following statement is true: for any finite family $\mathcal{F}$ of convex sets in~$\mathbb{R}^2$ such that the intersection of any $N$ or fewer members of~$\mathcal{F}$ contains at least one point of $S$, there is a point of $S$ common to all members of $\mathcal{F}$. We prove that the Helly numbers of \emph{exponential lattices} $\{\alpha^n \colon n \in \mathbb{N}_0\}^2$ are finite for every $\alpha>1$ and we determine their exact values in some instances. In particular, we obtain $H(\{2^n \colon n \in \mathbb{N}_0\}^2)=5$, solving a problem posed by Dillon (2021). For real numbers $\alpha, \beta > 1$, we also fully characterize exponential lattices $L(\alpha,\beta) = \{\alpha^n \colon n \in \mathbb{N}_0\} \times \{\beta^n \colon n \in \mathbb{N}_0\}$ with finite Helly numbers by showing that $H(L(\alpha,\beta))$ is finite if and only if $\log_\alpha(\beta)$ is rational.


title: Design and Characterization of Viscoelastic McKibben Actuators with Tunable Force-Velocity Curves

abstract: The McKibben pneumatic artificial muscle is a commonly studied soft robotic actuator, and its quasistatic force-length properties have been well characterized and modeled. However, its damping and force-velocity properties are less well studied. Understanding these properties will allow for more robust dynamic modeling of soft robotic systems. The force-velocity response of these actuators is of particular interest because these actuators are often used as hardware models of skeletal muscles for bioinspired robots, and this force-velocity relationship is fundamental to muscle physiology. In this work, we investigated the force-velocity response of McKibben actuators and the ability to tune this response through the use of viscoelastic polymer sheaths. These viscoelastic McKibben actuators (VMAs) were characterized using iso-velocity experiments inspired by skeletal muscle physiology tests. A simplified 1D model of the actuators was developed to connect the shape of the force-velocity curve to the material parameters of the actuator and sheaths. Using these viscoelastic materials, we were able to modulate the shape and magnitude of the actuators' force-velocity curves, and using the developed model, these changes were connected back to the material properties of the sheaths.


title: SHUNIT: Style Harmonization for Unpaired Image-to-Image Translation

abstract: We propose a novel solution for unpaired image-to-image (I2I) translation. To translate complex images with a wide range of objects to a different domain, recent approaches often use the object annotations to perform per-class source-to-target style mapping. However, there remains a point for us to exploit in the I2I. An object in each class consists of multiple components, and all the sub-object components have different characteristics. For example, a car in CAR class consists of a car body, tires, windows and head and tail lamps, etc., and they should be handled separately for realistic I2I translation. The simplest solution to the problem will be to use more detailed annotations with sub-object component annotations than the simple object annotations, but it is not possible. The key idea of this paper is to bypass the sub-object component annotations by leveraging the original style of the input image because the original style will include the information about the characteristics of the sub-object components. Specifically, for each pixel, we use not only the per-class style gap between the source and target domains but also the pixel's original style to determine the target style of a pixel. To this end, we present Style Harmonization for unpaired I2I translation (SHUNIT). Our SHUNIT generates a new style by harmonizing the target domain style retrieved from a class memory and an original source image style. Instead of direct source-to-target style mapping, we aim for source and target styles harmonization. We validate our method with extensive experiments and achieve state-of-the-art performance on the latest benchmark sets. The source code is available online: https://github.com/bluejangbaljang/SHUNIT.


title: Non-Lorentzian Ka\v{c}-Moody Algebras

abstract: We investigate two dimensional (2d) quantum field theories which exhibit Non- Lorentzian Ka\v{c}-Moody (NLKM) algebras as their underlying symmetry. Our investigations encompass both 2d Galilean (speed of light $c \rightarrow \infty$) and Carrollian ($c \rightarrow 0$) CFTs with additional number of infinite non-Abelian currents, stemming from an isomorphism between the two algebras. We alternate between an intrinsic and a limiting analysis. Our NLKM algebra is constructed first through a contraction and then derived from an intrinsically Carrollian perspective. We then go on to use the symmetries to derive a Non-Lorentzian (NL) Sugawara construction and ultimately write down the NL equivalent of the Knizhnik Zamolodchikov equations. All of these are also derived from contractions, thus providing a robust cross-check of our analyses.