title: Asteroseismic modelling of the subgiant {\mu} Herculis using SONG data: lifting the degeneracy between age and model input parameters

abstract: We model the oscillations of the SONG target $\mu$ Herculis to estimate the parameters of the star. The $\ell$ = 1 mixed modes of $\mu$ Her provide strong constraints on stellar properties. The mass and age given by our asteroseismic modelling are 1.10$^{+0.11}_{-0.06}$ M$_{\odot}$ and 7.55$^{+0.96}_{-0.79}$ Gyr. The initial helium abundance is also constrained at around $Y_{\rm{init}}$ = 0.28, suggesting a ratio in the elements enrichment law ($\Delta Y/\Delta Z$) around 1.3, which is closed to the solar value. The mixing-length parameter converges to about 1.7, which is $\sim$ 10\% lower than the solar value and consistent with the results from hydrodynamic simulations. Our estimates of stellar mass and age agree very well with the previous modelling results with different input physics. Adding asteroseismic information makes these determinations less model-dependent than is typically the case when only surface information is available. Our studies of the model dependence (mass, initial helium and metallicity fractions, and the mixing length parameter) of the age determination indicate that accurate stellar ages ($\lesssim$ 10\%) can be expected from asteroseismic modelling for stars similar to $\mu$ Her. The $\ell$ = 1 bumped modes, which are sensitive to the mean density of the helium core, provide a useful 'clock' that provides additional constraints on its age.


title: On the Origin of Black Hole Spin in High-mass X-Ray Binaries

abstract: Black hole (BH) spins in low-mass X-ray binaries (LMXBs) cover a range of values that can be explained by accretion after BH birth. In contrast, the three BH spin measurements in high-mass X-ray binaries (HMXBs) show only values near the maximum and likely have a different origin connected to the BH stellar progenitor. We explore here two possible scenarios to explain the high spins of BHs in HMXBs: formation in binaries that undergo mass transfer (MT) during the main sequence (MS; Case-A MT), and very close binaries undergoing chemically homogeneous evolution (CHE). We find that both scenarios are able to produce high-spin BHs if internal angular momentum (AM) transport in the progenitor star after its MS evolution is not too strong (i.e., weak coupling between the stellar core and its envelope). If instead efficient AM transport is assumed, we find that the resulting BH spins are always too low with respect to observations. The Case-A MT model provides a good fit for the BH spins, the masses of the two components, and the final orbital periods for two of the three BHs in HMXBs with measured spins. For one of them, the mass predicted for the BH companion is significantly lower than observed, but this depends strongly on the assumed efficiency of mass transfer. The CHE models predict orbital periods that are too large for all three cases considered here. We expect the Case-A MT to be much more frequent at the metallicities involved, so we conclude that the Case-A MT scenario is preferred. Finally, we predict that the stellar companions of HMXBs formed through the Case-A MT have enhanced nitrogen surface abundances, which can be tested by future observations.


title: Non-Cooperative Systems Modeling Repulsive Interaction of Bose-Einstein Condesates in RN

abstract: Inspired by so many possible applications of this class of problems, we seek solution for non-cooperative elliptic systems of two Schrodinger equations. General conditions are assumed under the potentials, which produces convenient spectral properties on the elliptic operator concerned, hence the non-cooperation characterizes it as a strongly indefinite problem. An Abstract Linking Theorem developed previously by the first two authors is the main tool, since variational approach is applied. Furthermore, super and asymptotically quadratic nonlinearities are considered.


title: Wavelength agile multi-photon microscopy with a fiber amplified diode laser

abstract: Multi-photon microscopy is a powerful tool in biomolecular research. Less complex and more cost effective excitation light sources will make this technique accessible to a broader community. Especially semiconductor diode seeded fiber lasers have proven to be robust, low cost and easy to use. However, their wavelength tuning range is often limited, so only a limited number of fluorophores can be accessed. Therefore, different approaches have been proposed to extend the spectral coverage of these lasers. Recently, we showed that four-wave mixing (FWM) assisted stimulated Raman scattering (SRS) can be harnessed to red-shift high power pulses from 1064 nm to a narrowband output at 1122 nm and 1186 nm and therefore extend the number of accessible fluorophores. In this contribution, we show the applicability of all three wavelengths for multi-photon microscopy and analyze the performance.


title: New algorithms for solving stochastic games

abstract: Stochastic games are a classical model in game theory in which two opponents interact and the environment changes in response to the players' behavior. The central solution concepts for these games are the discounted values and the value, which represent what playing the game is worth to the players for different levels of impatience. In the present manuscript, we provide algorithms for computing exact expressions for the discounted values and for the value, which are polynomial in the number of pure stationary strategies of the players. This result considerably improves all the existing algorithms, including the most efficient one, due to Hansen, Kouck\'y, Lauritzen, Miltersen and Tsigaridas (STOC 2011).


title: Conformal vector fields and Yamabe solitons

abstract: In this paper, we use less topological restrictions and more geometric and analytic conditions to obtain some sufficient conditions on Yamabe solitons such that their metrics are Yamabe metrics, that is, metrics of constant scalar curvature. More precisely, we use properties of conformal vector fields to find several sufficient conditions on the soliton vector fields of Yamabe solitons under which their metrics are of Yamabe metrics.


title: Weighted EGZ Constant for p-groups of rank 2

abstract: Let $G$ be a finite abelian group of exponent $n$, written additively, and let $A$ be a subset of $\mathbb{Z}$. The constant $s_A(G)$ is defined as the smallest integer $\ell$ such that any sequence over $G$ of length at least $\ell$ has an $A$-weighted zero-sum of length $n$ and $\eta_A(G)$ defined as the smallest integer $\ell$ such that any sequence over $G$ of length at least $\ell$ has an $A$-weighted zero-sum of length at most $n$. Here we prove that, for $\alpha \geq \beta$, and $A=\left\{x\in\mathbb{N}\; : \; 1 \le a \le p^{\alpha} \; \mbox{ and }\; \gcd(a, p) = 1\right \}$, we have $s_{A}(\mathbb{Z}_{p^{\alpha}}\oplus \mathbb{Z}_{p^\beta}) = \eta_A(\mathbb{Z}_{p^{\alpha}}\oplus \mathbb{Z}_{p^\beta}) + p^{\alpha}-1 = p^{\alpha} + \alpha +\beta$ and classify all the extremal $A$-weighted zero-sum free sequences.


title: Machine Parameters and Projected Luminosity Performance of Proposed Future Colliders at CERN

abstract: In response to a request from the CERN Scientific Policy Committee (SPC), the machine parameters and expected luminosity performance for several proposed post-LHC collider projects at CERN are compiled: three types of hadron colliders (HL-LHC upgrade, FCC-hh and HE-LHC), a circular lepton collider (FCC-ee), a linear lepton collider (CLIC), and three options for lepton-hadron colliders (LHeC, HE-LHeC, and FCC-eh). Particular emphasis is put on availability, physics run time, and efficiency. The information contained in this document was presented at the SPC Meeting of September 2018. It will serve as one of the inputs to the 2019/20 Update of the European Strategy for Particle Physics.


title: Batalin-Vilkovisky structure on Hochschild cohomology with coefficients in the dual algebra

abstract: We prove that Hochschild cohomology with coefficients in $A^*=Hom_k(A,k)$ together with an $A$-structural map $\psi:A^* \otimes_A A^* \to A^*$ is a Batalin-Vilkovisky algebra. This applies to symmetric, Frobenius and monomial path algebras.


title: Bi-Directional Lattice Recurrent Neural Networks for Confidence Estimation

abstract: The standard approach to mitigate errors made by an automatic speech recognition system is to use confidence scores associated with each predicted word. In the simplest case, these scores are word posterior probabilities whilst more complex schemes utilise bi-directional recurrent neural network (BiRNN) models. A number of upstream and downstream applications, however, rely on confidence scores assigned not only to 1-best hypotheses but to all words found in confusion networks or lattices. These include but are not limited to speaker adaptation, semi-supervised training and information retrieval. Although word posteriors could be used in those applications as confidence scores, they are known to have reliability issues. To make improved confidence scores more generally available, this paper shows how BiRNNs can be extended from 1-best sequences to confusion network and lattice structures. Experiments are conducted using one of the Cambridge University submissions to the IARPA OpenKWS 2016 competition. The results show that confusion network and lattice-based BiRNNs can provide a significant improvement in confidence estimation.


title: Confidence Estimation and Deletion Prediction Using Bidirectional Recurrent Neural Networks

abstract: The standard approach to assess reliability of automatic speech transcriptions is through the use of confidence scores. If accurate, these scores provide a flexible mechanism to flag transcription errors for upstream and downstream applications. One challenging type of errors that recognisers make are deletions. These errors are not accounted for by the standard confidence estimation schemes and are hard to rectify in the upstream and downstream processing. High deletion rates are prominent in limited resource and highly mismatched training/testing conditions studied under IARPA Babel and Material programs. This paper looks at the use of bidirectional recurrent neural networks to yield confidence estimates in predicted as well as deleted words. Several simple schemes are examined for combination. To assess usefulness of this approach, the combined confidence score is examined for untranscribed data selection that favours transcriptions with lower deletion errors. Experiments are conducted using IARPA Babel/Material program languages.


title: 2MASS J10274572+0629104: The very short period young M6 Dwarf Binary System identified in K2 data

abstract: We report the identification of a very low mass new binary system 2MASS J10274572+0629104, based on Kepler K2 photometry and Gaia DR2 astrometry. The K2 light curve is consistent with a beat pattern of two periodic signals, and using Lomb-Scargle periodogram, we find two rotation periods of 0.2114 \pm 0.0002 day and 0.2199 \pm 0.0003 day.We conclude that these rotation periods arise from two stars with similar spectral types of M6, and have near equal luminosity. It is the first ultracool binary system to be identified based on beat patterns in the light curve. Near-infrared spectroscopy yields RV = -9.8 \pm 0.6 km s^-1, v sin i = 21.5 \pm 1.1 km s^-1, T_eff = 3110 \pm 40 K, and log g = 5.2 \pm 0.2. The motions are consistent with a young age, as are the rotation periods, but the source does not appear to be part of any known moving group. Furthermore, we detected three strong white light flares in the K2 light curve, with estimated total bolometric (UV/optical/IR) energies to be 2.6 X 10^33, 5.0 X 10^33, and 3.5 X 10^33 erg respectively.


