title: Galaxy wide radio induced feedback in a radio quiet quasar

abstract: We report the discovery of a radio quiet type 2 quasar (SDSS J165315.06+234943.0 nicknamed the "Beetle" at z=0.103) with unambiguous evidence for active galactic nucleus (AGN) radio induced feedback acting across a total extension of ~46 kpc and up to ~26 kpc from the AGN. To the best of our knowledge, this is the first radio quiet system where radio induced feedback has been securely identified at >>several kpc from the AGN. Turbulent gas is also found far from the radio axis, ~25 kpc in the perpendicular direction. We propose a scenario in which the radio structures have perforated the interstellar medium of the galaxy and escaped into the circumgalactic medium. While advancing, they have interacted with in-situ gas modifying its properties. Our results show that jets of modest power can be the dominant feedback mechanism acting across huge volumes in radio quiet systems, including highly accreting luminous AGN, where radiative mode feedback may be expected.


title: Anomalous gauge interactions in photon collisions at the LHC and the FCC

abstract: The forward proton detectors recently installed and operating at the LHC open the possibility to observe photon collisions with high precision, providing a novel window on physics beyond the Standard Model. We review recent simulations and theoretical developments about the measurement of anomalous ${\gamma}{\gamma}{\gamma}{\gamma}$ and $Z{\gamma}{\gamma}{\gamma}$ interactions. The searches for these anomalous gauge interactions are expected to set bounds on a wide range of particles including generic electroweak particles, neutral particles with dimension-5 coupling to gauge bosons, polarizable dark particles, and are typically complementary from new physics searches in other channels.


title: Wang and Yau's Quasi-Local Energy for an Extreme Kerr Spacetime

abstract: There exist constant radial surfaces, $\mathcal{S}$, that may not be globally embeddable in $\mathbb{R}^3$ for Kerr spacetimes with $a>\sqrt{3}M/2$. To compute the Brown and York (B-Y) quasi-local energy (QLE), one must isometrically embed $\mathcal{S}$ into $\mathbb{R}^3$. On the other hand, the Wang and Yau (W-Y) QLE embeds $\mathcal{S}$ into Minkowski space. In this paper, we examine the W-Y QLE for surfaces that may or may not be globally embeddable in $\mathbb{R}^3$. We show that their energy functional, $E[\tau]$, has a critical point at $\tau=0$ for all constant radial surfaces in $t=constant$ hypersurfaces using Boyer-Lindquist coordinates. For $\tau=0$, the W-Y QLE reduces to the B-Y QLE. To examine the W-Y QLE in these cases, we write the functional explicitly in terms of $\tau$ under the assumption that $\tau$ is only a function of $\theta$. We then use a Fourier expansion of $\tau\left(\theta\right)$ to explore the values of $E[\tau\left(\theta\right)]$ in the space of coefficients. From our analysis, we discovered an open region of complex values for $E[\tau\left(\theta\right)]$. We also study the physical properties of the smallest real value of $E[\tau\left(\theta\right)]$, which lies on the boundary separating real and complex energies.


title: Elastic scattering of a quark from a color field: longitudinal momentum exchange

abstract: Perturbative QCD in the small Bjorken $x$ limit can be formulated as an effective theory known as the Color Glass Condensate (CGC) formalism. The CGC formalism takes into account the dynamics of large gluon densities at small $x$ and has been successfully applied to Deep Inelastic Scattering (DIS) and particle production in high energy hadronic and nuclear collisions in the small $x$ kinematic region. The effective degrees of of freedom in CGC are Wilson lines which enter in the effective quark (and gluon) propagators and re-sum multiple soft scatterings from the small $x$ gluon field of the target. It is however known that the CGC effective theory breaks down when one probes the moderately large $x$ (high $p_t$) kinematics where collinear factorization and DGLAP evolution of parton distribution functions should be the right framework. Here we propose a general framework which may allow one to eventually unify the two approaches and to calculate pQCD cross sections in both small and large Bjorken $x$ regions. We take the first step towards this goal by deriving an expression for the quark propagator in a background field which includes scatterings from both small and large $x$ modes of the gluon field of the target. We describe how this quark propagator can be used to calculate QCD structure functions $F_2$ and $F_L$ at all $x$ and thus generalize the dipole model of DIS. We outline this approach can also be used to extend the so-called hybrid approach to particle production in the forward rapidity region of high energy hadronic and nuclear collisions to all $x$ and $p_t$ regions and speculate on how one may apply the same techniques to extend the McLerran-Venugopalan effective action used in high energy heavy ion collisions to include high $p_t$ physics.


title: Detecting and monitoring foodborne illness outbreaks: Twitter communications and the 2015 U.S. Salmonella outbreak linked to imported cucumbers

abstract: This research uses Twitter, as a social media device, to track communications related to the 2015 U.S. foodborne illness outbreak linked to Salmonella in imported cucumbers from Mexico. The relevant Twitter data are analyzed in light of the timeline of the official announcements made by the Centers for Disease Control and Prevention (CDC). The largest number of registered tweets is associated with the period immediately following the CDC initial announcement and the official release of the first recall of cucumbers.


title: The Intrinsic Far-UV Spectrum of B1422+231

abstract: We present new spectroscopy of the z=3.62 gravitationally lensed quasar B1422+117 from the Gemini North GMOS integral field spectrograph. We observe significant differential magnifications between the broad emission lines and the continuum, as well as across the velocity structure of the Lyman-alpha line. We take advantage of this differential microlensing to algebraically decompose the quasar spectrum into the absorbed broad emission line and absorbed continuum components. We use the latter to derive the intrinsic Ly-alpha forest absorption spectrum. The proximity effect is clearly detected, with a proximity zone edge of 8.6-17.3 Mpc from the quasar, implying (perhaps intermittent) activity over at least 28 Myrs. The Ly-alpha line profile exhibits a blue excess that is inconsistent with a symmetric fit to the unabsorbed red side. This has important implications for the use of this fitting technique in estimating the absorbed blue Ly-alpha wings of Gunn-Peterson trough quasars.


title: Global regularity for a family of 3D models of the axisymmetric Navier-Stokes equations

abstract: We consider a family of 3D models for the axi-symmetric incompressible Navier-Stokes equations. The models are derived by changing the strength of the convection terms in the axisymmetric Navier-Stokes equations written using a set of transformed variables. We prove the global regularity of the family of models in the case that the strength of convection is slightly stronger than that of the original Navier-Stokes equations, which demonstrates the potential stabilizing effect of convection.


title: Memory matrix theory of the dc resistivity of a disordered antiferromagnetic metal with an effective composite operator

abstract: We perform the calculation of the dc resistivity as a function of temperature of the "strange-metal" state that emerges in the vicinity of a spin-density-wave phase transition in the presence of weak disorder. This scenario is relevant to the phenomenology of many important correlated materials, such as, e.g., the pnictides, the heavy-fermion compounds and the cuprates. To accomplish this task, we implement the memory-matrix approach that allows the calculation of the transport coefficients of the model beyond the quasiparticle paradigm. Our computation is also inspired by the $\epsilon=3-d$ expansion in a hot-spot model embedded in $d$-space dimensions recently put forth by Sur and Lee [Phys. Rev. B 91, 125136 (2015)], in which they find a new low-energy non-Fermi liquid fixed point that is perturbatively accessible near three dimensions. As a consequence, we are able to establish here the temperature and doping dependence of the electrical resistivity at intermediate temperatures of a two-dimensional disordered antiferromagnetic metallic model with a composite operator that couples the order-parameter fluctuations to the entire Fermi surface. We argue that our present theory provides a good basis in order to unify the experimental transport data, e.g., in the cuprates and the pnictide superconductors, within a wide range of doping regimes.


title: Evaluation of latent variances in Monte Carlo dose calculations with Varian TrueBeam photon phase-spaces used as a particle source

abstract: Purpose: To evaluate the latent variance (LV) of TrueBeam photon phase-space files (PSF) for open 10x10$cm^2$ and small fields. Methods: BEAMnrc/DOSXYZnrc was used to transport particles from Varian phase-space files (PSFA) through the secondary collimators. Transported particles were scored into another phase-space located under the jaws (PSFB), or transported further through the cone collimators and scored straight below, forming PSFC. PFSB were scored for 6MV-FFF, 6MV, 10MV-FFF, 10MV and 15MV beams with 10x10$cm^2$ field size, and PSFC were scored for 6MV beam under cones of 0.13, 0.25, 0.35, 1.0, 1.2, 1.5 and 4cm diameter. PSFB and PSFC were transported into a water phantom with particle recycling number ranging from 10 to 1000. For 10x10$cm^2$ fields 0.5x0.5x0.5cm^3 voxels were used to score the dose, whereas the dose was scored in 0.1x0.1x0.5cm^3 voxels for beams collimated with small cones. For small 0.25cm diameter cone-collimated 6MV beam, phantom voxel size varied as 0.02x0.02x0.5cm^3, 0.05x0.05x0.5cm^3 and 0.1x0.1x0.5cm^3. Dose variances were scored in all cases and LV evaluated as per Sempau et al. Results: For the 10x10cm^2 fields calculated LVs were greatest at the phantom surface and decreased with depth until reached a plateau at 5cm depth. LVs were 0.54%, 0.96%, 0.35%, 0.69% and 0.57% for the 6MV-FFF, 6MV, 10MV-FFF, 10MV and 15MV energies, respectively at 10cm depth. For the 6 MV phase-space collimated with cones of 0.13, 0.25, 0.35, 1.0cm diameter, the LVs calculated at 1.5cm depth were 75.6%, 25.4%, 17.6% and 8.0% respectively. Calculated LV for the 0.25cm cone-collimated 6MV beam were 61.2%, 40.7%, 22.5% in 0.02x0.02x0.5cm^3, 0.05x0.05x0.5cm^3 and 0.1x0.1x0.5cm^3 voxels respectively. Conclusion: Single PSF can be used to achieve sub-percent latent variance in open 10x10cm^2 field MC simulations, whereas more PSFs would have to be summed for small SRS fields.


title: Decoherence of two coupled singlet-triplet spin qubits

abstract: We study a pair of capacitively coupled singlet-triplet spin qubits. We characterize the two-qubit decoherence through two complementary measures, the decay time of coupled-qubit oscillations and the fidelity of entangled state preparation. We provide a quantitative map of their dependence on charge noise and field noise, and we highlight the magnetic field gradient across each singlet-triplet qubit as an effective tool to suppress decoherence due to charge noise.


