title: Non-equilibrium Photoionization and Hydrodynamic Simulations of Starburst-driven Outflows

abstract: Starburst-driven galactic outflows in star-forming galaxies have been observed to contain complex thermal structures and emission line features that are difficult to explain by adiabatic fluid models and plasmas in photoionization equilibrium (PIE) and collisional ionization equilibrium (CIE). We previously performed hydrodynamic simulations of starburst-driven outflows using the MAIHEM module for non-equilibrium atomic chemistry and radiative cooling functions in the hydrodynamics code FLASH, and calculated emission lines in combined CIE and PIE conditions. In the present study, we consider time-dependent non-equilibrium ionization (NEI) states produced by our MAIHEM simulations. Through extensive CLOUDY calculations made with the NEI states from our hydrodynamic simulations, we predict the UV and optical emission line profiles for starburst-driven outflows in time-evolving non-equilibrium photoionization conditions. Our hydrodynamic results demonstrate applications of non-equilibrium radiative cooling for H II regions in starburst galaxies hosting cool outflows.


title: Structure of IRAS 05168+3634 star-forming region

abstract: This study aims to determine the main physical parameters (N(H2) hydrogen column density and Td dust temperature) of the Interstellar medium, and their distribution in the extended star-forming region, which includes IRAS 05156+3643, 05162+3639, 05168+3634, 05177+3636, and 05184+3635 sources. We also provide a comparative analysis of the properties of the Interstellar medium and young stellar objects. Analysis of the results revealed that Interstellar medium forms relativity dense condensations around IRAS sources, which are interconnected by a filament structure. In general, in sub-regions Td varies from 11 to 24 K, and N(H2) - from 1.0 to 4.0 x 10^23 cm^(-2). The masses of the ISM vary from 1.7 x 10^4 to 2.1 x 10^5 Msol. All BGPSv2 objects identified in this star-forming region are located at the N(H2) maximum. The direction of the outflows, which were found in two sub-regions, IRAS 05168+3634 and 05184+3635, correlates well with the isodenses direction. The sub-regions with the highest N(H2) and Interstellar medium mass have the largest percentage of young stellar objects with Class I evolutionary stage. The wide spread of the evolutionary ages of stars in all sub-regions (from 10^5 to 10^7 years) suggests that the process of star formation in the considered region is sequential. In those sub-regions where the mass of the initial, parent molecular cloud is larger, this process is likely to proceed more actively. On the Gaia EDR3 database, it can be assumed that all sub-regions are embedded in the single molecular cloud and belong to the same star-forming region, which is located at a distance of about 1.9 kpc.


title: Model-Free Risk-Sensitive Reinforcement Learning

abstract: We extend temporal-difference (TD) learning in order to obtain risk-sensitive, model-free reinforcement learning algorithms. This extension can be regarded as modification of the Rescorla-Wagner rule, where the (sigmoidal) stimulus is taken to be either the event of over- or underestimating the TD target. As a result, one obtains a stochastic approximation rule for estimating the free energy from i.i.d. samples generated by a Gaussian distribution with unknown mean and variance. Since the Gaussian free energy is known to be a certainty-equivalent sensitive to the mean and the variance, the learning rule has applications in risk-sensitive decision-making.


title: Constraining the nucleon size with relativistic nuclear collisions

abstract: The notion of the "size" of nucleons and their constituents plays a pivotal role in the current paradigm of the formation and the fluctuations of the quark-gluon plasma produced in high-energy nuclear collision experiments. We report on state-of-the-art hydrodynamic results showing that the correlation between anisotropic flow, $v_n^2$, and the mean transverse momentum of hadrons, $[p_t]$, possesses a unique sensitivity to the nucleon size in off-central heavy-ion collisions. We argue that existing experimental measurements of this observable support a picture where the relevant length scale characterizing the colliding nucleons is of order 0.5 fm or smaller, and we discuss the broad implications of this finding for future global Bayesian analyses aimed at extracting initial-state and medium properties from nucleus-nucleus collision data, including $v_n^2$-$[p_t]$ correlations. Determinations of the nucleon size in heavy-ion collisions will provide a solid independent constraint on the initial state of small system collisions, and will establish a deep connection between collective flow data in nucleus-nucleus experiments and data on deep inelastic scattering on protons and nuclei.


title: Cooper pairs localization in tree-like networks of superconducting islands

abstract: We study inhomogeneous Cooper pairs distribution and localization effects in tree-like networks of superconducting islands coupled via Josephson weak links. Using a generalized Feynman's approach, reminiscent of the Bose-Hubbard model, we demonstrate that the Cooper pairs fraction which localizes on a specific network's island is limited by the network topology and, if present, by the repulsive interaction. These findings contribute to clarify the interplay between confinement effects induced by the network's topology and interaction and shed some light on recent experiments dealing with networks of Josephson junctions.


title: Estimating SARS-CoV-2 Seroprevalence

abstract: Governments and public health authorities use seroprevalence studies to guide responses to the COVID-19 pandemic. Seroprevalence surveys estimate the proportion of individuals who have detectable SARS-CoV-2 antibodies. However, serologic assays are prone to misclassification error, and non-probability sampling may induce selection bias. In this paper, nonparametric and parametric seroprevalence estimators are considered that address both challenges by leveraging validation data and assuming equal probabilities of sample inclusion within covariate-defined strata. Both estimators are shown to be consistent and asymptotically normal, and consistent variance estimators are derived. Simulation studies are presented comparing the estimators over a range of scenarios. The methods are used to estimate SARS-CoV-2 seroprevalence in New York City, Belgium, and North Carolina.


title: Emergent Mott insulators and non-Hermitian conservation laws in an interacting bosonic chain with noninteger filling and nonreciprocal hopping

abstract: We investigate the ground state and quantum dynamics of an interacting bosonic chain with the nonreciprocal hopping. In sharp contrast to its Hermitian counterpart, the ground state can support Mott insulators in systems with noninteger filling due to the competition between nonreciprocal hopping and the on-site interaction. For the quantum dynamics, conservation laws for non-Hermitian systems manifest a stark difference from their Hermitian counterpart. In particular, for any Hermitian operator that commutes with the Hamiltonian operator, its expectation value is guaranteed to be nonconserved in the non-Hermitian quantum dynamics. To systematically identify the non-Hermitian conservation law, we establish a generic approach for constructing the conserved quantities in non-Hermitian many-body quantum systems with completely real spectra, and illustrate it concretely by the system under study. The direct experimental observation of Mott insulators in systems with noninteger filling and non-Hermitian conservation laws can be performed by ultracold atoms in optical lattices with the engineered nonreciprocal hopping.

categories: cond-mat.mes-hall cond-mat.quant-gas cond-mat.str-el quant-ph


title: Efficient approximation of Jacobian matrices involving a non-uniform fast Fourier transform (NUFFT)

abstract: There is growing interest in learning Fourier domain sampling strategies (particularly for magnetic resonance imaging, MRI) using optimization approaches. For non-Cartesian sampling patterns, the system models typically involve non-uniform FFT (NUFFT) operations. Commonly used NUFFT algorithms contain frequency domain interpolation, which is not differentiable with respect to the sampling pattern, complicating the use of gradient methods. This paper describes an efficient and accurate approach for computing approximate gradients involving NUFFTs. Multiple numerical experiments validated the improved accuracy and efficiency of the proposed approximation. As an application to computational imaging, the NUFFT Jacobians were used to optimize non-Cartesian MRI sampling trajectories via data-driven stochastic optimization. Specifically, the sampling patterns were learned with respect to various model-based reconstruction algorithms, including quadratic regularized reconstruction and compressed sensing-based reconstruction. The proposed approach enables sampling optimization for image sizes that are infeasible with standard auto-differentiation methods due to memory limits. The synergistic acquisition and reconstruction lead to remarkably improved image quality. In fact, we show that model-based image reconstruction (MBIR) methods with suitably optimized imaging parameters can perform nearly as well as CNN-based methods.


title: Pressure profiles of distant Galaxy clusters with Planck-SPT data

abstract: We present a full set of numerical tools to extract Galaxy Cluster pressure profiles from the joint analysis of Planck and South Pole Telescope (SPT) observations. Pressure profiles are powerful tracers of the thermodynamic properties and the internal structure of the clusters. Tracing the pressure over the cosmic times allows to constraints the evolution of the cluster structure and the contribution of astrophysical phenomena. SPT and Planck are complementary to constrain the cluster structure at various spatial scales. The SPT cluster catalogue counts 677 cluster candidates up to redshift 1.7, it is a nearly mass limited sample, an ideal benchmark to test cluster evolution. We developed a pipeline to first separate the cluster signal from the background and foreground components and then jointly fit a parametric profile model on a combination of Planck and SPT data. We validate our algorithm on a sub-sample of six clusters, common to the SPT and the CHEX-MATE catalogues, comparing the results with the profiles obtained from X-ray observations with XMM-Newton.


title: Binary Formation in the Orion Nebula Cluster: Exploring the Sub-stellar Limit

abstract: We present results constraining the multiplicity of the very low mass stars and sub-stellar objects in the Orion Nebula Cluster (ONC). Our sample covers primary masses 0.012-0.1M$_{\odot}$ using archival Hubble Space Telescope data obtained with the Advanced Camera for Surveys using multiple filters. Studying the binary populations of clusters provides valuable constraints of how the birth environment affects binary formation and evolution. Prior surveys have shown that the binary populations of high-mass, high-density star clusters like the ONC may substantially differ from those in low-mass associations. Very low mass stellar and sub-stellar binaries at wide separations, $>$20AU, are statistically rare in the Galactic field and have been identified in stellar associations like Taurus-Auriga and Ophiuchus. They also may be susceptible to dynamical interactions, and their formation may be suppressed by feedback from on-going star formation. We implement a double point-spread function (PSF) fitting algorithm using empirical, position dependent PSF models to search for binary companions at projected separations $>$ 10 AU (25 mas). With this technique, we identify 7 very low mass binaries, 5 of which are new detections, resulting in a binary frequency of 12$^{+6}_{-3.2}$% over mass ratios of 0.5 - 1.0 and projected separations of 20 - 200 AU. We find an excess of very low mass binaries in the ONC compared to the Galactic field, with a probability of 10$^{-6}$ that the populations are statistically consistent. The sub-stellar population of the ONC may require further dynamical processing of the lowest binding energy binaries to resemble the field population.