title: Whole genome single nucleotide polymorphism genotyping of Staphylococcus aureus

abstract: Next-generation sequencing technology enables routine detection of bacterial pathogens for clinical diagnostics and genetic research. Whole genome sequencing has been of importance in the epidemiologic analysis of bacterial pathogens. However, few whole genome sequencing-based genotyping pipelines are available for practical applications. Here, we present the whole genome sequencing-based single nucleotide polymorphism (SNP) genotyping method and apply to the evolutionary analysis of methicillin-resistant Staphylococcus aureus. The SNP genotyping method calls genome variants using next-generation sequencing reads of whole genomes and calculates the pair-wise Jaccard distances of the genome variants. The method may reveal the high-resolution whole genome SNP profiles and the structural variants of different isolates of methicillin-resistant S. aureus (MRSA) and methicillin-susceptible S. aureus (MSSA) strains. The phylogenetic analysis of whole genomes and particular regions may monitor and track the evolution and the transmission dynamic of bacterial pathogens. The computer programs of the whole genome sequencing-based SNP genotyping method are available to the public at https://github.com/cyinbox/NGS.


title: Designing Informative Rating Systems: Evidence from an Online Labor Market

abstract: Platforms critically rely on rating systems to learn the quality of market participants. In practice, however, these ratings are often highly inflated, and therefore not very informative. In this paper, we first investigate whether the platform can obtain less inflated, more informative ratings by altering the meaning and relative importance of the levels in the rating system. Second, we seek a principled approach for the platform to make these choices in the design of the rating system. First, we analyze the results of a randomized controlled trial on an online labor market in which an additional question was added to the feedback form. Between treatment conditions, we vary the question phrasing and answer choices; in particular, the treatment conditions include several positive-skewed verbal rating scales with descriptive phrases or adjectives providing specific interpretation for each rating level. The online labor market test reveals that current inflationary norms can in fact be countered by re-anchoring the meaning of the levels of the rating system. In particular, the positive-skewed verbal rating scales yield rating distributions that significantly reduce rating inflation and are much more informative about seller quality. Second, we develop a model-based framework to compare and select among rating system designs, and apply this framework to the data obtained from the online labor market test. Our simulations demonstrate that our model-based framework for scale design and optimization can identify the most informative rating system and substantially improve the quality of information obtained over baseline designs. Overall, our study illustrates that rating systems that are informative in practice can be designed, and demonstrates how to design them in a principled manner.


title: BCL: A Cross-Platform Distributed Container Library

abstract: One-sided communication is a useful paradigm for irregular parallel applications, but most one-sided programming environments, including MPI's one-sided interface and PGAS programming languages, lack application level libraries to support these applications. We present the Berkeley Container Library, a set of generic, cross-platform, high-performance data structures for irregular applications, including queues, hash tables, Bloom filters and more. BCL is written in C++ using an internal DSL called the BCL Core that provides one-sided communication primitives such as remote get and remote put operations. The BCL Core has backends for MPI, OpenSHMEM, GASNet-EX, and UPC++, allowing BCL data structures to be used natively in programs written using any of these programming environments. Along with our internal DSL, we present the BCL ObjectContainer abstraction, which allows BCL data structures to transparently serialize complex data types while maintaining efficiency for primitive types. We also introduce the set of BCL data structures and evaluate their performance across a number of high-performance computing systems, demonstrating that BCL programs are competitive with hand-optimized code, even while hiding many of the underlying details of message aggregation, serialization, and synchronization.


title: Twistors, quartics, and del Pezzo fibrations

abstract: We investigate the structure of a variety of new Moishezon twistor spaces, by utilizing the pluri-half-anti-canonical map from the twistor spaces. Each of these twistor spaces is bimeromorphic to a double covering of a scroll of planes over a rational normal curve, and the branch divisor of the double cover is a cut of the scroll by a quartic hypersurface. In particular, the double covering has a pencil of Del Pezzo surfaces of degree two. Correspondingly, the twistor spaces have a pencil of rational surfaces with big anti-canonical class. The base locus of the last pencil is a cycle of rational curves, and it is an anti-canonical curve on smooth members of the pencil. These twistor spaces are naturally classified into four types according to the type of singularities of the branch divisor, or equivalently, those of the Del Pezzo surfaces in the pencil. We also show that the quartic hypersurface satisfies a strong constraint and as a result the defining polynomial of the quartic hypersurface has to be of a specific form. Together with our previous result, the present result completes a classification of Moishezon twistor spaces whose half-anti-canonical system is a pencil. Twistor spaces whose half-anti-canonical system is larger than pencil have been understood for a long time before. In the opposite direction, no example is known of a Moishezon twistor space whose half-anti-canonical system is smaller than a pencil.


title: On travelling waves of the non linear Schr\"odinger equation escaping a potential well

abstract: In this paper we consider the NLS equation with focusing nonlinearities in the presence of a potential. We investigate the compact soliton motions that correspond to a free soliton escaping the well created by the potential. We exhibit the dynamical system driving the exiting trajectory and construct associated nonlinear dynamics for untrapped motions. We show that the nature of the potential/soliton is fundamental, and two regimes may exist: one where the tail of the potential is fat and dictates the motion, one where the tail is weak and the soliton self interacts with the potential defects, hence leading to different motions.


title: Bell-state correlations in current of local Fermi liquid

abstract: We study Bell-state correlations for quasiparticle pairs excited in nonlinear current through a double quantum dot in the Kondo regime. Exploiting the renormalized perturbation expansion in the residual interactions of the local Fermi liquid and Bell's inequality for cross correlation of spin currents through distinct conduction channels, we derive an asymptotically exact form of Bell's correlation for the double dot at low bias voltages. We find that pairs of quasiparticles and holes excited by the residual exchange interaction can violate Bell's inequality for the cross correlations of the spin currents.


title: Stress-Testing Neural Models of Natural Language Inference with Multiply-Quantified Sentences

abstract: Standard evaluations of deep learning models for semantics using naturalistic corpora are limited in what they can tell us about the fidelity of the learned representations, because the corpora rarely come with good measures of semantic complexity. To overcome this limitation, we present a method for generating data sets of multiply-quantified natural language inference (NLI) examples in which semantic complexity can be precisely characterized, and we use this method to show that a variety of common architectures for NLI inevitably fail to encode crucial information; only a model with forced lexical alignments avoids this damaging information loss.


title: Theory for Inverse Stochastic Resonance in Nature

abstract: The inverse stochastic resonance (ISR) phenomenon consists in an unexpected depression in the response of a system under external noise, e.g., as observed in the behavior of the mean-firing rate in some pacemaker neurons in the presence of moderate values of noise. A possible requirement for such behavior is the existence of a bistable regime in the behavior of these neurons. We here explore theoretically the possible emergence of this behavior in a general bistable system, and conclude on conditions the potential function which drives the dynamics must accomplish. We show that such an intriguing, and apparently widely observed, phenomenon ensues in the case of an asymmetric potential function when the high activity minimum state of the system is metastable with the largest basin of attraction and the low activity state is the global minimum with a smaller basin of attraction. We discuss on the relevance of such a picture to understand the ISR features and to predict its general appearance in other natural systems that share the requirements described here. Finally, we report another intriguing non-standard stochastic resonance in our system, which occurs in the absence of any weak signal input into the system and whose emergence can be explained, with the ISR, within our theoretical framework in this paper in terms of the shape of the potential function.


title: Chemically Accurate 0-0 Energies with not-so-Accurate Excited State Geometries

abstract: Using a series of increasingly refined wavefunction methods able to tackle electronic excited states, namely ADC(2), CC2, CCSD, CCSDR(3) and CC3, we investigate the interplay between geometries and 0-0 energies. We show that, due to a strong and nearly systematic error cancelation between the vertical transition and geometrical reorganization energies, CC2 and CCSD structures can be used to obtain chemically-accurate 0-0 energies, though the underlying geometries are rather far from the reference ones and would deliver significant errors for many chemical and physical properties. This indicates that obtaining 0-0 energies matching experiment does not demonstrate the quality of the geometrical parameters. In contrast, accurate computation of vertical excitation energies is mandatory in order to reach chemical accuracy. By determining CC3 total energies on CCSD structures, we model a large set of compounds (including radicals) and electronic transitions (including singlet-triplet excitations) and successfully reach chemical accuracy in a near systematic way. Indeed, for this particular set, our protocol delivers a mean absolute error as small as $0.032$ eV, chemical accuracy (error smaller than $1$ kcal.mol$^{-1}$ or $0.043$ eV) being obtained in 80\%\ of the cases. In only three cases the error exceeds $0.15$ eV which is of the order of the typical error provided by TD-DFT or second-order wavefunction methods for this particular property. The present composite protocol is therefore very effective despite the fact that the geometries may not be considered as very accurate.


title: Parameterized Complexity of Equitable Coloring

abstract: A graph on $n$ vertices is equitably $k$-colorable if it is $k$-colorable and every color is used either $\left\lfloor n/k \right\rfloor$ or $\left\lceil n/k \right\rceil$ times. Such a problem appears to be considerably harder than vertex coloring, being $\mathsf{NP\text{-}Complete}$ even for cographs and interval graphs. In this work, we prove that it is $\mathsf{W[1]\text{-}Hard}$ for block graphs and for disjoint union of split graphs when parameterized by the number of colors; and $\mathsf{W[1]\text{-}Hard}$ for $K_{1,4}$-free interval graphs when parameterized by treewidth, number of colors and maximum degree, generalizing a result by Fellows et al. (2014) through a much simpler reduction. Using a previous result due to Dominique de Werra (1985), we establish a dichotomy for the complexity of equitable coloring of chordal graphs based on the size of the largest induced star. Finally, we show that \textsc{equitable coloring} is $\mathsf{FPT}$ when parameterized by the treewidth of the complement graph.


title: Variability in Proto-Planetary Nebulae: V. Velocity and Light Curve Analyses of IRAS 17436+5003, 18095+2704, and 19475+3119

