title: M Subdwarf Research. II. Atmospheric Parameters and Kinematics

abstract: Applying the revised M subdwarf classification criteria discussed in Paper I to LAMOST DR7, combining the M subdwarf sample from Savcheva et al, a new M subdwarf sample was constructed for further study. The atmospheric parameters for each object were derived fitting with the PHOENIX grid, combining with Gaia DR2, the relationship between the gravity and metallicity were explored according to the locus both in the color-absolute magnitude diagram and the reduced proper motion diagram. Objects that have both the largest gravity and the lowest metallicity are located away from the main-sequence cloud and may be considered as the intrinsic M subdwarfs, which can be classified as luminosity class VI. Another group of objects whose spectra show typical M subdwarf characters have lower gravity and relatively moderate metal deficiency and occupy part of the ordinary M dwarf region in both diagrams. The Galactic U , V , W space velocity components and their dispersion show that the local Galactic halo population sampled in the solar neighborhood is represented by objects of high gravity and an inconspicuous bimodal metallicity distribution, with a fraction of prograde orbits. The other M subdwarfs seem to partly belong to the thick disk component with a significant fraction of thin disk moderately metal-poor objects intricately mixed with them. However, the selection effects, especially the favored anti-center direction of investigation in the LAMOST sub-sample, but also contamination by multiplicity and parameter coupling could play important roles and need to be further investigated.


title: Maximizing the Minimum and Maximum Forcing Numbers of Perfect Matchings of Graphs

abstract: Let $G$ be a simple graph with $2n$ vertices and a perfect matching. The forcing number $f(G,M)$ of a perfect matching $M$ of $G$ is the smallest cardinality of a subset of $M$ that is contained in no other perfect matching of $G$. Among all perfect matchings $M$ of $G$, the minimum and maximum values of $f(G,M)$ are called the minimum and maximum forcing numbers of $G$, denoted by $f(G)$ and $F(G)$, respectively. Then $f(G)\leq F(G)\leq n-1$. Che and Chen (2011) proposed an open problem: how to characterize the graphs $G$ with $f(G)=n-1$. Later they showed that for a bipartite graph $G$, $f(G)=n-1$ if and only if $G$ is a complete bipartite graph $K_{n,n}$. In this paper, we completely solve the problem of Che and Chen, and show that $f(G)=n-1$ if and only if $G$ is a complete multipartite graph or a graph obtained from complete bipartite graph $K_{n,n}$ by adding arbitrary edges in the same partite set. For all graphs $G$ with $F(G)=n-1$, we prove that the forcing spectrum of each such graph $G$ forms an integer interval by matching 2-switches and the minimum forcing numbers of all such graphs $G$ form an integer interval from $\lfloor\frac{n}{2}\rfloor$ to $n-1$.


title: Exploring Global Information for Session-based Recommendation

abstract: Session-based recommendation (SBR) is a challenging task, which aims at recommending items based on anonymous behavior sequences. Most existing SBR studies model the user preferences based only on the current session while neglecting the item-transition information from the other sessions, which suffer from the inability of modeling the complicated item-transition pattern. To address the limitations, we introduce global item-transition information to strength the modeling of the dynamic item-transition. For fully exploiting the global item-transition information, two ways of exploring global information for SBR are studied in this work. Specifically, we first propose a basic GNN-based framework (BGNN), which solely uses session-level item-transition information on session graph. Based on BGNN, we propose a novel approach, called Session-based Recommendation with Global Information (SRGI), which infers the user preferences via fully exploring global item-transitions over all sessions from two different perspectives: (i) Fusion-based Model (SRGI-FM), which recursively incorporates the neighbor embeddings of each node on global graph into the learning process of session level item representation; and (ii) Constrained-based Model (SRGI-CM), which treats the global-level item-transition information as a constraint to ensure the learned item embeddings are consistent with the global item-transition. Extensive experiments conducted on three popular benchmark datasets demonstrate that both SRGI-FM and SRGI-CM outperform the state-of-the-art methods consistently.


title: FLAVA: Find, Localize, Adjust and Verify to Annotate LiDAR-Based Point Clouds

abstract: Recent years have witnessed the rapid progress of perception algorithms on top of LiDAR, a widely adopted sensor for autonomous driving systems. These LiDAR-based solutions are typically data hungry, requiring a large amount of data to be labeled for training and evaluation. However, annotating this kind of data is very challenging due to the sparsity and irregularity of point clouds and more complex interaction involved in this procedure. To tackle this problem, we propose FLAVA, a systematic approach to minimizing human interaction in the annotation process. Specifically, we divide the annotation pipeline into four parts: find, localize, adjust and verify. In addition, we carefully design the UI for different stages of the annotation procedure, thus keeping the annotators to focus on the aspects that are most important to each stage. Furthermore, our system also greatly reduces the amount of interaction by introducing a light-weight yet effective mechanism to propagate the annotation results. Experimental results show that our method can remarkably accelerate the procedure and improve the annotation quality.


title: Escherization with Generalized Distance Functions Focusing on Local Structural Similarity

abstract: The Escherization problem involves finding a closed figure that tiles the plane that is most similar to a given goal figure. In Koizumi and Sugihara's formulation of the Escherization problem, the tile and goal figures are represented as $n$-point polygons where the similarity between them is measured based on the difference in the positions between the corresponding points. This paper presents alternative similarity measures (distance functions) suitable for this problem. The proposed distance functions focus on the similarity of local structures in several different manners. The designed distance functions are incorporated into a recently developed framework of the exhaustive search of the templates for the Escherization problem. Efficient exhaustive and incomplete search algorithms for the formulated problems are also developed to obtain results within a reasonable computation time. Experimental results showed that the proposed algorithms found satisfactory tile shapes for fairly complicated goal figures in a reasonable computation time.


title: Fourier transform decay of distributions in Hardy-Morrey spaces

abstract: In this paper we establish decay estimates for Fourier transform on Hardy-Morrey spaces and its localizable version. Our work include some aspects to these spaces linked up with pointwise Fourier estimates, in particular a natural approach on cancellation moment conditions. As application, we discuss the optimality for continuity of Fourier multipliers and pseudodifferential operators in Hardy-Morrey spaces.


title: Flight Sensor Data and Beamforming based Integrated UAV Tracking with Channel Estimation using Gaussian Process Regression

abstract: With explosively increasing demands for unmanned aerial vehicle (UAV) applications, reliable link acquisition for serving UAVs is required. Considering the dynamic characteristics of UAV, it is hugely challenging to persist a reliable link without beam misalignment. In this paper, we propose a flight sensor data and beamforming signal based integrated UAV tracking scheme to deal with this problem. The proposed scheme provides a compatible integrated system considering the practical specification of the flight sensor data and the beamforming pilot signal. The UAV position tracking is comprised of two steps: 1) UAV position prediction by the flight sensor data and 2) position update with the beamforming signal using Gaussian process regression (GPR) method, which is a nonparametric machine learning. The flight sensor data can assist ground station (GS) or UAV nodes in designing the precoding and the receive beamforming matrix with drastically reduced overheads. The beamforming signal can accomplish high beamforming gain to be maintained even when the flight sensor data is absent. Therefore, the proposed scheme can support the moving target continuously by utilizing these two signals. The simulation results are provided to confirm that the proposed scheme outperforms other conventional beam tracking schemes. We also derive 3-dimensional (3D) beamforming gain and spectral efficiency (SE) from the mean absolute error (MAE) of the angular value estimation, which can be used as beamforming performance metrics of the data transmission link in advance.


title: Wall slip in primitive chain network simulations of shear startup of entangled polymers and its effect on the shear stress undershoot

abstract: In some recent experiments on entangled polymers of stress growth in startup of fast shear flows an undershoot in the shear stress is observed following the overshoot, i.e., before approaching the steady state. Whereas tumbling of the entangled chain was proposed to be at its origin, here we investigate another possible cause for the stress undershoot, i.e., slippage at the interface between polymer and solid wall. To this end, we extend the primitive chain network model to include slip at the interface between entangled polymeric liquids and solid walls with grafted polymers. We determine the slip velocity at the wall, and the shear rate in the bulk, by imposing that the shear stress in the bulk polymers is equal to that resulting from the polymers grafted at the wall. After confirming that the predicted results for the steady state are reasonable, we examine the transient behavior. The simulations confirm that slippage weakens the magnitude of the stress overshoot, as reported earlier. The undershoot is also weakened, or even disappears, because of a reduced coherence in molecular tumbling. In other words, the disentanglement between grafted and bulk chains, occurring throughout the stress overshoot region, does not contribute to the stress undershoot.


title: Unruh-DeWitt Detector Differentiation of Black Holes and Exotic Compact Objects

abstract: We study the response of a static Unruh-DeWitt detector outside an exotic compact object (ECO) with a general reflective boundary condition in 3+1 dimensions. The horizonless ECO, whose boundary is extremely close to the would-be event horizon, acts as a black hole mimicker. We find that the response rate is notably distinct from the black hole case, even when the ECO boundary is perfectly absorbing. For a (partially) reflective ECO boundary, we find resonance structures in the response rate that depend on the different locations of the ECO boundary and those of the detector. We provide a detailed analysis in connection with the ECO's vacuum mode structure and transfer function.


title: Survey and Open Problems in Privacy Preserving Knowledge Graph: Merging, Query, Representation, Completion and Applications

abstract: Knowledge Graph (KG) has attracted more and more companies' attention for its ability to connect different types of data in meaningful ways and support rich data services. However, the data isolation problem limits the performance of KG and prevents its further development. That is, multiple parties have their own KGs but they cannot share with each other due to regulation or competition reasons. Therefore, how to conduct privacy preserving KG becomes an important research question to answer. That is, multiple parties conduct KG related tasks collaboratively on the basis of protecting the privacy of multiple KGs. To date, there is few work on solving the above KG isolation problem. In this paper, to fill this gap, we summarize the open problems for privacy preserving KG in data isolation setting and propose possible solutions for them. Specifically, we summarize the open problems in privacy preserving KG from four aspects, i.e., merging, query, representation, and completion. We present these problems in details and propose possible technical solutions for them. Moreover, we present three privacy preserving KG-aware applications and simply describe how can our proposed techniques be applied into these applications.