title: $t$-Structures on stable derivators and Grothendieck hearts

abstract: We prove that given any strong, stable derivator and a $t$-structure on its base triangulated category $\cal D$, the $t$-structure canonically lifts to all the (coherent) diagram categories and each incoherent diagram in the heart uniquely lifts to a coherent one. We use this to show that the $t$-structure being compactly generated implies that the coaisle is closed under directed homotopy colimit which in turns implies that the heart is an (Ab.$5$) Abelian category. If, moreover, $\cal D$ is a well generated algebraic or topological triangulated category, then the heart of any accessibly embedded (in particular, compactly generated) $t$-structure has a generator. As a consequence, it follows that the heart of any compactly generated $t$-structure of a well generated algebraic or topological triangulated category is a Grothendieck category.


title: On the geometry of some Equivariantly related manifolds

abstract: We provide a topological procedure to obtain geometric realizations of both classical and `exotic' $G$-manifolds, such as spheres, bundles over spheres and Kervaire manifolds. As an application, we apply the process known as Cheeger deformations to produce new metrics of both positive Ricci and almost non-negative curvature on such objects.


title: Smoothness and long time existence for solutions of the porous medium equation on manifolds with conical singularities

abstract: We study the porous medium equation on manifolds with conical singularities. Given strictly positive initial values, we show that the solution exists in the maximal $L^{q}$-regularity space for all times and is instantaneously smooth in space and time, where the maximal $L^{q}$-regularity is obtained in the sense of Mellin-Sobolev spaces. Moreover, we obtain precise information concerning the asymptotic behavior of the solution close to the singularity. Finally, we show the existence of generalized solutions for non-negative initial data.


title: A test of the high-eccentricity migration scenario for close-in planets

abstract: In the high-eccentricity migration (HEM) scenario, close-in planets reach the vicinity of the central star on high-eccentricity orbits that become circularized---with a concomitant decrease in the semimajor axis---through a tidal interaction with the star. Giant planets that arrive with periastron distances that are smaller than the Roche limit $a_\mathrm{R}$ lose their gaseous envelopes, resulting in an inner edge to the surviving planets' distribution. The observational evidence for this effect, while extensive, is nonetheless somewhat ambiguous because of the effect of tidal orbital decay. Here we consider another key prediction of the HEM scenario---the existence of a spatial eccentricity gradient near the location where the circularization time becomes comparable to the planet's age for typical parameters. Previous studies already found evidence for this gradient and demonstrated that its properties are consistent with the circularization process being dominated by tidal dissipation in the planet (encapsulated by the tidal quality factor $Q^\prime_\mathrm{p}$). Our work extends these treatments by constructing explicit model distributions for comparison with the data and by carrying out backward-in-time integrations using observed system parameters. We show that circularization generally occurs outside the distribution's inner edge (which defines the boundary of the so-called sub-Jovian desert) and that typically $Q^\prime_\mathrm{p}\approx10^6$ in the circularization zone (to within a factor of 3). We also find tentative evidence for an eccentricity gradient in lower-mass planets, indicating that formation through HEM may be relevant down to Neptune scales.


title: Cheap and near exact CASSCF with large active spaces

abstract: We use the recently-developed Heat-bath Configuration Interaction (HCI) algorithm as an efficient active-space solver to perform multi-configuration self-consistent field calculations (HCISCF) with large active spaces. We give a detailed derivation of the theory and show that difficulties associated with non-variationality of the HCI procedure can be overcome by making use of the Lagrangian formulation to calculate the HCI relaxed two body reduced density matrix. HCISCF is then used to study the electronic structure of butadiene, pentacene, and Fe-porphyrin. One of the most striking results of our work is that the converged active space orbitals obtained from HCISCF are relatively insensitive to the accuracy of the HCI calculation. This allows us to obtain nearly converged CASSCF energies with an estimated error of less than 1 mHa using the orbitals obtained from the HCISCF procedure in which the integral transformation is the dominant cost. For example, an HCISCF calculation on Fe-Porphyrin model complex with an active space of (44e, 44o) took only 412 seconds per iteration on a single node containing 28 cores, out of which 185 seconds were spent in the HCI calculation and the remaining 227 seconds were mainly used for integral transformation. Finally, we also show that active-space orbitals can be optimized using HCISCF to substantially speed up the convergence of the HCI energy to the Full CI limit because HCI is not invariant to unitary transformations within the active space.


title: Asymptotic Stability of the Landau-Lifshitz Equation

abstract: The Landau--Lifshitz equation describes the behaviour of magnetic domains in ferromagnetic structures. Recently such structures have been found to be favourable for storing digital data. Stability of magnetic domains is important for this. Consequently, asymptotic stability of the equilibrium points in the Landau--Lifshitz equation are established. A suitable Lyapunov function is presented.


title: Center and isochronous center of a class of quasi-analytic switching systems

abstract: In this paper, we study the integrability and linearization of a class of quadratic quasi-analytic switching systems. We improve an existing method to compute the focus values and periodic constants of quasi-analytic switching systems. In particular, with our method, we demonstrate that the dynamical behaviors of quasi-analytic switching systems are more complex than that of continuous quasi-analytic systems, by showing the existence of six and seven limit cycles in the neighborhood of the origin and infinity, respectively, in a quadratic quasi-analytic switching system. Moreover, explicit conditions are obtained for classifying the centers and isochronous centers of the system.


title: Cosmicflows-3: Cold Spot Repeller?

abstract: The three-dimensional gravitational velocity field within z~0.1 has been modeled with the Wiener filter methodology applied to the Cosmicflows-3 compilation of galaxy distances. The dominant features are a basin of attraction and two basins of repulsion. The major basin of attraction is an extension of the Shapley concentration of galaxies. One basin of repulsion, the Dipole Repeller, is located near the anti-apex of the cosmic microwave background dipole. The other basin of repulsion is in the proximate direction toward the 'Cold Spot' irregularity in the cosmic microwave background. It has been speculated that a vast void might contribute to the amplitude of the Cold Spot from the integrated Sachs-Wolfe effect.


title: Audio signal encryption using chaotic H\'enon map and lifting wavelet transforms

abstract: We propose a new audio signal encryption scheme based on the chaotic H\'enon map. The scheme mainly comprises two phases: one is the preprocessing stage where the audio signal is transformed into a data by the lifting wavelet scheme and the other in which the transformed data is encrypted by chaotic data set and hyperbolic functions. Furthermore, we use dynamic keys and consider the key space size to be large enough to resist any kind of cryptographic attacks. A statistical investigation is also made to test the security and the efficiency of the proposed scheme.


title: Objective Classes for Micro-Facial Expression Recognition

abstract: Micro-expressions are brief spontaneous facial expressions that appear on a face when a person conceals an emotion, making them different to normal facial expressions in subtlety and duration. Currently, emotion classes within the CASME II dataset are based on Action Units and self-reports, creating conflicts during machine learning training. We will show that classifying expressions using Action Units, instead of predicted emotion, removes the potential bias of human reporting. The proposed classes are tested using LBP-TOP, HOOF and HOG 3D feature descriptors. The experiments are evaluated on two benchmark FACS coded datasets: CASME II and SAMM. The best result achieves 86.35\% accuracy when classifying the proposed 5 classes on CASME II using HOG 3D, outperforming the result of the state-of-the-art 5-class emotional-based classification in CASME II. Results indicate that classification based on Action Units provides an objective method to improve micro-expression recognition.


title: Critical behavior of a chiral superfluid in a bipartite square lattice

abstract: We study the critical behavior of Bose-Einstein condensation in the second band of a bipartite optical square lattice in a renormalization group framework at one-loop order. Within our field theoretical representation of the system, we approximate the system as a two-component Bose gas in three dimensions. We demonstrate that the system is in a different universality class than the previously studied condensation in a frustrated triangular lattice due to an additional Umklapp scattering term, which stabilizes the chiral superfluid order at low temperatures. We derive the renormalization group flow of the system and show that this order persists in the low energy limit. Furthermore, the renormalization flow suggests that the phase transition from the thermal phase to the chiral superfluid state is first order.


title: Spark complexes on good effective orbifold atlases categorically

abstract: Good atlases are defined for effective orbifolds, and a spark complex is constructed on each good atlas. It is proved that this process is 2-functorial with compatible systems playing as morphisms between good atlases, and that the spark character 2-functor factors through this 2-functor.


title: Perturbed redshifts from N-body simulations

abstract: In order to keep pace with the increasing data quality of astronomical surveys the observed source redshift has to be modeled beyond the well-known Doppler contribution. In this letter I want to examine the gauge issue that is often glossed over when one assigns a perturbed redshift to simulated data generated with a Newtonian N-body code. A careful analysis reveals the presence of a correction term that has so far been neglected. It is roughly proportional to the observed length scale divided by the Hubble scale and therefore suppressed inside the horizon. However, on gigaparsec scales it can be comparable to the gravitational redshift and hence amounts to an important relativistic effect.


title: Electron Excitation of High Dipole Moment Molecules Reexamined

abstract: Emission from high-dipole moment molecules such as HCN allows determination of the density in molecular clouds, and is often considered to trace the "dense" gas available for star formation. We assess the importance of electron excitation in various environments. The ratio of the rate coefficients for electrons and H$_2$ molecules, $\simeq$10$^5$ for HCN, yields the requirements for electron excitation to be of practical importance if $n({\rm H}_2) \leq\ 10^{5.5} ~ \rm cm^{-3}$ and $X({\rm e}^-) \geq\ 10^{-5}$, where the numerical factors reflect critical values $n_{\rm{}c}({\rm H_2})$ and $X^*({\rm{}e}^-)$. This indicates that in regions where a large fraction of carbon is ionized, $X({\rm e}^-)$ will be large enough to make electron excitation significant. The situation is in general similar for other "high density tracers", including HCO$^+$, CN, and CS. But there are significant differences in the critical electron fractional abundance, $X^*({\rm e}^-)$, defined by the value required for equal effect from collisions with H$_2$ and e$^-$. Electron excitation is, for example, unimportant for CO and C$^+$. Electron excitation may be responsible for the surprisingly large spatial extent of the emission from dense gas tracers in some molecular clouds (Pety et al. 2017; Kauffmann, Goldsmith et al. 2017). The enhanced estimates for HCN abundances and HCN/CO and HCN/HCO$^+$ ratios observed in the nuclear regions of luminous galaxies may be in part a result of electron excitation of high dipole moment tracers. The importance of electron excitation will depend on detailed models of the chemistry, which may well be non-steady state and non-static.


title: Monotonic properties of the shift and penetration factors

abstract: We study derivatives of the shift and penetration factors of collision theory with respect to energy, angular momentum, and charge. Definitive results for the signs of these derivatives are found for the repulsive Coulomb case. In particular, we find that the derivative of the shift factor with respect to energy is positive for the repulsive Coulomb case, a long anticipated but heretofore unproven result. These results are closely connected to the properties of the sum of squares of the regular and irregular Coulomb functions; we also present investigations of this quantity.

categories: quant-ph math-ph math.CA math.MP nucl-th


title: A Robust Indoor Scene Recognition Method based on Sparse Representation

abstract: In this paper, we present a robust method for scene recognition, which leverages Convolutional Neural Networks (CNNs) features and Sparse Coding setting by creating a new representation of indoor scenes. Although CNNs highly benefited the fields of computer vision and pattern recognition, convolutional layers adjust weights on a global-approach, which might lead to losing important local details such as objects and small structures. Our proposed scene representation relies on both: global features that mostly refers to environment's structure, and local features that are sparsely combined to capture characteristics of common objects of a given scene. This new representation is based on fragments of the scene and leverages features extracted by CNNs. The experimental evaluation shows that the resulting representation outperforms previous scene recognition methods on Scene15 and MIT67 datasets, and performs competitively on SUN397, while being highly robust to perturbations in the input image such as noise and occlusion.


title: Extinction and the Dimming of KIC 8462852

abstract: To test alternative hypotheses for the behavior of KIC 8462852, we obtained measurements of the star over a wide wavelength range from the UV to the mid-infrared from October 2015 through December 2016, using Swift, Spitzer and at AstroLAB IRIS. The star faded in a manner similar to the long-term fading seen in Kepler data about 1400 days previously. The dimming rate for the entire period reported is 22.1 +\- 9.7 milli-mag/yr in the Swift wavebands, with amounts of 21.0 +\- 4.5 mmag in the groundbased B measurements, 14.0 +\- 4.5 mmag in V, and 13.0 +\- 4.5 in R, and a rate of 5.0 +\- 1.2 mmag/yr averaged over the two warm Spitzer bands. Although the dimming is small, it is seen at >= 3 sigma by three different observatories operating from the UV to the IR. The presence of long-term secular dimming means that previous SED models of the star based on photometric measurements taken years apart may not be accurate. We find that stellar models with T_{eff} = 7000 - 7100 K and A_V ~ 0.73 best fit the Swift data from UV to optical. These models also show no excess in the near-simultaneous Spitzer photometry at 3.6 and 4.5 microns, although a longer wavelength excess from a substantial debris disk is still possible (e.g., as around Fomalhaut). The wavelength dependence of the fading favors a relatively neutral color (i.e., R_V >= 5, but not flat across all the bands) compared with the extinction law for the general ISM (R_V = 3.1), suggesting that the dimming arises from circumstellar material.


title: Nondestructive testing of grating imperfections using grating-based X-ray phase-contrast imaging

abstract: We reported the usage of grating-based X-ray phase-contrast imaging in nondestructive testing of grating imperfections. It was found that electroplating flaws could be easily detected by conventional absorption signal, and in particular, we observed that the grating defects resulting from uneven ultraviolet exposure could be clearly discriminated with phase-contrast signal. The experimental results demonstrate that grating-based X-ray phase-contrast imaging, with a conventional low-brilliance X-ray source, a large field of view and a reasonable compact setup, which simultaneously yields phase- and attenuation-contrast signal of the sample, can be ready-to-use in fast nondestructive testing of various imperfections in gratings and other similar photoetching products.


title: First principle study on the transition of electronic, mechanic and piezoelectric property of hexagonal Boron Nitride nanotube

abstract: One-dimensional nanostructures such as nanowires and nanotubes are stimulating tremendous interest due to their structural, electronic and magnetic properties. In this study, the first principle calculation was performed to investigate on the size effect on the electronic, mechanic and piezoelectric properties of hexagonal-BN nanotube (BNNT). It demonstrates that the nanotube diameter can modulate the direct bandgap significantly of the zigzag-BNNT, increasing from 3.96eV to 4.44eV with the larger diameter from 10 to 16 unit cells in circular direction. Both armchair and zigzag BNNT exhibit large elastic modulus, implying promising applications in nanoscale surface engineering, tribology and nanomanufacturing/nanofabrication. Outstanding piezoelectricity is also observed with large piezoelectric coefficient of 0.084C/m2 for 10-zigazag-BNNT and 0.13 C/m2 for 16-zigzag-BNNT, exhibiting excellent potential in producing mechao-electric power generators.


title: Efficient barycentric point sampling on meshes

abstract: We present an easy-to-implement and efficient analytical inversion algorithm for the unbiased random sampling of a set of points on a triangle mesh whose surface density is specified by barycentric interpolation of non-negative per-vertex weights. The correctness of the inversion algorithm is verified via statistical tests, and we show that it is faster on average than rejection sampling.


title: The Study of the Bloch Transformed Fields Scattered by Locally Perturbed Periodic Surfaces

abstract: Scattering problems with locally perturbed periodic surfaces have been studied both theoretically and numerically in recent years. In this paper, we will discuss the regularity results of the Bloch transform of the total fields. The idea is inspired by Theorem a in \cite{Kirsc1993}, which considered how the total field depends on the wave numbers and the incident angles, with a family of plain incident fields and a smooth enough periodic surface. We will show that when the incident field satisfies some certain conditions, the Bloch transform of the total field depends analytically on the quasi-periodicities one the straight line $\R$ except for a countable number of points, while near such points, a square-root like singularity exists. We also give some examples to show that the conditions are satisfied by a large number of commonly used incident fields. This result also provides a probability to improve the numerical solution of this kind of problems, which is expected to be discussed in our future papers.


title: Model-theoretic Characterizations of Large Cardinals

abstract: We consider compactness characterizations of large cardinals. Based on results of Benda \cite{b-sccomp}, we study compactness for omitting types in various logics. In $\bL_{\kappa, \kappa}$, this allows us to characterize any large cardinal defined in terms of normal ultrafilters, and we also analyze second-order and sort logic. In particular, we give a compactness for omitting types characterization of huge cardinals, which have consistency strength beyond Vop\v{e}nka's Principle.


title: The Status and Initial Results of the MAJORANA DEMONSTRATOR Experiment

abstract: Neutrinoless double-beta decay searches play a major role in determining the nature of neutrinos, the existence of a lepton violating process, and the effective Majorana neutrino mass. The MAJORANA Collaboration assembled an array of high purity Ge detectors to search for neutrinoless double-beta decay in Ge-76. The MAJORANA DEMONSTRATOR is comprised of 44.1 kg (29.7 kg enriched in Ge-76) of Ge detectors divided between two modules contained in a low-background shield at the Sanford Underground Research Facility in Lead, South Dakota, USA. The initial goals of the DEMONSTRATOR are to establish the required background and scalability of a Ge-based next-generation ton-scale experiment. Following a commissioning run that started in 2015, the first detector module started low-background data production in early 2016. The second detector module was added in August 2016 to begin operation of the entire array. We discuss results of the initial physics runs, as well as the status and physics reach of the full MAJORANA DEMONSTRATOR experiment.


title: Excitation of coupled phononic frequency combs via two-mode parametric three-wave mixing

abstract: This paper builds on the recent demonstration of three-wave mixing based phononic frequency comb. Here, in this process, an intrinsic coupling between the drive and resonant frequency leads to a frequency comb of spacing corresponding to the separation between drive and resonant frequency. Now, in this paper, we experimentally demonstrate the possibility to further excite multiple frequency combs with the same external drive through its coupling with other identical devices. In addition, we also experimentally identify interesting features associated with such a frequency comb generation process.


title: BMS in higher space-time dimensions and Non-relativistic BMS

abstract: Master's thesis. We present a study of the BMS Group is higher space-time dimensions, and the extension of this group to non-relativistic systems.


title: Variations in the 6.2 $\mu$m emission profile in starburst-dominated galaxies: a signature of polycyclic aromatic nitrogen heterocycles (PANHs)?

abstract: Analyses of the polycyclic aromatic hydrocarbon (PAH) feature profiles, especially the 6.2 $\mu$m feature, could indicate the presence of nitrogen incorporated in their aromatic rings. In this work, 155 predominantly starburst-dominated galaxies (including HII regions and Seyferts, for example), extracted from the Spitzer/IRS ATLAS project (Hern\'an-Caballero & Hatziminaoglou 2011), have their 6.2 $\mu$m profiles fitted allowing their separation into the Peeters' A, B and C classes (Peeters et al. 2002). 67% of these galaxies were classified as class A, 31% were as class B and 2% as class C. Currently class A sources, corresponding to a central wavelength near 6.22 $\mu$m, seem only to be explained by polycyclic aromatic nitrogen heterocycles (PANH, Hudgins et al. 2005), whereas class B may represent a mix between PAHs and PANHs emissions or different PANH structures or ionization states. Therefore, these spectra suggest a significant presence of PANHs in the interstellar medium (ISM) of these galaxies that could be related to their starburst-dominated emission. These results also suggest that PANHs constitute another reservoir of nitrogen in the Universe, in addition to the nitrogen in the gas phase and ices of the ISM.


title: Weighing Uranus' moon Cressida with the $\eta$ ring

abstract: The $\eta$ ring is one of the narrow rings of Uranus, consisting of a dense core that is 1-2 km wide and a diffuse outer sheet spanning about 40 km. Its dense core lies just exterior to the 3:2 Inner Lindblad Resonance of the small moon Cressida. We fit the $\eta$ ring radius residuals and longitudes from a complete set of both ground-based and Voyager stellar and radio occultations of the Uranian rings spanning 1977-2002. We find variations in the radial position of the $\eta$ ring that are likely generated by this resonance, and take the form of a 3-lobed structure rotating at an angular rate equal to the mean motion of the moon Cressida. The amplitude of these radial oscillations is $0.667\pm0.113$ km, which is consistent with the expected shape due to the perturbations from Cressida. The magnitude of these variations provides the first measurement of the mass and density of the moon Cressida ($m=2.5\pm0.4\times10^{17}$ kg and $\rho=0.86\pm0.16$ g/cm$^3$) or, indeed, any of Uranus' small inner moons. A better grasp of inner Uranian satellite masses will provide another clue to the composition, dynamical stability, and history of Uranus' tightly packed system of small moons.


title: Active Preference Learning for Personalized Portfolio Construction

abstract: In financial asset management, choosing a portfolio requires balancing returns, risk, exposure, liquidity, volatility and other factors. These concerns are difficult to compare explicitly, with many asset managers using an intuitive or implicit sense of their interaction. We propose a mechanism for learning someone's sense of distinctness between portfolios with the goal of being able to identify portfolios which are predicted to perform well but are distinct from the perspective of the user. This identification occurs, e.g., in the context of Bayesian optimization of a backtested performance metric. Numerical experiments are presented which show the impact of personal beliefs in informing the development of a diverse and high-performing portfolio.


title: Three Photon Entanglement from Ortho-Positronium Revisited

abstract: Entanglement of three photons from the decay of ortho-positronium is re-analyzed. We use the full three body phase space to write down the entangled states classified according to the spin directions of the ortho-positronium. Even in the case when the spin is perpendicular to the decay plane, we find non-negligible phases entering the entangled state. This has not been noticed before. We advocate a fixed quantization axis of the spin for the sake of generality. A brief discussion of a three dimensional formalism for photons, including correlations, versus a two dimensional one is given.


title: Secure by default - the case of TLS

abstract: Default configuration of various software applications often neglects security objectives. We tested the default configuration of TLS in dozen web and application servers. The results show that "secure by default" principle should be adopted more broadly by developers and package maintainers. In addition, system administrators cannot rely blindly on default security options.


title: Leaf Counting with Deep Convolutional and Deconvolutional Networks

abstract: In this paper, we investigate the problem of counting rosette leaves from an RGB image, an important task in plant phenotyping. We propose a data-driven approach for this task generalized over different plant species and imaging setups. To accomplish this task, we use state-of-the-art deep learning architectures: a deconvolutional network for initial segmentation and a convolutional network for leaf counting. Evaluation is performed on the leaf counting challenge dataset at CVPPP-2017. Despite the small number of training samples in this dataset, as compared to typical deep learning image sets, we obtain satisfactory performance on segmenting leaves from the background as a whole and counting the number of leaves using simple data augmentation strategies. Comparative analysis is provided against methods evaluated on the previous competition datasets. Our framework achieves mean and standard deviation of absolute count difference of 1.62 and 2.30 averaged over all five test datasets.


title: Energy density and path-length dependence of the fractional momentum loss in heavy-ion collisions at $\sqrt{s_{\rm NN}}$ from 62.4 to 5020 GeV

abstract: In this work a study of the fractional momentum loss ($S_{\rm loss}$) as a function of the characteristic path-length ($L$) and the Bjorken energy density times the equilibration time ($\epsilon_{\rm Bj}\tau_{0}$) for heavy-ion collisions at different $\sqrt{s_{\rm NN}}$ is presented. The study has been conducted using inclusive charged particles from intermediate to large transverse momentum ($5<p_{\rm T}<20$ GeV/$c$). Within uncertainties and for all the transverse momentum values which were explored, the fractional momentum loss linearly increases with $({\epsilon_{\rm Bj}\tau_{0}})^{3/8}$$L$. The functional form of $S_{\rm loss}$ vs. $({\epsilon_{\rm Bj}\tau_{0}})^{3/8}$$L$ seems to be universal. Moreover, for identified charged hadrons a linear relationship between $S_{\rm loss}$ and $L$ is also observed. The behaviour of data could provide important information aimed to understand the parton energy loss mechanism in heavy-ion collisions and some insight into the expected effect for small systems.


title: Thermal broadening of the power spectra of laser-trapped particles in vacuum

abstract: We show that at low pressures the spectral widths of the power spectra of laser-trapped particles are nearly independent from pressures and, due to the nonlinearities of the trap, reflect the thermal distribution of particles. In the experiments with nanoparticles trapped in an optical lattice, we identify two distinct features of the widths. First, the widths along an optical lattice are much broader than those in the other directions. Second, the spectral widths are narrower for larger nanoparticles. We develop a theory of thermal broadening and show that the spectral widths normalized by the frequencies of the center-of-mass motion directly reveal the ratio of the thermal energy to the trap depth. The presented model provides a good understanding of the observed features. Our model holds also for smaller particles such as atoms and molecules and can be readily extended to the general case with a single-beam optical trap.


title: Reconstruction of a compact Riemannian manifold from the scattering data of internal sources

abstract: Given a smooth non-trapping compact manifold with strictly con- vex boundary, we consider an inverse problem of reconstructing the manifold from the scattering data initiated from internal sources. This data consist of the exit directions of geodesics that are emaneted from interior points of the manifold. We show that under certain generic assumption of the metric, one can reconstruct an isometric copy of the manifold from such scattering data measured on the boundary.


title: Magnetization in pristine graphene with Zeeman splitting and variable spin-orbit coupling

abstract: The aim of this work is to describe the spin magnetization of graphene with Rashba spinorbit coupling and Zeeman effect. It is shown that the magnetization depends critically on the spin-orbit coupling l that is controlled with an external electric field. In turn, by manipulating the density of charge carriers, it is shown that spin up and down Landau levels mix introducing jumps in the spin magnetization. Two magnetic oscillations phases are described that can be tunable through the applied external fields. The maximum and minimum of the oscillations can be alternated by taking into account how the energy levels are filled when the Rashba-spin-orbit coupling is turned on. The results obtained are of importance to design superlattices with variable spin-orbit coupling with different configurations in which spin oscillations and spin filters can be developed.


title: Generalized Paxos Made Byzantine (and Less Complex)

abstract: One of the most recent members of the Paxos family of protocols is Generalized Paxos. This variant of Paxos has the characteristic that it departs from the original specification of consensus, allowing for a weaker safety condition where different processes can have a different views on a sequence being agreed upon. However, much like the original Paxos counterpart, Generalized Paxos does not have a simple implementation. Furthermore, with the recent practical adoption of Byzantine fault tolerant protocols, it is timely and important to understand how Generalized Paxos can be implemented in the Byzantine model. In this paper, we make two main contributions. First, we provide a description of Generalized Paxos that is easier to understand, based on a simpler specification and the pseudocode for a solution that can be readily implemented. Second, we extend the protocol to the Byzantine fault model.


title: Learning the Enigma with Recurrent Neural Networks

abstract: Recurrent neural networks (RNNs) represent the state of the art in translation, image captioning, and speech recognition. They are also capable of learning algorithmic tasks such as long addition, copying, and sorting from a set of training examples. We demonstrate that RNNs can learn decryption algorithms -- the mappings from plaintext to ciphertext -- for three polyalphabetic ciphers (Vigen\`ere, Autokey, and Enigma). Most notably, we demonstrate that an RNN with a 3000-unit Long Short-Term Memory (LSTM) cell can learn the decryption function of the Enigma machine. We argue that our model learns efficient internal representations of these ciphers 1) by exploring activations of individual memory neurons and 2) by comparing memory usage across the three ciphers. To be clear, our work is not aimed at 'cracking' the Enigma cipher. However, we do show that our model can perform elementary cryptanalysis by running known-plaintext attacks on the Vigen\`ere and Autokey ciphers. Our results indicate that RNNs can learn algorithmic representations of black box polyalphabetic ciphers and that these representations are useful for cryptanalysis.


title: The particle in a box in PT quantum mechanics and an electromagnetic analog

abstract: In PT quantum mechanics a fundamental principle of quantum mechanics, that the Hamiltonian must be hermitian, is replaced by another set of requirements, including notably symmetry under PT, where P denotes parity and T denotes time reversal. Here we study the role of boundary conditions in PT quantum mechanics by constructing a simple model that is the PT symmetric analog of a particle in a box. The model has the usual particle in a box Hamiltonian but boundary conditions that respect PT symmetry rather than hermiticity. We find that for a broad class of PT-symmetric boundary conditions the model respects the condition of unbroken PT-symmetry, namely that the Hamiltonian and the symmetry operator PT have simultaneous eigenfunctions, implying that the energy eigenvalues are real. We also find that the Hamiltonian is self-adjoint under the PT inner product. Thus we obtain a simple soluble model that fulfils all the requirements of PT quantum mechanics. In the second part of this paper we formulate a variational principle for PT quantum mechanics that is the analog of the textbook Rayleigh-Ritz principle. Finally we consider electromagnetic analogs of the PT-symmetric particle in a box. We show that the isolated particle in a box may be realized as a Fabry-Perot cavity between an absorbing medium and its conjugate gain medium. Coupling the cavity to an external continuum of incoming and outgoing states turns the energy levels of the box into sharp resonances. Remarkably we find that the resonances have a Breit-Wigner lineshape in transmission and a Fano lineshape in reflection; by contrast in the corresponding hermitian case the lineshapes always have a Breit-Wigner form in both transmission and reflection.


title: Nonlinear responses from the interaction of two progressing waves at an interface

abstract: For scalar semilinear wave equations, we analyze the interaction of two (distorted) plane waves at an interface between media of different nonlinear properties. We show that new waves are generated from the nonlinear interactions, which might be responsible for the observed nonlinear effects in applications. Also, we show that the incident waves and the nonlinear responses determine the location of the interface and some information of the nonlinear properties of the media. In particular, for the case of a jump discontinuity at the interface, we can determine the magnitude of the jump.


title: Hamiltonian Maker-Breaker games on small graphs

abstract: We look at the unbiased Maker-Breaker Hamiltonicity game played on the edge set of a complete graph $K_n$, where Maker's goal is to claim a Hamiltonian cycle. First, we prove that, independent of who starts, Maker can win the game for $n = 8$ and $n = 9$. Then we use an inductive argument to show that, independent of who starts, Maker can win the game if and only if $n \geq 8$. This, in particular, resolves in the affirmative the long-standing conjecture of Papaioannou. We also study two standard positional games related to Hamiltonicity game. For Hamiltonian Path game, we show that Maker can claim a Hamiltonian path if and only if $n \geq 5$, independent of who starts. Next, we look at Fixed Hamiltonian Path game, where the goal of Maker is to claim a Hamiltonian path between two predetermined vertices. We prove that if Maker starts the game, he wins if and only if $n \geq 7$, and if Breaker starts, Maker wins if and only if $n \geq 8$. Using this result, we are able to improve the previously best upper bound on the smallest number of edges a graph on $n$ vertices can have, knowing that Maker can win the Maker-Breaker Hamiltonicity game played on its edges. To resolve the outcomes of the mentioned games on small (finite) boards, we devise algorithms for efficiently searching game trees and then obtain our results with the help of a computer.


title: The Expanding Approvals Rule: Improving Proportional Representation and Monotonicity

abstract: Proportional representation (PR) is often discussed in voting settings as a major desideratum. For the past century or so, it is common both in practice and in the academic literature to jump to single transferable vote (STV) as the solution for achieving PR. Some of the most prominent electoral reform movements around the globe are pushing for the adoption of STV. It has been termed a major open problem to design a voting rule that satisfies the same PR properties as STV and better monotonicity properties. In this paper, we first present a taxonomy of proportional representation axioms for general weak order preferences, some of which generalise and strengthen previously introduced concepts. We then present a rule called Expanding Approvals Rule (EAR) that satisfies properties stronger than the central PR axiom satisfied by STV, can handle indifferences in a convenient and computationally efficient manner, and also satisfies better candidate monotonicity properties. In view of this, our proposed rule seems to be a compelling solution for achieving proportional representation in voting settings.


title: Accurate parameter estimation for Bayesian Network Classifiers using Hierarchical Dirichlet Processes

abstract: This paper introduces a novel parameter estimation method for the probability tables of Bayesian network classifiers (BNCs), using hierarchical Dirichlet processes (HDPs). The main result of this paper is to show that improved parameter estimation allows BNCs to outperform leading learning methods such as Random Forest for both 0-1 loss and RMSE, albeit just on categorical datasets. As data assets become larger, entering the hyped world of "big", efficient accurate classification requires three main elements: (1) classifiers with low-bias that can capture the fine-detail of large datasets (2) out-of-core learners that can learn from data without having to hold it all in main memory and (3) models that can classify new data very efficiently. The latest Bayesian network classifiers (BNCs) satisfy these requirements. Their bias can be controlled easily by increasing the number of parents of the nodes in the graph. Their structure can be learned out of core with a limited number of passes over the data. However, as the bias is made lower to accurately model classification tasks, so is the accuracy of their parameters' estimates, as each parameter is estimated from ever decreasing quantities of data. In this paper, we introduce the use of Hierarchical Dirichlet Processes for accurate BNC parameter estimation. We conduct an extensive set of experiments on 68 standard datasets and demonstrate that our resulting classifiers perform very competitively with Random Forest in terms of prediction, while keeping the out-of-core capability and superior classification time.


title: Discovery of Striking Difference of Molecular-Emission-Line Richness in the Potential Proto-Binary System NGC 2264 CMM3

abstract: We have conducted an interferometric line survey in the 0.8 mm band toward the young high-mass protostar candidate NGC 2264 CMM3 with ALMA. CMM3 is resolved into the two continuum peaks, CMM3A and CMM3B, at an angular separation of 0.9". Thus, CMM3 is found to be a binary system candidate. We have detected molecular outflows associated with CMM3A and CMM3B each, indicating active star formation. In addition to the two peaks, six faint continuum peaks are detected around CMM3A and CMM3B, most of which are thought to be evolved low-mass protostars. CMM3A is found to be rich in molecular line emission including complex organic molecules such as HCOOCH3 and CH3OCH3. The emission of complex organic molecules is distributed within a compact region around the continuum peak of CMM3A. Hence, CMM3A apparently harbors a hot core. On the other hand, CMM3B is deficient in molecular line emission, although its continuum flux is almost comparable to that of CMM3A. Possible origins of the striking difference between CMM3A and CMM3B are discussed.


title: Learning to Blame: Localizing Novice Type Errors with Data-Driven Diagnosis

abstract: Localizing type errors is challenging in languages with global type inference, as the type checker must make assumptions about what the programmer intended to do. We introduce Nate, a data-driven approach to error localization based on supervised learning. Nate analyzes a large corpus of training data -- pairs of ill-typed programs and their "fixed" versions -- to automatically learn a model of where the error is most likely to be found. Given a new ill-typed program, Nate executes the model to generate a list of potential blame assignments ranked by likelihood. We evaluate Nate by comparing its precision to the state of the art on a set of over 5,000 ill-typed OCaml programs drawn from two instances of an introductory programming course. We show that when the top-ranked blame assignment is considered, Nate's data-driven model is able to correctly predict the exact sub-expression that should be changed 72% of the time, 28 points higher than OCaml and 16 points higher than the state-of-the-art SHErrLoc tool. Furthermore, Nate's accuracy surpasses 85% when we consider the top two locations and reaches 91% if we consider the top three.


title: Role of dimensional crossover on spin-orbit torque efficiency in magnetic insulator thin films

abstract: Magnetic insulators (MIs) attract tremendous interest for spintronic applications due to low Gilbert damping and absence of Ohmic loss. Magnetic order of MIs can be manipulated and even switched by spin-orbit torques (SOTs) generated through spin Hall effect and Rashba-Edelstein effect in heavy metal/MI bilayers. SOTs on MIs are more intriguing than magnetic metals since SOTs cannot be transferred to MIs through direct injection of electron spins. Understanding of SOTs on MIs remains elusive, especially how SOTs scale with the film thickness. Here, we observe the critical role of dimensionality on the SOT efficiency by systematically studying the MI layer thickness dependent SOT efficiency in tungsten/thulium iron garnet (W/TmIG) bilayers. We first show that the TmIG thin film evolves from two-dimensional to three-dimensional magnetic phase transitions as the thickness increases, due to the suppression of long-wavelength thermal fluctuation. Then, we report the significant enhancement of the measured SOT efficiency as the thickness increases. We attribute this effect to the increase of the magnetic moment density in concert with the suppression of thermal fluctuations. At last, we demonstrate the current-induced SOT switching in the W/TmIG bilayers with a TmIG thickness up to 15 nm. The switching current density is comparable with those of heavy metal/ferromagnetic metal cases. Our findings shed light on the understanding of SOTs in MIs, which is important for the future development of ultrathin MI-based low-power spintronics.


title: Haircutting Non-cash Collateral

abstract: Haircutting non-cash collateral has become a key element of the post-crisis reform of the shadow banking system and OTC derivatives markets. This article develops a parametric haircut model by expanding haircut definitions beyond the traditional value-at-risk measure and employing a double-exponential jump-diffusion model for collateral market risk. Haircuts are solved to target credit risk measurements, including probability of default, expected loss or unexpected loss criteria. Comparing to data-driven approach typically run on proxy data series, the model enables sensitivity analysis and stress test, captures market liquidity risk, allows idiosyncratic risk adjustments, and incorporates relevant market information. Computational results for main equities, securitization, and corporate bonds show potential for uses in collateral agreements, e.g. CSAs, and for regulatory capital calculations.


title: Fast Locality-Sensitive Hashing Frameworks for Approximate Near Neighbor Search

abstract: The Indyk-Motwani Locality-Sensitive Hashing (LSH) framework (STOC 1998) is a general technique for constructing a data structure to answer approximate near neighbor queries by using a distribution $\mathcal{H}$ over locality-sensitive hash functions that partition space. For a collection of $n$ points, after preprocessing, the query time is dominated by $O(n^{\rho} \log n)$ evaluations of hash functions from $\mathcal{H}$ and $O(n^{\rho})$ hash table lookups and distance computations where $\rho \in (0,1)$ is determined by the locality-sensitivity properties of $\mathcal{H}$. It follows from a recent result by Dahlgaard et al. (FOCS 2017) that the number of locality-sensitive hash functions can be reduced to $O(\log^2 n)$, leaving the query time to be dominated by $O(n^{\rho})$ distance computations and $O(n^{\rho} \log n)$ additional word-RAM operations. We state this result as a general framework and provide a simpler analysis showing that the number of lookups and distance computations closely match the Indyk-Motwani framework, making it a viable replacement in practice. Using ideas from another locality-sensitive hashing framework by Andoni and Indyk (SODA 2006) we are able to reduce the number of additional word-RAM operations to $O(n^\rho)$.


title: Semiparametric GARCH via Bayesian model averaging

abstract: As the dynamic structure of the financial markets is subject to dramatic changes, a model capable of providing consistently accurate volatility estimates must not make strong assumptions on how prices change over time. Most volatility models impose a particular parametric functional form that relates an observed price change to a volatility forecast (news impact function). We propose a new class of functional coefficient semiparametric volatility models where the news impact function is allowed to be any smooth function, and study its ability to estimate volatilities compared to the well known parametric proposals, in both a simulation study and an empirical study with real financial data. We estimate the news impact function using a Bayesian model averaging approach, implemented via a carefully developed Markov chain Monte Carlo (MCMC) sampling algorithm. Using simulations we show that our flexible semiparametric model is able to learn the shape of the news impact function from the observed data. When applied to real financial time series, our new model suggests that the news impact functions are significantly different in shapes for different asset types, but are similar for the assets of the same type.


title: Central Moment Analogues to Linear Optically Generated Cluster States

abstract: Two-mode squeezed states in the limit of small squeezing, Hong-Ou-Mandel interference and post selection on coincidence counts are some of the staples of linear quantum optics. We show that by using classical expectations on central moments of intensities, we can remove the requirement of small squeezing necessary for high fidelity coincidence detection. Utilizing existing techniques to probabilistically generate a cluster state, we construct a statistical analogue with deterministic generation at the cost of losing the ability to feed forward and requiring statistical averaging.


title: Analyzing Exoplanet Phase Curve Information Content: Toward Optimized Observing Strategies

abstract: Secondary eclipses and phase curves reveal information about the reflectivity and heat distribution in exoplanet atmospheres. The phase curve is composed of a combination of reflected, and thermally emitted light from the planet, and for circular orbits the phase curve peaks during the secondary eclipse or at an orbital phase of 0.5. Physical mechanisms have been discovered which shift the phase curve maximum of tidally locked close in planets to the right, or left, of the secondary eclipse. These mechanisms include cloud formations, and atmospheric superrotation, both of which serve to shift the thermally bright hot-spot, or highly reflective bright spot, of the atmosphere away from the sub-stellar point. Here we present a methodology for optimizing observing strategies for both secondary eclipses and phase curves with the goal of maximizing the information gained about the planetary atmosphere while minimizing the (assumed) continuous observation time. We show that we can increase the duty cycle of observations aimed at the measurements of phase curve characteristics (Amplitude, Phase offset) by up to $50\%$ for future platforms such as CHEOPS and JWST. We apply this methodology to the test cases of the Spitzer phase curve of 55-Cancri-e, which displays an eastward shift in its phase curve maximum, as well as model-generated observations of an ultra-short period planet observed with CHEOPS.


title: Hierarchical Multi-scale Attention Networks for Action Recognition

abstract: Recurrent Neural Networks (RNNs) have been widely used in natural language processing and computer vision. Among them, the Hierarchical Multi-scale RNN (HM-RNN), a kind of multi-scale hierarchical RNN proposed recently, can learn the hierarchical temporal structure from data automatically. In this paper, we extend the work to solve the computer vision task of action recognition. However, in sequence-to-sequence models like RNN, it is normally very hard to discover the relationships between inputs and outputs given static inputs. As a solution, attention mechanism could be applied to extract the relevant information from input thus facilitating the modeling of input-output relationships. Based on these considerations, we propose a novel attention network, namely Hierarchical Multi-scale Attention Network (HM-AN), by combining the HM-RNN and the attention mechanism and apply it to action recognition. A newly proposed gradient estimation method for stochastic neurons, namely Gumbel-softmax, is exploited to implement the temporal boundary detectors and the stochastic hard attention mechanism. To amealiate the negative effect of sensitive temperature of the Gumbel-softmax, an adaptive temperature training method is applied to better the system performance. The experimental results demonstrate the improved effect of HM-AN over LSTM with attention on the vision task. Through visualization of what have been learnt by the networks, it can be observed that both the attention regions of images and the hierarchical temporal structure can be captured by HM-AN.


title: Conditional Lower Bound for Subgraph Isomorphism with a Tree Pattern

abstract: The kTree problem is a special case of Subgraph Isomorphism where the pattern graph is a tree, that is, the input is an $n$-node graph $G$ and a $k$-node tree $T$, and the goal is to determine whether $G$ has a subgraph isomorphic to $T$. We provide evidence that this problem cannot be computed significantly faster than $2^{k} \textsf{poly}(n)$, which matches the fastest algorithm known for this problem by Koutis and Williams [ICALP 2009 and TALG 2016]. Specifically, we show that if kTree can be solved in time $(2-\varepsilon)^k \textsf{poly}(n)$ for some constant $\varepsilon>0$, then Set Cover with $n'$ elements and $m'$ sets can be solved in time $(2-\delta)^{n'} \textsf{poly}(m')$ for a constant $\delta(\varepsilon) > 0$, which would refute the Set Cover Conjecture by Cygan et al. [CCC 2012 and TALG 2016]. Our techniques yield a new algorithm for the p-Partial Cover problem, a parameterized version of Set Cover that requires covering at least $p$ elements (rather than all elements). Its running time is $(2+\varepsilon)^p (m')^{O(1/\varepsilon)}$ for any fixed $\varepsilon>0$, which improves the previous $2.597^p \textsf{poly}(m')$-time algorithm by Zehavi [ESA 2015]. Our running time is nearly optimal, as a $(2-\varepsilon')^p \textsf{poly}(m')$-time algorithm would refute the Set Cover Conjecture.


title: Sequential Decision Model for Inference and Prediction on Non-Uniform Hypergraphs with Application to Knot Matching from Computational Forestry

abstract: In this paper, we consider the knot matching problem arising in computational forestry. The knot matching problem is an important problem that needs to be solved to advance the state of the art in automatic strength prediction of lumber. We show that this problem can be formulated as a quadripartite matching problem and develop a sequential decision model that admits efficient parameter estimation along with a sequential Monte Carlo sampler on graph matching that can be utilized for rapid sampling of graph matching. We demonstrate the effectiveness of our methods on 30 manually annotated boards and present findings from various simulation studies to provide further evidence supporting the efficacy of our methods.


title: Exotic Bifurcations Inspired by Walking Droplet Dynamics

abstract: We identify two rather novel types of (compound) dynamical bifurcations generated primarily by interactions of an invariant attracting submanifold with stable and unstable manifolds of hyperbolic fixed points. These bifurcation types - inspired by recent investigations of mathematical models for walking droplet (pilot-wave) phenomena - are introduced and illustrated. Some of the one-parameter bifurcation types are analyzed in detail and extended from the plane to higher-dimensional spaces. A few applications to walking droplet dynamics are analyzed.


title: SUSY Global Fits

abstract: We present comprehensive global fits of supersymmetric (SUSY) models from the Global and Modular Beyond-the-Standard-Model Inference Tool (GAMBIT) collaboration, based on arXiv:1705.07935 and arXiv:1705.07917. We investigate several variants of the minimal supersymmetric standard model, a fully constrained version (CMSSM) with universal scalar ($m_0$), gaugino ($m_{1/2}$) and trilinear masses ($A_0$) at the gauge coupling unification scale, a similar model that is relaxed by adding an extra parameter for the soft Higgs masses (NUHM1), another where the soft Higgs masses are also split (NUHM2) and finally a weak scale MSSM7 model. We use the public GAMBIT global fitting framework and take into account all relevant data to reveal the regions of parameter space with the highest likelihood. Our results reveal that all models have very heavy scenarios that are well out of reach of the LHC, but will be probed by forthcoming dark matter experiments, as well as a stop-co-annihilation region which has better prospects for detection in collider experiments. The stau co-annihilation region is excluded from the CMSSM at $2 \sigma$ but is present in the NUHM1 and NUHM2 variants. Finally by relaxing constraints in the NUHM1, NUHM2 and MSSM7 we see additional regions appear: lighter chargino co-annihilation region, sbottom co-annihilation and $h/Z$ funnels.


title: Detection of the number of principal components by extended AIC-type method

abstract: Estimating the number of principal components is one of the fundamental problems in many scientific fields such as signal processing (or the spiked covariance model). In this paper, we first demonstrate that, for fixed $p$, any penalty term of the form $k'(p-k'/2+1/2)C_n$ may lead to an asymptotically consistent estimator under the condition that $C_n\to\infty$ and $C_n/n\to0$. We also extend our results to the case $n,p\to\infty$, with $p/n\to c>0$. In this case, for $k=o(n^{\frac{1}{3}})$, we first investigate the limiting laws for the leading eigenvalues of the sample covariance matrix $S_n$ under the condition that $\lambda_k>1+\sqrt{c}$. At low SNR, since the AIC tends to underestimate the number of signals $k$, the AIC should be re-defined in this case. As a natural extension of the AIC for fixed $p$, we propose the extended AIC (EAIC), i.e., the AIC-type method with tuning parameter $\gamma=\varphi(c)=1/2+\sqrt{1/c}-\log(1+\sqrt{c})/c$, and demonstrate that the EAIC-type method, i.e., the AIC-type method with tuning parameter $\gamma>\varphi(c)$, can select the number of signals $k$ consistently. In the following two cases, (1) $p$ fixed, $n\to\infty$, (2) $n,p\to\infty$ with $p/n\to 0$, if the AIC is defined as the degeneration of the EAIC in the case $n,p\to\infty$ with $p/n\to c>0$, i.e., $\gamma=\lim_{c\rightarrow 0+0}\varphi(c)=1$, then we have essentially demonstrated that, to achieve the consistency of the AIC-type method in the above two cases, $\gamma>1$ is required. Moreover, we show that the EAIC-type method is essentially tuning-free and outperforms the well-known KN estimator proposed in Kritchman and Nadler (2008) and the BCF estimator proposed in Bai, Choi and Fujikoshi (2018). Numerical studies indicate that the proposed method works well.


title: High-order harmonic generation from 2D periodic potentials in circularly and bichromatic circularly polarized laser fields

abstract: We studied the high-order harmonic generation (HHG) from 2D solid materials in circularly and bichromatic circularly polarized laser fields numerically by simulating the dynamics of single-active-electron processes in 2D periodic potentials. Contrary to the absence of HHG in the atomic case, circular HHGs below the bandgap with different helicities are produced from intraband transitions in solids with $C_4$ symmetry driven by circularly polarized lasers. Harmonics above the bandgap are elliptically polarized due to the interband transitions. High-order elliptically polarized harmonics can be generated efficiently by both co-rotating and counter-rotating bicircular mid-infrared lasers. The cutoff energy, ellipticity, phase, and intensity of the harmonics can be tuned by the control of the relative phase difference between the 1$\omega$ and 2$\omega$ fields in bicircularly polarized lasers, which can be utilized to image the structure of solids.


title: Spectral and Combinatorial Properties of Some Algebraically Defined Graphs

abstract: Let $k\ge 3$ be an integer, $q$ be a prime power, and $\mathbb{F}_q$ denote the field of $q$ elements. Let $f_i, g_i\in\mathbb{F}_q[X]$, $3\le i\le k$, such that $g_i(-X) = -\, g_i(X)$. We define a graph $S(k,q) = S(k,q;f_3,g_3,\cdots,f_k,g_k)$ as a graph with the vertex set $\mathbb{F}_q^k$ and edges defined as follows: vertices $a = (a_1,a_2,\ldots,a_k)$ and $b = (b_1,b_2,\ldots,b_k)$ are adjacent if $a_1\ne b_1$ and the following $k-2$ relations on their components hold: $$ b_i-a_i = g_i(b_1-a_1)f_i\Bigl(\frac{b_2-a_2}{b_1-a_1}\Bigr)\;,\quad 3\le i\le k. $$ We show that graphs $S(k,q)$ generalize several recently studied examples of regular expanders and can provide many new such examples.


title: The rainbow connection number of enhanced power graph

abstract: Let $G$ be a finite group, the enhanced power graph of $G$, denoted by $\Gamma_G^e$, is the graph with vertex set $G$ and two vertices $x,y$ are edge connected in $\Gamma_{G}^e$ if there exist $z\in G$ such that $x,y\in\langle z\rangle$. Let $\zeta$ be a edge-coloring of $\Gamma_G^e$. In this article, we calculate the rainbow connection number of the enhanced power graph $\Gamma_G^e$.


title: POEMMA: Probe Of Extreme Multi-Messenger Astrophysics

abstract: The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) mission is being designed to establish charged-particle astronomy with ultra-high energy cosmic rays (UHECRs) and to observe cosmogenic tau neutrinos (CTNs). The study of UHECRs and CTNs from space will yield orders-of-magnitude increase in statistics of observed UHECRs at the highest energies, and the observation of the cosmogenic flux of neutrinos for a range of UHECR models. These observations should solve the long-standing puzzle of the origin of the highest energy particles ever observed, providing a new window onto the most energetic environments and events in the Universe, while studying particle interactions well beyond accelerator energies. The discovery of CTNs will help solve the puzzle of the origin of UHECRs and begin a new field of Astroparticle Physics with the study of neutrino properties at ultra-high energies.


title: Classification of symmetry properties of waveguide modes in presence of gain/losses, anisotropy/bianisotropy, or continuous/discrete rotational symmetry

abstract: We study the symmetric properties of waveguide modes in presence of gain/losses, anisotropy/bianisotropy, or continuous/discrete rotational symmetry. We provide a comprehensive approach to identity the modal symmetry by constructing a $4\times4$ waveguide Hamiltonian and searching the symmetric operation in association with the corresponding waveguides. We classify the chiral/time reversal/parity/parity time/rotational symmetry for different waveguides, and provide the criterion for the aforementioned symmetry operations. Lastly, we provide examples to illustrate how the symmetry operations can be used to classify the modal properties from the symmetric relation between modal profiles of several different waveguides.


title: A wavelet frame coefficient total variational model for image restoration

abstract: In this paper, we propose a vector total variation (VTV) of feature image model for image restoration. The VTV imposes different smoothing powers on different features (e.g. edges and cartoons) based on choosing various regularization parameters. Thus, the model can simultaneously preserve edges and remove noises. Next, the existence of solution for the model is proved and the split Bregman algorithm is used to solve the model. At last, we use the wavelet filter banks to explicitly define the feature operator and present some experimental results to show its advantage over the related methods in both quality and efficiency.


title: New leaves of the tree: percolation analysis for cosmic web with discrete points

abstract: Percolation analysis has long been used to quantify the connectivity of the cosmic web. Most of the previous work is based on density fields on grids. By smoothing into fields, we lose information about galaxy properties like shape or luminosity. Lack of mathematical model also limits our understanding of percolation analysis. In order to overcome these difficulties, we have studied percolation analysis based on discrete points. Using a Friends-of-Friends (FoF) algorithm, we generate the S-bb relation, between the fractional mass of the largest connected group (S) and the FoF linking length (bb). We propose a new model, the Probability Cloud Cluster Expansion Theory (PCCET) to relate the S-bb relation with correlation functions. We show that the S-bb relation reflects a combination of all orders of correlation functions. Using N-body simulation, we find that the S-bb relation is robust against redshift distortion and incompleteness in observation. From the Bolshoi simulation, with Halo Abundance Matching (HAM), we have generated a mock galaxy catalogue. Good matching of the projected two-point correlation function with observation is confirmed. However, comparing the mock catalogue with the latest galaxy catalogue from SDSS DR12, we have found significant differences in their S-bb relations. This indicates that the mock galaxy catalogue cannot accurately retain higher order correlation functions than the two-point correlation function, which reveals the limit of HAM method. As a new measurement, S-bb relation is applicable to a wide range of data types, fast to compute, robust against redshift distortion and incompleteness, and it contains information of all orders of correlation function.


title: A Data-Driven Distributionally Robust Bound on the Expected Optimal Value of Uncertain Mixed 0-1 Linear Programming

abstract: This paper studies the expected optimal value of a mixed 0-1 programming problem with uncertain objective coefficients following a joint distribution. We assume that the true distribution is not known exactly, but a set of independent samples can be observed. Using the Wasserstein metric, we construct an ambiguity set centered at the empirical distribution from the observed samples and containing the true distribution with a high statistical guarantee. The problem of interest is to investigate the bound on the expected optimal value over the Wasserstein ambiguity set. Under standard assumptions, we reformulate the problem into a copositive program, which naturally leads to a tractable semidefinite-based approximation. We compare our approach with a moment-based approach from the literature on three applications. Numerical results illustrate the effectiveness of our approach.


title: A Mixed-Membership Model for Social Network Clustering

abstract: We propose a simple mixed membership model for social network clustering in this paper. A flexible function is adopted to measure affinities among a set of entities in a social network. The model not only allows each entity in the network to possess more than one membership, but also provides accurate statistical inference about network structure. We estimate the membership parameters using an MCMC algorithm. We evaluate the performance of the proposed algorithm by applying our model to two empirical social network data, the Zachary club data and the bottlenose dolphin network data. We also conduct some numerical studies based on synthetic networks for further assessing the effectiveness of our algorithm. In the end, some concluding remarks and future work are addressed briefly.


title: Spatial- and Frequency-Wideband Effects in Millimeter-Wave Massive MIMO Systems

abstract: When there are a large number of antennas in massive MIMO systems, the transmitted wideband signal will be sensitive to the physical propagation delay of electromagnetic waves across the large array aperture, which is called the spatial-wideband effect. In this scenario, transceiver design is different from most of the existing works, which presume that the bandwidth of the transmitted signals is not that wide, ignore the spatial-wideband effect, and only address the frequency selectivity. In this paper, we investigate spatial- and frequency-wideband effects, called dual-wideband effects, in massive MIMO systems from array signal processing point of view. Taking mmWave-band communications as an example, we describe the transmission process to address the dual-wideband effects. By exploiting the channel sparsity in the angle domain and the delay domain, we develop the efficient uplink and downlink channel estimation strategies that require much less amount of training overhead and cause no pilot contamination. Thanks to the array signal processing techniques, the proposed channel estimation is suitable for both TDD and FDD massive MIMO systems. Numerical examples demonstrate that the proposed transmission design for massive MIMO systems can effectively deal with the dual-wideband effects.


title: Critical role of water in defect aggregation and chemical degradation of perovskite solar cells

abstract: The chemical stability of methylammonium lead iodide (\ce{MAPbI3}) under humid conditions remains the primary challenge facing halide perovskite solar cells. We investigate defect processes in the water-intercalated iodide perovskite (\ce{MAPbI3}\_\ce{H2O}) and monohydrated phase (\ce{MAPbI3}$\cdot$\ce{H2O}) within a first-principles thermodynamic framework. We consider the formation energies of isolated and aggregated vacancy defects with different charge states under I-rich and I-poor conditions. It is found that a \ce{PbI2} (partial Schottky) vacancy complex can be formed readily, while the \ce{MAI} vacancy complex is difficult to form in the hydrous compounds. Vacancies in the hydrous phases create deep charge transition levels, indicating the degradation of halide perovskite upon exposure to moisture. Electronic structure analysis supports a novel mechanism of water-mediated vacancy-pair formation.


title: Reinforcement Mechanism Design for e-commerce

abstract: We study the problem of allocating impressions to sellers in e-commerce websites, such as Amazon, eBay or Taobao, aiming to maximize the total revenue generated by the platform. We employ a general framework of reinforcement mechanism design, which uses deep reinforcement learning to design efficient algorithms, taking the strategic behaviour of the sellers into account. Specifically, we model the impression allocation problem as a Markov decision process, where the states encode the history of impressions, prices, transactions and generated revenue and the actions are the possible impression allocations in each round. To tackle the problem of continuity and high-dimensionality of states and actions, we adopt the ideas of the DDPG algorithm to design an actor-critic policy gradient algorithm which takes advantage of the problem domain in order to achieve convergence and stability. We evaluate our proposed algorithm, coined IA(GRU), by comparing it against DDPG, as well as several natural heuristics, under different rationality models for the sellers - we assume that sellers follow well-known no-regret type strategies which may vary in their degree of sophistication. We find that IA(GRU) outperforms all algorithms in terms of the total revenue.


title: Influence of water intercalation and hydration on chemical decomposition and ion transport in methylammonium lead halide perovskites

abstract: The use of methylammonium (MA) lead halide perovskites \ce{CH3NH3PbX3} (X=I, Br, Cl) in perovskite solar cells (PSCs) has made great progress in performance efficiency during recent years. However, the rapid decomposition of \ce{MAPbI3} in humid environments hinders outdoor application of PSCs, and thus, a comprehensive understanding of the degradation mechanism is required. To do this, we investigate the effect of water intercalation and hydration of the decomposition and ion migration of \ce{CH3NH3PbX3} using first-principles calculations. We find that water interacts with \ce{PbX6} and MA through hydrogen bonding, and the former interaction enhances gradually, while the latter hardly changes when going from X=I to Br and to Cl. Thermodynamic calculations indicate that water exothermically intercalates into the perovskite, while the water intercalated and monohydrated compounds are stable with respect to decomposition. More importantly, the water intercalation greatly reduces the activation energies for vacancy-mediated ion migration, which become higher going from X=I to Br and to Cl. Our work indicates that hydration of halide perovskites must be avoided to prevent the degradation of PSCs upon moisture exposure.


title: Clique-based Method for Social Network Clustering

abstract: In this article, we develop a clique-based method for social network clustering. We introduce a new index to evaluate the quality of clustering results, and propose an efficient algorithm based on recursive bipartition to maximize an objective function of the proposed index. The optimization problem is NP-hard, so we approximate the semi-optimal solution via an implicitly restarted Lanczos method. One of the advantages of our algorithm is that the proposed index of each community in the clustering result is guaranteed to be higher than some predetermined threshold, $p$, which is completely controlled by users. We also account for the situation that $p$ is unknown. A statistical procedure of controlling both under-clustering and over-clustering errors simultaneously is carried out to select localized threshold for each subnetwork, such that the community detection accuracy is optimized. Accordingly, we propose a localized clustering algorithm based on binary tree structure. Finally, we exploit the stochastic blockmodels to conduct simulation studies and demonstrate the accuracy and efficiency of our algorithms, both numerically and graphically.


title: On the Hartshorne-Hirschowitz theorem

abstract: The Hartshorne--Hirschowitz theorem says that a generic union of lines in $\mathbb{P}^n$, $(n\geq 3)$, has good postulation. The proof of Hartshorne and Hirschowitz in the initial case $\mathbb{P}^3$ is difficult and so long, which is handled by a method of specialization via a smooth quadric surface with the property of having two rulings of skew lines. We provide a proof in the case $\mathbb{P}^3$ based on a new degeneration of disjoint lines via a plane $H\cong\mathbb{P}^2$, which we call $(2,s)$-cone configuration, that is a schematic union of $s$ intersecting lines passing through a single point $P$ together with the trace of an $s$-multiple point supported at $P$ on the double plane $2H$. In the first part of this paper, we discuss our degeneration inductive approach. We prove that a $(2,s)$-cone configuration is a degeneration of $s$ disjoint lines in $\mathbb{P}^3$, or more generally in $\mathbb{P}^n$. In the second part of the paper, we use this degeneration in an effective method to show that a generic union of lines in $\mathbb{P}^3$ imposes independent conditions on the linear system $|\mathbb{O}_{\mathbb{P}^3}(d)|$ of surfaces of given degree $d$. The basic motivation behind our degeneration approach is that it looks more systematic that gives some hope of extensions to the analogous problem in higher dimensional spaces, that is the postulation problem for $m$-dimensional planes in $\mathbb{P}^{2m+1}$.


title: The Young Substellar Companion ROXs 12 B: Near-Infrared Spectrum, System Architecture, and Spin-Orbit Misalignment

abstract: ROXs 12 (2MASS J16262803-2526477) is a young star hosting a directly imaged companion near the deuterium-burning limit. We present a suite of spectroscopic, imaging, and time-series observations to characterize the physical and environmental properties of this system. Moderate-resolution near-infrared spectroscopy of ROXs 12 B from Gemini-North/NIFS and Keck/OSIRIS reveals signatures of low surface gravity including weak alkali absorption lines and a triangular $H$-band pseudo-continuum shape. No signs of Pa$\beta$ emission are evident. As a population, however, we find that about half (46 $\pm$ 14\%) of young ($\lesssim$15 Myr) companions with masses $\lesssim$20 $M_\mathrm{Jup}$ possess actively accreting subdisks detected via Pa$\beta$ line emission, which represents a lower limit on the prevalence of circumplanetary disks in general as some are expected to be in a quiescent phase of accretion. The bolometric luminosity of the companion and age of the host star (6$^{+4}_{-2}$ Myr) imply a mass of 17.5 $\pm$ 1.5 $M_\mathrm{Jup}$ for ROXs 12 B based on hot-start evolutionary models. We identify a wide (5100 AU) tertiary companion to this system, 2MASS J16262774-2527247, which is heavily accreting and exhibits stochastic variability in its $K2$ light curve. By combining $v$sin$i_*$ measurements with rotation periods from $K2$, we constrain the line-of-sight inclinations of ROXs 12 A and 2MASS J16262774-2527247 and find that they are misaligned by 60$^{+7}_{-11}$$^{\circ}$. In addition, the orbital axis of ROXs 12 B is likely misaligned from the spin axis of its host star ROXs 12 A, suggesting that ROXs 12 B formed akin to fragmenting binary stars or in an equatorial disk that was torqued by the wide stellar tertiary.


title: Cooperative "folding transition" in the sequence space facilitates function-driven evolution of protein families

abstract: In the protein sequence space, natural proteins form clusters of families which are characterized by their unique native folds whereas the great majority of random polypeptides are neither clustered nor foldable to unique structures. Since a given polypeptide can be either foldable or unfoldable, a kind of "folding transition" is expected at the boundary of a protein family in the sequence space. By Monte Carlo simulations of a statistical mechanical model of protein sequence alignment that coherently incorporates both short-range and long-range interactions as well as variable-length insertions to reproduce the statistics of the multiple sequence alignment of a given protein family, we demonstrate the existence of such transition between natural-like sequences and random sequences in the sequence subspaces for 15 domain families of various folds. The transition was found to be highly cooperative and two-state-like. Furthermore, enforcing or suppressing consensus residues on a few of the well-conserved sites enhanced or diminished, respectively, the natural-like pattern formation over the entire sequence. In most families, the key sites included ligand binding sites. These results suggest some selective pressure on the key residues, such as ligand binding activity, may cooperatively facilitate the emergence of a protein family during evolution. From a more practical aspect, the present results highlight an essential role of long-range effects in precisely defining protein families, which are absent in conventional sequence models.


title: Broad pore lifetime distributions: A fundamental concept for cell electroporation

abstract: We describe a concept that has the potential to change how we think about the underlying mechanisms of cell membrane electroporation (EP). Prior experimental, theoretical and modeling have emphasized a single pore lifetime as adequate for particular conditions. Here we introduce a much more complex response: The rapid creation of many types of pore structures, of which some are traditional transient lipidic pores (TPs), but the great majority are complex pores (CPs) based on both lipids and other molecules or molecular segments. At the inner leaflet of the cell plasma membrane (PM) non-lipidic molecules come from the over-crowded cytoplasm. At the outer leaflet they originate from the extracellular medium and extracellular matrix. Some partially or fully insert into TPs during or shortly after TP formation, or bind to the membrane nearby. This process is complex, leading to mostly short-lived structures, with relatively few lasting for long times. We speculate that the characteristic pore lifetimes range from $\sim$100 ns to 1,000 s, based on implications from experiments. The frequency-of-occurrence probably falls off extremely rapidly with increasing lifetime, $\tau_{CP}$, which implies that most are inaccessible to traditional experimental methods. It also suggests that unexpected behavior can occur early in pulsing, vanishing before post-pulse observations begin.


title: Axion driven cosmic magneto-genesis during the QCD crossover

abstract: We propose a mechanism for the generation of a magnetic field in the early universe during the QCD crossover assuming that dark matter is made of axions. Thermoelectric fields arise at pressure gradients in the primordial plasma due to the difference in charge, energy density and equation of state between the quark and lepton components. The axion field is coupled to the EM field, so when its spatial gradient is misaligned with the thermoelectric field, an electric current is driven. Due to the finite resistivity of the plasma an electric field appears that is generally rotational. For a QCD axion mass consistent with observational constraints and a conventional efficiency for turbulent dynamo amplification --- driven by the same pressure gradients responsible for the thermoelectric fields --- a magnetic field is generated on subhorizon scales. After significant Alfv\'enic unwinding it reaches a present day strength of $B \sim 10^{-13 }$ G on a characteristic scale $L_B \sim $ 20 pc. The resulting combination of $BL_B^{1/2}$ is significantly stronger than in any astrophysical scenario, providing a clear test for the cosmological origin of the field through $\gamma$-ray observations of distant blazars. The amplitude of the pressure gradients may be inferred from the detection of concomitant gravitational waves, while several experiments are underway to confirm or rule out the existence of axions.


title: On the inevitability of the consistency operator

abstract: We examine recursive monotonic functions on the Lindenbaum algebra of $\mathsf{EA}$. We prove that no such function sends every consistent $\varphi$ to a sentence with deductive strength strictly between $\varphi$ and $(\varphi\wedge\mathsf{Con}(\varphi))$. We generalize this result to iterates of consistency into the effective transfinite. We then prove that for any recursive monotonic function $f$, if there is an iterate of $\mathsf{Con}$ that bounds $f$ everywhere, then $f$ must be somewhere equal to an iterate of $\mathsf{Con}$.


title: Design of Configurable Sequential Circuits in Quantum-dot Cellular Automata

abstract: Quantum-dot cellular automata (QCA) is a likely candidate for future low power nano-scale electronic devices. Sequential circuits in QCA attract more attention due to its numerous application in digital industry. On the other hand, configurable devices provide low device cost and efficient utilization of device area. Since the fundamental building block of any sequential logic circuit is flip flop, hence constructing configurable, multi-purpose QCA flip-flops are one of the prime importance of current research. This work proposes a design of configurable flip-flop (CFF) which is the first of its kind in QCA domain. The proposed flip-flop can be configured to D, T and JK flip-flop by configuring its control inputs. In addition, to make more efficient configurable flip-flop, a clock pulse generator (CPG) is designed which can trigger all types of edges (falling, rising and dual) of a clock. The same CFF design is used to realize an edge configurable (dual/rising/falling) flip- flop with the help of CPG. The biggest advantage of using edge configurable (dual/rising/falling) flip-flop is that it can be used in 9 different ways using the same single circuit. All the proposed designs are verified using QCADesigner simulator.


title: Miraculous cancellations for quantum $SL_2$

abstract: In earlier work, Helen Wong and the author discovered certain "miraculous cancellations" for the quantum trace map connecting the Kauffman bracket skein algebra of a surface to its quantum Teichmueller space, occurring when the quantum parameter $q$ is a root of unity. The current paper is devoted to giving a more representation theoretic interpretation of this phenomenon, in terms of the quantum group $U_q(sl_2)$ and its dual Hopf algebra $SL_2^q$.


title: Finite energy of generalized suitable weak solutions to the Navier-Stokes equations and Liouville-type theorems in two dimensional domains

abstract: Introducing a new notion of generalized suitable weak solutions, we first prove validity of the energy inequality for such a class of weak solutions to the Navier-Stokes equations in the whole space $\mathbb{R}^n$. Although we need certain growth condition on the pressure, we may treat the class even with infinite energy quantity except for the initial velocity. We next handle the equation for vorticity in 2D unbounded domains. Under a certain condition on the asymptotic behavior at infinity, we prove that the vorticity and its gradient of solutions are both globally square integrable. As their applications, Loiuville-type theorems are obtained.


title: Design of Adiabatic MTJ-CMOS Hybrid Circuits

abstract: Low-power designs are a necessity with the increasing demand of portable devices which are battery operated. In many of such devices the operational speed is not as important as battery life. Logic-in-memory structures using nano-devices and adiabatic designs are two methods to reduce the static and dynamic power consumption respectively. Magnetic tunnel junction (MTJ) is an emerging technology which has many advantages when used in logic-in-memory structures in conjunction with CMOS. In this paper, we introduce a novel adiabatic hybrid MTJ/CMOS structure which is used to design AND/NAND, XOR/XNOR and 1-bit full adder circuits. We simulate the designs using HSPICE with 32nm CMOS technology and compared it with a non-adiabatic hybrid MTJ/CMOS circuits. The proposed adiabatic MTJ/CMOS full adder design has more than 7 times lower power consumtion compared to the previous MTJ/CMOS full adder.


title: Fenchel Dual Gradient Methods for Distributed Convex Optimization over Time-varying Networks

abstract: In the large collection of existing distributed algorithms for convex multi-agent optimization, only a handful of them provide convergence rate guarantees on agent networks with time-varying topologies, which, however, restrict the problem to be unconstrained. Motivated by this, we develop a family of distributed Fenchel dual gradient methods for solving constrained, strongly convex but not necessarily smooth multi-agent optimization problems over time-varying undirected networks. The proposed algorithms are constructed based on the application of weighted gradient methods to the Fenchel dual of the multi-agent optimization problem, and can be implemented in a fully decentralized fashion. We show that the proposed algorithms drive all the agents to both primal and dual optimality asymptotically under a minimal connectivity condition and at sublinear rates under a standard connectivity condition. Finally, the competent convergence performance of the distributed Fenchel dual gradient methods is demonstrated via simulations.