abstract: We have obtained contemporaneous light, color, and radial velocity data for three proto-planetary nebulae (PPNe) over the years 2007 to 2015. The light and velocity curves of each show similar periods of pulsation, with photometric periods of 42 and 50 days for IRAS 17436+5003, 102 days for IRAS 18095+2704, and 35 days for IRAS 19475+3119. The light and velocity curves are complex with multiple periods and small, variable amplitudes. Nevertheless, at least over limited time intervals, we were able to identify dominant periods in the light, color, and velocity curves and compare the phasing of each. The color curves appear to peak with or slightly after the light curves while the radial velocity curves peak about a quarter of a cycle before the light curves. Similar results were found previously for two other PPNe, although for them the light and color appeared to be in phase. Thus it appears that PPNe are brightest when smallest and hottest. These phase results differ from those found for classical Cepheid variables, where the light and velocity differ by half a cycle, and are hottest at about average size and expanding. However, they do appear to have similar phasing to the larger amplitude pulsations seen in RV Tauri variables. Presently, few pulsation models exist for PPNe, and these do not fit the observations well, especially the longer periods observed. Model fits to these new light and velocity curves would allow masses to be determined for these post-AGB objects, and thereby provide important constraints to post-AGB stellar evolution models of low and intermediate-mass stars.


title: Spectral Method for Multiplexed Phase Retrieval and Application in Optical Imaging in Complex Media

abstract: We introduce a generalized version of phase retrieval called multiplexed phase retrieval. We want to recover the phase of amplitude-only measurements from linear combinations of them. This corresponds to the case in which multiple incoherent sources are sampled jointly, and one would like to recover their individual contributions. We show that a recent spectral method developed for phase retrieval can be generalized to this setting, and that its performance follows a phase transition behavior. We apply this new technique to light focusing at depth in a complex medium. Experimentally, although we only have access to the sum of the intensities on multiple targets, we are able to separately focus on each ones, thus opening potential applications in deep fluorescence imaging and light delivery


title: A conservative implicit scheme for steady state solutions of diatomic gas flow in all flow regimes

abstract: An implicit scheme for steady state solutions of diatomic gas flow is presented. The method solves the Rykov model equation in the finite volume discrete velocity method (DVM) framework, in which the translational and rotational degrees of freedom are taken into account. At the cell interface, a difference scheme of the model equation is used to construct a multiscale flux (similar to discrete unified gas-kinetic scheme (DUGKS)), so that the cell size is not constrained by the cell Knudsen (Kn) number. The physical local time step is implemented to preserve the multiscale property in the nonuniform-mesh case. The implicit macroscopic prediction technique is adopted to find a predicted equilibrium state at each time level and the implicit macroscopic governing equation is solved along with the implicit microscopic system. Furthermore, an efficient integral error compensation technique is applied, which makes the scheme conservative and allows more flexible discretization for particle velocity space. In the test cases, the unstructured velocity-space mesh is used, the present method is proved to be efficient and accurate.


title: Enforcing public data archiving policies in academic publishing: A study of ecology journals

abstract: To improve the quality and efficiency of research, groups within the scientific community seek to exploit the value of data sharing. Funders, institutions, and specialist organizations are developing and implementing strategies to encourage or mandate data sharing within and across disciplines, with varying degrees of success. Academic journals in ecology and evolution have adopted several types of public data archiving policies requiring authors to make data underlying scholarly manuscripts freely available. Yet anecdotes from the community and studies evaluating data availability suggest that these policies have not obtained the desired effects, both in terms of quantity and quality of available datasets. We conducted a qualitative, interview-based study with journal editorial staff and other stakeholders in the academic publishing process to examine how journals enforce data archiving policies. We specifically sought to establish who editors and other stakeholders perceive as responsible for ensuring data completeness and quality in the peer review process. Our analysis revealed little consensus with regard to how data archiving policies should be enforced and who should hold authors accountable for dataset submissions. Themes in interviewee responses included hopefulness that reviewers would take the initiative to review datasets and trust in authors to ensure the completeness and quality of their datasets. We highlight problematic aspects of these thematic responses and offer potential starting points for improvement of the public data archiving process.


title: Concurrent coupling of atomistic simulation and mesoscopic hydrodynamics for flows over soft multi-functional surfaces

abstract: We develop an efficient parallel multiscale method that bridges the atomistic and mesoscale regimes, from nanometer to micron and beyond, via concurrent coupling of atomistic simulation and mesoscopic dynamics. In particular, we combine an all-atom molecular dynamics (MD) description for specific atomistic details in the vicinity of the functional surface, with a dissipative particle dynamics (DPD) approach that captures mesoscopic hydrodynamics in the domain away from the functional surface. In order to achieve a seamless transition in dynamic properties we endow the MD simulation with a DPD thermostat, which is validated against experimental results by modeling water at different temperatures. We then validate the MD-DPD coupling method for transient Couette and Poiseuille flows, demonstrating that the concurrent MD-DPD coupling can resolve accurately the continuum-based analytical solutions. Subsequently, we simulate shear flows over polydimethylsiloxane (PDMS)-grafted surfaces (polymer brushes) for various grafting densities, and investigate the slip flow as a function of the shear stress. We verify that a "universal" power law exists for the sliplength, in agreement with published results. Having validated the MD-DPD coupling method, we simulate time-dependent flows past an endothelial glycocalyx layer (EGL) in a microchannel. Coupled simulation results elucidate the dynamics of EGL changing from an equilibrium state to a compressed state under shear by aligning the molecular structures along the shear direction. MD-DPD simulation results agree well with results of a single MD simulation, but with the former more than two orders of magnitude faster than the latter for system sizes above one micron.


title: DiracSolver: a tool for solving the Dirac Equation

abstract: Advantageous numerical methods for solving the Dirac equations are derived. They are based on different stochastic optimization techniques, namely the Genetic algorithms, the Particle Swarm Optimization and the Simulated Annealing method, their use of which is favored from intuitive, practical, and theoretical arguments. Towards this end, we optimize appropriate parametric expressions representing the radial Dirac wave functions by employing methods that minimize multi parametric expressions in several physical applications. As a concrete application, we calculate the small (bottom) and large (top) components of the Dirac wave function for a bound muon orbiting around a very heavy (complex) nuclear system (the $^{208}$Pb nucleus), but the new approach may effectively be applied in other complex atomic, nuclear and molecular systems.


title: Stochastic Optimal Control of Epidemic Processes in Networks

abstract: We approach the development of models and control strategies of susceptible-infected-susceptible (SIS) epidemic processes from the perspective of marked temporal point processes and stochastic optimal control of stochastic differential equations (SDEs) with jumps. In contrast to previous work, this novel perspective is particularly well-suited to make use of fine-grained data about disease outbreaks and lets us overcome the shortcomings of current control strategies. Our control strategy resorts to treatment intensities to determine who to treat and when to do so to minimize the amount of infected individuals over time. Preliminary experiments with synthetic data show that our control strategy consistently outperforms several alternatives. Looking into the future, we believe our methodology provides a promising step towards the development of practical data-driven control strategies of epidemic processes.


title: Scalable Laplacian K-modes

abstract: We advocate Laplacian K-modes for joint clustering and density mode finding, and propose a concave-convex relaxation of the problem, which yields a parallel algorithm that scales up to large datasets and high dimensions. We optimize a tight bound (auxiliary function) of our relaxation, which, at each iteration, amounts to computing an independent update for each cluster-assignment variable, with guaranteed convergence. Therefore, our bound optimizer can be trivially distributed for large-scale data sets. Furthermore, we show that the density modes can be obtained as byproducts of the assignment variables via simple maximum-value operations whose additional computational cost is linear in the number of data points. Our formulation does not need storing a full affinity matrix and computing its eigenvalue decomposition, neither does it perform expensive projection steps and Lagrangian-dual inner iterates for the simplex constraints of each point. Furthermore, unlike mean-shift, our density-mode estimation does not require inner-loop gradient-ascent iterates. It has a complexity independent of feature-space dimension, yields modes that are valid data points in the input set and is applicable to discrete domains as well as arbitrary kernels. We report comprehensive experiments over various data sets, which show that our algorithm yields very competitive performances in term of optimization quality (i.e., the value of the discrete-variable objective at convergence) and clustering accuracy.


title: Analytic Variable Exponent Hardy Spaces

abstract: We introduce a variable exponent version of the Hardy space of analytic functions on the unit disk, we show some properties of the space, and give an example of a variable exponent $p(\cdot)$ that satisfies the $\log$-H\"older condition such that $H^{p(\cdot)}\neq H^q$ for any constant exponent $1<q<\infty$. We also consider the variable exponent version of the Hardy space on the upper-half plane.


title: An immersed boundary method for fluid-structure interaction based on overlapping domain decomposition

abstract: We present a novel framework inspired by the Immersed Boundary Method for predicting the fluid-structure interaction of complex structures immersed in flows with moderate to high Reynolds numbers. The main novelties of the proposed fluid-structure interaction framework are 1) the use of elastodynamics equations for the structure, 2) the use of a high-order Navier-Stokes solver for the flow, and 3) the variational transfer (L2-projection) for coupling the solid and fluid subproblems. The dynamic behavior of a deformable structure is simulated in a finite element framework by adopting a fully implicit scheme for its temporal integration. It allows for mechanical constitutive laws including nonhomogeneous and fiber-reinforced materials. The Navier-Stokes equations for the incompressible flow are discretized with high-order finite differences which allow for the direct numerical simulation of laminar, transitional and turbulent flows. The structure and the flow solvers are coupled by using an L2-projection method for the transfer of velocities and forces between the fluid grid and the solid mesh. This strategy allows for the numerical solution of coupled large scale problems based on nonconforming structured and unstructured grids. The framework is validated with the Turek-Hron benchmark and a newly proposed benchmark modelling the flow-induced oscillation of an inert plate. A three-dimensional simulation of an elastic beam in transitional flow is provided to show the solver's capability of coping with anisotropic elastic structures immersed in complex fluid flow.


title: Uncertainty quantification in Eulerian-Lagrangian simulations of (point-)particle-laden flows with data-driven and empirical forcing models