title: Rational curves on K3 surfaces of small genus

abstract: Let $\mathfrak B_g$ denote the moduli space of primitively polarized $K3$ surfaces $(S,H)$ of genus $g$ over $\mathbb C$. It is well-known that $\mathfrak B_g$ is irreducible and that there are only finitely many rational curves in $|H|$ for any primitively polarized $K3$ surface $(S,H)$. So we can ask the question of finding the monodromy group of such curves. The case of $g=2$ essentially follows from the results of Harris \cite{Ha} to be the full symmetric group $S_{324}$, here we solve the case $g=3$ and $4$.


title: Impact of COVID-19 on Mobility and Electric Vehicle Charging Load

abstract: The COVID-19 pandemic has depressed overall mobility across the country. The changes seen reflect responses to new COVID-19 cases, local health guidelines, and seasonality, making the relationship between mobility and COVID-19 unique from region to region. This paper presents a data-driven case study of electric vehicle (EV) charging and mobility in the wake of COVID-19. The study shows that the number of EV charging sessions and total energy consumed per day dropped by 40% immediately after the arrival of the first COVID-19 case in Utah. By contrast, the energy consumed per charging session fell by just 8% over the same periods and the distribution of session start and end times remained consistent throughout the year. While EV mobility dropped more dramatically than total vehicle mobility during the first wave of COVID-19 cases, and returned more slowly, both returned to stable levels near their mean values by September 2020, despite a dramatic third wave in new infections.


title: A three-dimensional thermal model of the human cochlea for magnetic cochlear implant surgery

abstract: In traditional cochlear implant surgery, physical trauma may occur during electrode array insertion. Magnetic guidance of the electrode array has been proposed to mitigate this medical complication. After insertion, the guiding magnet attached to the tip of the electrode array must be detached via a heating process and removed. This heating process may, however, cause thermal trauma within the cochlea. In this study, a validated three-dimensional finite element heat transfer model of the human cochlea is applied to perform an intracochlear thermal analysis necessary to ensure the safety of the magnet removal phase. Specifically, the maximum safe input power density to detach the magnet is determined as a function of the boundary conditions, heating duration, cochlea size, implant electrode array radius and insertion depth, magnet size, and cochlear fluid. A dimensional analysis and numerical simulations reveal that the maximum safe input power density increases with increasing cochlea size and the radius of the electrode array, whereas it decreases with increasing electrode array insertion depth and magnet size. The best cochlear fluids from the thermal perspective are perilymph and a soap solution. Even for the worst case scenario in which the cochlear walls are assumed to be adiabatic except at the round window, the maximum safe input power density is larger than that required to melt 1 $\rm{mm^3}$ of paraffin bonding the magnet to the implant electrode array. By combining the outcome of this work with other aspects of the design of the magnetic insertion process, namely the magnetic guidance procedure and medical requirements, it will be possible to implement a thermally safe patient-specific surgical procedure.

categories: physics.app-ph physics.bio-ph physics.flu-dyn physics.med-ph


title: A Deep Search for Emission From "Rock Comet" (3200) Phaethon At 1 AU

abstract: We present a deep imaging and spectroscopic search for emission from (3200) Phaethon, a large near-Earth asteroid that appears to be the parent of the strong Geminid meteoroid stream, using the 4.3 m Lowell Discovery Telescope. Observations were conducted on 2017 December 14-18 when Phaethon passed only 0.07 au from the Earth. We determine the $3\sigma$ upper level of dust and CN production rates to be 0.007-0.2 $\mathrm{kg~s^{-1}}$ and $2.3\times10^{22}~\mathrm{molecule~s^{-1}}$ through narrowband imaging. A search in broadband images taken through the SDSS $r'$ filter shows no 100-m-class fragments in Phaethon's vicinity. A deeper, but star-contaminated search also shows no sign of fragments down to 15 m. Optical spectroscopy of Phaethon and comet C/2017 O1 (ASASSN) as comparison confirms the absence of cometary emission lines from Phaethon and yields $3\sigma$ upper levels of CN, C$_2$ and C$_3$ of $\sim10^{24}$-$10^{25} \mathrm{molecule~s^{-1}}$, 2 orders of magnitude higher than the CN constraint placed by narrowband imaging, due to the much narrower on-sky aperture of the spectrographic slit. We show that narrowband imaging could provide an efficient way to look for weak gas emission from near-extinct bodies near the Earth, though these observations require careful interpretation. Assuming Phaethon's behavior is unchanged, our analysis shows that the DESTINY$^+$ mission, currently planning to explore Phaethon in 2026, may not be able to directly detect a gas coma.


title: ConvTransformer: A Convolutional Transformer Network for Video Frame Synthesis

abstract: Deep Convolutional Neural Networks (CNNs) are powerful models that have achieved excellent performance on difficult computer vision tasks. Although CNNs perform well whenever large labeled training samples are available, they work badly on video frame synthesis due to objects deforming and moving, scene lighting changes, and cameras moving in video sequence. In this paper, we present a novel and general end-to-end architecture, called convolutional Transformer or ConvTransformer, for video frame sequence learning and video frame synthesis. The core ingredient of ConvTransformer is the proposed attention layer, i.e., multi-head convolutional self-attention layer, that learns the sequential dependence of video sequence. ConvTransformer uses an encoder, built upon multi-head convolutional self-attention layer, to encode the sequential dependence between the input frames, and then a decoder decodes the long-term dependence between the target synthesized frames and the input frames. Experiments on video future frame extrapolation task show ConvTransformer to be superior in quality while being more parallelizable to recent approaches built upon convolutional LSTM (ConvLSTM). To the best of our knowledge, this is the first time that ConvTransformer architecture is proposed and applied to video frame synthesis.


title: Gaussian Process Modeling $\it{Fermi}$-LAT $\gamma$-ray Blazar Variability: A Sample of Blazars with $\gamma$-ray Quasi-periodicities

abstract: Blazar variability may be driven by stochastic processes. On the other hand, quasi-periodic oscillation (QPO) behaviors are recently reported to be detected in $\it{Fermi}$-LAT data of blazars. However, the significances of these QPO signals given by traditional Fourier-like methods are still questioned. We analyze $\gamma$-ray light curves of the QPO blazars with two Gaussian process methods, CARMA and $\it{celerite}$, to examine the appropriateness of Gaussian processes for characterizing $\gamma$-ray light curves of blazars and the existence of the reported QPOs. We collect a sample of 27 blazars with possible $\gamma$-ray periodicity and generate their $\sim11$ years $\it{Fermi}$-LAT light curves. We apply the Gaussian process models to the $\gamma$-ray light curves, and build their intrinsic power spectral densities (PSDs). The results show that in general the $\gamma$-ray light curves can be characterized by CARMA and $\textit{celerite}$ models, indicating that $\gamma$-ray variabilities of blazars are essentially Gaussian processes. The resulting PSDs are generally the red noise shapes with slopes between $-0.6$ and $-1.7$. Possible evidence for the $\gamma$-ray QPOs in PKS 0537$-$441 and PG 1553$+$113 are found in the Gaussian process modelings.


title: A Survey on Deep Learning Based Point-Of-Interest (POI) Recommendations

abstract: Location-based Social Networks (LBSNs) enable users to socialize with friends and acquaintances by sharing their check-ins, opinions, photos, and reviews. Huge volume of data generated from LBSNs opens up a new avenue of research that gives birth to a new sub-field of recommendation systems, known as Point-of-Interest (POI) recommendation. A POI recommendation technique essentially exploits users' historical check-ins and other multi-modal information such as POI attributes and friendship network, to recommend the next set of POIs suitable for a user. A plethora of earlier works focused on traditional machine learning techniques by using hand-crafted features from the dataset. With the recent surge of deep learning research, we have witnessed a large variety of POI recommendation works utilizing different deep learning paradigms. These techniques largely vary in problem formulations, proposed techniques, used datasets, and features, etc. To the best of our knowledge, this work is the first comprehensive survey of all major deep learning-based POI recommendation works. Our work categorizes and critically analyzes the recent POI recommendation works based on different deep learning paradigms and other relevant features. This review can be considered a cookbook for researchers or practitioners working in the area of POI recommendation.


title: Targeted Self Supervision for Classification on a Small COVID-19 CT Scan Dataset

abstract: Traditionally, convolutional neural networks need large amounts of data labelled by humans to train. Self supervision has been proposed as a method of dealing with small amounts of labelled data. The aim of this study is to determine whether self supervision can increase classification performance on a small COVID-19 CT scan dataset. This study also aims to determine whether the proposed self supervision strategy, targeted self supervision, is a viable option for a COVID-19 imaging dataset. A total of 10 experiments are run comparing the classification performance of the proposed method of self supervision with different amounts of data. The experiments run with the proposed self supervision strategy perform significantly better than their non-self supervised counterparts. We get almost 8% increase in accuracy with full self supervision when compared to no self supervision. The results suggest that self supervision can improve classification performance on a small COVID-19 CT scan dataset. Code for targeted self supervision can be found at this link: https://github.com/Mewtwo/Targeted-Self-Supervision/tree/main/COVID-CT


title: MobileDepth: Efficient Monocular Depth Prediction on Mobile Devices

abstract: Depth prediction is fundamental for many useful applications on computer vision and robotic systems. On mobile phones, the performance of some useful applications such as augmented reality, autofocus and so on could be enhanced by accurate depth prediction. In this work, an efficient fully convolutional network architecture for depth prediction has been proposed, which uses RegNetY 06 as the encoder and split-concatenate shuffle blocks as decoder. At the same time, an appropriate combination of data augmentation, hyper-parameters and loss functions to efficiently train the lightweight network has been provided. Also, an Android application has been developed which can load CNN models to predict depth map by the monocular images captured from the mobile camera and evaluate the average latency and frame per second of the models. As a result, the network achieves 82.7% {\delta}1 accuracy on NYU Depth v2 dataset and at the same time, have only 62ms latency on ARM A76 CPUs so that it can predict the depth map from the mobile camera in real-time.