title: Local Compatibility Boundary Conditions for High-Order Accurate Finite-Difference Approximations of PDEs

abstract: We describe a new approach to derive numerical approximations of boundary conditions for high-order accurate finite-difference approximations. The approach, called the Local Compatibility Boundary Condition (LCBC) method, uses boundary conditions and compatibility boundary conditions derived from the governing equations, as well as interior and boundary grid values, to construct a local polynomial, whose degree matches the order of accuracy of the interior scheme, centered at each boundary point. The local polynomial is then used to derive a discrete formula for each ghost point in terms of the data. This approach leads to centered approximations that are generally more accurate and stable than one-sided approximations. Moreover, the stencil approximations are local since they do not couple to neighboring ghost-point values which can occur with traditional compatibility conditions. The local polynomial is derived using continuous operators and derivatives which enables the automatic construction of stencil approximations at different orders of accuracy. The LCBC method is developed here for problems governed by second-order partial differential equations, and it is verified for a wide range of sample problems, both time-dependent and time-independent, in two space dimensions and for schemes up to sixth-order accuracy.


title: A Unified View of Relational Deep Learning for Drug Pair Scoring

abstract: In recent years, numerous machine learning models which attempt to solve polypharmacy side effect identification, drug-drug interaction prediction and combination therapy design tasks have been proposed. Here, we present a unified theoretical view of relational machine learning models which can address these tasks. We provide fundamental definitions, compare existing model architectures and discuss performance metrics, datasets and evaluation protocols. In addition, we emphasize possible high impact applications and important future research directions in this domain.


title: Laser Manipulation of Spin-Exchange Interaction Between Alkaline-Earth Atoms in $^1$S$_0$ and $^3$P$_2$ States

abstract: Ultracold gases of fermionic alkaline-earth (like) atoms are hopeful candidates for the quantum simulation of many-body physics induced by magnetic impurities (e.g., the Kondo physics), because there are spin-exchange interactions (SEIs) between two atoms in the electronic ground ($^1$S$_0$) and metastable ($^3$P) state, respectively. Nevertheless, this SEI cannot be tuned via magnetic Feshbach resonance. In this work we propose three methods to control the SEI between one atom in the $^1$S$_0$ state and another atom in the $^3$P$_2$ states or $^3$P$_2$-$^3$P$_0$ dressed states, with one or two laser beams.These methods are based on the spin-dependent AC-Stark shifts of the $^3$P$_2$ states, or the $^3$P$_2$-$^3$P$_0$ Raman coupling. We show that due to the structure of alkaline-earth (like) atoms, the heating effects induced by the laser beams of our methods are very weak. For instance, for ultracold Yb atoms, AC-Stark-shift difference of variant spin states of the $^3$P$_2(F=3/2)$ level, or the strength of the $^3$P$_2$-$^3$P$_0$ Raman coupling, could be of the order of $(2\pi)$MHz, while the heating rate (photon scattering rate) is only of the order of Hz. As a result, the Feshbach resonances, with which one can efficiently control the SEI by changing the laser intensity, may be induced by the laser beams with low-enough heating rate, even if the scattering lengths of the bare inter-atomic interaction are so small that being comparable with the length scale associated with the van der Waals interaction.


title: Metric definition of quasiconformality and exceptional sets

abstract: We show that a homeomorphism of Euclidean space is quasiconformal if and only if at each point there exists a sequence of uncentered open sets with bounded eccentricity shrinking to that point whose images also have bounded eccentricity. This generalizes the metric definition of quasiconformality of Gehring that uses balls instead. We also study exceptional sets for this definition, in connection with sets that are negligible for extremal distances. We introduce the class of CNED sets, generalizing the classical notion of NED sets studied by Ahlfors--Beurling. A set $A$ is CNED if the conformal modulus of a curve family is not affected when one restricts to the subfamily intersecting the set $A$ at countably many points. We show as our main theorem that CNED sets are exceptional for the definition of quasiconformality.


title: Amplitude Structure of Optical Vortices Determines Annihilation Dynamics

abstract: We show that annihilation dynamics between oppositely charged optical vortex pairs can be manipulated by the initial size of the vortex cores, consistent with hydrodynamics. When sufficiently close together, vortices with strongly overlapped cores annihilate more quickly than vortices with smaller cores that must wait for diffraction to cause meaningful core overlap. Numerical simulations and experimental measurements for vortices with hyperbolic tangent cores of various initial sizes show that hydrodynamics governs their motion, and reveal distinct phases of vortex recombination; decreasing the core size of an annihilating pair can prevent the annihilation event.


title: Examples of non--effective rays at the boundary of the Mori cone of blow--ups of the plane

abstract: In this paper we prove that no multiple of the linear system of plane curves of degree $d\geq 4$ with a point of multiplicity $d-m$ (with $2 \leq m \leq d$) and $m(2d-m)$ simple general points is effective.


title: Map-Assisted Constellation Design for mmWave WDM with OAM in Short-Range LOS Environment

abstract: We consider a system that integrates positioning and single-user millimeter wave (mmWave) communication, where the communication part adopts wavelength division multiplexing (WDM) and orbital angular momentum (OAM). This paper addresses the multi-dimensional constellation design in shortrange line-of-sight (LOS) environment, with stable communication links. We propose a map-assisted method to quantify the system parameters based on positions and reduce real-time computing overhead. We explore the possibility of using a few patterns in the maps, and investigate its performance loss. We first investigate the features of OAM beams, and find that the link gain ratio between any two sub-channels remains unchanged at some postions. Then, we prove that a fixed constellation can be adopted for the positions where the link gain matrices are sufficiently close to be proportional. Moreover, we prove that the system can adopt a fixed power vector to generate a multidimensional constellation if the difference between fixed power vector and optimal power vector is small. Finally, we figure out that the constellation design for all receiver locations can be represented by a few constellation sets.


title: Reconstructing Nonlinear Dynamical Systems from Multi-Modal Time Series

abstract: Empirically observed time series in physics, biology, or medicine, are commonly generated by some underlying dynamical system (DS) which is the target of scientific interest. There is an increasing interest to harvest machine learning methods to reconstruct this latent DS in a data-driven, unsupervised way. In many areas of science it is common to sample time series observations from many data modalities simultaneously, e.g. electrophysiological and behavioral time series in a typical neuroscience experiment. However, current machine learning tools for reconstructing DSs usually focus on just one data modality. Here we propose a general framework for multi-modal data integration for the purpose of nonlinear DS reconstruction and the analysis of cross-modal relations. This framework is based on dynamically interpretable recurrent neural networks as general approximators of nonlinear DSs, coupled to sets of modality-specific decoder models from the class of generalized linear models. Both an expectation-maximization and a variational inference algorithm for model training are advanced and compared. We show on nonlinear DS benchmarks that our algorithms can efficiently compensate for too noisy or missing information in one data channel by exploiting other channels, and demonstrate on experimental neuroscience data how the algorithm learns to link different data domains to the underlying dynamics.


title: Non effective planar linear systems at the boundary of the Mori cone

abstract: In this paper we prove that certain linear systems (and all their multiples) of plane curves with general base points and zero--self intersection are empty, thus exhibiting further examples of rays at the boundary of the Mori cone of a general blow--up of the plane.


title: The ABALONE Photosensor

abstract: The ABALONE is a new type of photosensor produced by PhotonLab, Inc. with cost effective mass production, robustness and high performance. This modern technology provides sensitivity to visible and UV light, exceptional radio-purity and excellent detection performance in terms of intrinsic gain, afterpulsing rate, timing resolution and single-photon sensitivity. For these reasons, the ABALONE can have many fields of application, including particle physics experiments, such as DARWIN, and medical imaging. This new hybrid photosensor, that works as light intensifier, is based on the acceleration in vacuum of photoelectrons generated in a traditional photosensor cathode and guided towards a window of scintillating material that can be read from the outside through a silicon photomultiplier. In this work we present the simulation of the ABALONE and the results from operation at room temperature. The goal of the characterization is the evaluation of the gain, the response in time and the single photoelectron spectrum as a function of the electric field and the photoelectron emission angle. Details of future tests will be also discussed.


title: SZ3: A Modular Framework for Composing Prediction-Based Error-Bounded Lossy Compressors

abstract: Today's scientific simulations require a significant reduction of data volume because of extremely large amounts of data they produce and the limited I/O bandwidth and storage space. Error-bounded lossy compressor has been considered one of the most effective solutions to the above problem. In practice, however, the best-fit compression method often needs to be customized/optimized in particular because of diverse characteristics in different datasets and various user requirements on the compression quality and performance. In this paper, we develop a novel modular, composable compression framework (namely SZ3), which involves three significant contributions. (1) SZ3 features a modular abstraction for the prediction-based compression framework such that the new compression modules can be plugged in easily. (2) SZ3 supports multialgorithm predictors and can automatically select the best-fit predictor for each data block based on the designed error estimation criterion. (3) SZ3 allows users to easily compose different compression pipelines on demand, such that both compression quality and performance can be significantly improved for their specific datasets and requirements. (4) In addition, we evaluate several lossy compressors composed from SZ3 using the real-world datasets. Specifically, we leverage SZ3 to improve the compression quality and performance for different use-cases, including GAMESS quantum chemistry dataset and Advanced Photon Source (APS) instrument dataset. Experiments show that our customized compression pipelines lead to up to 20% improvement in compression ratios under the same data distortion compared with the state-of-the-art approaches.


title: OpenFWI: Large-Scale Multi-Structural Benchmark Datasets for Seismic Full Waveform Inversion