abstract: An uncertainty quantification framework is developed for Eulerian-Lagrangian models of particle-laden flows, where the fluid is modeled through a system of partial differential equations in the Eulerian frame and inertial particles are traced as points in the Lagrangian frame. The source of uncertainty in such problems is the particle forcing, which is determined empirically or computationally with high-fidelity methods (data-driven). The framework relies on the averaging of the deterministic governing equations with the stochastic forcing and allows for an estimation of the first and second moment of the quantities of interest. Via comparison with Monte Carlo simulations, it is demonstrated that the moment equations accurately predict the uncertainty for problems whose Eulerian dynamics are either governed by the linear advection equation or the compressible Euler equations. In areas of singular particle interfaces and shock singularities significant uncertainty is generated. An investigation into the effect of the numerical methods shows that low-dissipative higher-order methods are necessary to capture numerical singularities (shock discontinuities, singular source terms, particle clustering) with low diffusion in the propagation of uncertainty.


title: Attentive Filtering Networks for Audio Replay Attack Detection

abstract: An attacker may use a variety of techniques to fool an automatic speaker verification system into accepting them as a genuine user. Anti-spoofing methods meanwhile aim to make the system robust against such attacks. The ASVspoof 2017 Challenge focused specifically on replay attacks, with the intention of measuring the limits of replay attack detection as well as developing countermeasures against them. In this work, we propose our replay attacks detection system - Attentive Filtering Network, which is composed of an attention-based filtering mechanism that enhances feature representations in both the frequency and time domains, and a ResNet-based classifier. We show that the network enables us to visualize the automatically acquired feature representations that are helpful for spoofing detection. Attentive Filtering Network attains an evaluation EER of 8.99$\%$ on the ASVspoof 2017 Version 2.0 dataset. With system fusion, our best system further obtains a 30$\%$ relative improvement over the ASVspoof 2017 enhanced baseline system.


title: Cooperative Holistic Scene Understanding: Unifying 3D Object, Layout, and Camera Pose Estimation

abstract: Holistic 3D indoor scene understanding refers to jointly recovering the i) object bounding boxes, ii) room layout, and iii) camera pose, all in 3D. The existing methods either are ineffective or only tackle the problem partially. In this paper, we propose an end-to-end model that simultaneously solves all three tasks in real-time given only a single RGB image. The essence of the proposed method is to improve the prediction by i) parametrizing the targets (e.g., 3D boxes) instead of directly estimating the targets, and ii) cooperative training across different modules in contrast to training these modules individually. Specifically, we parametrize the 3D object bounding boxes by the predictions from several modules, i.e., 3D camera pose and object attributes. The proposed method provides two major advantages: i) The parametrization helps maintain the consistency between the 2D image and the 3D world, thus largely reducing the prediction variances in 3D coordinates. ii) Constraints can be imposed on the parametrization to train different modules simultaneously. We call these constraints "cooperative losses" as they enable the joint training and inference. We employ three cooperative losses for 3D bounding boxes, 2D projections, and physical constraints to estimate a geometrically consistent and physically plausible 3D scene. Experiments on the SUN RGB-D dataset shows that the proposed method significantly outperforms prior approaches on 3D object detection, 3D layout estimation, 3D camera pose estimation, and holistic scene understanding.


title: Characters for Projective Modules in the BGG Category O for General Linear Lie Superalgebras

abstract: We determine the Verma multiplicities and the characters of projective modules for atypical blocks in the BGG Category O for the general linear Lie superalgebras $\frak{gl}(2|2)$ and $\frak{gl}(3|1)$. We then explicitly determine the composition factor multiplcities of Verma modules in the atypicality 2 block of $\frak{gl}(2|2)$.


title: Minimal Surfaces and Generalized Einstein-Maxwell-dilaton Theory

abstract: We present novel classes of non-stationary solutions to the five-dimensional generalized Einstein-Maxwell-dilaton theory with cosmological constant, in which the Maxwell's filed and the cosmological constant couple to the dilaton field. In the first class of solutions, the two non-zero coupling constants are different while in the second class of solutions, the two coupling constants are equal to each other. We find consistent cosmological solutions with positive, negative or zero cosmological constant, where the cosmological constant depends on the value of one coupling constant in the theory. Moreover, we discuss the physical properties of the five-dimensional solutions and the uniqueness of the solutions in five dimensions by showing the solutions with different coupling constants, can't be uplifted to any Einstein-Maxwell theory in higher dimensions.


title: Enhancing quadratic optomechanical coupling via nonlinear medium and lasers

abstract: We propose a scheme to significantly increase quadratic optomechanical couplings of optomechanical systems with the help of a nonlinear medium and two driving lasers. The nonlinear medium is driven by one laser and the optical cavity mode is driven by a strong laser. We derive an effective Hamiltonian using squeezing transformation and rotating wave approximation. The effective quadratic optomechanical coupling strength can be larger than the decay rate of the cavity mode by adjusting the two optical driving fields. The thermal noise of squeezed cavity mode can be suppressed totally with the help of a squeezed vacuum field. Then, a driving field is applied to the mechanical mode. We investigate the equal-time second order correlations and find there are photon, phonon, and photon-phonon blockades even the original single-photon quadratic coupling is much smaller than the decay rate of the optical mode. In addition, the sub-Poissonian window of the two-time second order correlations can be controlled by the mechanical driving field. Finally, we show the squeezing and entanglement of the model could be tuned by the driving fields of the nonlinear medium and mechanical mode.


title: Wavelength-resolved neutron tomography for crystalline materials

abstract: Wavelength-resolved (WR) neutron transmission tomography is an emerging technique to characterize engineering materials. While tomographic reconstruction for amorphous samples is straightforward, it is challenging to reconstruct samples with single-crystal domains because the attenuation of the sample varies as a function of its orientation with respect to the incident beam due to Bragg scattering. In this paper, we present an algorithm that can reconstruct samples with single-crystal domains from WR neutron tomographic measurements. In particular, we use a model-based iterative reconstruction (MBIR) technique that reconstructs the volume by identifying and leaving out the regions of the measurement that are affected by Bragg scatter. We combine the output of the MBIR method with an algorithm that matches the reconstruction to the identified Bragg scatter to reconstruct a feature that corresponds to the local crystallography of the sample being measured. Using simulated data, we demonstrate how our algorithm can reconstruct materials with single-crystal domains, thereby adding a powerful new capability for WR neutron imaging instruments.


title: User Constrained Thumbnail Generation using Adaptive Convolutions

abstract: Thumbnails are widely used all over the world as a preview for digital images. In this work we propose a deep neural framework to generate thumbnails of any size and aspect ratio, even for unseen values during training, with high accuracy and precision. We use Global Context Aggregation (GCA) and a modified Region Proposal Network (RPN) with adaptive convolutions to generate thumbnails in real time. GCA is used to selectively attend and aggregate the global context information from the entire image while the RPN is used to predict candidate bounding boxes for the thumbnail image. Adaptive convolution eliminates the problem of generating thumbnails of various aspect ratios by using filter weights dynamically generated from the aspect ratio information. The experimental results indicate the superior performance of the proposed model over existing state-of-the-art techniques.


title: Signatures of Intra-binary shock emission in the black widow pulsar binary PSR J2241-5236

abstract: We report on high-energy properties of the black widow pulsar PSR J2241$-$5236 in the X-ray and the Fermi-LAT (GeV gamma-ray) bands. In the LAT band, the phase-averaged gamma-ray light curve shows orbital modulation below $\sim$1 GeV with a chance probability ($p$) monotonically decreasing with time to $p\sim 10^{-5}$. The peak of the light curve is near the superior conjunction of the pulsar (binary phase $\phi_{\rm B}\approx 0.25$). We attribute the modulation to the intra-binary shock (IBS) emission and search for IBS signatures in the archival X-ray data. We find that the X-ray spectral fit requires a non-thermal component, which implies a possible IBS origin of the X-rays. We discuss our observations in the context of IBS scenarios.


title: Volume pinching theorems for CAT(1) spaces

abstract: We examine volume pinching problems of CAT(1) spaces. We characterize a class of compact geodesically complete CAT(1) spaces of small specific volume. We prove a sphere theorem for compact CAT(1) homology manifolds of small volume. We also formulate a criterion of manifold recognition for homology manifolds on volume growths under an upper curvature bound.


title: A local instability mechanism of the Navier-Stokes flow with swirl on the no-slip flat boundary

abstract: Using numerical simulations of the axisymmetric Navier-Stokes equations with swirl on a no-slip flat boundary, Hsu-Notsu-Yoneda [J. Fluid Mech. 2016] observed the creation of a high-vorticity region on the boundary near the axis of symmetry. In this paper, using a differential geometric approach, we prove that such flows indeed have a destabilizing effect, which is formulated in terms of a lower bound on the $L^\infty$-norm of derivatives of the velocity field on the boundary.


title: Spin-orbital excitons and their potential condensation in pentavalent iridates

abstract: We investigate magnetic excitations in iridium insulators with pentavalent Ir$^{5+}$ ($5d^4$) ions with strong spin-orbit coupling. We obtain a microscopic model based on the local Ir$^{5+}$ multiplets involving $J=0$ (singlet), $J=1$ (triplet), and $J=2$ (quintet) spin-orbital states. We get effective interactions between these multiplets on square and face-centered-cubic (fcc) structures of magnetic ions in the layered-perovskites and the double-perovkites, in particular Ba$_2$YIrO$_6$. Further, we derive an effective spin-orbital Hamiltonian in terms of bond bosons and explore possible instabilities towards magnetic and quadrupole orderings. Additionally, we study charge excitations with help of the variational cluster perturbation theory and calculate the electronic charge gap as a function of hopping and Coulomb interactions. Based on both electronic and magnetic phase diagrams, we verify the possibility of excitonic magnetism due to condensation of spin-orbital excitons in Ir$^{5+}$ iridates.


title: Visual Attention Network for Low Dose CT

abstract: Noise and artifacts are intrinsic to low dose CT (LDCT) data acquisition, and will significantly affect the imaging performance. Perfect noise removal and image restoration is intractable in the context of LDCT due to the statistical and technical uncertainties. In this paper, we apply the generative adversarial network (GAN) framework with a visual attention mechanism to deal with this problem in a data-driven/machine learning fashion. Our main idea is to inject visual attention knowledge into the learning process of GAN to provide a powerful prior of the noise distribution. By doing this, both the generator and discriminator networks are empowered with visual attention information so they will not only pay special attention to noisy regions and surrounding structures but also explicitly assess the local consistency of the recovered regions. Our experiments qualitatively and quantitatively demonstrate the effectiveness of the proposed method with clinic CT images.