title: Action Duration Prediction for Segment-Level Alignment of Weakly-Labeled Videos

abstract: This paper focuses on weakly-supervised action alignment, where only the ordered sequence of video-level actions is available for training. We propose a novel Duration Network, which captures a short temporal window of the video and learns to predict the remaining duration of a given action at any point in time with a level of granularity based on the type of that action. Further, we introduce a Segment-Level Beam Search to obtain the best alignment, that maximizes our posterior probability. Segment-Level Beam Search efficiently aligns actions by considering only a selected set of frames that have more confident predictions. The experimental results show that our alignments for long videos are more robust than existing models. Moreover, the proposed method achieves state of the art results in certain cases on the popular Breakfast and Hollywood Extended datasets.


title: Estimating fast transient detection pipeline efficiencies at UTMOST via real-time injection of mock FRBs

abstract: Dedicated surveys using different detection pipelines are being carried out at multiple observatories to find more Fast Radio Bursts (FRBs). Understanding the efficiency of detection algorithms and the survey completeness function is important to enable unbiased estimation of the underlying FRB population properties. One method to achieve end-to-end testing of the system is by injecting mock FRBs in the live data-stream and searching for them blindly. Mock FRB injection is particularly effective for machine-learning-based classifiers, for which analytic characterisation is impractical. We describe a first-of-its-kind implementation of a real-time mock FRB injection system at the upgraded Molonglo Observatory Synthesis Telescope (UTMOST) and present our results for a set of 20,000 mock FRB injections. The injections have yielded clear insight into the detection efficiencies and have provided a survey completeness function for pulse width, fluence and DM. Mock FRBs are recovered with uniform efficiency over the full range of injected DMs, however the recovery fraction is found to be a strong function of the width and Signal-to-Noise (SNR). For low widths ($\lesssim 20$ ms) and high SNR ($\gtrsim$ 9) the recovery is highly effective with recovery fractions exceeding 90%. We find that the presence of radio frequency interference causes the recovered SNR values to be systematically lower by up to 20% compared to the injected values. We find that wider FRBs become increasingly hard to recover for the machine-learning-based classifier employed at UTMOST. We encourage other observatories to implement live injection set-ups for similar testing of their surveys.


title: Dynamics of Sunspot Shock Waves in the Chromosphere and Transition Region

abstract: We study the dynamics of shock waves observed in the umbra of a sunspot using the spectroscopic observations from the Interface Region Imaging Spectrometer (IRIS). The presence of the shock significantly deforms the shape of the spectral lines of Mg II , C II , and Si IV . We found that C II 1335.66 {\AA} and Si IV 1393.75 {\AA} show double-peaked profiles that change to a single peak later on. However, the Mg II h 2803.53 {\AA} line first shows flat-top profiles that change into double-peaked followed by the single peak. To study the shock dynamics, we isolate the shock component from the spectra by fitting two Gaussians. We find that the lifetime of the shock is largest in Mg II h 2803.53 {\AA} line. Moreover, the plasma motion shows both acceleration and deceleration phase of the shock. Yet, in C II 1335.66 {\AA} and Si IV 1393.75 {\AA}, only deceleration phase is observed. We observe a strong correlation between the largest blueshift of the shock and deceleration for all three spectral lines. We find a positive (negative) correlation between intensities contributed due to the shocks in Mg II and C II (Si IV ). This is suggestive that the shocks are first amplified in C II , followed by a decline in the height range corresponding to Si IV . These results may indicate the dissipation of shocks above the formation height of C II , and the shocks may have important roles in the dynamics of the upper chromosphere and transition region above sunspots.


title: Simulation of the orbit and spin period evolution of the double pulsars PSR J0737-3039 from their birth to coalescence induced by the gravitational wave radiation

abstract: The complete orbital and sin period evolutions of the double neutron star (NS) system PSR J0737-3039 are simulated from the birth to coalescence, which include the two observed radio pulsars classified as primary NS PSR J0737-3039A and companion NS PSR J0737-3039B. By employing the characteristic age of PSR J0737-3039B to constrain the true age of the double pulsar system, the initial orbital period and initial binary separation are obtained as 2.89 hrs and $1.44 \times 10^{6}$ km, respectively, and the coalescence age or the lifetime from the birth to merger of PSR J0737-3039 is obtained to be $1.38 \times 10^{8}$ yr. At the last minute of coalescence, corresponding to that the gravitational wave frequency changes from 20 Hz to 1180 Hz, we present the binary separation of PSR J0737-3039 to be from 442 km to 30 km, while the spin periods of PSR J0737-3039A and PSR J0737-3039B are 27.10 ms and 4.63 s, respectively. From the standard radio pulsar emission model, before the system merged, the primary NS could still be observed by the radio telescope, while the companion NS has crossed the death line in the pulsar magnetic-field versus period ($B-P$) diagram which is usually considered to cease the life as a pulsar. It is for the first time that the whole life evolutionary simulation of the orbit and spin periods for double NS system is presented, which provides the useful information for observing the primary NS at the coalescence stage.


title: Faint Quasars Live in the Same Number Density Environments as Lyman Break Galaxies at z ~ 4

abstract: Characterizing high-z quasar environments is key to understanding the co-evolution of quasars and the surrounding galaxies. To restrict their global picture, we statistically examine the g-dropout galaxy overdensity distribution around 570 faint quasar candidates at z ~ 4, based on the Hyper Suprime-Cam Subaru Strategic Program survey. We compare the overdensity significances of g-dropout galaxies around the quasars with those around g-dropout galaxies, and find no significant difference between their distributions. A total of 4 (22) out of the 570 faint quasars, 0.7_{-0.4}^{+0.4} (3.9_{-0.8}^{+0.8}) %, are found to be associated with the > 4 sigma overdense regions within an angular separation of 1.8 (3.0) arcmin, which is the typical size of protoclusters at this epoch. This is similar to the fraction of g-dropout galaxies associated with the > 4 sigma overdense regions. This result is consistent with our previous work that 1.3_{-0.9}^{+0.9} % and 2.0_{-1.1}^{+1.1} % of luminous quasars detected in the Sloan Digital Sky Survey exist in the > 4 sigma overdense regions within 1.8 and 3.0 arcmin separations, respectively. Therefore, we suggest that the galaxy number densities around quasars are independent of their luminosity, and most quasars do not preferentially appear in the richest protocluster regions at z ~ 4. The lack of an apparent positive correlation between the quasars and the protoclusters implies that: i) the gas-rich major merger rate is relatively low in the protocluster regions, ii) most high-z quasars may appear through secular processes, or iii) some dust-obscured quasars exist in the protocluster regions.


title: Detecting systematic anomalies affecting systems when inputs are stationary time series

abstract: We develop an anomaly-detection method when systematic anomalies, possibly statistically very similar to genuine inputs, are affecting control systems at the input and/or output stages. The method allows anomaly-free inputs (i.e., those before contamination) to originate from a wide class of random sequences, thus opening up possibilities for diverse applications. To illustrate how the method works on data, and how to interpret its results and make decisions, we analyze several actual time series, which are originally non-stationary but in the process of analysis are converted into stationary. As a further illustration, we provide a controlled experiment with anomaly-free inputs following an ARMA time series model under various contamination scenarios.


title: Deep unfolding-based output feedback control design for linear systems with input saturation

abstract: In this paper, we propose a deep unfolding-based framework for the output feedback control of systems with input saturation. Although saturation commonly arises in several practical control systems, there is still a scarce of effective design methodologies that can directly deal with the severe non-linearity of the saturation operator. In this paper, we aim to design an anti-windup controller for enlarging the region of stability of the closed-loop system by learning from the numerical simulations of the closed-loop system. The data-driven framework we propose in this paper is based on a deep-learning technique called Neural Ordinary Differential Equations. Within our framework, we first obtain a candidate controller by using the deep-learning technique, which is then tested by the existing theoretical results already established in the literature, thereby avoiding the computational challenge in the conventional design methodologies as well as theoretically guaranteeing the performance of the system. Our numerical simulation shows that the proposed framework can significantly outperform a conventional design methodology based on linear matrix inequalities.


title: Cooperative Activity Detection: Sourced and Unsourced Massive Random Access Paradigms

abstract: This paper investigates the issue of cooperative activity detection for grant-free random access in the sixth-generation (6G) cell-free wireless networks with sourced and unsourced paradigms. First, we propose a cooperative framework for solving the problem of device activity detection in sourced random access. In particular, multiple access points (APs) cooperatively detect the device activity via exchanging low-dimensional intermediate information with their neighbors. This is enabled by the proposed covariance-based algorithm via exploiting both the sparsity-promoting and similarity-promoting terms of the device state vectors among neighboring APs. A decentralized approximate separating approach is introduced based on the forward-backward splitting strategy for addressing the formulated problem. Then, the proposed activity detection algorithm is adopted as a decoder of cooperative unsourced random access, where the multiple APs cooperatively detect the list of transmitted messages regardless of the identity of the transmitting devices. Finally, we provide sufficient conditions on the step sizes that ensure the convergence of the proposed algorithm in the sense of Bregman divergence. Simulation results show that the proposed algorithm is efficient for addressing both sourced and unsourced massive random access problems, while requires a shorter signature sequence and accommodates a significantly larger number of active devices with a reasonable antenna array size, compared with the state-of-art algorithms.


title: Dual-gated graphene devices for near-field nano-imaging

abstract: Graphene-based heterostructures display a variety of phenomena that are strongly tunable by electrostatic local gates. Monolayer graphene (MLG) exhibits tunable surface plasmon polaritons, as revealed by scanning nano-infrared experiments. In bilayer graphene (BLG), an electronic gap is induced by a perpendicular displacement field. Gapped BLG is predicted to display unusual effects such as plasmon amplification and domain wall plasmons with significantly larger lifetime than MLG. Furthermore, a variety of correlated electronic phases highly sensitive to displacement fields have been observed in twisted graphene structures. However, applying perpendicular displacement fields in nano-infrared experiments has only recently become possible (Ref. 1). In this work, we fully characterize two approaches to realizing nano-optics compatible top-gates: bilayer $\text{MoS}_2$ and MLG. We perform nano-infrared imaging on both types of structures and evaluate their strengths and weaknesses. Our work paves the way for comprehensive near-field experiments of correlated phenomena and plasmonic effects in graphene-based heterostructures.