abstract: Full waveform inversion (FWI) is widely used in geophysics to reconstruct high-resolution velocity maps from seismic data. The recent success of data-driven FWI methods results in a rapidly increasing demand for open datasets to serve the geophysics community. We present OpenFWI, a collection of large-scale multi-structural benchmark datasets, to facilitate diversified, rigorous, and reproducible research on FWI. In particular, OpenFWI consists of 12 datasets (2.1TB in total) synthesized from multiple sources. It encompasses diverse domains in geophysics (interface, fault, CO2 reservoir, etc.), covers different geological subsurface structures (flat, curve, etc.), and contains various amounts of data samples (2K - 67K). It also includes a dataset for 3D FWI. Moreover, we use OpenFWI to perform benchmarking over four deep learning methods, covering both supervised and unsupervised learning regimes. Along with the benchmarks, we implement additional experiments, including physics-driven methods, complexity analysis, generalization study, uncertainty quantification, and so on, to sharpen our understanding of datasets and methods. The studies either provide valuable insights into the datasets and the performance, or uncover their current limitations. We hope OpenFWI supports prospective research on FWI and inspires future open-source efforts on AI for science. All datasets and related information can be accessed through our website at https://openfwi-lanl.github.io/


title: Causal Inference with Hidden Mediators

abstract: Proximal causal inference was recently proposed as a framework to identify causal effects from observational data in the presence of hidden confounders for which proxies are available. In this paper, we extend the proximal causal inference approach to settings where identification of causal effects hinges upon a set of mediators which are not observed, yet error prone proxies of the hidden mediators are measured. Specifically, (i) We establish causal hidden mediation analysis, which extends classical causal mediation analysis methods for identifying natural direct and indirect effects under no unmeasured confounding to a setting where the mediator of interest is hidden, but proxies of it are available. (ii) We establish hidden front-door criterion, which extends the classical front-door criterion to allow for hidden mediators for which proxies are available. (iii) We show that the identification of a certain causal effect called population intervention indirect effect remains possible with hidden mediators in settings where challenges in (i) and (ii) might co-exist. We view (i)-(iii) as important steps towards the practical application of front-door criteria and mediation analysis as mediators are almost always measured with error and thus, the most one can hope for in practice is that the measurements are at best proxies of mediating mechanisms. We propose identification approaches for the parameters of interest in our considered models. For the estimation aspect, we propose an influence function-based estimation method and provide an analysis for the robustness of the estimators.


title: Prediction of sunspot and plage coverage for Solar Cycle 25

abstract: Solar variability occurs over a broad range of spatial and temporal scales, from the Sun's brightening over its lifetime to the fluctuations commonly associated with magnetic activity over minutes to years. The latter activity includes most prominently the 11-year sunspot solar cycle and its modulations. Space weather events, in the form of solar flares, solar energetic particles, coronal mass ejections, and geomagnetic storms, have long been known to approximately follow the solar cycle occurring more frequently at solar maximum than solar minimum. These events can significantly impact our advanced technologies and critical infrastructures, making the prediction for the strength of future solar cycles particularly important. Several methods have been proposed to predict the strength of the next solar cycle, cycle 25, with results that are generally not always consistent. Most of these methods are based on the international sunspot number time series, or other indicators of solar activity. We present here a new approach that uses more than 100 years of measured fractional areas of the visible solar disk covered by sunspots and plages and an empirical relationship for each of these two indices of solar activity in even-odd cycles. We anticipate that cycle 25 will peak in 2024 and will last for about 12 years, slightly longer than cycle 24. We also found that, in terms of sunspot and plage areas coverage, the amplitude of cycle 25 will be substantially similar or slightly higher than cycle 24.


title: Mechanism of magnetic diode in artificial honeycomb lattice

abstract: Spin diode is important prerequisite to practical manifestation of spin electronics. Yet, a functioning magnetic diode at room temperature is still illusive. Here, we reveal diode-type phenomena due to magnetic charge mediated conduction in artificial honeycomb geometry, made of concave shape single domain permalloy element. We find that honeycomb lattice defies symmetry by populating vertices with low and high multiplicity magnetic charges, causing asymmetric magnetization, in applied current of opposite polarity. High multiplicity units create highly resistive network, thereby inhibiting magnetic charge dynamics propelled electrical conduction. However, practical realization of this effect requires modest demagnetization factor in constituting element. Concave structure fulfills the condition. Subsequently, magnetic diode behavior emerges across broad thermal range of $T$ = 40K - 300K. The finding is a departure from the prevailing notion of spin-charge interaction as the sole guiding principle behind spintronics. Consequently, a new vista, mediated by magnetic charge interaction, is envisaged for spintronic research


title: Decoupled coordinates for machine learning-based molecular fragment linking

abstract: Recent developments in machine-learning based molecular fragment linking have demonstrated the importance of informing the generation process with structural information specifying the relative orientation of the fragments to be linked. However, such structural information has not yet been provided in the form of a complete relative coordinate system. Mathematical details for a decoupled set of bond lengths, bond angles and torsion angles are elaborated and the coordinate system is demonstrated to be complete. Significant impact on the quality of the generated linkers is demonstrated numerically. The amount of reliable information within the different types of degrees of freedom is investigated. Ablation studies and an information-theoretical analysis are performed. The presented benefits suggest the application of a complete and decoupled relative coordinate system as a standard good practice in linker design.


title: Optomechanical Response of a Strongly Interacting Fermi Gas

abstract: We study a Fermi gas with strong, tunable interactions dispersively coupled to a high-finesse cavity. Upon probing the system along the cavity axis, we observe a strong optomechanical Kerr nonlinearity originating from the density response of the gas to the intracavity field and measure it as a function of interaction strength. We find that the zero-frequency density response function of the Fermi gas increases by a factor of two from the Bardeen-Cooper-Schrieffer to the Bose-Einstein condensate regime. The results are in quantitative agreement with a theory based on operator-product expansion, expressing the density response in terms of universal functions of the interactions, the contact and the internal energy of the gas. This provides an example of a driven-dissipative, strongly correlated system with a strong nonlinear response, opening up perspectives for the sensing of weak perturbations or inducing long-range interactions in Fermi gases.


title: Tutorial on Rational Rotation $C^*$--Algebras

abstract: The rotation algebra $\mathcal A_{\theta}$ is the universal $C^*$--algebra generated by unitary operators $U, V$ satisfying the commutation relation $UV = \omega V U$ where $\omega= e^{2\pi i \theta}.$ They are rational if $\theta = p/q$ with $1 \leq p \leq q-1,$ othewise irrational. Operators in these algebras relate to the quantum Hall effect \cite{boca,rammal,simon}, kicked quantum systems \cite{lawton1, wang}, and the spectacular solution of the Ten Martini problem \cite{avila}. Brabanter \cite{brabanter} and Yin \cite{yin} classified rational rotation $C^*$--algebras up to $*$-isomorphism. Stacey \cite{stacey} constructed their automorphism groups. They used methods known to experts: cocycles, crossed products, Dixmier-Douady classes, ergodic actions, K--theory, and Morita equivalence. This expository paper defines $\mathcal A_{p/q}$ as a $C^*$--algebra generated by two operators on a Hilbert space and uses linear algebra, Fourier series and the Gelfand-Naimark-Segal construction \cite{gelfand} to prove its universality. It then represents it as the algebra of sections of a matrix algebra bundle over a torus to compute its isomorphism class. The remarks section relates these concepts to general operator algebra theory. We write for mathematicians who are not $C^*$--algebra experts.


title: On an tangent equation by primes

abstract: In this paper we introduce a new diophantine equation with prime numbers. Let $[\, \cdot\,]$ be the floor function. We prove that when $1<c<\frac{23}{21}$ and $\theta>1$ is a fixed, then every sufficiently large positive integer $N$ can be represented in the form \begin{equation*} N=\big[p^c_1\tan^\theta(\log p_1)\big]+ \big[p^c_2\tan^\theta(\log p_2)\big]+ \big[p^c_3\tan^\theta(\log p_3)\big]\,, \end{equation*} where $p_1,\, p_2,\, p_3$ are prime numbers. We also establish an asymptotic formula for the number of such representations.


title: Charmonium-like resonances in coupled $D \bar D$-$D_s \bar D_s$ scattering

abstract: Charmonium-like resonances and bound states with isospin zero and $J^{PC}=0^{++},~1^{--},~2^{++},~3^{--}$ are extracted on the lattice. Coupled $D\bar D$ and $D_s\bar D_s$ scattering suggests three charmonium-like states with $J^{PC}=0^{++}$ in addition to $\chi_{c0}(1P)$: a so far unobserved $D\bar D$ bound state just below threshold, a conventional resonance likely related to $\chi_{c0}(3860)/\chi_{c0}(2P)$ and a narrow resonance just below the $D_s\bar D_s$ threshold with a large coupling to $D_s\bar D_s$ likely related to $X(3915)/\chi_{c0}(3930)$. One-channel $D\bar D$ scattering renders resonances and bound states with $J^{PC}= 1^{--},~2^{++},~3^{--}$ related to the observed conventional charmonia. Lattice QCD ensembles from the CLS consortium with $m_{\pi}\simeq 280$ MeV are utilized.


title: Apples and Oranges: Comparing black holes in X-ray binaries and gravitational-wave sources

abstract: The component black holes (BHs) observed in gravitational-wave (GW) binary black hole (BBH) events tend to be more massive and slower spinning than those observed in black hole X-ray binaries (BH-XRBs). Without modeling their evolutionary histories, we investigate whether these apparent tensions in the BH populations can be explained by GW observational selection effects alone. We find that this is indeed the case for the discrepancy between BH masses in BBHs and the observed high-mass X-ray binaries (HMXBs), when we account for statistical uncertainty from the small sample size of just three HMXBs. On the other hand, the BHs in observed low-mass X-ray binaries (LMXBs) are significantly lighter than the astrophysical BBH population, but this may just be due to a correlation between component masses in a binary system. Given their light stellar companions, we expect light BHs in LMXBs. The observed spins in HMXBs and LMXBs, however, are in tension with the inferred BBH spin distribution at the $>99.9\%$ level. We discuss possible scenarios behind the significantly larger spins in observed BH-XRBs. One possibility is that a small subpopulation (conservatively $<30\%$) of BBHs have rapidly spinning primary components, indicating that they may have followed a similar evolutionary pathway to the observed HMXBs. In LMXBs, it has been suggested that BHs can spin up by accretion. If LMXB natal spins follow the BBH spin distribution, we find LMXBs must gain an average dimensionless spin of $0.47^{+0.10}_{-0.11}$, but if their natal spins follow the observed HMXB spins, the average spinup must be $<0.03$.