title: On the 2D Dirac oscillator in the presence of vector and scalar potentials in the cosmic string spacetime in the context of spin and pseudospin symmetries

abstract: The Dirac equation with both scalar and vector couplings describing the dynamics of a two-dimensional Dirac oscillator in the cosmic string spacetime is considered. We derive the Dirac-Pauli equation and solve it in the limit of the spin and the pseudo-spin symmetries. We analyze the presence of cylindrical symmetric scalar potentials which allows us to provide analytic solutions for the resultant field equation. By using an appropriate ansatz, we find that the radial equation is a biconfluent Heun-like differential equation. The solution of this equation provides us with more than one expression for the energy eigenvalues of the oscillator. We investigate these energies and find that there is a quantum condition between them. We study this condition in detail and find that it requires the fixation of one of the physical parameters involved in the problem. Expressions for the energy of the oscillator are obtained for some values of the quantum number $n$. Some particular cases which lead to known physical systems are also addressed.


title: Cross-over from trion-hole to exciton-polaron in n-doped semiconductor quantum wells

abstract: We present a theoretical study of photo-absorption in n-doped two-dimensional (2D) and quasi-2D semiconductors that takes into account the interaction of the photocreated exciton with Fermi-sea (FS) electrons through (i) Pauli blocking, (ii) Coulomb screening, and (iii) excitation of FS electron-hole pairs---that we here restrict to one. The system we tackle is thus made of one exciton plus zero or one FS electron-hole pair. At low doping, the system ground state is predominantly made of a "trion-hole"---a trion (two opposite-spin electrons plus a valence hole) weakly bound to a FS hole---with a small exciton component. As the trion is poorly coupled to photon, the intensity of the lowest absorption peak is weak; it increases with doping, thanks to the growing exciton component, due to a larger coupling between 2-particle and 4-particle states. Under a further doping increase, the trion-hole complex is less bound because of Pauli blocking by FS electrons, and its energy increases. The lower peak then becomes predominantly due to an exciton dressed by FS electron-hole pairs, that is, an exciton-polaron. As a result, the absorption spectra of $n$-doped semiconductor quantum wells show two prominent peaks, the nature of the lowest peak turning from trion-hole to exciton-polaron under a doping increase. Our work also nails down the physical mechanism behind the increase with doping of the energy separation between the trion-hole peak and the exciton-polaron peak, even before the anti-crossing, as experimentally observed.


title: Birds of a Feather Flock Together? A Study of Developers' Flocking and Migration Behavior in GitHub and Stack Overflow

abstract: Interactions between individuals and their participation in community activities are governed by how individuals identify themselves with their peers. We want to investigate such behavior for developers while they are learning and contributing on socially collaborative environments, specifically code hosting sites and question/answer sites. In this study, we investigate the following questions about advocates, developers who can be identified as well-rounded community contributors and active learners. Do advocates flock together in a community? How do flocks of advocates migrate within a community? Do these flocks of advocates migrate beyond a single community? To understand such behavior, we collected 12,578 common advocates across a code hosting site - GitHub and a question/answering site - Stack Overflow. These advocates were involved in 1,549 projects on GitHub and were actively asking 114,569 questions and responding 408,858 answers and 1,001,125 comments on Stack Overflow. We performed an in-depth empirical analysis using social networks to find the flocks of advocates and their migratory pattern on GitHub, Stack Overflow, and across both communities. We found that 7.5% of the advocates create flocks on GitHub and 8.7% on Stack Overflow. Further, these flocks of advocates migrate on an average of 5 times on GitHub and 2 times on Stack Overflow. In particular, advocates in flocks of two migrate more frequently than larger flocks. However, this migration behavior was only common within a single community. Our findings indicate that advocates' flocking and migration behavior differs substantially from the ones found in other social environments. This suggests a need to investigate the factors that demotivate the flocking and migration behavior of advocates and ways to enhance and integrate support for such behavior in collaborative software tools.


title: On the number of total prime factors of an odd perfect number

abstract: Let $N$ be an odd perfect number. Let $\omega(N)$ be the number of distinct prime factors of $N$ and let $\Omega(N)$ be the total number of prime factors of $N$. We prove that if $(3,N)=1$, then $ \frac{302}{113}\omega - \frac{286}{113} \leq \Omega. $ If $3|N$, then $\frac{66}{25}\omega-5\leq\Omega.$ This is an improvement on similar prior results by the author which was an improvement of a result of Ochem and Rao. We also establish new lower bounds on $\omega(N)$ in terms of the smallest prime factor of $N$ and establish new lower bounds on $N$ in terms of its smallest prime factor.


title: Annual parallax measurements of a semiregular variable star SV Pegasus with VERA

abstract: Many studies have shown that there are clear sequences in the period-luminosity relationship (PLR) for Mira variables and semiregular variables (SRVs) in the Large Magellanic Cloud (LMC). To investigate the PLR for SRVs in our galaxy, we examined the annual parallax measurement and conducted K'-band photometric monitoring of an SRV star SV Pegasus (SV Peg). We measured the position change of the associating H$_2$O maser spots by phase-referencing VLBI observations with VERA at 22 GHz, spanning approximately 3 yr, and detected an annual parallax of $\pi = 3.00 \pm 0.06$ mas, corresponding to a distance of $D=333 \pm 7$ pc. This result is in good agreement with the Hipparcos parallax and improves the accuracy of the distance from 35 $\%$ to 2 $\%$. However, the GAIA DR2 catalog gave a parallax of $\pi=1.12\pm0.28$ mas for SV Peg. This indicates that the GAIA result might be blurred by the effect of the stellar size because the estimated stellar radius was $\sim 5$ mas, which is comparable to the parallax. We obtained a K'-band mean magnitude of $m_{K'} = -0.48$ mag and a period of $P=177$ days from our photometric monitoring with a 1-m telescope. Using the trigonometric distance, we derived an absolute magnitude of $M_{K'}=-8.09 \pm 0.05$ mag. This result shows that the position of SV Peg in the PLR falls on the C' sequence found in the PLR in the LMC, which is similar to other SRVs in our galaxy.


title: Theory of chiral effects in magnetic textures with spin-orbit coupling

abstract: We present a theoretical study of two-dimensional spatially and temporally varying magnetic textures in the presence of spin-orbit coupling (SOC) of both the Rashba and Dresselhaus types. We show that the effective gauge field due to these SOCs, contributes to the dissipative and reactive spin torques in exactly the same way as in electromagnetism. Our calculations reveal that Rashba (Dresselhaus) SOC induces a chiral dissipation in interfacial (bulk) inversion asymmetric magnetic materials. Furthermore, we show that in addition to chiral dissipation $\alpha_c$, these SOCs also produce a chiral renormalization of the gyromagnetic ratio $\tilde{\gamma}_c$, and show that the latter is intrinsically linked to the former via a simple relation $\alpha_c = (\tau/\tau_{\rm ex}) \tilde{\gamma}_c$, where $\tau_{\rm ex}$ and $\tau$ are the exchange time and the electron relaxation time, respectively. Finally, we propose a theoretical scheme based on the Scattering theory to calculate and investigate the properties of damping in chiral magnets. Our findings should in principle provide a guide for material engineering of effects related to chiral dynamics in magnetic textures with SOC.


title: Connecting the Dots: Identifying Network Structure via Graph Signal Processing

abstract: Network topology inference is a prominent problem in Network Science. Most graph signal processing (GSP) efforts to date assume that the underlying network is known, and then analyze how the graph's algebraic and spectral characteristics impact the properties of the graph signals of interest. Such an assumption is often untenable beyond applications dealing with e.g., directly observable social and infrastructure networks; and typically adopted graph construction schemes are largely informal, distinctly lacking an element of validation. This tutorial offers an overview of graph learning methods developed to bridge the aforementioned gap, by using information available from graph signals to infer the underlying graph topology. Fairly mature statistical approaches are surveyed first, where correlation analysis takes center stage along with its connections to covariance selection and high-dimensional regression for learning Gaussian graphical models. Recent GSP-based network inference frameworks are also described, which postulate that the network exists as a latent underlying structure, and that observations are generated as a result of a network process defined in such a graph. A number of arguably more nascent topics are also briefly outlined, including inference of dynamic networks, nonlinear models of pairwise interaction, as well as extensions to directed graphs and their relation to causal inference. All in all, this paper introduces readers to challenges and opportunities for signal processing research in emerging topic areas at the crossroads of modeling, prediction, and control of complex behavior arising in networked systems that evolve over time.


title: Compression Performance of Grayscale-based Image Encryption for Encryption-then-Compression Systems

abstract: This paper considers a new grayscale-based image encryption for Encryption-then-Compression (EtC) systems with JPEG compression. Firstly, generation methods of grayscale-based images are discussed in terms of the selection of color space. In addition, a new JPEG quantization table for the grayscale-based images is proposed to provide a better compression performance. Moreover, the quality of both images uploaded to Social Network Services (SNS) and downloaded from SNS, are discussed and evaluated. In the experiments, encrypted images are compressed using various compression parameters and quantization tables, and uploaded to Twitter and Facebook. The results proved that the selection of color space and the proposed quantization table can improve the compression performances of not only uploaded images but also downloaded ones.


title: Symbiotic Radio: A New Communication Paradigm for Passive Internet-of-Things

abstract: In this paper, a novel technique, called symbiotic radio (SR), is proposed for passive Internet-of-Things (IoT), in which a backscatter device (BD) is integrated with a primary transmission. The primary transmitter is designed to assist the primary and BD transmissions, and the primary receiver decodes the information from the primary transmitter as well as the BD. We consider a multiple-input single-output (MISO) SR and the symbol period for BD transmission is designed to be either the same as or much longer than that of the primary system, resulting in parasitic or commensal relationship between the primary and BD transmissions. We first derive the achievable rates for the primary system and the BD transmission. Then, we formulate two transmit beamforming optimization problems, i.e., the weighted sum-rate maximization problem and the transmit power minimization problem, and solve these non-convex problems by applying semi-definite relaxation technique. In addition, a novel transmit beamforming structure is proposed to reduce the computational complexity of the solutions. Simulation results show that when the BD transmission rate is properly designed, the proposed SR not only enables the opportunistic transmission for the BD via energy-efficient passive backscattering, but also enhances the achievable rate of the primary system by properly exploiting the additional signal path from the BD.