title: Super-clustering of consecutive numbers in $p$-shifted random permutations

abstract: Let $A^{(n)}_{l;k}\subset S_n$ denote the event that the set of $l$ consecutive numbers $\{k,k+1,\cdots, k+l-1\}$ appear in a set of $l$ consecutive positions. Let $p=\{p_j\}_{j=1}^\infty$ be a distribution on $\mathbb{N}$ with $p_j>0$. Let $P_n$ denote the probability measure on $S_n$ corresponding to the $p$-shifted random permutation. Our main result, under the additional assumption that $\{p_j\}_{j=1}^\infty$ is non-increasing, is that $$ \begin{aligned} &\lim_{l\to\infty}\lim_{n\to\infty}P_n(A^{(n )}_{l,k})=\big(\prod_{j=1}^{k-1}\sum_{i=1}^jp_i\big) \big(\prod_{j=1}^\infty\sum_{i=1}^jp_i\big), \end{aligned} $$ and that if $\lim_{n\to\infty}\min(k_n,n-k_n)=\infty$, then $$ \begin{aligned} &\lim_{l\to\infty}\lim_{n\to\infty}P_n(A^{(n )}_{l,k_n})= \big(\prod_{j=1}^\infty\sum_{i=1}^jp_i\big)^2. \end{aligned} $$ In particular these limits are positive if and only if $\sum_{j=1}^\infty jp_j<\infty$. We say that super-clustering occurs when the limits are positive. We also give a new characterization of the class of $p$-shifted probability distributions on $S_\infty$.


title: The partition function of the four-vertex model in inhomogeneous external field and trace statistics

abstract: The exactly solvable four-vertex model with the fixed boundary conditions in the presence of inhomogeneous linearly growing external field is considered. The partition function of the model is calculated and represented in the determinantal form. The established connection with the boxed plane partitions allows us to calculate the generating function of plane partitions with the fixed sums of their diagonals. The obtained results are another example of the connection of integrable models with the enumerative combinatorics.

categories: cond-mat.stat-mech math-ph math.MP nlin.SI


title: Discriminative Localized Sparse Representations for Breast Cancer Screening

abstract: Breast cancer is the most common cancer among women both in developed and developing countries. Early detection and diagnosis of breast cancer may reduce its mortality and improve the quality of life. Computer-aided detection (CADx) and computer-aided diagnosis (CAD) techniques have shown promise for reducing the burden of human expert reading and improve the accuracy and reproducibility of results. Sparse analysis techniques have produced relevant results for representing and recognizing imaging patterns. In this work we propose a method for Label Consistent Spatially Localized Ensemble Sparse Analysis (LC-SLESA). In this work we apply dictionary learning to our block based sparse analysis method to classify breast lesions as benign or malignant. The performance of our method in conjunction with LC-KSVD dictionary learning is evaluated using 10-, 20-, and 30-fold cross validation on the MIAS dataset. Our results indicate that the proposed sparse analyses may be a useful component for breast cancer screening applications.


title: CLIPPER: A Graph-Theoretic Framework for Robust Data Association

abstract: We present CLIPPER (Consistent LInking, Pruning, and Pairwise Error Rectification), a framework for robust data association in the presence of noise and outliers. We formulate the problem in a graph-theoretic framework using the notion of geometric consistency. State-of-the-art techniques that use this framework utilize either combinatorial optimization techniques that do not scale well to large-sized problems, or use heuristic approximations that yield low accuracy in high-noise, high-outlier regimes. In contrast, CLIPPER uses a relaxation of the combinatorial problem and returns solutions that are guaranteed to correspond to the optima of the original problem. Low time complexity is achieved with an efficient projected gradient ascent approach. Experiments indicate that CLIPPER maintains a consistently low runtime of 15 ms where exact methods can require up to 24 s at their peak, even on small-sized problems with 200 associations. When evaluated on noisy point cloud registration problems, CLIPPER achieves 100% precision and 98% recall in 90% outlier regimes while competing algorithms begin degrading by 70% outliers. In an instance of associating noisy points of the Stanford Bunny with 990 outlier associations and only 10 inlier associations, CLIPPER successfully returns 8 inlier associations with 100% precision in 138 ms. Code is available at https://mit-acl.github.io/clipper.


title: Vacancies and dopants in two-dimensional tin monoxide: An ab initio study

abstract: Layered tin monoxide (SnO) offers an exciting two-dimensional (2D) semiconducting system with great technological potential for next-generation electronics and photocatalytic applications. Using a combination of first-principles simulations and strain field analysis, this study investigates the structural dynamics of oxygen (O) vacancies in monolayer SnO and their functionalization by complementary lightweight dopants, namely C, Si, N, P, S, F, Cl, H and H$_{2}$. Our results show that O vacancies are the dominant native defect under Sn-rich growth conditions with active diffusion characteristics that are comparable to that of graphene vacancies. Depending on the choice of substitutional species and its concentration within the material, significant opportunities are revealed in the doped-SnO system for facilitating $n$/$p$-type tendencies, work function reduction, and metallization of the monolayer. N and F dopants are found to demonstrate superior mechanical compatibility with the host lattice, with F being especially likely to take part in substitution and lead to degenerately doped phases with high open-air stability. The findings reported here suggest that post-growth filling of O vacancies in Sn-rich conditions presents a viable channel for doping 2D tin monoxide, opening up new avenues in harnessing defect-engineered SnO nanostructures for emergent technologies.


title: Exchange interactions and magnetic force theorem

abstract: We critically reexamine the problem of interatomic exchange interactions, which describe the total energy change caused by infinitesimal rotations of spins near some equilibrium state. For the small variations, such interactions can be always related to the response function. However, the form of this relation can depend on additional approximations. Particularly, the commonly used magnetic force theorem (MFT) prescribes the linear relation between the exchange interactions and the response function, while the exact theory requires this dependence to be inverse. We explore the origin and consequences of these differences in the definition for the wide class of materials: ferromagnetic Ni, antiferromagnetic NiO, half-metallic CrO2, multiferroic HoMnO3, and layered magnets CrCl3 and CrI3. While in most of these cases, MFT produces quite reasonable results and can be rigorously justifies in the long wavelength and strong-coupling limits, the exact formulation appears to be more consistent, especially in dealing with two important issues, which typically arise in the theory of exchange interactions: (i) the treatment of the ligand states, and (ii) the choice of the suitable variable for the description of infinitesimal rotations of spins. Both issues can be efficiently resolved by employing the ideas of adiabatic spin dynamics supplemented with the exact expression for the exchange interactions. Particularly, we propose a simple "downfolding" procedure for the elimination of the ligand spins by transferring their effect to the interaction parameters between the localized spins. Furthermore, we argue that the rotations of spin moments are more suitable for the description of low-energy excitations, while the rotations of the whole magnetization matrix cause much stronger perturbation in the system of spins.


title: Virtues of Patience in Strategic Queuing Systems