title: Causal versus Marginal Shapley Values for Robotic Lever Manipulation Controlled using Deep Reinforcement Learning

abstract: We investigate the effect of including domain knowledge about a robotic system's causal relations when generating explanations. To this end, we compare two methods from explainable artificial intelligence, the popular KernelSHAP and the recent causal SHAP, on a deep neural network trained using deep reinforcement learning on the task of controlling a lever using a robotic manipulator. A primary disadvantage of KernelSHAP is that its explanations represent only the features' direct effects on a model's output, not considering the indirect effects a feature can have on the output by affecting other features. Causal SHAP uses a partial causal ordering to alter KernelSHAP's sampling procedure to incorporate these indirect effects. This partial causal ordering defines the causal relations between the features, and we specify this using domain knowledge about the lever control task. We show that enabling an explanation method to account for indirect effects and incorporating some domain knowledge can lead to explanations that better agree with human intuition. This is especially favorable for a real-world robotics task, where there is considerable causality at play, and in addition, the required domain knowledge is often handily available.


title: Geometry of the Gaussian graphical model of the cycle

abstract: We prove a conjecture due to Sturmfels and Uhler concerning the degree of the projective variety associated to the Gaussian graphical model of the cycle. We involve new methods based on the intersection theory in the space of complete quadrics.


title: Source-Level Bitwise Branching for Temporal Verification

abstract: There is increasing interest in applying verification tools to programs that have bitvector operations. SMT solvers, which serve as a foundation for these tools, have thus increased support for bitvector reasoning through bit-blasting and linear arithmetic approximations. Still, verification tools are limited on termination and LTL verification of bitvector programs. In this work, we show that similar linear arithmetic approximation of bitvector operations can be done at the source level through transformations. Specifically, we introduce new paths that over-approximate bitvector operations with linear conditions/constraints, increasing branching but allowing us to better exploit the well-developed integer reasoning and interpolation of verification tools. We present two sets of rules, namely rewriting rules and weakening rules, that can be implemented as bitwise branching of program transformation, the branching path can facilitate verification tools widen verification tasks over bitvector programs. Our experiment shows this exploitation of integer reasoning and interpolation enables competitive termination verification of bitvector programs and leads to the first effective technique for LTL verification of bitvector programs.


title: Effective weak and vague convergence of measures on the real line

abstract: We expand our effective framework for weak convergence of measures on the real line by showing that effective convergence in the Prokhorov metric is equivalent to effective weak convergence. In addition, we establish a framework for the study of the effective theory of vague convergence of measures. We introduce a uniform notion and a non-uniform notion of vague convergence, and we show that both these notions are equivalent. However, limits under effective vague convergence may not be computable even when they are finite. We give an example of a finite incomputable effective vague limit measure, and we provide a necessary and sufficient condition so that effective vague convergence produces a computable limit. Finally, we determine a sufficient condition for which effective weak and vague convergence of measures coincide. As a corollary, we obtain an effective version of the equivalence between classical weak and vague convergence of sequences of probability measures.


title: Complete singular collineations and quadrics

abstract: We construct wonderful compactifications of the spaces of linear maps, and symmetric linear maps of a given rank as blow-ups of secant varieties of Segre and Veronese varieties. Furthermore, we investigate their birational geometry and their relations with some spaces of degree two stable maps.


title: New Wormhole Solutions in a Viable $f(R)$ Gravity Model

abstract: Traversable wormhole solutions in General Relativity require exotic matter sources that violate the null energy condition (NEC), and such behavior maybe avoided in modified gravity. In this study, we analyze the energy conditions for static, spherically symmetric traversable Morris-Thorne wormholes in a recently proposed viable $f(R)$ gravity model. We numerically analyze solutions considering both constant and variable redshift functions, and present wormhole space-times respecting the NEC, supported by a phantom energy-like equation of state for the source. Moreover, we analyze the stability of the space-times using the generalized Tolman-Oppenheimer-Volkov equation. We demonstrate the effects of certain parameters in the $f(R)$ model in determining energy condition violations, and establish that stable wormholes can be formulated only at the expense of violating the NEC.


title: An investigation of factors affecting high-precision Sr isotope analyses (87Sr/86Sr and {\delta}88/86Sr) by MC-ICP-MS

abstract: The abundance and fractionation of the stable strontium (Sr) isotope system are being increasingly utilized to move forward our understanding in geological and cosmological processes. Two analytical techniques are commonly used to measure stable Sr isotopes: 1) double-spike thermal ionization mass spectrometry (DS-TIMS) and 2) Zr-doped sample-standard bracketing multi-collector inductively coupled plasma mass spectrometry (Zr-doped SSB via MC-ICP-MS). Relative to DS-TIMS, Zr-doped SSB via MC-ICP-MS allows simultaneous determinations of both 87Sr/86Sr and 88Sr/86Sr ratios, increasing measurement efficiency and sample throughput. However, this technique is currently associated with greater uncertainties in measurement precision and accuracy. In this study, we evaluated potential factors that can affect the quality of Sr isotope measurements during Zr-doped SSB. Our tests show that incomplete Sr recovery during chromatographic separation, mismatches of Sr and Zr concentrations and acid molarity between sample and bracketing standard, and cation contamination could all affect the precision and accuracy of Sr isotope measurements. We present evidence that, with updated preparation procedures and diligent concentration checks, a long-term reproducibility (2{\sigma}SD: 87Sr/86Sr = +/-0.000015 and {\delta}88/86Sr = +/-0.03 permil) comparable to that of DS-TIMS is achievable when using the Zr-doped SSB method via MC-ICP-MC.


title: idSTLPy: A Python Toolbox for Active Perception and Control