title: Dynamic Assortment Optimization with Changing Contextual Information

abstract: In this paper, we study the dynamic assortment optimization problem under a finite selling season of length $T$. At each time period, the seller offers an arriving customer an assortment of substitutable products under a cardinality constraint, and the customer makes the purchase among offered products according to a discrete choice model. Most existing work associates each product with a real-valued fixed mean utility and assumes a multinomial logit choice (MNL) model. In many practical applications, feature/contexutal information of products is readily available. In this paper, we incorporate the feature information by assuming a linear relationship between the mean utility and the feature. In addition, we allow the feature information of products to change over time so that the underlying choice model can also be non-stationary. To solve the dynamic assortment optimization under this changing contextual MNL model, we need to simultaneously learn the underlying unknown coefficient and makes the decision on the assortment. To this end, we develop an upper confidence bound (UCB) based policy and establish the regret bound on the order of $\widetilde O(d\sqrt{T})$, where $d$ is the dimension of the feature and $\widetilde O$ suppresses logarithmic dependence. We further established the lower bound $\Omega(d\sqrt{T}/K)$ where $K$ is the cardinality constraint of an offered assortment, which is usually small. When $K$ is a constant, our policy is optimal up to logarithmic factors. In the exploitation phase of the UCB algorithm, we need to solve a combinatorial optimization for assortment optimization based on the learned information. We further develop an approximation algorithm and an efficient greedy heuristic. The effectiveness of the proposed policy is further demonstrated by our numerical studies.


title: Special value formula for the twisted triple product $L$-function and an application to the restricted $L^2$-norm problem

abstract: We establish explicit Ichino's formulae for the central values of the triple product $L$-functions with emphasis on the calculations for the real place. The key ingredient for our computations is Proposition 6.8 which generalizes a result of Michel-Venkatesh. As an application we prove the optimal upper bound of a sum of restricted $L^2$-norms of the $L^2$-normalized newforms on certain quadratic extensions with prime level and bounded spectral parameter following the methods of Blomer.


title: Gluing maps and cobordism maps for sutured monopole Floer homology

abstract: In this paper we construct gluing maps and cobordism maps for sutured monopole Floer homology.


title: Formal Verification of Neural Network Controlled Autonomous Systems

abstract: In this paper, we consider the problem of formally verifying the safety of an autonomous robot equipped with a Neural Network (NN) controller that processes LiDAR images to produce control actions. Given a workspace that is characterized by a set of polytopic obstacles, our objective is to compute the set of safe initial conditions such that a robot trajectory starting from these initial conditions is guaranteed to avoid the obstacles. Our approach is to construct a finite state abstraction of the system and use standard reachability analysis over the finite state abstraction to compute the set of the safe initial states. The first technical problem in computing the finite state abstraction is to mathematically model the imaging function that maps the robot position to the LiDAR image. To that end, we introduce the notion of imaging-adapted sets as partitions of the workspace in which the imaging function is guaranteed to be affine. We develop a polynomial-time algorithm to partition the workspace into imaging-adapted sets along with computing the corresponding affine imaging functions. Given this workspace partitioning, a discrete-time linear dynamics of the robot, and a pre-trained NN controller with Rectified Linear Unit (ReLU) nonlinearity, the second technical challenge is to analyze the behavior of the neural network. To that end, we utilize a Satisfiability Modulo Convex (SMC) encoding to enumerate all the possible segments of different ReLUs. SMC solvers then use a Boolean satisfiability solver and a convex programming solver and decompose the problem into smaller subproblems. To accelerate this process, we develop a pre-processing algorithm that could rapidly prune the space feasible ReLU segments. Finally, we demonstrate the efficiency of the proposed algorithms using numerical simulations with increasing complexity of the neural network controller.


title: Self-cleaning of hydrophobic rough surfaces by coalescence-induced wetting transition

abstract: The superhydrophobic leaves of a lotus plant and other natural surfaces with self-cleaning function have been studied intensively for the development of artificial biomimetic surfaces. Surface roughness generated by hierarchical structures is a crucial property required for superhydrophobicity and self-cleaning. Here, we demonstrate a novel self-cleaning mechanism of textured surfaces attributed to a spontaneous coalescence-induced wetting transition. We focus on the wetting transition as it represents a new mechanism, which can explain why droplets on rough surfaces are able to change from the highly adhesive Wenzel state to the low-adhesion Cassie-Baxter state and achieve self-cleaning. In particular, we perform many-body dissipative particle dynamics simulations of liquid droplets sitting on mechanically textured substrates. We quantitatively investigate the wetting behavior of an isolated droplet as well as coalescence of droplets for both Cassie-Baxter and Wenzel states. Our simulation results reveal that droplets in the Cassie-Baxter state have much lower contact angle hysteresis and smaller hydrodynamic resistance than droplets in the Wenzel state. When small neighboring droplets coalesce into bigger ones on textured hydrophobic substrates, we observe a spontaneous wetting transition from a Wenzel state to a Cassie-Baxter state, which is powered by the surface energy released upon coalescence of the droplets. For superhydrophobic surfaces, the released surface energy may be sufficient to cause a jumping motion of droplets off the surface, in which case adding one more droplet to coalescence may increase the jumping velocity by one order of magnitude. When multiple droplets are involved, we find that the spatial distribution of liquid components in the coalesced droplet can be controlled by properly designing the overall arrangement of droplets and the distance between them.


title: Angular momentum evolution of galaxies: the perspective of hydrodynamical simulations

abstract: Until a decade ago, galaxy formation simulations were unable to simultaneously reproduce the observed angular momentum (AM) of galaxy disks and bulges. Improvements in the interstellar medium and stellar feedback modelling, together with advances in computational capabilities, have allowed the current generation of cosmological galaxy formation simulations to reproduce the diversity of AM and morphology that is observed in local galaxies. In this review I discuss where we currently stand in this area from the perspective of hydrodynamical simulations, specifically how galaxies gain their AM, and the effect galaxy mergers and gas accretion have on this process. I discuss results which suggest that a revision of the classical theory of disk formation is needed, and finish by discussing what the current challenges are.


title: Physics Guided RNNs for Modeling Dynamical Systems: A Case Study in Simulating Lake Temperature Profiles

abstract: This paper proposes a physics-guided recurrent neural network model (PGRNN) that combines RNNs and physics-based models to leverage their complementary strengths and improve the modeling of physical processes. Specifically, we show that a PGRNN can improve prediction accuracy over that of physical models, while generating outputs consistent with physical laws, and achieving good generalizability. Standard RNNs, even when producing superior prediction accuracy, often produce physically inconsistent results and lack generalizability. We further enhance this approach by using a pre-training method that leverages the simulated data from a physics-based model to address the scarcity of observed data. The PGRNN has the flexibility to incorporate additional physical constraints and we incorporate a density-depth relationship. Both enhancements further improve PGRNN performance. Although we present and evaluate this methodology in the context of modeling the dynamics of temperature in lakes, it is applicable more widely to a range of scientific and engineering disciplines where mechanistic (also known as process-based) models are used, e.g., power engineering, climate science, materials science, computational chemistry, and biomedicine.


title: A framework for automated anomaly detection in high frequency water-quality data from in situ sensors

abstract: River water-quality monitoring is increasingly conducted using automated in situ sensors, enabling timelier identification of unexpected values. However, anomalies caused by technical issues confound these data, while the volume and velocity of data prevent manual detection. We present a framework for automated anomaly detection in high-frequency water-quality data from in situ sensors, using turbidity, conductivity and river level data. After identifying end-user needs and defining anomalies, we ranked their importance and selected suitable detection methods. High priority anomalies included sudden isolated spikes and level shifts, most of which were classified correctly by regression-based methods such as autoregressive integrated moving average models. However, using other water-quality variables as covariates reduced performance due to complex relationships among variables. Classification of drift and periods of anomalously low or high variability improved when we applied replaced anomalous measurements with forecasts, but this inflated false positive rates. Feature-based methods also performed well on high priority anomalies, but were also less proficient at detecting lower priority anomalies, resulting in high false negative rates. Unlike regression-based methods, all feature-based methods produced low false positive rates, but did not and require training or optimization. Rule-based methods successfully detected impossible values and missing observations. Thus, we recommend using a combination of methods to improve anomaly detection performance, whilst minimizing false detection rates. Furthermore, our framework emphasizes the importance of communication between end-users and analysts for optimal outcomes with respect to both detection performance and end-user needs. Our framework is applicable to other types of high frequency time-series data and anomaly detection applications.


title: Lagrangian densities of hypergraph cycles

abstract: The Lagrangian density of an $r$-uniform hypergraph $F$ is $r!$ multiplying the supremum of the Lagrangians of all $F$-free $r$-uniform hypergraphs. For an $r$-graph $H$ with $t$ vertices, it is clear that $\pi_{\lambda}(H)\ge r!\lambda{(K_{t-1}^r)}$. We say that an $r$-unform hypergraph $H$ with $t$ vertices is perfect if $\pi_{\lambda}(H)= r!\lambda{(K_{t-1}^r)}$. A theorem of Motzkin-Straus implies that all $2$-uniform graphs are perfect. It is interesting to explore what kind of hypergraphs are perfect. A hypergraph is linear if any 2 edges have at most 1 vertex in common. We propose the following conjecture: (1) For $r\ge 3$, there exists $n$ such that a linear $r$-unofrm hypergraph with at least $n$ vertices is perfect. (2) For $r\ge 3$, there exists $n$ such that if $G, H$ are perfect $r$-uniform hypergraphs with at least $n$ vertices, then $G\bigsqcup H$ is perfect. Regarding this conjecture, we obtain a partial result: Let $S_{2,t}=\{123,124,125,126,...,12(t+2)\}$. (An earlier result of Sidorenko states that $S_{2,t}$ is perfect \cite{Sidorenko-89}.) Let $H$ be a perfect $3$-graph with $s$ vertices. Then $F=S_{2,t}\bigsqcup H$ is perfect if $s\geq 3$ and $t\geq 3$.