abstract: We consider the problem of selfish agents in discrete-time queuing systems, where competitive queues try to get their packets served. In this model, a queue gets to send a packet each step to one of the servers, which will attempt to serve the oldest arriving packet, and unprocessed packets are returned to each queue. We model this as a repeated game where queues compete for the capacity of the servers, but where the state of the game evolves as the length of each queue varies, resulting in a highly dependent random process. Earlier work by the authors [EC'20] shows that with no-regret learners, the system needs twice the capacity as would be required in the coordinated setting to ensure queue lengths remain stable despite the selfish behavior of the queues. In this paper, we demonstrate that this way of evaluating outcomes is myopic: if more patient queues choose strategies that selfishly maximize their long-run success rate, stability can be ensured with just $\frac{e}{e-1}\approx 1.58$ times extra capacity, better than what is possible assuming the no-regret property. As these systems induce highly dependent processes, our analysis draws heavily on techniques from probability theory. Though these systems are random under any fixed policies by the queues, we show that, surprisingly, these systems have deterministic and explicit asymptotic behavior. We show that the asymptotic growth rates of queues can be written as a ratio of a submodular and modular function, which provides significant game-theoretic properties. Our equilibrium analysis then relies on a novel deformation argument towards a more analyzable solution that differs significantly from previous price of anarchy results. While the intermediate points will not be equilibria, this analytic structure will ensure that this deformation is monotonic along this continuous path.


title: 591 high velocity stars in the Galactic halo selected from LAMOST DR7 and Gaia DR2

abstract: In this paper, we report 591 high velocity star candidates (HiVelSCs) selected from over 10 million spectra of the data release seven (DR7) of the Large Sky Area Multi-object Fiber Spectroscopic Telescope and the second Gaia data release, with three-dimensional velocities in the Galactic rest-frame larger than 445 km/s. We show that at least 43 HiVelSCs are unbound to the Galaxy with escape probabilities larger than 50%, and this number decreases to eight if the possible parallax zero-point error is corrected. Most of these HiVelSCs are metal-poor and slightly alpha-enhanced inner halo stars. Only 14% of them have [Fe/H] > -1, which may be the metal-rich "in situ" stars in halo formed in the initial collapse of the Milky Way or metal-rich stars formed in the disk or bulge but kinematically heated. The low ratio of 14% implies that the bulk of stellar halo was formed from the accretion and tidal disruption of satellite galaxies. In addition, HiVelSCs on the retrograde orbits have a slightly lower metallicities on average compared with these on the prograde orbits, meanwhile metal-poor HiVelSCs with [Fe/H] < -1 have an even faster mean retrograde velocity compared with metal-rich HiVelSCs. To investigate the origins of HiVelSCs, we perform orbit integrations and divide them into four types, i.e., hypervelocity stars, hyper-runaway stars, runaway stars and fast halo stars. A catalog for these 591 HiVelSCs, including radial velocities, atmospheric parameters, Gaia astrometric parameters, spatial positions, and velocities, etc., is available at http://paperdata.chinavo.org/LYB/lamostdr7_gaiadr2_hvs_591.csv.


title: A note on a question of Markman

abstract: Let $\mathcal{F}$ be a vector bundle on a complex projective algebraic variety $X$. If $\mathcal{F}$ deforms along a first order deformation of $X$, its Mukai vector remains of Hodge type along this deformation. We prove an analogous statement for all polyvector fields, not only for those in $H^1(X,T_X)$ corresponding to deformations of the complex structure. This answers a question of Markman. We also explore a Lie theoretic analogue of the statement above.


title: Collaborative Storytelling with Large-scale Neural Language Models

abstract: Storytelling plays a central role in human socializing and entertainment. However, much of the research on automatic storytelling generation assumes that stories will be generated by an agent without any human interaction. In this paper, we introduce the task of collaborative storytelling, where an artificial intelligence agent and a person collaborate to create a unique story by taking turns adding to it. We present a collaborative storytelling system which works with a human storyteller to create a story by generating new utterances based on the story so far. We constructed the storytelling system by tuning a publicly-available large scale language model on a dataset of writing prompts and their accompanying fictional works. We identify generating sufficiently human-like utterances to be an important technical issue and propose a sample-and-rank approach to improve utterance quality. Quantitative evaluation shows that our approach outperforms a baseline, and we present qualitative evaluation of our system's capabilities.


title: Categorical models of Linear Logic with fixed points of formulas

abstract: We develop a denotational semantics of muLL, a version of propositional Linear Logic with least and greatest fixed points extending David Baelde's propositional muMALL with exponentials. Our general categorical setting is based on the notion of Seely category and on strong functors acting on them. We exhibit two simple instances of this setting. In the first one, which is based on the category of sets and relations, least and greatest fixed points are interpreted in the same way. In the second one, based on a category of sets equipped with a notion of totality (non-uniform totality spaces) and relations preserving them, least and greatest fixed points have distinct interpretations. This latter model shows that muLL enjoys a denotational form of normalization of proofs.


title: Inelastic collision dynamics of a single cold ion immersed in a Bose-Einstein condensate

abstract: We investigate inelastic collision dynamics of a single cold ion in a Bose-Einstein condensate. We observe rapid ion-atom-atom three-body recombination leading to formation of weakly bound molecular ions followed by secondary two-body molecule-atom collisions quenching the rovibrational states towards deeper binding energies. In contrast to previous studies exploiting hybrid ion traps, we work in an effectively field-free environment and generate a free low-energy ionic impurity directly from the atomic ensemble via Rydberg excitation and ionization. This allows us to implement an energy-resolved field-dissociation technique to trace the relaxation dynamics of the recombination products. Our observations are in good agreement with numerical simulations based on Langevin capture dynamics and provide complementary means to study stability and reaction dynamics of ionic impurities in ultracold quantum gases.


title: Portable Magnetic Particle Spectrometer (MPS) for Future Rapid and Wash-free Bioassays

abstract: Nowadays, there is an increasing demand for more accessible routine diagnostics for patients with respect to high accuracy, ease of use, and low cost. However, the quantitative and high accuracy bioassays in large hospitals and laboratories usually require trained technicians and equipment that is both bulky and expensive. In addition, the multi-step bioassays and long turnaround time could severely affect the disease surveillance and control especially in pandemics such as influenza and COVID-19. In view of this, a portable, quantitative bioassay device will be valuable in regions with scarce medical resources and help relieve burden on local healthcare systems. Herein, we introduce the MagiCoil diagnostic device, an inexpensive, portable, quantitative and rapid bioassay platform based on magnetic particle spectrometer (MPS) technique. MPS detects the dynamic magnetic responses of magnetic nanoparticles (MNPs) and uses the harmonics from oscillating MNPs as metrics for sensitive and quantitative bioassays. This device does not require trained technicians to operate and employs a fully automatic, one-step, wash-free assay with user friendly smartphone interface. Using a streptavidin-biotin binding system as a model, we show that the detection limit of the current portable device for streptavidin is 64 nM (equal to 5.12 pmole). In addition, this MPS technique is very versatile and allows for the detection of different diseases just by changing the surface modifications on MNPs.


title: Charged current semi-inclusive deeply inelastic scattering at the Electron-Ion Collider

abstract: We present a systematic calculation of the charged current semi-inclusive deeply inelastic scattering process at leading order twist-3 level in the parton model. We consider the general form of the negatively charged beam scattering off the polarized target which has spin-1. The calculations are carried out by applying the collinear expansion where multiple gluon scattering is taken into account and gauge links are obtained automatically. We first present the general form of the differential cross section in terms of the structure functions by kinematic analysis and then present the structure functions in terms of the the gauge invariant parton distribution functions up to twist-3 level. Considering the angle modulations and polarizations of the cross section, we calculate the complete azimuthal asymmetries in the charged current semi-inclusive deeply inelastic scattering process. The charge asymmetries are also considered in this paper with the introduction of the definitions of the $plus$ and $minus$ cross sections.


title: On a cylinder freely floating in oblique waves

abstract: The coupled motion is investigated for a mechanical system consisting of water and a body freely floating in it. Water occupies either a half-space or a layer of constant depth into which an infinitely long surface-piercing cylinder is immersed, thus allowing us to study the so-called oblique waves. Under the assumption that the motion is of small amplitude near equilibrium and describes time-harmonic oscillations, the phenomenon's linear setting reduces to a spectral problem with the radian frequency as the spectral parameter. If the radiation condition holds, then the total energy is finite and the equipartition of kinetic and potential energy holds for the whole system. On this basis, it is proved that no wave modes are trapped under some restrictions on their frequencies; in the case when a symmetric cylinder has two immersed parts restrictions are imposed on the type of mode as well.


title: PIFE-PIC: Parallel Immersed-Finite-Element Particle-In-Cell For 3-D Kinetic Simulations of Plasma-Material Interactions

abstract: This paper presents a recently developed particle simulation code package PIFE-PIC, which is a novel three-dimensional (3-D) Parallel Immersed-Finite-Element (IFE) Particle-in-Cell (PIC) simulation model for particle simulations of plasma-material interactions. This framework is based on the recently developed non-homogeneous electrostatic IFE-PIC algorithm, which is designed to handle complex plasma-material interface conditions associated with irregular geometries using a Cartesian-mesh-based PIC. Three-dimensional domain decomposition is utilized for both the electrostatic field solver with IFE and the particle operations in PIC to distribute the computation among multiple processors. A simulation of the orbital-motion-limited (OML) sheath of a dielectric sphere immersed in a stationary plasma is carried out to validate PIFE-PIC and profile the parallel performance of the code package. Furthermore, a large-scale simulation of plasma charging at a lunar crater containing 2 million PIC cells (10 million FE/IFE cells) and about 520 million particles, running for 20,000 PIC steps in about 109 wall-clock hours, is presented to demonstrate the high-performance computing capability of PIFE-PIC.


title: An Investigation of Commercial Iron Oxide Nanoparticles: Advanced Structural and Magnetic Properties Characterization

abstract: Magnetic nanoparticles (MNPs) have been extensively used as tiny heating sources in magnetic hyperthermia therapy, contrast agents in magnetic resonance imaging (MRI), tracers in magnetic particle imaging (MPI), carriers for drug/gene delivery, etc. There have emerged many magnetic nanoparticle/microbeads suppliers since the last decade, such as Ocean NanoTech, Nanoprobes, US Research Nanomaterials, Miltenyi Biotec, micromod Partikeltechnologie GmbH, and nanoComposix, etc. In this paper, we report the physical and magnetic characterizations on iron oxide nanoparticle products from Ocean NanoTech. Standard characterization tools such as Vibrating-Sample Magnetometer (VSM), X-Ray Diffraction (XRD), Dynamic Light Scattering (DLS), Transmission Electron Microscopy (TEM), and Zeta Potential Analyzer are used to provide magnetic nanoparticle customers and researchers with an overview of these iron oxide nanoparticle products. In addition, the dynamic magnetic responses of these iron oxide nanoparticles in aqueous solutions are investigated under low and high frequency alternating magnetic fields, giving a standardized operating procedure for characterizing the MNPs from Ocean NanoTech, thereby yielding the best of magnetic nanoparticles for different applications.


title: Sequential Targeting: an incremental learning approach for data imbalance in text classification

abstract: Classification tasks require a balanced distribution of data to ensure the learner to be trained to generalize over all classes. In real-world datasets, however, the number of instances vary substantially among classes. This typically leads to a learner that promotes bias towards the majority group due to its dominating property. Therefore, methods to handle imbalanced datasets are crucial for alleviating distributional skews and fully utilizing the under-represented data, especially in text classification. While addressing the imbalance in text data, most methods utilize sampling methods on the numerical representation of the data, which limits its efficiency on how effective the representation is. We propose a novel training method, Sequential Targeting(ST), independent of the effectiveness of the representation method, which enforces an incremental learning setting by splitting the data into mutually exclusive subsets and training the learner adaptively. To address problems that arise within incremental learning, we apply elastic weight consolidation. We demonstrate the effectiveness of our method through experiments on simulated benchmark datasets (IMDB) and data collected from NAVER.


title: DoDNet: Learning to segment multi-organ and tumors from multiple partially labeled datasets

abstract: Due to the intensive cost of labor and expertise in annotating 3D medical images at a voxel level, most benchmark datasets are equipped with the annotations of only one type of organs and/or tumors, resulting in the so-called partially labeling issue. To address this, we propose a dynamic on-demand network (DoDNet) that learns to segment multiple organs and tumors on partially labeled datasets. DoDNet consists of a shared encoder-decoder architecture, a task encoding module, a controller for generating dynamic convolution filters, and a single but dynamic segmentation head. The information of the current segmentation task is encoded as a task-aware prior to tell the model what the task is expected to solve. Different from existing approaches which fix kernels after training, the kernels in dynamic head are generated adaptively by the controller, conditioned on both input image and assigned task. Thus, DoDNet is able to segment multiple organs and tumors, as done by multiple networks or a multi-head network, in a much efficient and flexible manner. We have created a large-scale partially labeled dataset, termed MOTS, and demonstrated the superior performance of our DoDNet over other competitors on seven organ and tumor segmentation tasks. We also transferred the weights pre-trained on MOTS to a downstream multi-organ segmentation task and achieved state-of-the-art performance. This study provides a general 3D medical image segmentation model that has been pre-trained on a large-scale partially labelled dataset and can be extended (after fine-tuning) to downstream volumetric medical data segmentation tasks. The dataset and code areavailableat: https://git.io/DoDNet


title: Optimizing Approximate Leave-one-out Cross-validation to Tune Hyperparameters

abstract: For a large class of regularized models, leave-one-out cross-validation can be efficiently estimated with an approximate leave-one-out formula (ALO). We consider the problem of adjusting hyperparameters so as to optimize ALO. We derive efficient formulas to compute the gradient and hessian of ALO and show how to apply a second-order optimizer to find hyperparameters. We demonstrate the usefulness of the proposed approach by finding hyperparameters for regularized logistic regression and ridge regression on various real-world data sets.


title: Certified Monotonic Neural Networks

abstract: Learning monotonic models with respect to a subset of the inputs is a desirable feature to effectively address the fairness, interpretability, and generalization issues in practice. Existing methods for learning monotonic neural networks either require specifically designed model structures to ensure monotonicity, which can be too restrictive/complicated, or enforce monotonicity by adjusting the learning process, which cannot provably guarantee the learned model is monotonic on selected features. In this work, we propose to certify the monotonicity of the general piece-wise linear neural networks by solving a mixed integer linear programming problem.This provides a new general approach for learning monotonic neural networks with arbitrary model structures. Our method allows us to train neural networks with heuristic monotonicity regularizations, and we can gradually increase the regularization magnitude until the learned network is certified monotonic. Compared to prior works, our approach does not require human-designed constraints on the weight space and also yields more accurate approximation. Empirical studies on various datasets demonstrate the efficiency of our approach over the state-of-the-art methods, such as Deep Lattice Networks.


title: Simulation of 2D ballistic deposition of porous nanostructured thin-films

abstract: A "two-dimensional ballistic deposition" (2D-BD) code has been developed to study the geometric effects in ballistic deposition of thin-film growth. Circular discs are used as depositing specie to understand the shadowing effects during the evolution of a thin-film. We carried out the 2D-BD simulations for the angles of deposition $20^0$-$80^0$ in steps of $10^0$. Standard deviations $1^0$, $2^0$, $4^0$, $6^0$ and $10^0$ are used for each angle of deposition with disc size of $1.5 \overset{\circ}{A}$ to understand its effect on the microstructure of the thin-films. Angle of growth, porosity and surface roughness properties have been studied for the afore-mentioned angles of deposition and their standard deviations. Ballistic deposition simulations with the discs of different sizes have been carried out to understand the effect of size in ballistic deposition. The results from this code are compared with the available theoretical and experimental results. The code is used to simulate a collimated glancing angle deposition (C-GLAD) experiment. We obtain a good qualitative match for various features of the deposits.


title: Goldblatt-Thomason Theorems for Modal Intuitionistic Logics

abstract: We prove a Goldblatt-Thomason theorem for dialgebraic intuitionistic logics, and instantiate it to Goldblatt-Thomason theorems for a wide variety of modal intuitionistic logics from the literature.


title: Discovery of an M-type Companion to the Herbig Ae Star V1787 Ori

abstract: The intermediate mass Herbig Ae star V1787 Ori is a member of the L1641 star-forming region in the Orion A molecular cloud. We report the detection of an M-type companion to V1787 Ori at a projected separation of 6.66" (corresponding to 2577 au), from the analysis of VLT/NACO adaptive optics $K_s$-band image. Using astrometric data from Gaia DR2, we show that V1787 Ori A and B share similar distance ($d$ $\sim$ 387 pc) and proper motion, indicating that they are physically associated. We estimate the spectral type of V1787 Ori B to be M5 $\pm$ 2 from color--spectral type calibration tables and template matching using SpeX spectral library. By fitting PARSEC models in the Pan-STARRS color-magnitude diagram, we find that V1787 Ori B has an age of 8.1$^{+1.7}_{-1.5}$ Myr and a mass of 0.39$^{+0.02}_{-0.05}$ $M_\odot$. We show that V1787 Ori is a pre-main sequence wide binary system with a mass ratio of 0.23. Such a low mass ratio system is rarely identified in Herbig Ae/Be binary systems. We conclude this work with a discussion on possible mechanisms for the formation of V1787 Ori wide binary system.


title: Complexity Controlled Generative Adversarial Networks

abstract: One of the issues faced in training Generative Adversarial Nets (GANs) and their variants is the problem of mode collapse, wherein the training stability in terms of the generative loss increases as more training data is used. In this paper, we propose an alternative architecture via the Low-Complexity Neural Network (LCNN), which attempts to learn models with low complexity. The motivation is that controlling model complexity leads to models that do not overfit the training data. We incorporate the LCNN loss function for GANs, Deep Convolutional GANs (DCGANs) and Spectral Normalized GANs (SNGANs), in order to develop hybrid architectures called the LCNN-GAN, LCNN-DCGAN and LCNN-SNGAN respectively. On various large benchmark image datasets, we show that the use of our proposed models results in stable training while avoiding the problem of mode collapse, resulting in better training stability. We also show how the learning behavior can be controlled by a hyperparameter in the LCNN functional, which also provides an improved inception score.


title: Designer spin order in diradical nanographenes

abstract: The magnetic properties of carbon materials are at present the focus of an intense research effort in physics, chemistry and materials science due to their potential applications in spintronics and quantum computations. Although the presence of spins in open-shell nanographenes has been recently confirmed, the ability to control magnetic coupling sign has remained elusive, but the most desirable. Here, we demonstrate an effective approach of engineering magnetic ground states in atomically precise open-shell bipartite/nonbipartite nanographenes using combined scanning probe techniques and mean-field Hubbard model calculations. The magnetic coupling sign between two spins has been controlled via breaking bipartite lattice symmetry of nanographenes. In addition, the exchange-interaction strength between two spins has been widely tuned by finely tailoring their spin density overlap, realizing a large exchange-interaction strength of 42 meV. Our demonstrated method provides ample opportunities for designer above-room-temperature magnetic phases and functionalities in graphene nanomaterials.


title: A global universality of two-layer neural networks with ReLU activations

abstract: In the present study, we investigate a universality of neural networks, which concerns a density of the set of two-layer neural networks in a function spaces. There are many works that handle the convergence over compact sets. In the present paper, we consider a global convergence by introducing a norm suitably, so that our results will be uniform over any compact set.


title: Invariance of plurigenera fails in positive and mixed characteristic

abstract: We construct smooth families of elliptic surface pairs with terminal singularities over a DVR of positive or mixed characteristic $(X,B)\to \mathrm{Spec}R$, such that $P_m(X_k,B_k)>P_m(X_K,B_K)$ for all sufficiently divisible $m>0$. In particular, this shows that invariance of all sufficiently divisible plurigenera does not follow from the MMP and Abundance Conjectures.


title: StressNet: Deep Learning to Predict Stress With Fracture Propagation in Brittle Materials

abstract: Catastrophic failure in brittle materials is often due to the rapid growth and coalescence of cracks aided by high internal stresses. Hence, accurate prediction of maximum internal stress is critical to predicting time to failure and improving the fracture resistance and reliability of materials. Existing high-fidelity methods, such as the Finite-Discrete Element Model (FDEM), are limited by their high computational cost. Therefore, to reduce computational cost while preserving accuracy, a novel deep learning model, "StressNet," is proposed to predict the entire sequence of maximum internal stress based on fracture propagation and the initial stress data. More specifically, the Temporal Independent Convolutional Neural Network (TI-CNN) is designed to capture the spatial features of fractures like fracture path and spall regions, and the Bidirectional Long Short-term Memory (Bi-LSTM) Network is adapted to capture the temporal features. By fusing these features, the evolution in time of the maximum internal stress can be accurately predicted. Moreover, an adaptive loss function is designed by dynamically integrating the Mean Squared Error (MSE) and the Mean Absolute Percentage Error (MAPE), to reflect the fluctuations in maximum internal stress. After training, the proposed model is able to compute accurate multi-step predictions of maximum internal stress in approximately 20 seconds, as compared to the FDEM run time of 4 hours, with an average MAPE of 2% relative to test data.


title: Cosheaf Theory and Overconvergent Verdier Duality for Rigid Analytic Spaces

abstract: This paper develops aspects of cosheaf theory on rigid analytic spaces, and demonstrates a sheaf-cosheaf Verdier duality equivalence theorem for overconvergent sheaves on separated, paracompact spaces, analogous to Jacob Lurie's treatment of Verdier duality for locally compact, Hausdorff topological spaces.


title: An ac-susceptibility study of magnetic relaxation phenomena in the antiskyrmion hosting tetragonal Mn-Pt(Pd)-Sn system

abstract: Here, we report an exhaustive study of the frequency-dependent ac-magnetic susceptibility of the $D_{2d}$ symmetric Heusler system Mn-Pt(Pd)-Sn that hosts antiskyrmions over a wide temperature range. Magnetic relaxation studies using Cole-Cole formalism reveal a Debye-type relaxation with a nearly negligible distribution in relaxation times. In contrast to the archetypical skyrmion hosts, the high Curie temperature ($ T_C $) of the present system ensures shorter switching times, and, correspondingly, higher frequencies are required to probe the relaxation dynamics. We find a non-monotonic variation in the characteristic relaxation time with distinct maxima at the phase boundaries \textit{via} helical $\longrightarrow$ antiskyrmion $\longrightarrow$ field-polarized states, indicating slower magnetization dynamics over the region of phase coexistence. The temperature-dependent relaxation time across different phases is of the order of $ 10^{-5} - 10^{-4} $ s and follows the well-known Arrhenius law with reasonable values of the energy barriers. The present study concerning the magnetization dynamics in the antiskyrmion host tetragonal Heusler system is an important contribution towards the basic understanding of the dynamical aspects of antiskyrmions for their potential applications.


title: Novel approach for the study of coherent elastic neutrino-nucleus scattering

abstract: We propose the use of isotopically highly enriched detectors for the precise study of coherent-elastic neutrino-nucleus scattering (CEvNS). CEvNS has been measured for the first time in CsI and recently confirmed with a liquid argon detector. It is expected that several new experimental setups will measure this process with increasing accuracy. Taking Ge detectors as a working example, we demonstrate that a combination of different isotopes is an excellent option to do precision neutrino physics with CEvNS, test Standard Model predictions, and probe new physics scenarios. Experiments based on this new idea can make simultaneous differential CEvNS measurements with detectors of different isotopic composition. Particular combination of observables could be used to cancel systematic errors. While many applications are possible, we illustrate the idea with three examples: testing the dominant quadratic dependence on the number of neutrons, $N$, that is predicted by the theoretical models; constraining the average neutron root mean square (rms) radius; and testing the weak mixing angle and the sensitivity to new physics. In all three cases we find that the extra sensitivity provided by this method will potentially allow high-precision robust measurements with CEvNS and particularly, will resolve the characteristic degeneracies appearing in new physics scenarios.


title: Efficient Conditional Pre-training for Transfer Learning

abstract: Almost all the state-of-the-art neural networks for computer vision tasks are trained by (1) pre-training on a large-scale dataset and (2) finetuning on the target dataset. This strategy helps reduce dependence on the target dataset and improves convergence rate and generalization on the target task. Although pre-training on large-scale datasets is very useful, its foremost disadvantage is high training cost. To address this, we propose efficient filtering methods to select relevant subsets from the pre-training dataset. Additionally, we discover that lowering image resolutions in the pre-training step offers a great trade-off between cost and performance. We validate our techniques by pre-training on ImageNet in both the unsupervised and supervised settings and finetuning on a diverse collection of target datasets and tasks. Our proposed methods drastically reduce pre-training cost and provide strong performance boosts. Finally, we improve standard ImageNet pre-training by 1-3% by tuning available models on our subsets and pre-training on a dataset filtered from a larger scale dataset.


title: Deep Snapshot HDR Imaging Using Multi-Exposure Color Filter Array

abstract: In this paper, we propose a deep snapshot high dynamic range (HDR) imaging framework that can effectively reconstruct an HDR image from the RAW data captured using a multi-exposure color filter array (ME-CFA), which consists of a mosaic pattern of RGB filters with different exposure levels. To effectively learn the HDR image reconstruction network, we introduce the idea of luminance normalization that simultaneously enables effective loss computation and input data normalization by considering relative local contrasts in the "normalized-by-luminance" HDR domain. This idea makes it possible to equally handle the errors in both bright and dark areas regardless of absolute luminance levels, which significantly improves the visual image quality in a tone-mapped domain. Experimental results using two public HDR image datasets demonstrate that our framework outperforms other snapshot methods and produces high-quality HDR images with fewer visual artifacts.


title: One Shot Learning for Speech Separation

abstract: Despite the recent success of speech separation models, they fail to separate sources properly while facing different sets of people or noisy environments. To tackle this problem, we proposed to apply meta-learning to the speech separation task. We aimed to find a meta-initialization model, which can quickly adapt to new speakers by seeing only one mixture generated by those people. In this paper, we use model-agnostic meta-learning(MAML) algorithm and almost no inner loop(ANIL) algorithm in Conv-TasNet to achieve this goal. The experiment results show that our model can adapt not only to a new set of speakers but also noisy environments. Furthermore, we found out that the encoder and decoder serve as the feature-reuse layers, while the separator is the task-specific module.


title: Photon surfaces in less symmetric spacetimes

abstract: We investigate photon surfaces and their stability in a less symmetric spacetime, a general static warped product with a warping function acting on a Riemannian submanifold of codimension two. We find a one-dimensional pseudopotential that gives photon surfaces as its extrema regardless of the spatial symmetry of the submanifold. The maxima and minima correspond to unstable and stable photon surfaces, respectively. It is analogous to the potential giving null circular orbits in a spherically symmetric spacetime. We also see that photon surfaces indeed exist for the spacetimes which are solutions to the Einstein equation. The parameter values for which the photon surfaces exist are specified. As we show finally, the pseudopotential arises due to the separability of the null geodesic equation, and the separability comes from the existence of a Killing tensor in the spacetime. The result leads to the conclusion that photon surfaces may exist even in a less symmetric spacetime if the spacetime admits a Killing tensor.


title: GAN based ball screw drive picture database enlargement for failure classification

abstract: The lack of reliable large datasets is one of the biggest difficulties of using modern machine learning methods in the field of failure detection in the manufacturing industry. In order to develop the function of failure classification for ball screw surface, sufficient image data of surface failures is necessary. When training a neural network model based on a small dataset, the trained model may lack the generalization ability and may perform poorly in practice. The main goal of this paper is to generate synthetic images based on the generative adversarial network (GAN) to enlarge the image dataset of ball screw surface failures. Pitting failure and rust failure are two possible failure types on ball screw surface chosen in this paper to represent the surface failure classes. The quality and diversity of generated images are evaluated afterwards using qualitative methods including expert observation, t-SNE visualization and the quantitative method of FID score. To verify whether the GAN based generated images can increase failure classification performance, the real image dataset was augmented and replaced by GAN based generated images to do the classification task. The authors successfully created GAN based images of ball screw surface failures which showed positive effect on classification test performance.


title: First observation of magnetic flux rope inside electron diffusion region

abstract: Magnrtic flux ropes (MFRs) play a crucial role during magnetic reconnection. They are believed to be primarily generated by tearing mode instabilities in the electron diffusion region (EDR). However, they have never been observed inside the EDR. Here, we present the first observation of an MFR inside an EDR. The bifurcated non-force-free MFR, with a width of 27.5de in the L direction and 4.8de in the N direction, was moving away from the X-line. Inside the MFR, strong energy dissipation was detected. The MFR can modulate the electric field in the EDR. We reconstructed magnetic topology of the electron-scale MFR. Our study promotes understanding of MRFs' initial state and its role in electron-scale processes during magnetic reconnection.


title: The strange case of the transient HBL blazar 4FGL J1544.3-0649

abstract: We present a multifrequency study of the transient $\gamma$-ray source 4FGL J1544.3-0649, a blazar that exhibited a remarkable behaviour raising from the state of an anonymous mid-intensity radio source, never detected at high energies, to that of one of the brightest extreme blazars in the X-ray and $\gamma$-ray sky. Our analysis shows that the averaged $\gamma$-ray spectrum is well described by a powerlaw with a photon index of $1.87\pm0.04$, while the flux above 100 MeV is $(8.0\pm0.9)\times10^{-9}\:{\rm photon\: cm^{-2}\: s^{-1}}$, which increases during the active state of the source. The X-ray flux and spectral slope are both highly variable, with the highest 2-10 keV flux reaching $(1.28\pm0.05)\times10^{-10}\:{\rm erg\:cm^{-2}\:s^{-1}}$. On several observations the X-ray spectrum hardened to the point implying as SED peak moving to energies larger than 10 keV. As in many extreme blazars the broadband spectral energy distribution can be described by a homogeneous one-zone synchrotron-self-Compton leptonic model. We briefly discuss the potential implications for high-energy multi-messenger astrophysics in case the dual behaviour shown by 4FGL J1544.3-0649 does not represent an isolated case, but rather a manifestation of a so far unnoticed relatively common phenomenon.


title: Energy Supply for Heating the Slow Solar Wind Observed by Parker Solar Probe between 0.17 and 0.7 au

abstract: Energy supply sources for the heating process in the slow solar wind remain unknown. The Parker Solar Probe (PSP) mission provides a good opportunity to study this issue. Recently, PSP observations have found that the slow solar wind experiences stronger heating inside 0.24 au. Here for the first time we measure in the slow solar wind the radial gradient of the low-frequency breaks on the magnetic trace power spectra and evaluate the associated energy supply rate. We find that the energy supply rate is consistent with the observed perpendicular heating rate calculated based on the gradient of the magnetic moment. Based on this finding, one could explain why the slow solar wind is strongly heated inside 0.25 au but expands nearly adiabatically outside 0.25 au. This finding supports the concept that the energy added from the energy-containing range is transferred by an energy cascade process to the dissipation range, and then dissipates to heat the slow solar wind. The related issues for further study are discussed.


title: Shuffle and Learn: Minimizing Mutual Information for Unsupervised Hashing

abstract: Unsupervised binary representation allows fast data retrieval without any annotations, enabling practical application like fast person re-identification and multimedia retrieval. It is argued that conflicts in binary space are one of the major barriers to high-performance unsupervised hashing as current methods failed to capture the precise code conflicts in the full domain. A novel relaxation method called Shuffle and Learn is proposed to tackle code conflicts in the unsupervised hash. Approximated derivatives for joint probability and the gradients for the binary layer are introduced to bridge the update from the hash to the input. Proof on $\epsilon$-Convergence of joint probability with approximated derivatives is provided to guarantee the preciseness on update applied on the mutual information. The proposed algorithm is carried out with iterative global updates to minimize mutual information, diverging the code before regular unsupervised optimization. Experiments suggest that the proposed method can relax the code optimization from local optimum and help to generate binary representations that are more discriminative and informative without any annotations. Performance benchmarks on image retrieval with the unsupervised binary code are conducted on three open datasets, and the model achieves state-of-the-art accuracy on image retrieval task for all those datasets. Datasets and reproducible code are provided.


title: Reconstruct Kaplan--Meier Estimator as M-estimator and Its Confidence Band

abstract: The Kaplan--Meier (KM) estimator, which provides a nonparametric estimate of a survival function for time-to-event data, has wide application in clinical studies, engineering, economics and other fields. The theoretical properties of the KM estimator including its consistency and asymptotic distribution have been extensively studied. We reconstruct the KM estimator as an M-estimator by maximizing a quadratic M-function based on concordance, which can be computed using the expectation--maximization (EM) algorithm. It is shown that the convergent point of the EM algorithm coincides with the traditional KM estimator, offering a new interpretation of the KM estimator as an M-estimator. Theoretical properties including the large-sample variance and limiting distribution of the KM estimator are established using M-estimation theory. Simulations and application on two real datasets demonstrate that the proposed M-estimator is exactly equivalent to the KM estimator, while the confidence interval and band can be derived as well.


title: Compressive Shack-Hartmann Wavefront Sensor based on Deep Neural Networks

abstract: The Shack-Hartmann wavefront sensor is widely used to measure aberrations induced by atmospheric turbulence in adaptive optics systems. However if there exists strong atmospheric turbulence or the brightness of guide stars is low, the accuracy of wavefront measurements will be affected. In this paper, we propose a compressive Shack-Hartmann wavefront sensing method. Instead of reconstructing wavefronts with slope measurements of all sub-apertures, our method reconstructs wavefronts with slope measurements of sub-apertures which have spot images with high signal to noise ratio. Besides, we further propose to use a deep neural network to accelerate wavefront reconstruction speed. During the training stage of the deep neural network, we propose to add a drop-out layer to simulate the compressive sensing process, which could increase development speed of our method. After training, the compressive Shack-Hartmann wavefront sensing method can reconstruct wavefronts in high spatial resolution with slope measurements from only a small amount of sub-apertures. We integrate the straightforward compressive Shack-Hartmann wavefront sensing method with image deconvolution algorithm to develop a high-order image restoration method. We use images restored by the high-order image restoration method to test the performance of our the compressive Shack-Hartmann wavefront sensing method. The results show that our method can improve the accuracy of wavefront measurements and is suitable for real-time applications.

categories: astro-ph.IM astro-ph.GA astro-ph.SR cs.CV


title: A Stationary Kyle Setup: Microfounding propagator models

abstract: We provide an economically sound micro-foundation to linear price impact models, by deriving them as the equilibrium of a suitable agent-based system. Our setup generalizes the well-known Kyle model, by dropping the assumption of a terminal time at which fundamental information is revealed so to describe a stationary market, while retaining agents' rationality and asymmetric information. We investigate the stationary equilibrium for arbitrary Gaussian noise trades and fundamental information, and show that the setup is compatible with universal price diffusion at small times, and non-universal mean-reversion at time scales at which fluctuations in fundamentals decay. Our model provides a testable relation between volatility of prices, magnitude of fluctuations in fundamentals and level of volume traded in the market.


title: Point Spread Function Estimation for Wide Field Small Aperture Telescopes with Deep Neural Networks and Calibration Data

abstract: The point spread function (PSF) reflects states of a telescope and plays an important role in development of data processing methods, such as PSF based astrometry, photometry and image restoration. However, for wide field small aperture telescopes (WFSATs), estimating PSF in any position of the whole field of view is hard, because aberrations induced by the optical system are quite complex and the signal to noise ratio of star images is often too low for PSF estimation. In this paper, we further develop our deep neural network (DNN) based PSF modelling method and show its applications in PSF estimation. During the telescope alignment and testing stage, our method collects system calibration data through modification of optical elements within engineering tolerances (tilting and decentering). Then we use these data to train a DNN (Tel--Net). After training, the Tel--Net can estimate PSF in any field of view from several discretely sampled star images. We use both simulated and experimental data to test performance of our method. The results show that the Tel--Net can successfully reconstruct PSFs of WFSATs of any states and in any positions of the FoV. Its results are significantly more precise than results obtained by the compared classic method - Inverse Distance Weight (IDW) interpolation. Our method provides foundations for developing of deep neural network based data processing methods for WFSATs, which require strong prior information of PSFs.


title: Isotropic Scaling Features Measured Locally in the Solar Wind Turbulence with Stationary Background Field

abstract: The scaling anisotropy is crucial to interpret the nonlinear interactions in solar wind turbulence. Previous observations provide diverse results and the structure function analyses are also reported to be an approach to investigate the scaling anisotropy based on a local magnetic field. However, the determination of the sampling angle with respect to the local background magnetic field requires that the observed time series for the average are time stationary. Whether or not this required time stationarity is compatible with the measurements has not been investigated. Here we utilize the second-order structure function method to study the scaling anisotropy with a time-stationary background field.We analyze 88 fast solar wind intervals each with time durations larger than 2 days measured by WIND spacecraft in the period 2005-2018. We calculate the local magnetic field as the average of the time series B whose time stationarity is fulfilled by our criterion phi smaller than 10 deg (phi is the angle between the two averaged magnetic fields after cutting B into two halves). We find for the first time the isotropic scaling feature of the magnetic-trace structure functions with scaling indices -0.63 pm 0.08 and -0.70 pm 0.04, respectively, with the local magnetic field parallel and perpendicular to the solar wind velocity directions. The scaling for the velocity-trace structure functions is also isotropic and the indices are -0.47 pm 0.10 and -0.51 pm 0.09. We also find that with increasing phi threshold to 90 deg, the scaling index of the magnetic-trace structure function in the parallel direction decreases to -0.81, while the rms of the instantaneous angle between magnetic field and solar wind velocity increases up to 45deg at the timescale 150 s, indicating a mix of perpendicular measurements into parallel ones at large scales.


title: Joint Trajectory and Power Allocation Design for Secure Artificial Noise aided UAV Communications

abstract: This paper investigates an average secrecy rate (ASR) maximization problem for an unmanned aerial vehicle (UAV) enabled wireless communication system, wherein a UAV is employed to deliver confidential information to a ground destination in the presence of a terrestrial passive eavesdropper. By employing an artificial noise (AN) injection based secure two-phase transmission protocol, we aim at jointly optimizing the UAV's trajectory, network transmission power, and AN power allocation over a given time horizon to enhance the ASR performance. Specifically, we divide the original non-convex problem into four subproblems, and propose a successive convex approximation based efficient iterative algorithm to solve it suboptimally with guaranteed convergence. Simulation results demonstrate significant security advantages of our designed scheme over other known benchmarks, particularly for stringent flight durations.


title: Perturbation analysis of baroclinic torque in low-Mach-number flows

abstract: In this paper, we propse a series expansion of the baroclinic torque in low-Mach-number flows, so that the accuracy and universality of any buoyancy term could be examined analytically, and new types of buoyancy terms could be constructed and validated. We first demonstrate that the purpose of introducing a buoyancy term is to approximate the baroclinic torque, and straightforwardly the accuracy of any buoyancy term could be measured by the deviation of its curl from the baroclinic torque. Then a regular perturbation method is introduced for the elliptic equation of the hydrodynamic pressure in low-Mach-number flows, resulting in a sequence of Poisson equations, whose solutions lead to the series representation of the baroclinic torque and the new types of buoyancy terms. With the error definition of buoyancy terms and the series representation of the baroclinic torque, the classical gravitational and centrifugal buoyancy term, as well as some other previously proposed buoyancy terms are revisited. Finally, numerical simulations confirm that, with a decreasing density variation or an increasing order of our newly proposed buoyancy term, the simplified equations with one of the new types of buoyancy terms can converge to the original low-Mach-number equations.


title: A phenomenological theory of the optical magnetization reversal

abstract: All-optical switching of the magnetization in magnetic nanostructures by femtosecond circularly polarized laser pulses has been demonstrated in several systems. We present a Landau-Lifshitz-Lambda (LLL) model which describes the magnetization dynamics using three density states: two ferromagnetic grounds states and an excited optical state. One of the ferromagnetic ground states is optically excited by circularly polarized light to a spin reversed state, which is then "Coulomb collapsed" to the magnetization reversed ground state. The time evolution of the optically excited states is described by a Lindblad master equation, in which the optical excitation is introduced via the Hamiltonian. Dissipation terms are introduced via Lindblad operators. The LLL model combines the precessional motion of the magnetization described by the Landau-Lifshitz theory, with the response of the three level Lambda system. The optical excitation lasts for the duration of the laser pulse and the system relaxes at a fast rate due to the electron-electron interaction. We study the solution of the eigenvalues problem of the optical equation of motion for the magnetization and identify a coherent and incoherent regimes and derive an LLL model that can be integrated with existing micromagnetic codes to describe optical excitation of magnetic materials.


title: On mixed Hodge-Riemann relations for translation-invariant valuations and Aleksandrov-Fenchel inequalities

abstract: A version of the Hodge-Riemann relations for valuations was recently conjectured and proved in several special cases by the first-named author. The Lefschetz operator considered there arises as either the product or the convolution with the mixed volume of several Euclidean balls. Here we prove that in (co-)degree one the Hodge-Riemann relations persist if the balls are replaced by several different (centrally symmetric) convex bodies with smooth boundary with positive Gauss curvature. While these mixed Hodge-Riemann relations for the convolution directly imply the Aleksandrov-Fenchel inequality, they yield for the dual operation of the product a new inequality. This new inequality strengthens classical consequences of the Aleksandrov-Fenchel inequality for lower dimensional convex bodies and generalizes some of the geometric inequalities recently discovered by S. Alesker


title: SIMF: Single-Instruction Multiple-Flush Mechanism for Processor Temporal Isolation

abstract: Microarchitectural timing attacks are a type of information leakage attack, which exploit the time-shared microarchitectural components, such as caches, translation look-aside buffers (TLBs), branch prediction unit (BPU), and speculative execution, in modern processors to leak critical information from a victim process or thread. To mitigate such attacks, the mechanism for flushing the on-core state is extensively used by operating-system-level solutions, since on-core state is too expensive to partition. In these systems, the flushing operations are implemented in software (using cache maintenance instructions), which severely limit the efficiency of timing attack protection. To bridge this gap, we propose specialized hardware support, a single-instruction multiple-flush (SIMF) mechanism to flush the core-level state, which consists of L1 caches, BPU, TLBs, and register file. We demonstrate SIMF by implementing it as an ISA extension, i.e., flushx instruction, in scalar in-order RISC-V processor. The resultant processor is prototyped on Xilinx ZCU102 FPGA and validated with state-of-art seL4 microkernel, Linux kernel in multi-core scenarios, and a cache timing attack. Our evaluation shows that SIMF significantly alleviates the overhead of flushing by more than a factor of two in execution time and reduces dynamic instruction count by orders-of-magnitude.