abstract: This paper describes a Python toolbox for active perception and control synthesis of probabilistic signal temporal logic (PrSTL) formulas of switched linear systems with additive Gaussian disturbances and measurement noises. We implement a counterexample-guided synthesis strategy that combines Bounded Model Checking, linear programming, and sampling-based motion planning techniques. We illustrate our approach and the toolbox throughout the paper with a motion planning example for a vehicle with noisy localization. The code is available at \url{https://codeocean.com/capsule/0013534/tree}.


title: A sharp isoperimetric inequality in metric measure spaces with non-negative Ricci curvature

abstract: We prove a sharp dimension-free isoperimetric inequality, involving the volume entropy, in non-compact metric measure spaces with non-negative synthetic Ricci curvature.


title: The outermost stellar halo of NGC 5128 (Centaurus A): Radial structure

abstract: The extended stellar halos of galaxies contain important clues for investigating their assembly history and evolution. We investigate the resolved stellar content and the extended halo of NGC 5128 as a function of galactocentric distance. We used HST images to resolve individual red giant branch (RGB) stars in 28 independent pointings. Star counts from deep VI color-magnitude diagrams reaching at least 1.5 mag below the tip of the RGB are used to derive the surface density distribution of the halo. The contamination by Milky Way stars is assessed with a new control field, with models, and by combining optical and near-IR photometry. We present a new calibration of the WFC3 F606W+F814W photometry to the ground-based VI photometric system. The photometry shows that the stellar halo of NGC 5128 is dominated by old RGB stars that are present in all fields. The V-band surface brightness changes from 23 to 32 mag/arcsec$^2$ between 8.3 kpc from the galaxy center to our outermost halo fields located 140 kpc away from the center along the major axis and 92 kpc along the minor axis. Within ~30 kpc, we also find evidence for a 2-3 Gyr old population traced by bright asymptotic giant branch stars. This population contributes only up to 10% in total stellar mass if it is 2 Gyr old, but a larger fraction of 30-40% is required if its age is 3 Gyr. The stellar surface density profile is well fit by a r$^{1/4}$ curve or a power-law $\sim r^{-3.1}$ over the full radial range, with no obvious break in the slope, but with large field-to-field scatter. The ellipticity measured from integrated-light photometry in the inner parts, $e=(b/a)=0.77$, flattens to $e=0.54 \pm 0.02$ beyond 30 kpc. Considering the flattening of the outer halo, the projection of the elliptical isophote on the semimajor axis for our most distant field reaches nearly 30 effective radii. [abridged]


title: Muonphilic Dark Matter explanation of gamma-ray galactic center excess: a comprehensive analysis

abstract: The Galactic center gamma-ray excess (GCE) is a long-standing unsolved problem. One of candidate solutions, the dark matter (DM) annihilation, has been recently tested with other astrophysical observations, such as AMS-02 electron-positron spectra, Fermi Dwarf spheroidal galaxies gamma-ray data, and so on. By assuming that the DM particles annihilate purely into a normal charged fermion pair, Di Mauro and Winkle (2021) claimed that only a muon-pair is compatible with the null detection of all the corresponding astrophysical measurements and can explain GCE simultaneously. On the other hand, a muonphilic DM model may also lead to a signal in the recent Fermilab muon $g-2$ measurement or be constrained by the latest PandaX-4T limit. In this work, we comprehensively study interactions between DM and muon, including various combinations of DM and mediator spins. In agreement with GCE (not only $2\mu$ but also $4\mu$ final states), we test these interactions against all the thermal DM constraints. Our results show that only the parameter space near the resonance region of mediator can explain GCE and relic density simultaneously, and larger parameter spaces are still allowed if other poorly-known systematic uncertainties are included. Regardless of the DM spin, only the interactions with the spin-0 mediator can explain the recent muon $g-2$ excess on top of GCE, relic density, and other DM and mediator constraints.


title: Variational Inference with Holder Bounds

abstract: The recent introduction of thermodynamic integration techniques has provided a new framework for understanding and improving variational inference (VI). In this work, we present a careful analysis of the thermodynamic variational objective (TVO), bridging the gap between existing variational objectives and shedding new insights to advance the field. In particular, we elucidate how the TVO naturally connects the three key variational schemes, namely the importance-weighted VI, Renyi-VI, and MCMC-VI, which subsumes most VI objectives employed in practice. To explain the performance gap between theory and practice, we reveal how the pathological geometry of thermodynamic curves negatively affects TVO. By generalizing the integration path from the geometric mean to the weighted Holder mean, we extend the theory of TVO and identify new opportunities for improving VI. This motivates our new VI objectives, named the Holder bounds, which flatten the thermodynamic curves and promise to achieve a one-step approximation of the exact marginal log-likelihood. A comprehensive discussion on the choices of numerical estimators is provided. We present strong empirical evidence on both synthetic and real-world datasets to support our claims.


title: Active Tuning of Resonant Lattice Kerker Effect

abstract: The Kerker effect has been generalized in nanophotonics and meta-optics, and has recently been of great interest by relating to various fascinating functionalities such as scattering management and perfect transmission, reflection or absorption. One of the most interesting generalizations is the resonant lattice Kerker effect in periodic nanostructures. However, its active tuning has not been explored yet. Here, we report, for the first time, the active control of the resonant lattice Kerker effect in periodic Ge2Se2Te5 nanodisks. By changing the crystalline fraction, we show that the electric dipole lattice resonance (ED-LR), the magnetic dipole resonance (MDR), and thus the resonant lattice Kerker effect are all red-shifted. We therefore realize the transition from the ED-LR to the resonant lattice Kerker effect, which enables multilevel tuning of reflection, transmission and absorption with modulation depths above 86%. Taking advantage of the MDR redshifts, we also observe broadband and multilevel tuning of transmission with modulation depth of 87% over a broadband range of 588 nm. Our work establishes a new path for designing high-performance active nanophotonic devices.


title: Finite-Time Consensus Learning for Decentralized Optimization with Nonlinear Gossiping

abstract: Distributed learning has become an integral tool for scaling up machine learning and addressing the growing need for data privacy. Although more robust to the network topology, decentralized learning schemes have not gained the same level of popularity as their centralized counterparts for being less competitive performance-wise. In this work, we attribute this issue to the lack of synchronization among decentralized learning workers, showing both empirically and theoretically that the convergence rate is tied to the synchronization level among the workers. Such motivated, we present a novel decentralized learning framework based on nonlinear gossiping (NGO), that enjoys an appealing finite-time consensus property to achieve better synchronization. We provide a careful analysis of its convergence and discuss its merits for modern distributed optimization applications, such as deep neural networks. Our analysis on how communication delay and randomized chats affect learning further enables the derivation of practical variants that accommodate asynchronous and randomized communications. To validate the effectiveness of our proposal, we benchmark NGO against competing solutions through an extensive set of tests, with encouraging results reported.


title: Intra- and interband excitations induced residue decay of the Bose polaron in a one-dimensional double-well

abstract: We investigate the polaronic properties of a single impurity immersed in a weakly interacting bosonic environment confined within a one-dimensional double-well potential using an exact diagonalization approach. We find that an increase of the impurity-bath coupling results in a vanishing residue, signifying the occurrence of the polaron orthogonality catastrophe. Asymptotic configurations of the systems' ground state wave function in the strongly interacting regime are obtained by means of a Schmidt decomposition, which in turn accounts for the observed orthogonality catastrophe of the polaron. We exemplify that depending on the repulsion of the Bose gas, three distinct residue behaviors appear with respect to the impurity-bath coupling. These residue regimes are characterized by two critical values of the bosonic repulsion and originate from the interplay between the intra- and the interband excitations of the impurity. Moreover, they can be clearly distinguished in the corresponding species reduced density matrices with the latter revealing a phase separation on either the one- or the two-body level. The impact of the interspecies mass-imbalance on the impurity's excitation processes is appreciated yielding an interaction shift of the residue regions. Our results explicate the interplay of intra- and interband excitation processes for the polaron generation in multiwell traps and for designing specific polaron entangled states motivating their exposure in current experiments.


title: Quantum tangent kernel

abstract: Quantum kernel method is one of the key approaches to quantum machine learning, which has the advantages that it does not require optimization and has theoretical simplicity. By virtue of these properties, several experimental demonstrations and discussions of the potential advantages have been developed so far. However, as is the case in classical machine learning, not all quantum machine learning models could be regarded as kernel methods. In this work, we explore a quantum machine learning model with a deep parameterized quantum circuit and aim to go beyond the conventional quantum kernel method. In this case, the representation power and performance are expected to be enhanced, while the training process might be a bottleneck because of the barren plateaus issue. However, we find that parameters of a deep enough quantum circuit do not move much from its initial values during training, allowing first-order expansion with respect to the parameters. This behavior is similar to the neural tangent kernel in the classical literatures, and such a deep variational quantum machine learning can be described by another emergent kernel, quantum tangent kernel. Numerical simulations show that the proposed quantum tangent kernel outperforms the conventional quantum kernel method for an ansatz-generated dataset. This work provides a new direction beyond the conventional quantum kernel method and explores potential power of quantum machine learning with deep parameterized quantum circuits.


title: Distinguishing the aftershock sequences of the DPRK5 and DPRK6 underground nuclear tests

abstract: The rate of aftershocks in the sequence initiated by the DPRK underground tests has been increasing since January 2021. In total, 22 reliable aftershocks were detected between January 13 and October 1, 2021. Their characteristics are similar to the aftershocks in one of two clusters: 1) the fifth DPRK (DPRK5) test (mb(IDC)=5.09) conducted on September 9, 2016, which induced the first DPRK aftershock in the sequence detected at 1:50:48 UTC on September 11, 2016; 2) the sixth DPRK (DPRK6) explosion (mb(IDC)=6.07), which generates its aftershock sequence with characteristics significantly different from the aftershocks in the DPRK5 sequence. The length, intensity, and alternating character of these sequences suggest specific mechanisms of energy release likely associated with the interaction of the damaged zones of the DPRK5 and DPRK6 and the collapse of their cavities with progressive propagation of the collapsing chimneys to the free surface. According to the depth estimates based on the moment tensor modelling, the DPRK5 and DPRK6 were conducted at practically the same depths. The difference in magnitudes suggests that their damaged zones differ by a factor of 2 or more. The first aftershock of the DPRK6 (mb(IDC=4.12) 8.5 minutes after the test is evidence of the cavity collapse and creation of a chimney, which did not reach the surface. The activity in 2021 indicates that the chimney collapse is not finished yet. One can expect more aftershocks in the near future, likely ended with the chimney reaching the free surface.


title: Higher partial wave contamination in finite-volume formulae for 1-to-2 transitions

abstract: It is common practice to truncate the finite-volume formula for $K\to\pi\pi$, and other one-to-two transitions, to only include the lowest partial wave, as in the original derivation by Lellouch and L\"uscher. However, as the precision of lattice calculations increases, it may become important to asses the systematic effect of this approximation. With this motivation, we compare the $S$-wave-only ($\ell=0$) results with those truncated at the next lowest value of angular momentum.


title: A data-based parametrization of parton distribution functions

abstract: Since the first determination of a structure function many decades ago, all methodologies used to determine structure functions or parton distribution functions (PDFs) have employed a common prefactor as part of the parametrization. The NNPDF collaboration pioneered the use of neural networks to overcome the inherent bias of constraining the space of solution with a fixed functional form while still keeping the same common prefactor as a preprocessing. Over the years various, increasingly sophisticated, techniques have been introduced to counter the effect of the prefactor on the PDF determination. In this paper we present a methodology to remove the prefactor entirely, thereby significantly simplifying the methodology, without a loss of efficiency and finding good agreement with previous results.


title: A spatio-temporal analogue of the Omori-Utsu law of aftershock sequences

abstract: A spatio-temporal version of the well-known Omori-Utsu law of aftershock sequences is proposed. This 'diffusive Omori-Utsu law' satisfies a nonlinear partial differential equation (PDE). A similarity reduction is obtained that reduces the PDE to an ordinary differential equation (ODE). A nonzero constant solution of this ODE leads to the usual Omori-Utsu law. An exact and explicit similarity solution is found that corresponds to the original Omori law. An initial value problem for the 'diffusive Omori-Utsu law' is also considered, and whose spatio-temporal dynamics are described by bounding functions that satisfy nonlinear, but linearisable, PDEs. Numerical results are also provided.


title: Novel theoretical considerations of the coefficient of earth pressure at rest

abstract: One of the well-known and widely applied parameters in soil mechanics is the so-called coefficient of earth pressure at rest (denoted by a K0). Field and laboratory investigations show that this parameter is correlated with the friction angle and the overconsolidation ratio of the soil deposit. For normally consolidated clays, the expression developed by Jaky is considered to hold good. The extension of Jaky's equation for overconsolidated soils is built mainly on empirical observations; and there is lack of a theoretical framework that satisfactorily explains the relationship. In this paper, these relationships are investigated from the nature of energy dissipation of soil aggregates. For the same, the cyclic stress-dilatancy relationship proposed by the author is employed. Then, the stress ratio that maximizes entropy is derived; and this stress ratio is shown to be reasonably close to the coefficient earth pressure ratio proposed by Jaky. An exponential relationship between K0 and OCR is derived; and for a special case, the exponent is shown to be the sine of the critical state friction angle. Finally, some generalizations are given and aspects which need further investigation are highlighted.


title: Microwave frequency dissemination systems as sensitive and low-cost interferometers for earthquake detection on commercially deployed fiber cables

abstract: We experimentally demonstrate a microwave frequency dissemination system operating as a sensitive interferometric sensor of seismic waves on commercially deployed fiber networks in Attika, Greece. Efficient detection of seismic waves from distant epicenters (>400km)


title: Reactive transport experiments of coupled carbonation and serpentinization in a natural serpentinite. Implication for hydrogen production and carbon geological storage

abstract: Serpentinization and carbonation of ultramafic formations is a ubiquitous phenomenon, which deeply influences the biogeochemical cycles of water, hydrogen, carbon... while supporting the particular biosphere around the oceanic hydrothermal vents. Carbonation of peridotites and other mafic and ultramafic rocks is also a hot topic in the current energy landscape as the engineered sequestration of mineral CO2 in these formations could help reduce the atmospheric emissions and cope with climate change. In this study, we present two reactive percolation experiments performed on a natural serpentinite dredged from the ultraslow South-West Indian Oceanic Ridge. The serpentinite cores (length 3-4 cm and dia. 5.6 mm) were subjected for about 10 days to the continuous injection of a NaHCO3-saturated brine at respectively 160{\deg}C and 280{\deg}C. Petrographic and petrophysical results as well as outlet fluid compositions were compared to numerical batch simulations performed with the PHREEQC open software allowing to reconstruct the mineralogical evolution of both cores. The most striking observation is the fast and dramatic decrease of the permeability for both experiments principally due to the precipitation of carbonates. On the contrary, serpentine was found to be less impacting as it precipitates in low-flow zones, out of the main percolation paths. In total, about 5.6% of the total injected CO2 was retained in the core, at 280{\deg}C. In the same time, hydrogen was consistently produced with a total recovered H2 corresponding to 0.8% of the maximum H2 possible. The global behavior of the cores is interpreted as the result from an interplay between interacting spatio-temporal lengthscales controlled by the Damkohler number.


title: Experimental study and reaction path modeling of the carbonation of natural serpentinites

abstract: Mineralization of carbon dioxide is often seen as an attractive alternative to classical Carbon Capture and Storage (CCS) technologies, allowing the sequestration of \ce{CO2} as a solid mineral with no risk of aquifer contamination or leakage back to the atmosphere. While olivine and pyroxenes are known to easily and quickly react with dissolved \ce{CO2}, fresh peridotites are quite rare and ultramafic rocks usually contain significant amounts of serpentine, which presents a lower reactivity. The purpose of this study was then to analyze the reactivity of two natural rocks: a partially serpentinized lherzolite and a fully altered serpentinite. Results confirm that serpentine is much slower to react and gets altered only if the activity of \ce{CO2} is high enough and if all olivine and pyroxenes have already been consumed. Resulting carbonates are mostly Mg-rich calcite or Mg-depleted dolomite with the occurrence of eitelite (\ce{Na2Mg(CO3)2} in the case of high Na activities. The carbonation of these serpentinite was however associated in some cases with a heavy precipitation of hazardous asbestiform chrysotile, which could be a potential threat for engineered carbonation processes.


title: Use of low-fidelity models with machine-learning error correction for well placement optimization

abstract: Well placement optimization is commonly performed using population-based global stochastic search algorithms. These optimizations are computationally expensive due to the large number of multiphase flow simulations that must be conducted. In this work, we present an optimization framework in which these simulations are performed with low-fidelity (LF) models. These LF models are constructed from the underlying high-fidelity (HF) geomodel using a global transmissibility upscaling procedure. Tree-based machine-learning methods, specifically random forest and light gradient boosting machine, are applied to estimate the error in objective function value (in this case net present value, NPV) associated with the LF models. In the offline (preprocessing) step, preliminary optimizations are performed using LF models, and a clustering procedure is applied to select a representative set of 100--150 well configurations to use for training. HF simulation is then performed for these configurations, and the tree-based models are trained using an appropriate set of features. In the online (runtime) step, optimization with LF models, with the machine-learning correction, is conducted. Differential evolution is used for all optimizations. Results are presented for two example cases involving the placement of vertical wells in 3D bimodal channelized geomodels. We compare the performance of our procedure to optimization using HF models. In the first case, 25 optimization runs are performed with both approaches. Our method provides an overall speedup factor of 46 relative to optimization using HF models, with the best-case NPV within 1% of the HF result. In the second case fewer HF optimization runs are conducted (consistent with actual practice), and the overall speedup factor with our approach is about 8. In this case, the best-case NPV from our procedure exceeds the HF result by 3.8%


title: Mobility-Dependent and Mobility-Compensated Effective Reproduction Number of COVID-19 Viral Variants: New Metric for Infectivity Evaluation

abstract: During epidemics, estimation of the effective reproduction number (ERN) associated with infectious disease is a challenging topic for policy development and medical resource management. There is still an open question about the dominant factors to characterize in corona virus disease 2019 (COVID-19), although recent studies based on nonlinear regression with machine learning suggested mobility. The emergence of new viral variants is common in widespread pandemics. However, understanding the potential ERN of new variants is required for policy revision, including lockdown constraints. In this study, we proposed time-averaged mobility at transit stations as a surrogate to correlate with ERN using data from three prefectures in Japan. The latency and duration to average over the mobility were 6-8 days and 6-7, respectively (R2 was 0.109-0.512 in Tokyo, 0.365-0.607 in Osaka, and 0.317-0.631 in Aichi). The same linear correlation was confirmed in Singapore and London. The mobility-adjusted ERN of the alpha variant was 15%-30%, and was 20%-40% higher than the standard type in Osaka, Aichi, and London. Similarly, the ERN of the delta variant was 20%-40% higher than that of the standard type in Osaka and Aichi. The proposed metric would be useful for proper evaluation of the infectivity of different variants in terms of ERN.


title: The Gompertz Growth of COVID-19 Outbreaks is Caused by Super-Spreaders

abstract: In individual SARS-CoV-2 outbreaks, the count of confirmed cases and deaths follow a Gompertz growth function for locations of very different sizes. This lack of dependence on region size leads us to hypothesize that virus spread depends on universal properties of the network of social interactions. We test this hypothesis by simulating the propagation of a virus on networks of different topologies. Our main finding is that Gompertz growth observed for early outbreaks occurs only for a scale-free network, in which nodes with many more neighbors than average are common. These nodes that have very many neighbors are infected early in the outbreak and then spread the infection very rapidly. When these nodes are no longer infectious, the remaining nodes that have most neighbors take over and continue to spread the infection. In this way, the rate of spread is fastest at the very start and slows down immediately. Geometrically it is seen that the "surface" of the epidemic, the number of susceptible nodes in contact with the infected nodes, starts to rapidly decrease very early in the epidemic and as soon as the larger nodes have been infected. In our simulation, the speed and impact of an outbreak depend on three parameters: the average number of contacts each node makes, the probability of being infected by a neighbor, and the probability of recovery. Intelligent interventions to reduce the impact of future outbreaks need to focus on these critical parameters in order to minimize economic and social collateral damage.


title: Unitarily inequivalent quantum cosmological bouncing models

abstract: By quantising the background as well as the perturbations in a simple one fluid model, we show that there exists an ambiguity in the choice of relevant variables, potentially leading to incompatible observational physical predictions. In a classical or quantum inflationary background, the exact same canonical transformations lead to unique predictions, so the ambiguity we put forward demands a semiclassical background with a sufficiently strong departure from classical evolution. The latter condition happens to be satisfied in bouncing scenarios, which may thus be having predictability issues. Inflationary models could evade such a problem because of the monotonic behavior of their scale factor; they do, however, initiate from a singular state which bouncing scenarios aim at solving.


title: Using Graph-Theoretic Machine Learning to Predict Human Driver Behavior

abstract: Studies have shown that autonomous vehicles (AVs) behave conservatively in a traffic environment composed of human drivers and do not adapt to local conditions and socio-cultural norms. It is known that socially aware AVs can be designed if there exists a mechanism to understand the behaviors of human drivers. We present an approach that leverages machine learning to predict, the behaviors of human drivers. This is similar to how humans implicitly interpret the behaviors of drivers on the road, by only observing the trajectories of their vehicles. We use graph-theoretic tools to extract driver behavior features from the trajectories and machine learning to obtain a computational mapping between the extracted trajectory of a vehicle in traffic and the driver behaviors. Compared to prior approaches in this domain, we prove that our method is robust, general, and extendable to broad-ranging applications such as autonomous navigation. We evaluate our approach on real-world traffic datasets captured in the U.S., India, China, and Singapore, as well as in simulation.


title: The stable rank of $\mathbb{Z}[x]$ is $3$

abstract: Grunewald, Mennicke and Vaserstein proved that the Bass stable rank of $\mathbb{Z}[x]$, the ring of the univariate polynomials over $\mathbb{Z}$, is $3$. This note addresses minor errors found in their proof. Using their method, we show in addition that the unimodular row $(3, x + 1, x^2 + 16)$ is not stable.


title: Consistent Estimation for PCA and Sparse Regression with Oblivious Outliers

abstract: We develop machinery to design efficiently computable and consistent estimators, achieving estimation error approaching zero as the number of observations grows, when facing an oblivious adversary that may corrupt responses in all but an $\alpha$ fraction of the samples. As concrete examples, we investigate two problems: sparse regression and principal component analysis (PCA). For sparse regression, we achieve consistency for optimal sample size $n\gtrsim (k\log d)/\alpha^2$ and optimal error rate $O(\sqrt{(k\log d)/(n\cdot \alpha^2)})$ where $n$ is the number of observations, $d$ is the number of dimensions and $k$ is the sparsity of the parameter vector, allowing the fraction of inliers to be inverse-polynomial in the number of samples. Prior to this work, no estimator was known to be consistent when the fraction of inliers $\alpha$ is $o(1/\log \log n)$, even for (non-spherical) Gaussian design matrices. Results holding under weak design assumptions and in the presence of such general noise have only been shown in dense setting (i.e., general linear regression) very recently by d'Orsi et al. [dNS21]. In the context of PCA, we attain optimal error guarantees under broad spikiness assumptions on the parameter matrix (usually used in matrix completion). Previous works could obtain non-trivial guarantees only under the assumptions that the measurement noise corresponding to the inliers is polynomially small in $n$ (e.g., Gaussian with variance $1/n^2$). To devise our estimators, we equip the Huber loss with non-smooth regularizers such as the $\ell_1$ norm or the nuclear norm, and extend d'Orsi et al.'s approach [dNS21] in a novel way to analyze the loss function. Our machinery appears to be easily applicable to a wide range of estimation problems.


title: An Empirical Comparison of the Quadratic Sieve Factoring Algorithm and the Pollard Rho Factoring Algorithm

abstract: One of the most significant challenges on cryptography today is the problem of factoring large integers since there are no algorithms that can factor in polynomial time, and factoring large numbers more than some limits(200 digits) remain difficult. The security of the current cryptosystems depends on the hardness of factoring large public keys. In this work, we want to implement two existing factoring algorithms - pollard-rho and quadratic sieve - and compare their performance. In addition, we want to analyze how close is the theoretical time complexity of both algorithms compared to their actual time complexity and how bit length of numbers can affect quadratic sieve's performance. Finally, we verify whether the quadratic sieve would do better than pollard-rho for factoring numbers smaller than 80 bits.


title: A Modification of the Random Cutting Model

abstract: We propose a modification to the random destruction of graphs: Given a finite network with a distinguished set of sources and targets, remove (cut) vertices at random, discarding components that do not contain a source node. We investigate the number of cuts required until all targets are removed, and the size of the remaining graph. This model interpolates between the random cutting model going back to Meir and Moon and site percolation. We prove several general results, including that the size of the remaining graph is a tight family of random variables for compatible sequences of expander-type graphs, and determine limiting distributions for binary caterpillar trees and complete binary trees.


title: Isomonodromic deformations along a stratum of the coalescence locus

abstract: We consider deformations of a differential system with Poincare' rank 1 at infinity and Fuchsian singularity at zero along a stratum of a coalescence locus. We give necessary and sufficient conditions for the deformation to be strongly isomonodromic, both as an explicit Pfaffian system (integrable deformation) and as a non linear system of PDEs on the residue matrix A at the Fuchsian singularity. This construction is complementary to that of [13]. For the specific system here considered, the results generalize those of [26], by giving up the generic conditions, and those of [3], by giving up the Lidskii generic assumption. The importance of the case here considered originates form its applications in the study of strata of Dubrovin-Frobenius manifolds and F-manifolds.


title: Consensus-based Optimization and Ensemble Kalman Inversion for Global Optimization Problems with Constraints

abstract: We introduce a practical method for incorporating equality and inequality constraints in global optimization methods based on stochastic interacting particle systems, specifically consensus-based optimization (CBO) and ensemble Kalman inversion (EKI). Unlike other approaches in the literature, the method we propose does not constrain the dynamics to the feasible region of the state space at all times; the particles evolve in the full space, but are attracted towards the feasible set by means of a penalization term added to the objective function and, in the case of CBO, an additional relaxation drift. We study the properties of the method through the associated mean-field Fokker--Planck equation and demonstrate its performance in numerical experiments on several test problems.


title: Shifted edge labeled tableaux and localizations

abstract: We prove cases of a conjectural rule of H. Yadav, A. Yong, and the author for structure coefficients of the D. Anderson-W. Fulton ring. In particular, we give a combinatorial description for certain localization coefficients of this ring, which is related to the equivariant cohomology of isotropic Grassmannians.


title: Stein Variational Probabilistic Roadmaps

abstract: Efficient and reliable generation of global path plans are necessary for safe execution and deployment of autonomous systems. In order to generate planning graphs which adequately resolve the topology of a given environment, many sampling-based motion planners resort to coarse, heuristically-driven strategies which often fail to generalize to new and varied surroundings. Further, many of these approaches are not designed to contend with partial-observability. We posit that such uncertainty in environment geometry can, in fact, help drive the sampling process in generating feasible, and probabilistically-safe planning graphs. We propose a method for Probabilistic Roadmaps which relies on particle-based Variational Inference to efficiently cover the posterior distribution over feasible regions in configuration space. Our approach, Stein Variational Probabilistic Roadmap (SV-PRM), results in sample-efficient generation of planning-graphs and large improvements over traditional sampling approaches. We demonstrate the approach on a variety of challenging planning problems, including real-world probabilistic occupancy maps and high-dof manipulation problems common in robotics.


title: An equidistribution involving invisible inversions

abstract: We provide two explicit bijections demonstrating that, among permutations, the number of invisible inversions is equidistributed with the number of occurrences of the vincular pattern 13-2 after sorting the set of runs.


title: The discrepancy of unsatisfiable matrices and a lower bound for the Koml\'{o}s conjecture constant

abstract: We construct simple, explicit matrices with columns having unit $\ell^2$ norm and discrepancy approaching $1 + \sqrt{2} \approx 2.414$. This number gives a lower bound, the strongest known as far as we are aware, on the constant appearing in the Koml\'{o}s conjecture. The "unsatisfiable matrices" giving this bound are built by scaling the entries of clause-variable matrices of certain unsatisfiable Boolean formulas. We show that, for a given formula, such a scaling maximizing a lower bound on the discrepancy may be computed with a convex second-order cone program. Using a dual certificate for this program, we show that our lower bound is optimal among those using unsatisfiable matrices built from formulas admitting read-once resolution proofs of unsatisfiability. We also conjecture that a generalization of this certificate shows that our bound is optimal among all bounds using unsatisfiable matrices.


title: Approximating Invertible Maps by Recovery Channels: Optimality and an Application to Non-Markovian Dynamics

abstract: We investigate the problem of reversing quantum dynamics, specifically via optimal Petz recovery maps. We focus on typical decoherence channels, such as dephasing, depolarizing and amplitude damping. We illustrate how well a physically implementable recovery map simulates an inverse evolution. We extend this idea to explore the use of recovery maps as an approximation of inverse maps, and apply it in the context of non-Markovian dynamics. We show how this strategy attenuates non-Markovian effects, such as the backflow of information.


title: Post-Minkowskian Radial Action from Soft Limits and Velocity Cuts

abstract: We consider gravitational massive scalar-scalar scattering from unitarity and demonstrate how intermediate soft graviton behavior and the concept of extracting classical physics from localization of integrands on velocity cuts devise an efficient extraction scheme for computing the classical post-Minkowskian radial action perturbatively. We demonstrate the computational efficiency by deriving the scattering amplitudes in the probe regime to the fifth post-Minkowskian order in arbitrary dimensions.


title: Safe, efficient and socially-compatible decision of automated vehicles: a case study of unsignalized intersection driving

abstract: Safe and smooth interacting with other vehicles is one of the ultimate goals of driving automation. However, recent reports of demonstrative deployments of automated vehicles (AVs) indicate that AVs are still difficult to meet the expectation of other interacting drivers, which leads to several AV accidents involving human-driven vehicles (HVs). This is most likely due to the lack of understanding about the dynamic interaction process, especially about the human drivers. By investigating the causes of 4,300 video clips of traffic accidents, we find that the limited dynamic visual field of drivers is one leading factor in inter-vehicle interaction accidents, especially in those involving trucks. A game-theoretic decision algorithm considering social compatibility is proposed to handle the interaction with a human-driven truck at an unsignalized intersection. Starting from a probabilistic model for the visual field characteristics of truck drivers, social fitness and reciprocal altruism in the decision are incorporated in the game payoff design. Human-in-the-loop experiments are carried out, in which 24 subjects are invited to drive and interact with AVs deployed with the proposed algorithm and two comparison algorithms. Totally 207 cases of intersection interactions are obtained and analyzed, which shows that the proposed decision-making algorithm can not only improve both safety and time efficiency, but also make AV decisions more in line with the expectation of interacting human drivers. These findings can help inform the design of automated driving decision algorithms, to ensure that AVs can be safely and efficiently integrated into the human-dominated traffic.


title: Gradient descent globally solves average-case non-resonant physical design problems

abstract: Optimization problems occurring in a wide variety of physical design problems, including but not limited to optical engineering, quantum control, structural engineering, involve minimization of a simple cost function of the state of the system (e.g. the optical fields, the quantum state) while being constrained by the physics of the system. The physics constraints often makes such problems non-convex and thus only locally solvable, leaving open the question of finding the globally optimal design. In this paper, I consider design problems whose physics is described by bi-affine equality constraints, and show that under assumptions on the stability of these constraints and the physical system being non-resonant, gradient descent globally solves a typical physical design problem. The key technical contributions of this paper are (i) outline a criteria that ensure the convergence of gradient descent to an approximate global optima in the limit of large problem sizes, and (ii) use random matrix theory to outline ensembles of physically motivated problems which, on an average, satisfy this convergence criteria.


title: Barrier States Embedded Iterative Dynamic Game for Robust and Safe Trajectory Optimization

abstract: Considering uncertainties and disturbances is an important, yet challenging, step in successful decision making. The problem becomes more challenging in safety-constrained environments. In this paper, we propose a robust and safe trajectory optimization algorithm through solving a constrained min-max optimal control problem. The proposed method leverages a game theoretic differential dynamic programming approach with barrier states to handle parametric and non-parametric uncertainties in safety-critical control systems. Barrier states are embedded into the differential game's dynamics and cost to portray the constrained environment in a higher dimensional state space and certify the safety of the optimized trajectory. Moreover, to find a convergent optimal solution, we propose to perform line-search in a Stackleberg (leader-follower) game fashion instead of picking a constant learning rate. The proposed algorithm is evaluated on a velocity-constrained inverted pendulum model in a moderate and high parametric uncertainties to show its efficacy in such a comprehensible system. The algorithm is subsequently implemented on a quadrotor in a windy environment in which sinusoidal wind turbulences applied in all directions.


title: Gettering in polySi/SiOx passivating contacts enables Si-based tandem solar cells with high thermal and contamination resilience

abstract: Multijunction solar cells in a tandem configuration could further lower the costs of electricity if crystalline Si (c-Si) is used as bottom cell. However, for direct monolithic integration on c-Si, only a restricted number of top and bottom cell architectures are compatible, due to either epitaxy or high temperature constraints, where the interface between subcells is subject to a trade-off between transmittance, electrical interconnection, and bottom cell degradation. Using polySi/SiOx passivating contacts for Si, this degradation can be largely circumvented by tuning the polySi/SiOx stacks to promote gettering of contaminants admitted into the Si bottom cell during the top cell synthesis. Applying this concept to the low-cost top cell chalcogenides Cu2ZnSnS4 (CZTS), CuGaSe2 (CGSe) and AgInGaSe2 (AIGSe), fabricated under harsh S or Se atmospheres above 550 {\deg}C, we show that increasing the heavily-doped polySi layer thickness from 40 to up to 400 nm prevents a reduction in Si carrier lifetime by one order of magnitude, with final lifetimes above 500 us uniformly across areas up to 20 cm2. In all cases, the increased resilience was correlated with a 99.9% reduction in contaminant concentration in the c-Si bulk, provided by the thick polySi layer, which acts as a buried gettering layer in the tandem structure without compromising the Si passivation quality. The Si resilience decreased as AIGSe > CGSe > CZTS, in accordance with the measured Cu contamination profiles and higher annealing temperatures. An efficiency of up to 7% was achieved for a CZTS/Si tandem, where the Si bottom cell is no longer the limiting factor.


title: Uniqueness of boundary tangent cones for $2$-dimensional area-minimizing currents

abstract: In this paper we show that, if $T$ is an area-minimizing $2$-dimensional integral current with $\partial T = Q [\![ \Gamma ]\!]$, where $\Gamma$ is a $C^{1,\alpha}$ curve for $\alpha>0$ and $Q$ an arbitrary integer, then $T$ has a unique tangent cone at every boundary point, with a polynomial convergence rate. The proof is a simple reduction to the case $Q=1$, studied by Hirsch and Marini.


title: Nuclear two point correlation functions on a quantum-computer

abstract: The calculation of dynamic response functions is expected to be an early application benefiting from rapidly developing quantum hardware resources. The ability to calculate real-time quantities of strongly-correlated quantum systems is one of the most exciting applications that can easily reach beyond the capabilities of traditional classical hardware. Response functions of fermionic systems at moderate momenta and energies corresponding roughly to the Fermi energy of the system are a potential early application because the relevant operators are nearly local and the energies can be resolved in moderately short real time, reducing the spatial resolution and gate depth required. This is particularly the case in quasielastic electron and neutrino scattering from nuclei, a topic of great interest in the nuclear and particle physics communities and directly related to experiments designed to probe neutrino properties. In this work we use current quantum hardware and error mitigation protocols to calculate response functions for a highly simplified nuclear model through calculations of a 2-point real time correlation function for a modified Fermi-Hubbard model in two dimensions with three distinguishable nucleons on four lattice sites.


title: Turing patterns resulting from a Sturm-Liouville problem

abstract: Pattern formation in reaction-diffusion systems where the diffusion terms correspond to a Sturm-Liouville problem are studied. These correspond to a problem where the diffusion coefficient depends on the spatial variable: $\nabla \cdot \left( \mathcal{D} ( {\mathbf x} ) \nabla {\mathbf u} \right)$. We found that the conditions for Turing instability are the same as in the case of homogeneous diffusion but the nonlinear analysis must be generalized to consider general orthogonal eigenfunctions instead of the standard Fourier approach. The particular case $\mathcal{D} (x)= 1-x^2$, where solutions are linear combinations of Legendre polynomials, is studied in detail. From the developed general nonlinear analysis, conditions for producing stripes and spots are obtained, which are numerically verified using the Schaneknberg system. Unlike to the case with homogeneous diffusion, and due to the properties of the Legendre polynomials, stripped and spotted patterns with variable wavelength are produced, and a change from stripes to spots is predicted when the wavelength increases. The patterns obtained can model biological systems where stripes or spots accumulate close to the boundaries and the theory developed here can be applied to study Turing patterns associated to other eigenfunctions related with Sturm-Liouville problems.


title: Multiple low-energy excitons and optical response of $d^0$ double perovskite Ba$_2$ScTaO$_6$

abstract: Large bandgap insulators are considered promising for applications such as photocatalysts, dielectric resonators and interference filters. Based on synchrotron X-ray diffraction, diffuse reflectance measurement and density functional theory, we report the crystal structure, optical response, and electronic properties of the synthesized $d^0$ double perovskite Ba$_2$ScTaO$_6$. In contrast to earlier prediction, the electronic bandgap is found to be large, $\sim 4.66$ eV. The optical response is characterized by the presence of multiple exciton modes extending up to the visible range. A detailed investigation of the direct gap excitons based on the Elliot formula is presented. Density functional theory based investigation of the electronic properties within generalized gradient approximation severely underestimates the electronic gap. To reach a quantitative agreement, we consider different available flavors of the modified-Becke-Johnson exchange-correlation potential and discuss their effects on the electronic and optical properties.


title: A QoE Model in Point Cloud Video Streaming

abstract: Point cloud video has been widely used by augmented reality (AR) and virtual reality (VR) applications as it allows users to have an immersive experience of six degrees of freedom (6DoFs). Yet there is still a lack of research on quality of experience (QoE) model of point cloud video streaming, which cannot provide optimization metric for streaming systems. Besides, position and color information contained in each pixel of point cloud video, and viewport distance effect caused by 6DoFs viewing procedure make the traditional objective quality evaluation metric cannot be directly used in point cloud video streaming system. In this paper we first analyze the subjective and objective factors related to QoE model. Then an experimental system to simulate point cloud video streaming is setup and detailed subjective quality evaluation experiments are carried out. Based on collected mean opinion score (MOS) data, we propose a QoE model for point cloud video streaming. We also verify the model by actual subjective scoring, and the results show that the proposed QoE model can accurately reflect users' visual perception. We also make the experimental database public to promote the QoE research of point cloud video streaming.


title: Quantum trajectories, interference, and state localisation in dephasing assisted quantum transport

abstract: Dephased quantum transport of excitations occurs when energetic fluctuations in a system are sufficient to suppress the built-up of coherent amplitudes. While this has been extensively studied in many different systems, a unified and comprehensive understanding of quantum assisted transport via on-site and coupling-induced dephasing processes is lacking. The aim of the present work is to present a simple and unified understanding of the role of these two key dephasing processes in dephasing assisted transport. Our work explicitly links continuous dephasing to classical and quantum transitions. We present a natural quantum trajectories explanation of how different coupling and dephasing terms alter the diffusion rate of excitations and how this is impacted by the onset of Anderson localized eigenstates. Our results provide insight in understanding quantum transport in molecular semiconductors, artificial lattices and quantum features of excitonic solids.


title: Numerical Approximation in CFD Problems Using Physics Informed Machine Learning

abstract: The thesis focuses on various techniques to find an alternate approximation method that could be universally used for a wide range of CFD problems but with low computational cost and low runtime. Various techniques have been explored within the field of machine learning to gauge the utility in fulfilling the core ambition. Steady advection diffusion problem has been used as the test case to understand the level of complexity up to which a method can provide solution. Ultimately, the focus stays over physics informed machine learning techniques where solving differential equations is possible without any training with computed data. The prevalent methods by I.E. Lagaris et.al. and M. Raissi et.al are explored thoroughly. The prevalent methods cannot solve advection dominant problems. A physics informed method, called as Distributed Physics Informed Neural Network (DPINN), is proposed to solve advection dominant problems. It increases the lexibility and capability of older methods by splitting the domain and introducing other physics-based constraints as mean squared loss terms. Various experiments are done to explore the end to end possibilities with the method. Parametric study is also done to understand the behavior of the method to different tunable parameters. The method is tested over steady advection-diffusion problems and unsteady square pulse problems. Very accurate results are recorded. Extreme learning machine (ELM) is a very fast neural network algorithm at the cost of tunable parameters. The ELM based variant of the proposed model is tested over the advection-diffusion problem. ELM makes the complex optimization simpler and Since the method is non-iterative, the solution is recorded in a single shot. The ELM based variant seems to work better than the simple DPINN method. Simultaneously scope for various development in future are hinted throughout the thesis.


title: Abnormal extremals of left-invariant sub-Finsler quasimetrics on four-dimensional Lie groups with three-dimensional generating distributions

abstract: We find three-dimensional subspaces of four-dimensional connected Lie algebras, generating these algebras, and abnormal extremals on connected Lie groups with these Lie algebras and with left-invariant sub-Finsler quasimetrics defined by seminorms on such subspaces. In terms of the structure constants of Lie algebras and dual seminorms, we establish a criterion for the strong abnormality of these extremals.


title: Fabrication and characterization of suspended microstructures of ultrananocrystalline diamond

abstract: The fabrication of various suspended microstructures made from ultra-nanocrystalline diamond was performed without sacrificial layer. Arrays of bridges and cantilevers of various dimensions were successfully fabricated by a well-controlled process. Mechanical characterizations and finite element calculations allowed us to estimate both the Young's modulus as well as the residual stress built in the ultra-nanocrystalline diamond film.


title: The geometry of mixed-Euclidean metrics on symmetric positive definite matrices

abstract: Several Riemannian metrics and families of Riemannian metrics were defined on the manifold of Symmetric Positive Definite (SPD) matrices. Firstly, we formalize a common general process to define families of metrics: the principle of deformed metrics. We relate the recently introduced family of alpha-Procrustes metrics to the general class of mean kernel metrics by providing a sufficient condition under which elements of the former belongs to the latter. Secondly, we focus on the principle of balanced bilinear forms that we recently introduced. We give a new sufficient condition under which the balanced bilinear form is a metric. It allows us to introduce the Mixed-Euclidean (ME) metrics which generalize the Mixed-Power-Euclidean (MPE) metrics. We unveal their link with the (u, v)-divergences and the ($\alpha$, $\beta$)-divergences of information geometry and we provide an explicit formula of the Riemann curvature tensor. We show that the sectional curvature of all ME metrics can take negative values and we show experimentally that the sectional curvature of all MPE metrics but the log-Euclidean, power-Euclidean and power-affine metrics can take positive values.