title: A Large-scale Study of Social Media Sources in News Articles

abstract: In this study, we closely look at the use of social media contents as source or reference in the U.S. news media. Specifically, we examine about 60 thousand news articles published within the 5 years period of 2013-2017 by 153 U.S. media outlets and analyze use of social media content as source compared to other sources. We designed a social media source extraction algorithm and investigated the extent and nature of social media source usage across different news topics. Our results show that uses of social media content in news almost doubled in five years. Unreliable media outlets rely on social media more than the mainstream media. Both mainstream and unreliable sites prefer Twitter to Facebook as a source of information.


title: Anisotropy of Dirac cones and van Hove singularity in an organic Dirac fermion system

abstract: We propose an experimental method to examine the in-plane anisotropy of electronic structure in layered conductors. In the method, we measure the interlayer magnetoresistance as a function of in-plane magnetic field orientation. We applied it to an organic Dirac fermion system a-(BEDT-TTF)2I3 to experimentally determine the orientation of the anisotropic Dirac cones. It is concluded that the long axis of the elliptic constant-energy contours of the Dirac cone is tilted by approximately -30 deg from the crystalline a-axis to b-axis under hydrostatic pressures. Additionally, we observed a signature of van Hove singularity (which is a saddle point of the band dispersion) at 30-40 K above or below the Dirac point. The ridgeline of the saddle point is estimated as almost parallel to the crystalline b-axis.


title: A new gap for complete hypersurfaces with constant mean curvature in space forms

abstract: Let $M$ be an $n$-dimensional closed hypersurface with constant mean curvature and constant scalar curvature in an unit sphere. Denote by $H$ and $S$ the mean curvature and the squared length of the second fundamental form respectively. We prove that if $S > \alpha (n, H)$, where $n\geq 4$ and $H\neq 0$, then $S > \alpha (n, H) + B_n\frac{n H^2}{n - 1}$. Here \[ \alpha (n, H) = n + \frac{n^3}{2 (n - 1)} H^2 - \frac{n (n - 2)}{2 (n - 1)}\sqrt{n^2 H^4 + 4 (n - 1) H^2}, \] $B_n=\frac{1}{5}$ for $4\leq n \leq 20$, and $B_n=\frac{49}{250}$ for $n>20$. Moreover, we obtain a gap theorem for complete hypersurfaces with constant mean curvature and constant scalar curvature in space forms.


title: Physical and Chemical Conditions of the Protostellar Envelope and the Protoplanetary Disk in HL Tau

abstract: We report our SMA observations of the Class \RN{1}-\RN{2} protostar HL Tau in the $^{13}$CO (2--1), C$^{18}$O (2--1), SO(5$_6$--4$_5$), and the 1.3 mm dust-continuum emission and our analyses of the ALMA long baseline data of the HCO$^{+}$ (1--0) emission. The 1.3 mm continuum emission observed with the SMA shows compact ($\sim$0$\farcs$8 $\times$ 0$\farcs$5) and extended ($\sim$6$\farcs$5 $\times$ 4$\farcs$3) components, tracing the protoplanetary disk and the protostellar envelope, respectively. The $^{13}$CO, C$^{18}$O, and HCO$^+$ show a compact ($\sim$ 200 AU) component at velocities higher than 3 km s$^{-1}$ from the systemic velocity and an extended ($\sim$ 1000 AU) component at the lower velocities. The high-velocity component traces the Keplerian rotating disk, and the low-velocity component traces the infalling envelope. The HCO$^+$ high-velocity component is fitted with a Keplerian disk model with a central stellar mass of 1.4 $M_{\odot}$. The radial intensity profiles of the $^{13}$CO and C$^{18}$O along the disk major axis are fitted with a disk+envelope model, and the gas masses of the disk and envelope are estimated to be $2\mbox{--}40\times10^{-4}$ $M_{\odot}$ and $2.9\times10^{-3}$ $M_{\odot}$, respectively. The disk dust mass has been estimated to be $1\mbox{--}3 \times 10^{-3}$ $M_{\odot}$ in the literature. Thus, our estimated disk gas mass suggests that the gas-to-dust mass ratio in the disk is $<$10, a factor of ten lower than the estimated ratio in the envelope. We discuss the possible gas depletion or CO depletion in the planet-forming candidate HL Tau in the context of disk and envelope evolution.


title: Efficient Humanoid Contact Planning using Learned Centroidal Dynamics Prediction

abstract: Humanoid robots dynamically navigate an environment by interacting with it via contact wrenches exerted at intermittent contact poses. Therefore, it is important to consider dynamics when planning a contact sequence. Traditional contact planning approaches assume a quasi-static balance criterion to reduce the computational challenges of selecting a contact sequence over a rough terrain. This however limits the applicability of the approach when dynamic motions are required, such as when walking down a steep slope or crossing a wide gap. Recent methods overcome this limitation with the help of efficient mixed integer convex programming solvers capable of synthesizing dynamic contact sequences. Nevertheless, its exponential-time complexity limits its applicability to short time horizon contact sequences within small environments. In this paper, we go beyond current approaches by learning a prediction of the dynamic evolution of the robot centroidal momenta, which can then be used for quickly generating dynamically robust contact sequences for robots with arms and legs using a search-based contact planner. We demonstrate the efficiency and quality of the results of the proposed approach in a set of dynamically challenging scenarios.


title: GraphIE: A Graph-Based Framework for Information Extraction

abstract: Most modern Information Extraction (IE) systems are implemented as sequential taggers and only model local dependencies. Non-local and non-sequential context is, however, a valuable source of information to improve predictions. In this paper, we introduce GraphIE, a framework that operates over a graph representing a broad set of dependencies between textual units (i.e. words or sentences). The algorithm propagates information between connected nodes through graph convolutions, generating a richer representation that can be exploited to improve word-level predictions. Evaluation on three different tasks --- namely textual, social media and visual information extraction --- shows that GraphIE consistently outperforms the state-of-the-art sequence tagging model by a significant margin.


title: Provably Accelerated Randomized Gossip Algorithms

abstract: In this work we present novel provably accelerated gossip algorithms for solving the average consensus problem. The proposed protocols are inspired from the recently developed accelerated variants of the randomized Kaczmarz method - a popular method for solving linear systems. In each gossip iteration all nodes of the network update their values but only a pair of them exchange their private information. Numerical experiments on popular wireless sensor networks showing the benefits of our protocols are also presented.


title: Local-in-time Solvability and Space Analyticity for the Navier-Stokes Equations with BMO-type Initial Data

abstract: It is proved that there exists a local-in-time solution $u\in C([0,T),bmo(\mathbb{R}^d)^d)$ of the Navier-Stokes equations such that every $u(t)$ has an analytic extension on a complex domain whose size only depends on $t$ (and increases with $t$) and the external force $f$, assuming only that the initial velocity $u_0$ is a local $BMO$ function. Our method for proving is a combination and refinement of the work by Gruji\'c and Kukavica [13], Guberovi\'c [15] and Kozono et al. [18]. One challenging step is the estimation of the heat and Stokes semigroups from $BMO$-type spaces to $L^\infty$; a result itself of independent interest. We also apply the idea to the analyticity of vorticity with the assistance of Calder\'on-Zygmund theory.


title: Updated determination of $N^*$ resonance parameters using a unitary, multichannel formalism

abstract: Results are presented for an updated multichannel energy-dependent partial-wave analysis of $\pi N$ scattering. Our earlier work incorporated single-energy amplitudes for $\pi N \rightarrow \pi N$, $\gamma N \rightarrow \pi N$, $\pi N \rightarrow \pi\pi N$, $\pi N \rightarrow \eta N$, and $\pi N \rightarrow K \Lambda$. The present work incorporates new single-energy solutions for $\gamma p \rightarrow \eta p$ up to a c.m.\ energy of 1990~MeV, $\gamma p \rightarrow K^+ \Lambda$ up to a c.m.\ energy of 2230~MeV, and $\gamma n \rightarrow \eta n$ up to a c.m.\ energy of 1885~MeV, as well as updated single-energy solutions for $\pi N \rightarrow \eta N$, $\pi N \rightarrow K \Lambda$, and $\gamma N \rightarrow \pi N$. In this paper, we present and discuss the resonance parameters obtained from a combined fit of all these single-energy amplitudes. Our determined energy-dependent amplitudes provide an excellent description of the corresponding measured observables.


title: Multirobot Coordination with Counting Temporal Logics

abstract: In many multirobot applications, planning trajectories in a way to guarantee that the collective behavior of the robots satisfies a certain high-level specification is crucial. Motivated by this problem, we introduce counting temporal logics---formal languages that enable concise expression of multirobot task specifications over possibly infinite horizons. We first introduce a general logic called counting linear temporal logic plus (cLTL+), and propose an optimization-based method that generates individual trajectories such that satisfaction of a given cLTL+ formula is guaranteed when these trajectories are synchronously executed. We then introduce a fragment of cLTL+, called counting linear temporal logic (cLTL), and show that a solution to planning problem with cLTL constraints can be obtained more efficiently if all robots have identical dynamics. In the second part of the paper, we relax the synchrony assumption and discuss how to generate trajectories that can be asynchronously executed, while preserving the satisfaction of the desired cLTL+ specification. In particular, we show that when the asynchrony between robots is bounded, the method presented in this paper can be modified to generate robust trajectories. We demonstrate these ideas with an experiment and provide numerical results that showcase the scalability of the method.


title: Attention-based sequence-to-sequence model for speech recognition: development of state-of-the-art system on LibriSpeech and its application to non-native English

abstract: Recent research has shown that attention-based sequence-to-sequence models such as Listen, Attend, and Spell (LAS) yield comparable results to state-of-the-art ASR systems on various tasks. In this paper, we describe the development of such a system and demonstrate its performance on two tasks: first we achieve a new state-of-the-art word error rate of 3.43% on the test clean subset of LibriSpeech English data; second on non-native English speech, including both read speech and spontaneous speech, we obtain very competitive results compared to a conventional system built with the most updated Kaldi recipe.


title: Quantum coherence control at near 1000 K

abstract: Quantum coherence control usually requires extremely low temperature environments. Even for spins in diamond, a remarkable exception, the coherence signal is lost as temperature approaches 700 K. Here we demonstrate quantum coherence control of the electron spins of nitrogen-vacancy centers in nanodiamonds at temperatures near 1000 K. The scheme is based on initialization and readout of the spins at room temperature and control at high temperature, which is enabled by pulse laser heating and rapid diffusion cooling of nanodiamonds on amorphous carbon films. Using high-temperature spin control, we observe the magnetic phase transition of a single nickel nanoparticle at about 615 K. This work enables nano-thermometry and nano-magnetometry in the high-temperature regime.


title: Two-loop perturbative corrections to the constrained effective potential in thermal QCD

abstract: In this paper, we compute the constrained QCD effective potential up to two-loop order with finite quark mass and chemical potential. We present the explicit calculations by using the double line notation and analytical expressions for massless quarks are obtained in terms of the Bernoulli polynomials or Polyakov loops. Our results explicitly show that the constrained QCD effective potential is independent on the gauge fixing parameter. In addition, as compared to the massless case, the constrained QCD effective potential with massive quarks develops a completely new term which is only absent when the background field vanishes. Furthermore, we discuss the relation between the one- and two-loop constrained effective potential. The surprisingly simple proportionality that exists in the pure gauge theories, however, is in general no longer true when fermions are taken into account. On the other hand, for high baryon density $\mu_B$ and low temperature $T$, in the massless limit, we do also find a similar proportionality between the one- and two-loop fermionic contributions in the constrained effective potential up to ${\cal O}(T/\mu_B)$.


title: Towards End-to-End Code-Switching Speech Recognition

abstract: Code-switching speech recognition has attracted an increasing interest recently, but the need for expert linguistic knowledge has always been a big issue. End-to-end automatic speech recognition (ASR) simplifies the building of ASR systems considerably by predicting graphemes or characters directly from acoustic input. In the mean time, the need of expert linguistic knowledge is also eliminated, which makes it an attractive choice for code-switching ASR. This paper presents a hybrid CTC-Attention based end-to-end Mandarin-English code-switching (CS) speech recognition system and studies the effect of hybrid CTC-Attention based models, different modeling units, the inclusion of language identification and different decoding strategies on the task of code-switching ASR. On the SEAME corpus, our system achieves a mixed error rate (MER) of 34.24%.


title: Roles of photospheric motions and flux emergence in the major solar eruption on 2017 September 6

abstract: We study the magnetic field evolution in the active region (AR) 12673 that produced the largest solar flare in the past decade on 2017 September 6. Fast flux emergence is one of the most prominent features of this AR. We calculate the magnetic helicity from photospheric tangential flows that shear and braid field lines (shear-helicity), and from normal flows that advect twisted magnetic flux into the corona (emergence-helicity), respectively. Our results show that the emergence-helicity accumulated in the corona is $-1.6\times10^{43}~Mx^2$ before the major eruption, while the shear-helicity accumulated in the corona is $-6\times10^{43}~Mx^2$, which contributes about 79\% of the total helicity. The shear-helicity flux is dominant throughout the overall investigated emergence phase. Our results imply that the emerged fields initially contain relatively low helicity. Much more helicity is built up by shearing and converging flows acting on preexisted and emerging flux. Shearing motions are getting stronger with the flux emergence, and especially on both sides of the polarity inversion line of the core field region. The evolution of the vertical currents shows that most of the intense currents do not appear initially with the emergence of the flux, which implies that most of the emerging flux is probably not strongly current-carrying. The helical magnetic fields (flux rope) in the core field region are probably formed by long-term photospheric motions. The shearing and converging motions are continuously generated driven by the flux emergence. AR 12673 is a representative as photospheric motions contribute most of the nonpotentiality in the AR with vigorous flux emergence.


title: Sharp inequalities for the numerical radius of block operator matrices

abstract: In this paper, we present several sharp upper bounds for the numerical radii of the diagonal and off-diagonal parts of the $2\times2$ block operator matrix $\begin{bmatrix}A&B\\ C&D\end{bmatrix}$. Among extensions of some results of Kittaneh et al., it is shown that if $T=\begin{bmatrix}A&0\\ 0&D\end{bmatrix}$, and $f$ and $g$ are non-negative continuous functions on $[0,\infty)$ such that $f(t)g(t)=t\,\,(t\geq 0)$, then for all nonnegative nondecreasing convex functions $h$ on $[0,\infty)$ , we obtain that \begin{align*}h\left(w^r(T)\right)\leq \max\left(\left\|\frac{1}{p}h\left(f^{pr}(\left|A\right|)\right)+ \frac{1}{q}h\left(g^{qr}(\left|A^*\right|)\right)\right\|, \left\|\frac{1}{p}h\left(f^{pr}(\left|D\right|)\right)+ \frac{1}{q}h\left(g^{qr}(\left|D^*\right|)\right)\right\|\right), \end{align*} where $p, q>1$ with $\frac{1}{p}+\frac{1}{q}=1$ and $r\min(p,q)\geq 2$.


title: Sample size considerations for comparing dynamic treatment regimens in a sequential multiple-assignment randomized trial with a continuous longitudinal outcome

abstract: Clinicians and researchers alike are increasingly interested in how best to personalize interventions. A dynamic treatment regimen (DTR) is a sequence of pre-specified decision rules which can be used to guide the delivery of a sequence of treatments or interventions that are tailored to the changing needs of the individual. The sequential multiple-assignment randomized trial (SMART) is a research tool which allows for the construction of effective DTRs. We derive easy-to-use formulae for computing the total sample size for three common two-stage SMART designs in which the primary aim is to compare mean end-of-study outcomes for two embedded DTRs which recommend different first-stage treatments. The formulae are derived in the context of a regression model which leverages information from a longitudinal outcome collected over the entire study. We show that the sample size formula for a SMART can be written as the product of the sample size formula for a standard two-arm randomized trial, a deflation factor that accounts for the increased statistical efficiency resulting from a longitudinal analysis, and an inflation factor that accounts for the design of a SMART. The SMART design inflation factor is typically a function of the anticipated probability of response to first-stage treatment. We review modeling and estimation for DTR effect analyses using a longitudinal outcome from a SMART, as well as the estimation of standard errors. We also present estimators for the covariance matrix for a variety of common working correlation structures. Methods are motivated using the ENGAGE study, a SMART aimed at developing a DTR for increasing motivation to attend treatments among alcohol- and cocaine-dependent patients.


title: Magnetohydrodynamic Simulation of the X9.3 Flare on 2017 September 6: Evolving Magnetic Topology

abstract: Three-dimensional magnetic topology is crucial to understanding the explosive release of magnetic energy in the corona during solar flares. Much attention has been given to the pre-flare magnetic topology to identify candidate sites of magnetic reconnection, yet it is unclear how the magnetic reconnection and its attendant topological changes shape the eruptive structure and how the topology evolves during the eruption. Here we employed a realistic, data-constrained magnetohydrodynamic simulation to study the evolving magnetic topology for an X9.3 eruptive flare that occurred on 2017 September 6. The simulation successfully reproduces the eruptive features and processes in unprecedented detail. The numerical results reveal that the pre-flare corona contains multiple twisted flux systems with different connections, and during the eruption, these twisted fluxes form a coherent flux rope through tether-cutting-like magnetic reconnection below the rope. Topological analysis shows that the rising flux rope is wrapped by a quasi-separatrix layer, which intersects itself below the rope, forming a topological structure known as hyperbolic flux tube, where a current sheet develops, triggering the reconnection. By mapping footpoints of the newly-reconnected field lines, we are able to reproduce both the spatial location and, for the first time, the temporal separation of the observed flare ribbons, as well as the dynamic boundary of the flux rope's feet. Futhermore, the temporal profile of the total reconnection flux is comparable to the soft X-ray light curve. Such a sophisticated characterization of the evolving magnetic topology provides important insight into the eventual understanding and forecast of solar eruptions.


title: Novel Strongly Spin-Orbit Coupled Quantum Dimer Magnet: Yb$_2$Si$_2$O$_7$

abstract: The quantum dimer magnet (QDM) is the canonical example of quantum magnetism. The QDM state consists of entangled nearest-neighbor spin dimers and often exhibits a field-induced triplon Bose-Einstein condensate (BEC) phase. We report on a new QDM in the strongly spin-orbit coupled, distorted honeycomb-lattice material Yb$_2$Si$_2$O$_7$. Our single crystal neutron scattering, specific heat, and ultrasound velocity measurements reveal a gapped singlet ground state at zero field with sharp, dispersive excitations. We find a field-induced magnetically ordered phase reminiscent of a BEC phase, with exceptionally low critical fields of $H_{c1} \sim 0.4$ T and $H_{c2} \sim 1.4$ T. Using inelastic neutron scattering in an applied magnetic field we observe a Goldstone mode (gapless to within $\delta E$ = 0.037 meV) that persists throughout the entire field-induced magnetically ordered phase, suggestive of the spontaneous breaking of U(1) symmetry expected for a triplon BEC. However, in contrast to other well-known cases of this phase, the high-field ($\mu$$_0$$H\geq1.2$T) part of the phase diagram in Yb$_2$Si$_2$O$_7$ is interrupted by an unusual regime signaled by a change in the field dependence of the ultrasound velocity and magnetization, as well as the disappearance of a sharp anomaly in the specific heat. These measurements raise the question of how anisotropy in strongly spin-orbit coupled materials modifies the field induced phases of QDMs.


title: Attentive Neural Network for Named Entity Recognition in Vietnamese

abstract: We propose an attentive neural network for the task of named entity recognition in Vietnamese. The proposed attentive neural model makes use of character-based language models and word embeddings to encode words as vector representations. A neural network architecture of encoder, attention, and decoder layers is then utilized to encode knowledge of input sentences and to label entity tags. The experimental results show that the proposed attentive neural network achieves the state-of-the-art results on the benchmark named entity recognition datasets in Vietnamese in comparison to both hand-crafted features based models and neural models.